From daa9c93755067fe6e9ba5fdc53fd2f8cd6e7034b Mon Sep 17 00:00:00 2001 From: Juan Escobar <escj@aero.obs-mip.fr> Date: Mon, 18 Sep 2017 09:47:24 +0200 Subject: [PATCH] Juan 18/09/2017: prandtl.f90, use epsilon XMNH_TINY_12 for R*4 --- src/MNH/prandtl.f90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/MNH/prandtl.f90 b/src/MNH/prandtl.f90 index a66ecfb7b..b21884eca 100644 --- a/src/MNH/prandtl.f90 +++ b/src/MNH/prandtl.f90 @@ -193,6 +193,7 @@ END MODULE MODI_PRANDTL !! change of YCOMMENT !! 2012-02 Y. Seity, add possibility to run with reversed !! vertical levels +!! 2017-09 J.Escobar, use epsilon XMNH_TINY_12 for R*4 !! -------------------------------------------------------------------------- ! !* 0. DECLARATIONS @@ -366,7 +367,7 @@ END DO ! DO JSV=1,ISV ZW2=SIGN(1.,PREDS1(:,:,:,JSV)) - PREDS1(:,:,:,JSV)= ZW2(:,:,:) * MAX(1.E-30, ZW2(:,:,:)*PREDS1(:,:,:,JSV)) + PREDS1(:,:,:,JSV)= ZW2(:,:,:) * MAX(XMNH_TINY_12, ZW2(:,:,:)*PREDS1(:,:,:,JSV)) END DO ! !--------------------------------------------------------------------------- -- GitLab