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
kt=len(pct_clay) k=0 while k < kt: s_sat, theta, epsi, taug, tauw, film,psi,chb=calc_Kaff_soil(pct_clay[k], 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, BT, alphaV,factw) Kaff_o2g_full=Kaff_o2g_full*o2scal k1_o2_full=k1_o2_full/o2scal Ncell=10 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 dammfls.append(dammf) sufls.append(suf) supfls.append(supf)