From ba8704d6f22a3bb9b0bb320ec9457a9d50f81324 Mon Sep 17 00:00:00 2001 From: Christine Lac <christine.lac@meteo.fr> Date: Fri, 12 May 2017 14:16:31 +0200 Subject: [PATCH] Maud:avoid negative values for SVS --- src/MNH/shallow_mf_pack.f90 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/MNH/shallow_mf_pack.f90 b/src/MNH/shallow_mf_pack.f90 index 421a54d6b..a91a1d972 100644 --- a/src/MNH/shallow_mf_pack.f90 +++ b/src/MNH/shallow_mf_pack.f90 @@ -114,6 +114,7 @@ END MODULE MODI_SHALLOW_MF_PACK !! V.Masson 09/2010 !! Modification R. Honnert 07/2012 : introduction of vertical wind !! for the height of the thermal +!! M. Leriche 02/2017 : avoid negative values for sv tendencies !! -------------------------------------------------------------------------- ! !* 0. DECLARATIONS @@ -365,8 +366,8 @@ PRVS(:,:,:) = PRVS(:,:,:) +MYM( & DO JSV=1,ISV IF (LNOMIXLG .AND. JSV >= NSV_LGBEG .AND. JSV<= NSV_LGEND) CYCLE - PRSVS(:,:,:,JSV) = PRSVS(:,:,:,JSV) + & - PRHODJ(:,:,:)*ZDSVDT(:,:,:,JSV) + PRSVS(:,:,:,JSV) = MAX((PRSVS(:,:,:,JSV) + & + PRHODJ(:,:,:)*ZDSVDT(:,:,:,JSV)),XSVMIN(JSV)) END DO !!! 7. call to MesoNH budgets -- GitLab