From ef0b871077b6d84c119fbd1178356b79283044ee Mon Sep 17 00:00:00 2001 From: Liang-Jun Zhu Date: Wed, 29 Nov 2023 15:07:04 +0800 Subject: [PATCH] Merge src/carbon_zhang2.f90 of master branch --- src/carbon_zhang2.f90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/carbon_zhang2.f90 b/src/carbon_zhang2.f90 index 4baedce..dd8f2dd 100644 --- a/src/carbon_zhang2.f90 +++ b/src/carbon_zhang2.f90 @@ -326,7 +326,8 @@ subroutine carbon_zhang2 OX = 0. !OX = 1 - (0.9* sol_z(k,j)/1000.) / (sol_z(k,j)/1000.+ exp(1.50-3.99*sol_z(k,j)/1000.)) !OX = 1 - (0.8* sol_z(k,j)) / (sol_z(k,j)+ exp(1.50-3.99*sol_z(k,j))) - OX=1.-0.8*((sol_z(kk,j)+sol_z(kk-1,j))/2)/(((sol_z(kk,j)+sol_z(kk-1,j))/2)+EXP(18.40961-0.023683632*((sol_z(kk,j)+sol_z(kk-1,j))/2))) + OX=1.-0.8*((sol_z(kk,j)+sol_z(kk-1,j))/2)/ & + & (((sol_z(kk,j)+sol_z(kk-1,j))/2)+EXP(18.40961-0.023683632*((sol_z(kk,j)+sol_z(kk-1,j))/2))) !! compute combined factor CS = 0. CS=MIN(10.,SQRT(cdg*sut)*0.9*OX*X1)