From a8f8dd2a2f168f7ddaf4292f28c9b426a48725b9 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Fri, 19 Aug 2022 12:12:08 +0200
Subject: [PATCH] Quentin 19/08/2022: bugfix bit-repro for SV variables (for
MesoNH ARMLES test case)
---
src/common/turb/mode_prandtl.F90 | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/common/turb/mode_prandtl.F90 b/src/common/turb/mode_prandtl.F90
index d636ed778..f43bb30ab 100644
--- a/src/common/turb/mode_prandtl.F90
+++ b/src/common/turb/mode_prandtl.F90
@@ -349,7 +349,7 @@ DO JK=1,D%NKT
ELSE
ZW2(JIJ,JK)=1.
END IF
- PREDTH1(JIJ,JK)= ZW2(JIJ,JK) * MAX(1.E-30, ZW2(JIJ,JK)*PREDTH1(JIJ,JK))
+ PREDTH1(JIJ,JK)= ZW2(JIJ,JK) * MAX(CST%XMNH_TINY_12, ZW2(JIJ,JK)*PREDTH1(JIJ,JK))
!
IF (KRR /= 0) THEN ! moist case
IF(PREDR1(JIJ,JK) < 0.) THEN
@@ -357,7 +357,7 @@ DO JK=1,D%NKT
ELSE
ZW2(JIJ,JK)=1.
END IF
- PREDR1(JIJ,JK)= ZW2(JIJ,JK) * MAX(1.E-30, ZW2(JIJ,JK)*PREDR1(JIJ,JK))
+ PREDR1(JIJ,JK)= ZW2(JIJ,JK) * MAX(CST%XMNH_TINY_12, ZW2(JIJ,JK)*PREDR1(JIJ,JK))
END IF
ENDDO
ENDDO
@@ -381,7 +381,7 @@ DO JSV=1,KSV
ELSE
ZW2(JIJ,JK)=1.
END IF
- PREDS1(JIJ,JK,JSV)= ZW2(JIJ,JK) * MAX(1.E-30, ZW2(JIJ,JK)*PREDS1(JIJ,JK,JSV))
+ PREDS1(JIJ,JK,JSV)= ZW2(JIJ,JK) * MAX(CST%XMNH_TINY_12, ZW2(JIJ,JK)*PREDS1(JIJ,JK,JSV))
ENDDO
ENDDO
ENDDO
--
GitLab