Initialize Plant Variables
Subroutine Inplnt (cupinp.f)
This subroutine reads in plant-related input data.
READS IN:
IBIDIR 0: no bidirectional radiation fluxes,
1: bidirectional radiation fluxes
ITASSL 0: assuming top canopy layer has no tassels
1: assuming tassels present
IPP 1: only reflectance in principle plane is
calculated
NVZENP Number of view zenith angles in principle plane
VWZENP View zenith angles
NOZENV
VIEWZN
NOAZMV
VIEWAZ
ITOT Number of leaf angle classes
JMAX Maximum number of canopy layers (pres. 19)
JMIN Minimum number of canopy layers (pres. 2)
DFMIN Thinnest leaf layer in terms of LAI units
CLUMP Clumping factor for canopy structure.
Set to 1 for no clumping.
KMAX Maximum number of wavlengths
EMIS Thermal emissivity of the leaves
EMISOL Thermal emissivity of the soil
FACTIR Reduction factor for emissivity of sky in
atmospheric window
RSOIL(K) Soil Reflectance: (1): visible,
(2): near infrared, (3): thermal
RLLIVE(K) Reflectance of live veg at each wavelength
TLLIVE(K) Transmittance of live vegetation
RLDEAD(K) Reflectance of dead vegetation
TLDEAD(K) Transmittance of dead vegetation
A-CVISR Coefficients of parabola describing hemispherical
reflectance
A-CNIRR Coefficients of parabola describing hemispherical
reflectance
A-CVIST Coefficients of parabola describing hemispherical
transmittance
A-CNIRT Coefficients of parabola describing hemispherical
transmittance
ISPHER 1: Leaf angle distribution is spherical
0: Other inputs are used
NALPHA Number of leaf angle classes input in GR(ALPHA)
GR(NALPHA) Fraction of leaf area input at each angle alpha.
IMUNU Index controlling input leaf inclination angle
distribution. May have values 0, 1, 2.
1: meu & neu input directly to calc leaf angle dist
2: mean and std. dev. of leaf angle dist. input
XMEU One parameter of the beta distribution for the
leaf inclination angle distribution.
imunu=1: value of meu
imunu=2: mean angle of leaf angle dist(0=hor leaves)
XNEU Second parameter of the beta distribution for the
leaf inclination angle distribution.
imunu=1: value of neu
imunu=2: std. dev. of leaf angle dist.
Canopy Type meu neu (See Subroutine LAD)
planophile 2.770 1.172
erectophile 1.172 2.770
plagiophile 3.326 3.326
extremophile 0.433 0.433
uniform 1.000 1.000
spherical 1.101 1.930
IMUNUA Index controlling input leaf azimuth angle
distribution.
0: Leaves symmetric about azimuth (meu=neu=1.0)
1: Mean and std. dev. of azimuth angle dist
XMEUAZ One parameter of the beta distribution for the
leaf azimuthal angle distribution.
XNEUAZ Second parameter of the beta distribution for the
leaf azimuthal angle distribution.
BETA0 Azimuth angle of minimum leaf area in azimuthal
beta distribution.
NBETA No. of azimuth angle classes (50)
ROWSPC Row spacing (m)
PLSPC Plant spacing (m)
DMAX Maximum width of the leaf (m)
ILAIUP Code for updating canopy structure data,
(=3 for reading in).
PLTPM2 Number of plants per square meter
REFHTW Reference height for wind (= temp and vp)
ZDHCR =.2, critical layer height/canopy height
AMFULL Canopy velocity profile coefficient (1.5)
Z0DH Canopy roughness length divided by canopy height.
DISPDH Canopy displacement height divided by canopy height.
NLABCY Number of layers above the canopy
ZABCPY Height of layers above the canopy. Top one = upper
boundary condition.(m)
NLBCPY Number of layers below the canopy
ZBCPY Height of layers below the canopy (m)
RCUT20 Light response characteristics of stomata
RSMIN Light response characteristics of stomata
ANSTOM Light response characteristics of stomata
RADN Light response characteristics of stomata
TRSOPT Temp. response characteristics of stomata
TRSMAX Temp. response characteristics of stomata
TRSMIN Temp. response characteristics of stomata
RSEXP Temp. response characteristics of stomata
RSM Temp. response characteristics of stomata
PSI1 Leaf water potential parameter for stomata
PSI2 Leaf water potential parameter for stomata
CIDCAI If = -1, use CIDCA in code, otherwise use
CIDCA=CIDCAI
IC3C4 Code for C3 plant (0), C4 plant (1)
0: C3 Farquhar
1: C4 Corn Berry model w/ Index for stom
2: C4 Big Blue Stem Berry model w/ Index for stom
3: C3 Soybean Berry model w/ Index for stom
10: C3 Farquhar
11: C4 Corn Farquhar model stom use D1 & D2
12: C4 Big Blue Stem Farquhar model stom use D1 & D2
13: C3 Soybean Farquhar model stom use D1 & D2
CAIR Ambient atmospheric CO2 concentration (ppm)
RASTOM Ratio of stomatal resistance on top to
bottom of leaf
CUCOND Cuticilar conductivity to water vapor
(mol/m2/sec)
D1 VPD for start of stomatal closure (mb)
BKV Coefficient for vapor pressure deficit effect on
stomata.
CONMIN Minimum stomatal conductance (moles m-2 s-1)
RXCHAM Boundary layer resistance of one leaf side in
chamber for measuring P.S. (s/m)
NOTEMP Variable to override the temperature
dependence of photosynthesis
(=0, normal temperature dependence)
(=1, temperature = TFIX)
TFIX Fixed temperature of photosynthesis if the
temperature dependence is overriden.
RROOT Hydraulic root resistance (m4kg-1s-1)
FRSTEM Fraction of intercepted precip. that runs
down the stem
FRWTMX Max. fraction of leaf area wetted by precip.
PINTMX Maximum layer of water (mm) that can be held
on a wet leaf in a uniform film
NODYIN Number of days for which daily input values
are listed
LDAYIN Julian Day number of daily values to follow
XTLAI Total leaf area index at day LDAYIN
XFRLV Fraction of green vegetation at day LDAYIN
XHT Height of canopy at day LDAYIN
XAROT Rooting depth coefficient at day LDAYIN
XZLDH Relative height of the lowest green leaf at
day LDAYIN
XZMDH Relative height of the maximum leaf area density
at day LDAYIN
TSOLBC Soil temperature at soil boundary condition
at day LDAYIN
WSOLBC Volumetric water content at soil boundary condition
at day LDAYIN
COMPUTES:
ALEAF(K) Absorption of vegetation
RLFHEM
TLFHEM
These are all calculated again in cupidg2!!
RLFDIF [620-660]
TLFDIF
FR(i) Fraction of leaf area in each I'th leaf
inclination angle class including all azimuth angles.
A spherical distribution is assumed.
XINTZ extin. coef. at different zenith averaged over azimuth
XINTV Extinc. coeff at different viewing directions
RSFAC Factor to adjust leaf boundary layer resistance
for stomatal conductance or resistance ratio
These values are interpolated for each day of the run from the day values which have been input:
TLAI(IDAY) Total leaf area index at day IDAY (from XTLAI)
FRLIVE(IDAY) Fraction of LAI which is live vegetation on day IDAY
(from XFRLV)
HT(IDAY) Height of canopy at day IDAY (from XHT)
AROTIN(IDAY) Rooting depth coefficient at day IDAY (from XAROT)
ZLDH1(IDAY) Relative height of the lowest green leaf at
day IDAY (from XZLDH)
ZMDH1(IDAY) Relative height of the highest green leaf at
day IDAY (from XZMDH)
TSBC(IDAY) Soil temperature at soil boundary condition
at day IDAY (from TSOLBC)
WSBC(IDAY) Volumetric water content at soil boundary condition
at day IDAY (from WSOLBC)
Also, many coeffients for the photosynthesis equations are set.
Return to cupidg2