def calc_mean_func(pct_claym, pct_sandm):
s_sat, theta, epsi, taug, tauw, film, psi, chb = calc_Kaff_soil(
pct_claym, pct_sandm)
factw = s_sat**(1. / df)
Kaff_o2g_full, Kaff_o2g, k2, k1_o2, k1_o2_full, kappa_tops = calc_Kaff_O2(
s_sat, theta, epsi, taug, tauw, film, DZ, Ncello, BT, alphaV, factw)
Kaff_o2g_full = Kaff_o2g_full * o2scal
k1_o2_full = k1_o2_full / o2scal
Ncell = Ncello
k1_s, Kaff_s, Kaff_s_0 = calc_Kaff_SC(s_sat, theta, epsi, taug, tauw, film,
DZ, Ncell, alphaV)
#damm model
k2 = k2 * k2f
Vmax = BT * k2
dammf = damm_flx(O2, S, factw, Kaff_s * 1.e0, Kaff_o2g, kappa_tops, Vmax)
damfmax = np.max(dammf)
dammf = dammf / damfmax
#su model
suf = su_flx(O2, S, factw, BT, k2 * 1.e0, k1_o2, k1_s, kappa_tops)
sufmax = np.max(suf)
suf = suf / sufmax
#supeca model
ms = 0.0
Km_s = 2.0
supf = supmic_flx(O2, S, factw, BT, ms, Km_s, Kaff_o2g, Kaff_s, kappa_tops,
k2)
supfmax = np.max(supf)
supf = supf / supfmax
return dammf, suf, supf
pct_clay = 19.
kt = len(S)
k = 0
k = 0
ms = [0.0, 50., 100., 200.]
supfmax0 = 0.
while k < kt:
s_sat, theta, epsi, taug, tauw, film, psi, chb = calc_Kaff_soil(
pct_clay, pct_sand)
factw = s_sat**(1. / df)
Kaff_o2g_full, Kaff_o2g, k2, k1_o2, k1_o2_full, kappa_tops = calc_Kaff_O2(
s_sat, theta, epsi, taug, tauw, film, DZ, Ncello[0], BT, alphaV, factw)
Ncell = 10
k1_s, Kaff_s, Kaff_s_0 = calc_Kaff_SC(s_sat, theta, epsi, taug, tauw, film,
DZ, Ncell, alphaV)
#supeca model
k2 = k2 * k2f
Km_s = 2.0
supf = supmic_flx(O2, S[2], factw, BT, ms[k], Km_s, Kaff_o2g, Kaff_s,
kappa_tops, k2)
supfmax = np.max(supf)
if k == 0:
supfmax0 = supfmax
supf = supf / supfmax0
supflsMs.append(supf)
Kaff_o2g = Kaff_o2g * 10.
supf = supmic_flx(O2, S[2], factw, BT, ms[k], Km_s, Kaff_o2g, Kaff_s,
kappa_tops, k2)
supfmax = np.max(supf)
supf = supf / supfmax0