diff --git a/src/MNH/shallow_mf_pack.f90 b/src/MNH/shallow_mf_pack.f90 index 421a54d6b3598337e2c498090c7bedbd62f2d385..a91a1d9722befa1e1af18c90a00d0707c8783ca8 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