From d006e513ef1d1eaea727aeeb08c3004066f30b68 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Tue, 2 Jan 2024 11:55:05 +0100
Subject: [PATCH] E. Jezequel 02/01/2024: bugfix ZPSI3 needed in case of KRR=0
 in turb

---
 src/common/turb/mode_turb_ver.F90 | 9 +++++++--
 1 file changed, 7 insertions(+), 2 deletions(-)

diff --git a/src/common/turb/mode_turb_ver.F90 b/src/common/turb/mode_turb_ver.F90
index d82a09460..8908cdb51 100644
--- a/src/common/turb/mode_turb_ver.F90
+++ b/src/common/turb/mode_turb_ver.F90
@@ -461,8 +461,13 @@ ENDIF
 GUSERV = KRR/=0
 !
 CALL PHI3(D,CSTURB,TURBN,ZREDTH1,ZREDR1,ZRED2TH3,ZRED2R3,ZRED2THR3,TURBN%CTURBDIM,GUSERV,ZPHI3)
-IF(KRR/=0) &
-CALL PSI3(D,CSTURB,TURBN,ZREDR1,ZREDTH1,ZRED2R3,ZRED2TH3,ZRED2THR3,TURBN%CTURBDIM,GUSERV,ZPSI3)
+IF(KRR/=0) THEN
+  CALL PSI3(D,CSTURB,TURBN,ZREDR1,ZREDTH1,ZRED2R3,ZRED2TH3,ZRED2THR3,TURBN%CTURBDIM,GUSERV,ZPSI3)
+ELSE
+  !$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+  ZPSI3(IIJB:IIJE,1:IKT) = 1.
+  !$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
+END IF
 !
 ! Prandtl numbers for scalars
 !
-- 
GitLab