Quentin 13/01/2023: bugfix diag with sea SALT (by. P.Tulet 02/12/2022 mail)

src/MNH/phys_paramn.f90

@@ -1028,9 +1028,15 @@ IF( CDCONV == 'KAFR' .OR. CSCONV == 'KAFR' ) THEN
         ZNDST(:,:,:,JSV)   = XN0MIN(IMODEIDX)
       ENDDO
       !
+    IF (CPROGRAM == "MESONH") THEN
       DO JSV=NSV_DSTBEG,NSV_DSTEND
         ZSVDST(:,:,:,JSV-NSV_DSTBEG+1) = XRSVS(:,:,:,JSV) * XTSTEP / XRHODJ(:,:,:) 
       ENDDO
+    ELSE
+      DO JSV=NSV_DSTBEG,NSV_DSTEND
+        ZSVDST(:,:,:,JSV-NSV_DSTBEG+1) = XSVT(:,:,:,JSV)
+      ENDDO
+    ENDIF
       CALL PPP2DUST(ZSVDST(IIB:IIE,IJB:IJE,IKB:IKE,:), XRHODREF(IIB:IIE,IJB:IJE,IKB:IKE),&
               PSIG3D=ZSIGDST(IIB:IIE,IJB:IJE,IKB:IKE,:), PRG3D=ZRGDST(IIB:IIE,IJB:IJE,IKB:IKE,:),   &
               PN3D=ZNDST(IIB:IIE,IJB:IJE,IKB:IKE,:))
@@ -1055,9 +1061,15 @@ IF( CDCONV == 'KAFR' .OR. CSCONV == 'KAFR' ) THEN
         ZNSLT(:,:,:,JSV)   = XN0MIN_SLT(IMODEIDX)
       ENDDO
       !
+    IF (CPROGRAM == "MESONH") THEN
       DO JSV=NSV_SLTBEG,NSV_SLTEND
         ZSVSLT(:,:,:,JSV-NSV_SLTBEG+1) = XRSVS(:,:,:,JSV) * XTSTEP / XRHODJ(:,:,:) 
       ENDDO
+    ELSE
+      DO JSV=NSV_SLTBEG,NSV_SLTEND
+        ZSVSLT(:,:,:,JSV-NSV_SLTBEG+1) = XSVT(:,:,:,JSV)
+      ENDDO
+    END IF
       CALL PPP2SALT(ZSVSLT(IIB:IIE,IJB:IJE,IKB:IKE,:), XRHODREF(IIB:IIE,IJB:IJE,IKB:IKE),&
               PSIG3D=ZSIGSLT(IIB:IIE,IJB:IJE,IKB:IKE,:), PRG3D=ZRGSLT(IIB:IIE,IJB:IJE,IKB:IKE,:),   &
               PN3D=ZNSLT(IIB:IIE,IJB:IJE,IKB:IKE,:))