Soline B 6/06/2019: minor correction on sea salt: deposition and sea salt order

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

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