diff --git a/src/MNH/phys_paramn.f90 b/src/MNH/phys_paramn.f90 index 2411666072a4afe4840ebd7fa0738139875ebb78..c01f484c51e41dbe8029d077e254174719206c58 100644 --- a/src/MNH/phys_paramn.f90 +++ b/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,:))