Skip to content
Snippets Groups Projects
Commit 48cb2a5f authored by WAUTELET Philippe's avatar WAUTELET Philippe
Browse files

Philippe 20/06/2023: turb: fix from ZSOLVER/turb.f90 for GZD

parent 6b79ffe8
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment