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

Quentin 19/08/2022: bugfix bit-repro for SV variables (for MesoNH ARMLES test case)

parent 08b69d4a
No related branches found
No related tags found
No related merge requests found
...@@ -349,7 +349,7 @@ DO JK=1,D%NKT ...@@ -349,7 +349,7 @@ DO JK=1,D%NKT
ELSE ELSE
ZW2(JIJ,JK)=1. ZW2(JIJ,JK)=1.
END IF END IF
PREDTH1(JIJ,JK)= ZW2(JIJ,JK) * MAX(1.E-30, ZW2(JIJ,JK)*PREDTH1(JIJ,JK)) PREDTH1(JIJ,JK)= ZW2(JIJ,JK) * MAX(CST%XMNH_TINY_12, ZW2(JIJ,JK)*PREDTH1(JIJ,JK))
! !
IF (KRR /= 0) THEN ! moist case IF (KRR /= 0) THEN ! moist case
IF(PREDR1(JIJ,JK) < 0.) THEN IF(PREDR1(JIJ,JK) < 0.) THEN
...@@ -357,7 +357,7 @@ DO JK=1,D%NKT ...@@ -357,7 +357,7 @@ DO JK=1,D%NKT
ELSE ELSE
ZW2(JIJ,JK)=1. ZW2(JIJ,JK)=1.
END IF END IF
PREDR1(JIJ,JK)= ZW2(JIJ,JK) * MAX(1.E-30, ZW2(JIJ,JK)*PREDR1(JIJ,JK)) PREDR1(JIJ,JK)= ZW2(JIJ,JK) * MAX(CST%XMNH_TINY_12, ZW2(JIJ,JK)*PREDR1(JIJ,JK))
END IF END IF
ENDDO ENDDO
ENDDO ENDDO
...@@ -381,7 +381,7 @@ DO JSV=1,KSV ...@@ -381,7 +381,7 @@ DO JSV=1,KSV
ELSE ELSE
ZW2(JIJ,JK)=1. ZW2(JIJ,JK)=1.
END IF END IF
PREDS1(JIJ,JK,JSV)= ZW2(JIJ,JK) * MAX(1.E-30, ZW2(JIJ,JK)*PREDS1(JIJ,JK,JSV)) PREDS1(JIJ,JK,JSV)= ZW2(JIJ,JK) * MAX(CST%XMNH_TINY_12, ZW2(JIJ,JK)*PREDS1(JIJ,JK,JSV))
ENDDO ENDDO
ENDDO ENDDO
ENDDO ENDDO
......
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