Skip to content
Snippets Groups Projects
Commit d006e513 authored by RODIER Quentin's avatar RODIER Quentin
Browse files

E. Jezequel 02/01/2024: bugfix ZPSI3 needed in case of KRR=0 in turb

parent 5f2c6f12
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/common/turb/mode_turb_ver.F90
+ 7
− 2
View file @ d006e513
...@@ -461,8 +461,13 @@ ENDIF ...@@ -461,8 +461,13 @@ ENDIF
GUSERV = KRR/=0 GUSERV = KRR/=0
! !
CALL PHI3(D,CSTURB,TURBN,ZREDTH1,ZREDR1,ZRED2TH3,ZRED2R3,ZRED2THR3,TURBN%CTURBDIM,GUSERV,ZPHI3) CALL PHI3(D,CSTURB,TURBN,ZREDTH1,ZREDR1,ZRED2TH3,ZRED2R3,ZRED2THR3,TURBN%CTURBDIM,GUSERV,ZPHI3)
IF(KRR/=0) & IF(KRR/=0) THEN
CALL PSI3(D,CSTURB,TURBN,ZREDR1,ZREDTH1,ZRED2R3,ZRED2TH3,ZRED2THR3,TURBN%CTURBDIM,GUSERV,ZPSI3) 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 ! Prandtl numbers for scalars
! !
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme