diff --git a/src/SURFEX/coupling_sltn.F90 b/src/SURFEX/coupling_sltn.F90
index 284fc093527a26c6342ca61c0ac9e7b190aa9b7b..4437552ac7edaedd31c20e0e3e905fb0299b2abc 100644
--- a/src/SURFEX/coupling_sltn.F90
+++ b/src/SURFEX/coupling_sltn.F90
@@ -297,21 +297,21 @@ DO JN = 1, JPMODE_SLT
   IF (LVARSIG_SLT) THEN ! cas 3 moment
 
     PSFSLT(:,2+(JN-1)*3) = ZSFSLT_MDE(:,JORDER_SLT(JN)) &
-                            * ((SLT%XEMISRADIUS_SLT(JORDER_SLT(JN))**3) &
-                            * EXP(4.5 * LOG(SLT%XEMISSIG_SLT(JORDER_SLT(JN)))**2)) &
+                            * ((SLT%XEMISRADIUS_SLT(JN)**3) &
+                            * EXP(4.5 * LOG(SLT%XEMISSIG_SLT(JN))**2)) &
                             * ZCONVERTFACM3_SLT
 
   ELSEIF (LRGFIX_SLT) THEN ! cas 1 moment
     PSFSLT(:,JN) =  ZSFSLT_MDE(:,JORDER_SLT(JN)) &
-                      * (SLT%XEMISRADIUS_SLT(JORDER_SLT(JN))**3) &
-                      * EXP(4.5 * LOG(SLT%XEMISSIG_SLT(JORDER_SLT(JN)))**2) &
+                      * (SLT%XEMISRADIUS_SLT(JN)**3) &
+                      * EXP(4.5 * LOG(SLT%XEMISSIG_SLT(JN))**2) &
                       * ZCONVERTFACM3_SLT
 
   ELSE ! cas 2 moments
 
     PSFSLT(:,2+(JN-1)*2) = ZSFSLT_MDE(:,JORDER_SLT(JN)) &
-                            * ((SLT%XEMISRADIUS_SLT(JORDER_SLT(JN))**3) &
-                            * EXP(4.5 * LOG(SLT%XEMISSIG_SLT(JORDER_SLT(JN)))**2)) &
+                            * ((SLT%XEMISRADIUS_SLT(JN)**3) &
+                            * EXP(4.5 * LOG(SLT%XEMISSIG_SLT(JN))**2)) &
                             * ZCONVERTFACM3_SLT
 ! -- PIERRE / MARINE SSA - MODIF --
   END IF
diff --git a/src/SURFEX/mode_dslt_surf.F90 b/src/SURFEX/mode_dslt_surf.F90
index c5e6d2fab1f7f7ac0145f4115f16d9f1a321e119..408f179f1c2fea06682802fac0b1df65a26c53c1 100644
--- a/src/SURFEX/mode_dslt_surf.F90
+++ b/src/SURFEX/mode_dslt_surf.F90
@@ -301,7 +301,7 @@ DO JN=1,SIZE(KM0)
     ZSIGMA(:) = PEMISSIG(JN)
 
     ZM(:,J0) = ZM(:,J3) /               &
-              ((PEMISRADIUS(KM3(JN))**3)*EXP(4.5 * LOG(ZSIGMA(:))**2))  
+              ((PEMISRADIUS(JN)**3)*EXP(4.5 * LOG(ZSIGMA(:))**2))  
 
   ELSE ! compute M6 from M0, M3 and SIGMA