From 40ef9ddf0a59939a5e28f6e697e2ffc790262b6e Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Fri, 13 Jan 2023 15:46:21 +0100
Subject: [PATCH] Quentin 13/01/2023: bugfix diag with sea SALT (by. P.Tulet
02/12/2022 mail)
---
src/MNH/phys_paramn.f90 | 12 ++++++++++++
1 file changed, 12 insertions(+)
diff --git a/src/MNH/phys_paramn.f90 b/src/MNH/phys_paramn.f90
index 241166607..c01f484c5 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,:))
--
GitLab