ea = rh * ea_str #vapor pressure (umol h2O/mol air)
#correct for leaf temperatures using leaf height
t_diff = 18 - 0.4 * ht
tl = surtemp_dm[xxx][time] + t_diff
z = 0.2
#---------------Photosynthesis Function---------------#
#alter this line of code for when implementing different photosynthesis functions
wue, nue, A, E, cs, ci, gsw, gs, gbw, gb, gm, cc, dd = photo(
tk_25, ekc, eko, etau, ev, ej, toptv, toptj,
na_dm[xxx][time], qeff, PAR, tl, ea, chl, ij, kc25,
ko25, o, ca, rh, m, a, frnr, flnr, ra, jm, g0, b,
dia, u, q, vwc_min, vwc_max, vwc_dm[xxx][time], z)
#test to make sure wue and nue are positive at not 'nan'
if wue[0] == -999 and nue[0] == -999:
continue
wue_d += [wue[0]]
nue_d += [nue[0]]
gsw_d += [gsw[0]]
A_d += [A[0]]
E_d += [E[0]]
vpd_d += [dd[0]]
if xxx == 0:
pa_v, 3528
) #saturation vapor pressure of air (Pa-->umol h20/mol air)
ea = rh * ea_str #vapor pressure (umol h2O/mol air)
#correct for leaf temperatures using leaf height
t_diff = 18 - 0.4 * ht[0][iii]
tl = t_dm + t_diff
#---------------Photosynthesis Function---------------#
#alter this line of code for when implementing different photosynthesis functions
wue, nue, A, E, cs, ci, gsw, gs, gbw, gb, gm, cc, dd = photo(
tk_25, ekc, eko, etau, ev, ej, toptv, toptj, na[0][iii], qeff,
PAR, tl, ea, chl[0][iii], ij, kc25, ko25, o, ca, rh, m, a,
frnr, flnr, ra, jm, g0, b, dia[0][iii], u, q, vwc_min, vwc_max,
vwc_dm)
#test to make sure wue and nue are positive at not 'nan'
if wue[0] == -999 and nue[0] == -999:
continue
wue_d += [wue[0]]
nue_d += [nue[0]]
gsw_d += [gsw[0]]
A_d += [A[0]]
E_d += [E[0]]
vpd_d += [dd[0]]
if 'nm' in dict_params[ii]:
nm = nm_c
s = s_c #nm and s are dependent variables
if 'm' in dict_params[ii]:
m = m_c
if 'chl' in dict_params[ii]:
chl = chl_c
if 'tl' in dict_params[ii]:
tl = tl_c
if 'vwc' in dict_params[ii]:
vwc = vwc_c
if 'ij' in dict_params[ii]:
ij = ij_c
wue, nue, A, gsw, E, na = photo(s, nm, tl, ea, chl, crc, rub_max,
ij, vwc, kc25, ko25, o, tau25, ca,
rh, m, a, frnr, flnr, ra, j_b,
j_m_max, q, vwc_min, vwc_max, b)
if wue == -999 and nue == -999:
continue
#---------------Test for Low NUE Values---------------#
# if any(nue<15):
# break
#---------------Make Array of Values for Each Meadow---------------#
if i + 1 < 33:
if iii == 0:
wue_f_per1 += [wue]