diff --git a/src/common/turb/mode_turb_ver.F90 b/src/common/turb/mode_turb_ver.F90
index d82a0946034a69fb93fb24a903621ef40c773d39..8908cdb51383e16c9face55a6ebf417b59e4bf34 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
 !