Commit 9f0198f5 authored by Markus Pahlow's avatar Markus Pahlow

increase the threshold for division by 0 to 1e-10

parent 613287a5
......@@ -291,7 +291,7 @@ CONTAINS
REAL(dp), INTENT(INOUT) :: zN
REAL(dp) :: FNmax, fnF, fvF
FNmax = phy%F0N*felim*phy%fTNF(bct)
fnF = phy%sVPh/(SQRT(MAX(FNmax, 1E-30_dp)*phy%Q0P/phy%QN) + phy%sVPh) ! (E54)
fnF = phy%sVPh/(SQRT(MAX(FNmax, 1E-10_dp)*phy%Q0P/phy%QN) + phy%sVPh) ! (E54)
fvF = MAX(phy%QsN/phy%QN - phy%zF*(phy%QN - phy%Q0N), 0E0_dp)
phy%FixN = fvF*fnF*(1 - phy%Q0P/phy%QP)*FNmax
IF (phy%VN.LT.phy%FixN) THEN ! fix N2
......@@ -318,7 +318,7 @@ CONTAINS
REAL(dp), INTENT(INOUT) :: RC
REAL(dp), INTENT(IN) :: phin, zN
REAL(dp) :: paravg
phy%tch = phy%theta/MAX(phy%fC, 1E-30_dp) ! (A1)
phy%tch = phy%theta/MAX(phy%fC, 1E-10_dp) ! (A1)
IF (phy%PARb.GT.0E0_dp) THEN
phy%SI = physi(phy, phin, phy%PARs, phy%PARb, 1E0_dp)
ELSE
......
Markdown is supported
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