From b661b18447835710a7f3a3ef78e60ead741c3625 Mon Sep 17 00:00:00 2001
From: Gaelle Tanguy <gaelle.tanguy@meteo.fr>
Date: Tue, 22 Mar 2016 14:27:38 +0000
Subject: [PATCH] Gaelle 22/03/2016 : restore min value of ZSIGMA
---
src/MNH/condensation.f90 | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/src/MNH/condensation.f90 b/src/MNH/condensation.f90
index 13c9bdcd2..9c0c714eb 100644
--- a/src/MNH/condensation.f90
+++ b/src/MNH/condensation.f90
@@ -116,6 +116,7 @@ END MODULE MODI_CONDENSATION
!! Safety limitation to .99*Pressure for saturation vapour pressure
!! 2012-02 Y. Seity, add possibility to run with reversed vertical levels
!! 2015 C.Lac Change min value of ZSIGMA to be in agreement with AROME
+!! 2016 G.Delautier Restore min value of ZSIGMA (instability)
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
@@ -321,8 +322,8 @@ DO JK=IKTB,IKTE
ZA*ZA*ZDRW*ZDRW - 2.*ZA*ZB*ZDRW*ZDTL + ZB*ZB*ZDTL*ZDTL) + &
ZSIG_CONV * ZSIG_CONV ) )
END IF
-! ZSIGMA= MAX( 1.E-10, ZSIGMA )
- ZSIGMA= MAX( 1.E-12, ZSIGMA )
+ ZSIGMA= MAX( 1.E-10, ZSIGMA )
+! ZSIGMA= MAX( 1.E-12, ZSIGMA )
! normalized saturation deficit
ZQ1 = ZSBAR/ZSIGMA
--
GitLab