variables.list (for E6.save.f) rhino ~/bosque 7 July 01
Listed in a modified collating order
numbers before letters
case does not matter except that, when two letters of different
case appear, capital letter goes first
Thus, we get
I0, I00, IL, Ilavg,....
This is not the standard collating order (numbers, then caps., then
lower case), but it makes searching easier and helps catch errors
in which two variables meant to be different only differ in
case, which usually makes no distinction in Fortran
In MAIN
Variables that are not arrays (single variables)
character*2
ampm (-) = AM vs. PM designator in weather records
character*3
test (-) = 1st 3 bytes of keyboard input to check for keyword
character*8
date0 (-) = date for beginning of simulation (mm-dd-yy format)
date0wea (-) = date at which weather record begins (mm-dd-yy)
datef (-) = date for end of simulation (mm-dd-yy format)
character*17
sval (-) = new numerical value of any parameter being changed
character*20
chaval (-) = keyword and numerical value for changing a parameter
character*24
lpenfile (-) = name of file that has previously calculated
array of light-penetration probabilities
newlpenfile = similar, if this is first time they're calculated
weafile = name of weather data file
integer
iac (-) = accumulator of number of days in preceding months
iconv (-) = number of iteration of eair, Tair in which convergence
occurred
iday (d) = for a weather record, the number of day within month
iflag (-) = (not used)
ihr (h) = for a weather record, the hour of the day
ihrcorr (h) = ihr, corrected for 12 h addition in afternoon
index1 (-) = weather record to retrieve when interpolating data
indexx (-) = variable used in debugging, to pick out one sky angle
ipensave (-) = switch: 1 = save newly calculated array of
light-penetration probabilities, 0 = don't
ipre (-) = switch: 1 = light-penetration probability arrays
exist in a file; use them; 0 = they don't exist
iRHpct (-) = air relative humidity, in %, in weather record
iter (-) = iteration in solution for eair, Tair in canopy
itoffset (min) = minutes offset of simulation times from an
even multiple of 10 min. (so we can match sapflow records)
jd (d) = Julian day of current weather record
jd0 (d) = Julian day at which to begin simulation
jd0wea (d) = Julian day at which the weather-date file begins
jdf (d) = Julian day at which to end simulation
minute0 (min) = minute after midnight at which to begin each day's
simulation
minutef (min) = minute after midnight at which to end each day's
simulation
month (mo) = for a weather record, the month of the year
nbins (-) = number of bins to use in histogram of leaf
PAR irradiance, IL
ncha (-) = ID number of the parameter to be changed
nd (-) = number of steps of size 'stepsize' to take through canopy
while accumulating total length of path within tree crowns
ndmax (-) = max. number of steps (prevents one from wasting time
looking outside all known trees)
nphi (-) = number of azimuth angles to use in sampling locations
within central tree
nradius-) = number of radii from center to use in sampling locations
within central tree
ntheta-) = number of zenith angles to use in sampling locations
within central tree
ntimes (-) = number of the current simulation interval
ntimetot (-) = final value of ntimes at end of simulation
ntrees (-) = number of trees around the central tree that we will
describe
logical
BBswitch (-) = switch: 1 = use Ball-Berry humidity response,
0 = do not
real
A (mol m-2 [10 min]-1) = CO2 assimilation rate of leaf under
consideration
AAtree (mol s-1) = accumulator for CO2 assim. rate of entire tree
AAtree10 (mol s-1) = same, when gs is increased by 10%)
aax (m) = a distance scale for sampling inside central tree's
ellipsoid
Agrand (mol) = accumulator for CO2 assimlation of central tree
integrated over the entire simualation time
Agrand10 (mol) = same, when gs has been increased 10% at each time
Aleaf (m2) = total leaf area on central tree
alph (m) = a distance scale for sampling inside central tree
aNIR (-) = absorptivity of leaf for NIR
aPAR (-) = absorptivity of leaf for PAR
arg (-) = temporary variable in several locations
avgI (mol m-2 s-1) = average PAR irradiance over all loc. in tree
bax (m) = a distance scale for sampling inside central tree's
ellipsoid
bBB (mol m-2 s-1) = intercept in Ball-Berry model of gs
bet (m) = intermediate in calc. of arbitary tree's projected
canopy length along sun direction
BTIR (W m-2K-1) = derivative w/r to T of QTIR flux density from leaf
c (J-1 K m2 s) = factor in computing numerical derivatives of
stand latent and sens. heat fluxes w/r to eair and Tair
Ca (Pa) = partial pressure of CO2 in ambient air (at ref. height
above canopy)
CAlocal (-) = control coefficient for gs over A of tree, for the
current simulation interval
CElocal (-) = control coefficient for gs over E of tree, for the
current simulation interval
Ci (Pa) = partial pressure of CO2 inside leaf
cosalphas (-) = cosine of solar elevation angle
costhetaD (-) = cosine of zenith angle of sky direction being
sampled while computing diffuse skylight penetration
costhetape (-) = dot product of direction vectors, in computation
of distance that direct solar beam travels inside tree canopy
Cpair (J mol-1 K-1) = heat capacity of air
Cs (Pa) = partial pressure of CO2 at surface of leaf, beneath
leaf boundary layer
D00 (micromol m-2 s-1) = radiant flux density (not radiance) in
diffuse skylight at top of canopy
D000 (micromol m-2 s-1) = initial estimate of D00 (which see)
D0PAR (micromol m-2 s-1) = same as D00, but interpolated to
this specific time from weather data
Daa (Pa) = air-to-air vapor pressure deficit above canopy
deair (Pa) = increment in eair within canopy for this iteration
denom (Pa K) = intermediate in calculation of numerical derivatives
for eair, Tair iterations 4 and 5
derrordeair (-) = derivative of error in eair of canopy with
respect to increment in eair
dleaf (m) = characteristic linear dimension of leaf
dpc (m) = distance of a point from center of selected tree
dpen (-) = penetration probability of direct beam to canopy locaiton
drow (m) = (mean) distance between tree rows
Ds (radians) = solar declination N of equator
dt (m) = distance along path from chosen point in tree to top of
canopy
dTair (oC or K) = increment in canopy air temperature during
iterative solution
dTerrordTair (-) = derivative of error in Tair of canopy w/r to
increment in Tair
dtime (h) = cumulative (decimal) time from beginning of year
dtree (m) = (mean) distance between columns of trees
dummy (-) = dummy argument to function Ecalc (copious arguments
are passed in common blocks, instead)
dx (m) = x-component of step along direct beam path while checking
for interception by trees
dy (m) = same, for y-component
dz (m) = same, for z-component
E (mol m-2 s-1) = transpiration rate of chosen leaf
E1, E2, E3, E4, E5 (mol s-1) = transpiration rate of central tree,
for respective iterations in eair and Tair
E1stand, E2stand, E3stand, E4stand, E5stand (mol m-2 s-1) = same,
per groudn area of stand, taking central tree as avg.tree
eair (Pa) = partial pressure of w.v. inside canopy
eair0 (Pa) = partial pressure of w.v. above canopy, at ref. height
eair1,eair2,eair3,eair4,eair5 (Pa) = eair in canopy, at successive
iterations for solution
eairtol (Pa) = limit of acceptable error in solving for eair
Egrand10Lh (l h-1) = transpiration rate of central tree, integrated
over entire simulation time, after gs is increased 10% at each time
EgrandLh (l h-1) = same, at correct gs
ELaavg (mol m-2 s-1) = average transpiration rate per unit leaf
area, in whole canopy
ELaavg10 (mol m-2 s-1) = same, after gs is increased 10%
enlocations (-) = number of discrete locations sampled in central tree
epsleaf - not used yet
Epe (mol s-1 Pa-1) = derivative of tree transpiration rate w/r to
increment in eair
EpT (mol s-1 K-1) = derivative of tree transpiration rate w/r to
increment in Tair
error (Pa) = error in consistency of eair with transport eqn.
error1, error2, error3, error4, error5 (Pa) = error in eair,
at successive iterations of solution
Estandlastavg (mol m-2 s-1) = stand tranpsiration rate, averaged
over previous two simulation intervals
Estandm1 (mol m-2 s-1) = stand transpir. rate, at immediately
preceding time interval
Estandm2 (mol m-2 s-1) = same, two intervals ago
Etree (mol s-1) = transpiration rate of central tree
Etree10 (mol s-1) = same, with gs increased by 10% on all leaves
Etree10Lh (l h-1) = same as Etree10, converted to l h-1
EtreeLh (l h-1) = transpiration rate of central tree, in l h-1
Fcrown (-) = fraction of total solar flux intercepted by central
tree over its ground area, integrated over whole simulation time
fd (m-1) = foliage density inside crowns of trees
fraci (-) = intermediate in constructing zenith angles to sample
within central tree
gammaa - not used now
gbair (mol m-2 s-1) = conductance of leaf boundary layer, molar units
gbcan (mol m-2 s-1) = same, for canopy b. l.
gcan (mol m-2 s-1) = total conductance for w.v. of all leaves in canopy
gtotCO2 (mol m-2 s-1) = conductance of stomata plus leaf b. l. for CO2
gtottree (mol m-2 s-1) =accumulator for conductance of all leaves on
tree
gtotw (mol m-2 s-1) = conductance of stomata plus leaf b. l., for w.v.
h (radians) = hour angle for solar azimuth
H1stand, H2stand, H3stand, H4stand, H5stand (mol m-2 s-1) = sensible
heat flux of stand per ground area, for respective iterations of
solution for eair and Tair in canopy
ha (-) = notional (default) average relative humidity at leaf surface
hcan (m) = max. height of canopy (among all trees)
Hpe (W m-2 Pa-1) = derivative of stand sens. heat flux w/r to increment
in eair
HpT (W m-2 K-1) = derivative of stand sens. heat flux w/r to increment
in Tair
hs0 (-) = mean rel. humidity at leaf surface; neg. value means
model should use Ball-Berry model
Hstandlastavg (mol m-2 s-1) = sens. heat flux of stand per ground
area, averaged over previous two time intervals
Hstandm1 (mol m-2 s-1) = sens. heat flux of stand at previous time
interval
Hstandm2 (mol m-2 s-1) = sens. heat flux of stand, two time
intervals ago
Htree (mol s-1) = transpir. rate of central tree
htree (m) - height of individual tree, when finding max. in canopy,
hcan - has same name as transpir. rate but use is finished
before Htree is calc'd
Hvap (J mol-1) = heat of vap. of water
I00 (micromol m-2 s-1) = radiant flux density in PAR of direct
solar beam, normal to its direc. of propag.
I00start (micromol m-2 s-1) = initial est. for I00
I0PAR (micromol m-2 s-1) = I00, after interpolation from weather data
Iground (micromol m-2) = sum over all simul. intervals of
total PAR flux density incident per ground area of stand
ILavg (micromol m-2 s-1) = average PAR irradiance on all leaves in
central tree
Isum (micromol m-2 s-1) = sum of PAR irradiance over all canopy
locations
Kc25 (Pa) = Michaelis constant for binding of CO2 to Rubisco, at ref. T
of 25oC
Ko25 (Pa) = same, for binding of O2
kpen (-) = average cosine of angle between leaves and sun
kTIR - not used now
lat (deg.) = latitude of site, N of equator
mBB (-) = slope in Ball-Berry model
Mcrown (micromol h-1) = total light intercepted by all leaves in
central tree
Mground (micromol h-1) = total light intercepted by flat ground
overrad (radians-1) = conversion factor from deg. to radians
p (Pa mol-1 m2 s) = factor in calculating numerical derivatives
in solution for eair, Tair
Pair (Pa) = total pressure of air
Pclear (-) = fraction of time in simulation interval that sky is
clear (estimated)
phican (radians) = currently selected azimuth within central tree
phiD (radians) = currently selected azimuth within sky for calc. of
diffuse skylight penetration
phis (radians) = solar azmiuth
pi (-) = math. constant, 3.141592...
PPFD (mol m-2 s-1) = PAR irradiance on leaf
Q00 (mol CO2 [mol photons]-1) = initial quantum yield of CO2
assimilation, at saturating CO2
Qrad (W m-2) = total radiative energy gain per area of leaf
Qsky (W m-2) = TIR energy flux density from sky
Qsol (MJ m-2 h-1) = solar energy flux density
Qveg (W m-2) = TIR energy flux density from surface at T of vegetation
ra (s m-1) - not used now (leaf b. l. resistance)
radmax (m) = max. radius from location in central tree to which to
search for other tree crowns
RdovrA (-) = dark respir.of leaf as fraction of saturated A at
Tleaf = recent mean T for last 2 weeks
rds (m) = radius of canopy along solar direction
RdTmean (mol m-2 s-1) = leaf respir. rate at Tleaf = Tmean
remain (-) = remainder in converting decimal time to closest
interval in weather data
rhoair (mol m-3) = molar density of air
rhoovrC1 (mol m-3) = proportionality between gbcan and windspeed
sigma (W m-2 K-4) = Stefan-Boltzmann constant
sinalphas (-) = sine of solar elevation angle
sinthetaD (-) = sine of solar zenith angle currently being sampled
in calc. of diffuse skylight penetration
st (m) = cumulative path of sun through crowns of all trees in stand
along current solar direction
stepsize (m) = increment in path search along solar beam direction
stmax (m) = max. value of st, at which search quits (penetration
probability is very low)
sumI (mol m-2 s-1) = sum of PAR irradiance over all leaves over all
simulation intervals
sumpen (-) = sum of fractions of diffuse skylight penetrating to
canopy location
Tair (oC) = air temperature inside canopy
Tair0 (oC) = air temperature above canopy, at ref. height
Tair1,Tair2,Tair3,Tair4,Tair5 = successive iterated solutions for Tair
Tairtol (oC or K) = tolerance for error in Tair
Terror (oC or K) = final error in Tair at last iteration needed
Terror1, Terror2, Terror3, Terror4, Terror5 = error in Tair at
respective iteration
thetaPS (-) = convexity parameter, for sharpness of transition from
light-limited to light-sat'd CO2 assimilation
TLavg (oC) = leaf temperature, averaged over all leaves in central tree
Tleafval (oC) = leaf temperature, renamed so it can be passed in common
TLtree (oC) = sum of (leaf T) x (leaf conductance) over all
locations in tree
Tmean (oC) = mean leaf T during photoperiod over last 2 weeks
(T to which leaf respiration has acclimated)
ToF (oF) = Fahrenheit air T in current weather record
tol (mol m-2 s-1) = tolerance for error in solving for gs
u (m s-1) = windspeed
umph (mi h-1) = windspeed, in English units
Vcmax25 (micromol m-2 s-1, immed. converted to mol m-2 s-1) = photosyn.
capacity per leaf area (light- and CO2-saturated), at 25oC
VPDavg (Pa) = average vapor-pressure deficit from leaf to air, over
whole canopy
wt (-) = fraction of leaf area at this location in shade only
wt1 (-)= weighting function for previous hourly weather record, in
interpolation of weather variables
wt2 (-) = same, for next hourly record
x (radians) = sine of solar azimuth; also (m), x-coordinate of
current point in central tree
xaux (m) = x-coordinate of center of central tree
xnew (m) = x-coordinate of current position in whole canopy during
search for interception by crowns
xo (m) = x-coordinate of currently-sampled point in central tree
yaux (m) = y-coordinate of center of central tree
ynew (m) = y-coordinate of current position in whole canopy during
search for interception by crowns
yo (m) = y-coordinate of currently-sampled point in central tree
zaux (m) = z-coordinate of center of central tree
znew (m) = z-coordinate of current position in whole canopy during
search for interception by crowns
zo (m) = z-coordinate of currently-sampled point in central tree
Arrays
character*3
ID(12) (-) = 3-letter keyword to recognize parameter that user wants
to change
integer
itime(15000) (min) = time of simulation interval, measured from
midnight on first day of weather record
ndm(12) (d) = number of days in given month
real
ae(3) (Pa) = absolute value of error in solving for eair in canopy
aeT(3) (oC) = absolute value of error in solving for Tair in canopy
alphas(10000) (radians) = solar elevation angle for each simulation
interval
alphaswea(10000) (radians) = solar elevation angle for each weather
record
atree(100) (m) = horizontal full axis of ellipsoid for given tree
btree(100) (m) = vertical full axis of ellipsoid for given tree
costheta(20) (-) = cosine of zenith angle of point being sampled in
central tree
D0PARwea(10000) (micromol m-2 s-1) = diffuse skylight flux density in
PAR, estimated for given weather record
Daawea(10000) (Pa) = air-to-air vapor pressure deficit for given
weather record
eairwea(10000)
I0PARwea(10000) (micromol m-2 s-1) = direct-beam flux density in
PAR, estimated for given weather record
lpenstore(20,20,20) (-) = penetration probability for direct beam
to all canopy locations, for current simulation interval
PARwea(10000) (micromol m-2 s-1) = total PAR irradiance on flat
ground area, for each weather record
Pclearwea(10000) (-) = estimated frac. of time that sky is clear
(direct beam is present), for given weather record
phi(20) (radians) = sky azimuths being sampled for diffuse skylight
phic(100) (radians) = azimuthal inclination of ellipsoid of each tree
phiswea(10000) (radians) = solar azimuth for each weather record
Ppendiffuse(20,20,20) (-) = total penetration fraction of diffuse
skylight, for each location in central tree
r(20,20) (m) = array of radii being sampled in central tree, for
each choice of theta and phi
sintheta(20) (-) = sine of zenith angle of point being sampled in
central tree
Tairwea(10000) (oC) = air temperature above canopy, for each weather
record
thetac(100) (radians) = zenith inclination of ellipsoid of each tree
timerel(10000) (h) = time of current simul. interval, after start
of simulation
uwea(10000) (m s-1) = windspeed for each weather record
xtree(100) (m) = x-coordinate of center of each tree
ytree(100) (m) = y-coordinate of center of each tree
ztree(100) (m) = z-coordinate of center of each tree
do-loop indices
i
i2
iline (-) = number of an individual record in weather-data file
iphi
iradius
itheta
itree
j
nt
ntimes
subroutines
BBsolve
CiACssolve
control
Ecalc
esat
esatpr
julian
Tleaf (function)
In BBsolve
Variables that are not arrays (single variables)
integer
nruns (-) = number of binary searches to do to narrow down the
solution for gs
logical
BBswitch (-) = switch: 1 = use Ball-Berry humidity response,
0 = don't
real
A (mol m-2 s-1) = CO2 assim. rate of leaf
Cs (Pa) = partial pressure of CO2 at leaf surface, beneath
leaf boundary layer
E (mol m-2 s-1) = transpiration rate of leaf
eair (Pa) = partial pressure of water vapor in ambient air
around leaf
eL (Pa) = saturated partial pressure of w.v., at leaf temperature
gbair (mol m-2 s-1) = conductance of leaf boundary layer, molar units
gs0 (mol m-2 s-1) = initial estimate of gs
gshi (mol m-2 s-1) = current upper limit in search for gs
gslo (mol m-2 s-1) = current lower limit in search for gs
gsmid (mol m-2 s-1) = midpoint of search for gs
gsmin (mol m-2 s-1) = minimum value of gs allowed physically
gsminus (mol m-2 s-1) = current lower limit in search for gs,
during binary search
gsplus (mol m-2 s-1) = current upper limit in search for gs,
during binary search
gtotCO2 (mol m-2 s-1) = conductance of stomata plus leaf b. l.
for CO2
gtotw (mol m-2 s-1) = conductance of stomata plus leaf b. l.,
for w.v.
ha (-) = relative humidity of air
hs (-) = relative humidity at leaf surface, beneath leaf
boundary layer
hs0 (-) = notional (average) value of hs
IBBhi (mol m-2 s-1) = Ball-Berry index for gshi
IBBlo (mol m-2 s-1) = Ball-Berry index for gslo
IBBmid (mol m-2 s-1) = Ball-Berry index for gsmid
mBB (-) = slope in Ball-Berry model
Pair (Pa) = total pressure of air
PPFD (mol m-2 s-1) = PAR irradiance on leaf
Qrad (W m-2) = total radiative energy gain per area of leaf
Tair (oC) = air temperature inside canopy
testhi (mol m-2 s-1) = error in solution for gs at gshi
testlo (mol m-2 s-1) = error in solution for gs at gslo
testmid (mol m-2 s-1) = error in solution for gs at gsmid
TL (oC) = leaf temperature
tol (mol m-2 s-1) = tolerance for error in solving for gs
do-loop indices
n
Arrays
(NONE)
In Tleaf
Variables that are not arrays (single variables)
real
A (mol m-2 s-1) = CO2 assim. rate of leaf
Bcc (W m-2 K-1) = heat transfer coefficient for leaf
Bnet (W m-2 K-1) = derivative of total energy balance w/r to Tleaf
BTIR (W m-2K-1) = derivative w/r to T of QTIR flux density from leaf
Cs (Pa) = partial pressure of CO2 at surface of leaf, beneath
leaf boundary layer
dTleaf (oC or K) = increment in Tleaf
E (mol m-2 s-1) = transpiration rate of chosen leaf
eair (Pa) = partial pressure of water vapor in ambient air
around leaf
eL (Pa) = saturated partial pressure of w.v., at leaf temperature
ep (Pa K-1) = derivative of eL w/r to Tleaf
gbair (mol m-2 s-1) = conductance of leaf boundary layer, molar units
gtotCO2 (mol m-2 s-1) = conductance of stomata plus leaf b. l.
for CO2
gtotw (mol m-2 s-1) = conductance of stomata plus leaf b. l.,
for w.v.
Hvap (J mol-1) = heat of vap. of water
Pair (Pa) = total pressure of air
PPFD (mol m-2 s-1) = PAR irradiance on leaf
Qmcc (W m-2) = rate of sensible heat loss from leaf
QmE (W m-2) = rate of latent heat loss from leaf
Qmrad (W m-2) = rate of TIR radiant heat loss from leaf
Qnet (W m-2) = net energy flux imbalance of leaf
Qrad (W m-2) = total radiative energy gain per area of leaf
sigma (W m-2 K-4) = Stefan-Boltzmann constant
Tair (oC) = air temperature inside canopy
TL (oC) = leaf temperature
Tleaf (oC) = leaf temperature
Arrays
(NONE)
In CiACssolve
Variables that are not arrays (single variables)
integer
iflag - not used currently
iter (-) = iteration in solution for eair, Tair in canopy
real
A (mol m-2 [10 min]-1) = CO2 assimilation rate of leaf under
consideration
ALL (mol m-2 s-1) = light-limited rate of CO2 assimilation
Amax (mol m-2 s-1) = light-saturated rate of CO2 assimilation
c0,c1,c2,c3,c4 (monster units) = coefficients in quartic solution
for Ci
Ca (Pa) = partial pressure of CO2 in ambient air (at ref. height
above canopy)
Ci (Pa) = partial pressure of CO2 inside leaf
Cs (Pa) = partial pressure of CO2 at surface of leaf, beneath
leaf boundary layer
dCi (Pa) = increment in solution for Ci
denom (J mol-1 K) = intermediate in calc. of temperature activation
functions for CO2 assimilation
dfac (K2) = similar
E (mol m-2 s-1) = transpiration rate of chosen leaf
eair (Pa) = partial pressure of w.v. inside canopy
F (monster units) = function whose root is sought to solve for Ci
G (Pa) = shorthand for gamma
gamma (Pa) = CO2 compensation partial pressure
gbair (mol m-2 s-1) = conductance of leaf boundary layer, molar units
gt (mol m-2 s-1) = shorthand for gtotCO2
gt2 (mol2 m-4 s-2) = shorthand for square of gtotCO2
gtotCO2 (mol m-2 s-1) = conductance of stomata plus leaf b. l.
for CO2
gtotw (mol m-2 s-1) = conductance of stomata plus leaf b. l.,
for w.v.
IL (mol m-2 s-1) = PAR irradiance on leaf
K (Pa) = shorthand for Kco
Kc (Pa) = Michaelis constant for binding of CO2 to Rubisco, at
current leaf temperature
Kc25 (Pa) = Michaelis constant for binding of CO2 to Rubisco,
at ref. T of 25oC
Kco (Pa) = effective Michaelis constant for binding of CO2 to
Rubisco, in presence of current O2 partial pressure
Ko (Pa) = Michaelis constant for binding of O2 to Rubisco, at
Ko25 (Pa) = Michaelis constant for binding of O2 to Rubisco,
at ref. T of 25oC
P (Pa) = shorthand for Pair
P2 (Pa2) = shorthand for square of Pair
Pair (Pa) = total pressure of air
PPFD (mol m-2 s-1) = PAR irradiance on leaf
Q00 (mol CO2 [mol photons]-1) = initial quantum yield of CO2
assimilation, at saturating CO2
QI (mol CO2 s-1) = shorthand for Q00*IL
Qrad (W m-2) = total radiative energy gain per area of leaf
Rd (mol m-2 s-1) = dark respiration rate of leaf
RdTmean (mol m-2 s-1) = leaf respir. rate at Tleaf = Tmean
T (-) = shorthand for thetaPS
Tair (oC) = air temperature inside canopy
thetaPS (-) = convexity parameter, for sharpness of transition from
light-limited to light-sat'd CO2 assimilation
TL (oC) = leaf temperature
Tmean (oC) = mean leaf T during photoperiod over last 2 weeks
(T to which leaf respiration has acclimated)
V (mol m-2 s-1) = shorthand for Vcmax
Vcmax (micromol m-2 s-1, immed. converted to mol m-2 s-1) = photosyn.
capacity per leaf area (light- and CO2-saturated)
Vcmax25 (micromol m-2 s-1, immed. converted to mol m-2 s-1) =
photosyn. capacity per leaf area (light- and CO2-saturated),
at 25oC
Arrays
(NONE)
In ESAT
real
ESAT
T
In ESATPR
real
ESATPR
T
In F
real
c0,c1,c2,c3,c4
Ci (Pa) = partial pressure of CO2 inside leaf
In julian
character*8
date (-) = current date, in mm-dd-yy format
integer
iac (-) = accumulator of number of days in preceding months
nday (d) = number of day within current month
nday1 (10 d) = 1st digit in nday
nday2 (10 d) = 2nd digit in nday
ndm(12) (d) = number of days in each month
nmonth (months) = number of current month in year
nmonth1 (10 months) = 1st digit of nmonth
nmonth2 (10 months) = 2nd digit of nmonth
do-loop indices
i
In Ecalc
Variables that are not arrays (single variables)
integer
ipensave (-) = switch: 1 = save newly calculated array of
light-penetration probabilities, 0 = don't
ipre (-) = switch: 1 = light-penetration probability arrays
exist in a file; use them; 0 = they don't exist
iter (-) = iteration in solution for eair, Tair in canopy
itest - not used currently (for debugging)
nbins (-) = number of bins to use in histogram of leaf
PAR irradiance, IL
nd (-) = number of steps of size 'stepsize' to take through canopy
while accumulating total length of path within tree crowns
ndmax (-) = max. number of steps (prevents one from wasting time
looking outside all known trees)
nphi (-) = number of azimuth angles to use in sampling locations
within central tree
nradius-) = number of radii from center to use in sampling locations
within central tree
ntheta-) = number of zenith angles to use in sampling locations
within central tree
ntimes (-) = number of the current simulation interval
ntrees (-) = number of trees around the central tree that we will
describe
real
A (mol m-2 [10 min]-1) = CO2 assimilation rate of leaf under
consideration
AAtree (mol s-1) = accumulator for CO2 assim. rate of entire tree
AAtree10 (mol s-1) = same, when gs is increased by 10%)
Aleaf (m2) = total leaf area on central tree
aNIR (-) = absorptivity of leaf for NIR
aPAR (-) = absorptivity of leaf for PAR
Asum (mol m-2 s-1) = sum of leaf assim. rate over different
bins of PAR irradiance
Asum10 (mol m-2 s-1) = same, for gs increased by 10%
Bcc (W m-2 K-1) = heat transfer coefficient for leaf
bet (m) = intermediate in calc. of arbitary tree's projected
canopy length along sun direction
cosalphas (-) = cosine of solar elevation angle
costhetape (-) = dot product of direction vectors, in computation
of distance that direct solar beam travels inside tree canopy
Cs (Pa) = partial pressure of CO2 at surface of leaf, beneath
leaf boundary layer
D0PAR (micromol m-2 s-1) = same as D00, but interpolated to
this specific time from weather data
dI (micromol m-2 s-1) = increment in PAR irradiance, between bins
dpc (m) = distance of a point from center of selected tree
dt (m) = distance along path from chosen point in tree to top of
canopy
dx (m) = x-component of step along direct beam path while checking
for interception by trees
dy (m) = same, for y-component
dz (m) = same, for z-component
E (mol m-2 s-1) = transpiration rate of chosen leaf
eair (Pa) = partial pressure of w.v. inside canopy
Esum (mol m-2 s-1) = sum of leaf transpir. rate over different
bins of PAR irradiance
Esum10 (mol m-2 s-1) = same, for gs increased by 10%
Etree (mol s-1) = transpiration rate of central tree
Etree10 (mol s-1) = same, with gs increased by 10% on all leaves
fbin (-) = fraction of total irrad. range per bin
fd (m-1) = foliage density inside crowns of trees
gbair (mol m-2 s-1) = conductance of leaf boundary layer, molar units
gsminus (mol m-2 s-1) = current lower limit in search for gs,
during binary search
gsplus (mol m-2 s-1) = current upper limit in search for gs,
during binary search
gstest (mol m-2 s-1) = test value of gs, 10% higher than true sol.
gtot (mol m-2 s-1) = conductance of stomata plus leaf b. l.,
for w.v.
gtotCO2 (mol m-2 s-1) = conductance of stomata plus leaf b. l.
for CO2
gtotsum (mol m-2 s-1) = accumulator for total leaf conductance
for this location in canopy
gtottree (mol m-2 s-1) =accumulator for conductance of all leaves on
tree
gtotw (mol m-2 s-1) = conductance of stomata plus leaf b. l.,
for w.v. - only passed to keep common block intact
hcan (m) = max. height of canopy (among all trees)
hs (-) = relative humidity at leaf surface, beneath leaf
boundary layer
Hsum (W m-2) = sensible heat flux from leaf at this location
Htree (mol s-1) = transpir. rate of central tree
I0PAR (micromol m-2 s-1) = I00, after interpolation from weather data
Iadd (micromol m-2 s-1) = average PAR irradiance for leaves at this
location in canopy
IL (micromol m-2 s-1) = PAR irradiance on leaf, in particular bin
Isum (micromol m-2 s-1) = PAR irradiance sum for leaf at this
location
kpen (-) = average cosine of angle between leaves and sun
lpen (-) = probability of direct-beam penetration to this
location in canopy at this simulation interval
Pair (Pa) = total pressure of air
Pclear (-) = fraction of time in simulation interval that sky is
clear (estimated)
phican (radians) = currently selected azimuth within central tree
PPFD (mol m-2 s-1) = PAR irradiance on leaf
Qmcc (W m-2) = rate of sensible heat loss from leaf
Qrad (W m-2) = total radiative energy gain per area of leaf
Qsky (W m-2) = TIR energy flux density from sky
QSW (W m-2) = total shortwave (PAR + NIR) energy flux density
absorbed by this leaf
QTIR (W m-2) = total TIR energy flux density absorbed by this leaf
Qveg (W m-2) = TIR energy flux density from surface at T of
vegetation
ra - not used anymore
radius (m) = radial distance in central tree, at which leaf
performance is being sampled
rds (m) = radius of canopy along solar direction
rhoair (mol m-3) = molar density of air
sigma (W m-2 K-4) = Stefan-Boltzmann constant
sinalphas (-) = sine of solar elevation angle
Skyview (-) = fraction of diffuse skylight reaching this location
in central tree
st (m) = cumulative path of sun through crowns of all trees in stand
along current solar direction
stepsize (m) = increment in path search along solar beam direction
stmax (m) = max. value of st, at which search quits (penetration
probability is very low)
sumI (mol m-2 s-1) = sum of PAR irradiance over all leaves over all
simulation intervals
Tair (oC) = air temperature inside canopy
Tleafval (oC) = leaf temperature, renamed so it can be passed
in common
TLsum (W m-2) = accumulator for average (leaf T) x (leaf conduc-
tance) at this location
TLtree (oC) = sum of (leaf T) x (leaf conductance) over all
locations in tree
Tsky (oC) = effective radiative temperature of sky
Tveg (oC) = effective radiative temperature of vegetation
wt (-) = fraction of leaf area at this location in shade only
x (radians) = sine of solar azimuth; also (m), x-coordinate of
current point in central tree
xnew (m) = x-coordinate of current position in whole canopy during
search for interception by crowns
xo (m) = x-coordinate of currently-sampled point in central tree
y = y-coordinate of current point in central tree
ynew (m) = y-coordinate of current position in whole canopy during
search for interception by crowns
yo (m) = y-coordinate of currently-sampled point in central tree
z = z-coordinate of current point in central tree
znew (m) = z-coordinate of current position in whole canopy during
search for interception by crowns
zo (m) = z-coordinate of currently-sampled point in central tree
Arrays
real
atree(100) (m) = horizontal full axis of ellipsoid for given tree
btree(100) (m) = vertical full axis of ellipsoid for given tree
costheta(20) (-) = cosine of zenith angle of point being sampled in
central tree
lpenstore(20,20,20) (-) = penetration probability for direct beam
to all canopy locations, for current simulation interval
phi(20) (radians) = sky azimuths being sampled for diffuse skylight
phic(100) (radians) = azimuthal inclination of ellipsoid of each tree
Ppendiffuse(20,20,20) (-) = total penetration fraction of diffuse
skylight, for each location in central tree
r(20,20) (m) = array of radii being sampled in central tree, for
each choice of theta and phi
sintheta(20) (-) = sine of zenith angle of point being sampled in
central tree
thetac(100) (radians) = zenith inclination of ellipsoid of each tree
xtree(100) (m) = x-coordinate of center of each tree
ytree(100) (m) = y-coordinate of center of each tree
ztree(100) (m) = z-coordinate of center of each tree
do-loop indices
i - steps in distance along direct beam through the canopy
iphi - loop over azimuthal angles for positions sampled in
central tree
iradius - loop over radii for positions sampled in central tree
itheta - loop over zenith angles for positions sampled in
central tree
nb - loop over bins of PAR irradiance
nt - loop over all trees in stand
In control
real
deair (Pa) = increment in eair for attempt at solution
eair (Pa) = previous iteration's value of w.v. partial pressure
in canopy
eairtest (Pa) = provisional value of eair
fractcorr (-) = fraction of computed correction deair allowed, after
constraint of 2000 Pa as max. step
In Tcontrol
real
dTair (oC or K) = increment in Tair for attempt at solution
Tair (oC) = previous iteration's value of air temperature inside canopy
Tairtest (oC) = provisional value of Tair