Commit 0ff28e31 authored by Markus Pahlow's avatar Markus Pahlow

prevent spurious high tch under extreme Fe limitation

parent 3abbae33
......@@ -341,9 +341,9 @@ CONTAINS
PARb = MAX(phy%I0, phy%PARb)
depfrac = LOG(phy%PARs/PARb)
paravg = (phy%PARs - PARb)/depfrac
aim = phy%alpha*paravg/phy%mu0
phy%tch = 1E0_dp/phy%zC + (1E0_dp - lambertw((1E0_dp + RC/(phy%daylen*phy%mu0))&
*EXP(1E0_dp + MIN(aim/phy%zC, 2E2_dp)), 0, 0))/aim
aim = MIN(phy%alpha*paravg/phy%mu0, 7.e2*phy%zC)
phy%tch = 1.E0_dp/phy%zC + (1E0_dp - lambertw((1.E0_dp + RC/(phy%daylen*phy%mu0))&
*EXP(1E0_dp + aim/phy%zC), 0, 0))/aim
phy%SI = physi(phy, phin, phy%PARs, PARb, depfrac)
phy%theta = phy%tch*phy%fC
RC = RC*depfrac/phin ! for calculation of phy%A in SR grow
......
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