From 48cb2a5fecd3dcf95ef37c6c95cada8bfca4921b Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Tue, 20 Jun 2023 09:59:41 +0200
Subject: [PATCH] Philippe 20/06/2023: turb: fix from ZSOLVER/turb.f90 for GZD
---
src/MNH/turb.f90 | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/src/MNH/turb.f90 b/src/MNH/turb.f90
index af75646a8..fce209f97 100644
--- a/src/MNH/turb.f90
+++ b/src/MNH/turb.f90
@@ -2334,6 +2334,7 @@ IF (.NOT. ORMC01) THEN
ZD=ZALPHA*(PZZ(JI,JJ,KKTE+1)-PZZ(JI,JJ,JK))
IF ( (PLM(JI,JJ,JK)>ZD) .AND. GZD ) THEN
PLM(JI,JJ,JK)=ZD
+ ELSE
GZD = .FALSE.
ENDIF
END DO
@@ -2343,6 +2344,7 @@ IF (.NOT. ORMC01) THEN
ZD = ZALPHA * ( 0.5 * ( PZZ(JI, JJ, JK) + PZZ(JI, JJ, JK+KKL) ) - PZZ(JI, JJ, KKB) ) * PDIRCOSZW(JI, JJ)
IF ( ( PLM(JI,JJ,JK) > ZD ) .AND. GZD ) THEN
PLM(JI,JJ,JK) = ZD
+ ELSE
GZD = .FALSE.
ENDIF
END DO
@@ -2685,6 +2687,7 @@ IF (.NOT. ORMC01) THEN
ZD=ZALPHA*(PZZ(JI,JJ,KKTE+1)-PZZ(JI,JJ,JK))
IF ( (PLM(JI,JJ,JK)>ZD) .AND. GZD ) THEN
PLM(JI,JJ,JK)=ZD
+ ELSE
GZD = .FALSE.
ENDIF
END DO
@@ -2695,6 +2698,7 @@ IF (.NOT. ORMC01) THEN
*PDIRCOSZW(JI,JJ)
IF ( (PLM(JI,JJ,JK)>ZD) .AND. GZD ) THEN
PLM(JI,JJ,JK)=ZD
+ ELSE
GZD = .FALSE.
ENDIF
END DO
--
GitLab