LIST OF ALL COMMON STATEMENTS IN XCUPID5B
Name Description
RAD1 Leaf, soil and layer optical properties:
EMIS Leaf thermal emissivity
EMISOL Soil thermal emissivity
RSOIL(3) Soil reflectance for VIS, NIR and thermal
RLEAF(3) Leaf reflectance for VIS, NIR and thermal
TLEAF(3) Leaf transmissivity for VIS, NIR, thermal
ALEAF(3) Leaf absorptivity for VIS, NIR and thermal
EXPDIF Non-interception factor for diffuse (one layer)
RLAYR(3) Layer reflectance for VIS, NIR and thermal
TLAYR(3) Layer transmissivity for VIS, NIR, thermal
(INPLNT,RADIAT,LFEBAL)
-----------------------------------------------------------------------------
RAD2 Some radiation transfer results from RADIAT
DSTRAD(3,10,20) Total radiation into layer and angle class
DSTNET(10,20) Net total radiation into layer and angle class
FRAREA((10,20) Fraction of LAI in angle class in sunlit LAI of layer
TEMPLF(10,20) Leaf temperature
TSOIL(MH) Soil surface temperature
(RADIAT,STOMA,PROFL2,PHOTO1,ROOTEX)
-----------------------------------------------------------------------------
RAD3 Radiation input variables
RADTOP(3,MH) Radiation onto a horizontal plane above the canopy
FBEAM1(3,MH) Direct beam fraction on a horizontal plane above canopy
COSZEN(MH) cosine of solar zenith angle
ZENANG(MH) solar zenith angle
HFDAY length of a half day (hr)
RATIOD Ratio of actual/potential radiation during the day
RATION Ratio of actual/potential radiation in the night
RATIO(MH) Ratio of actual/potential radiation for each hour
(SKYIR,RADIN4,PENMAN)
-----------------------------------------------------------------------------
RAD4 Radiation fluxes into layers
D(3,20) Downward diffuse flux into layer (fraction of RADTOP)
U(3,20) Upward diffuse flux into layer (fraction of RADTOP)
BMFLX(3,20) Direct beam flux density into layer
RNET(20) Net total radiation flux density into a layer
RNDIV(20) Net total radiation flux density lost in a layer
TSFC Soil surface temperature
RNLAM(3,20) Net total radiation into a layer by wavelength
(RADIAT,PROFL2,PENMAN)
----------------------------------------------------------------------------
RAD5 Source terms for diffuse radiation
SOURUP(3,MH) Source of upward diff. rad. in a layer (fr.of RADTOP)
SOURDN(3,MH) Source of downward diff. rad. in a layer (fr. of " )
(RADIAT)
-----------------------------------------------------------------------------
MISC1 Constants and print out code numbers
PI 3.14159
PID180 3.14159/180 (converts from degrees into radians)
SIGMA 5.67 E-08
IWRITE(9,99) Code numbers for print out
(INFIX,INPLNT,RADIAT,LFEBAL,PENMAN,SIMPET)
-----------------------------------------------------------------------------
MISC2 Mostly layer and angle class structure data
ITOT Number of leaf angle classes
ITOTP1 ITOT + 1
JTOT Number of canopy layers
FR(10) Fraction of LAI of angle classes in layer
CT(10) Fraction of perpend. beam flux inc. at angle
TOTLAI Total leaf area index of canopy
DF Leaf area index of one layer
DT Length of time step for a single computation
(INPLNT,RADIAT,WINDFAC,HEIGHT,RBOUND,STOMA,LFEBAL,STRESS,HITE2,PROFL2,
PHOTO1,INTERC,DRIPS,ROOTEX,PENMAN)
-----------------------------------------------------------------------------
MISC3 Mostly canopy structure input data
H Height of canopy (m)
COVER NOT USED ANY MORE!
ZDH(20) Height of layer divided by canopy height
ROWSPC Row spacing (m)
ZOSOIL Depth of soil layer
ZLDH Fraction of canopy height below which no leaves present
ZMDH Fraction of canopy height with maximum leaf density
(INPLNT,INSOIL,WINDFAC,HEIGHT,HITE2,PROFL2,PENMAN)
-----------------------------------------------------------------------------
MISC4 Whole model layer structure
Z(50) Heights of the top of layers
ZMID(50) Heights of the layer midpoints
JZCPY Index number of lowest canopy layer
JZCPY1 Index number of layer just below canopy
JZSFC Index number of layer just above soil surface
JZSFM1 Index number of layer above layer JZSFC
JZBOT Index number of lowest soil layer (boundary)
JZCRIT # of layer in canopy where infl. of soil on wind prof.
JZBM1 Index number of layer above layer JZBOT
(RADIAT,HITE2,PROFL2,ROOTEX)
-----------------------------------------------------------------------------
MISC5 Physical layer increments
DELTZ(50) Layer height increments between layer midpoints (m)
(HITE2,PROFL2)
-----------------------------------------------------------------------------
TIME Temporal input data
MONTH Month of the year
JDAY Day of the month
IYEAR Calendar year
ICUMDY Cumulative day of the year
TIMLOC(MH) Local time in decimal fractions of hours
(DATE,DECLIN,ZENITH,PROFL2,)
-----------------------------------------------------------------------------
ASTRON Geometric data for location
EQTM Equation of time for the sun
DECL Solar declination angle
SINDEC sine of solar declination angle
COSDEC cosine of solar declination angle
DECMAX sine of the maximum solar declination angle
SINLAT sine of the latitude
COSLAT cosine of the latitude
TANLAT tangent of the latitude
DLONG time difference of solar time from local time (hr)
(INFIX,DECLIN,ZENITH)
-----------------------------------------------------------------------------
WIND1 Data for calculation of wind speed profiles
FWIND(20) Wind speed in canopy layer below JZCRIT ????
WIND(MH) Wind speed at ref. height above the canopy
DMAX Maximum size of leaf (m)
REFHTW Reference height for wind speed
Z0 Roughness length of canopy (m)
DISP Canopy displacement height (m)
AM Factor for wind profiles in canopy
ZCRIT Height of the lowest leaf (m)
(INPLNT,WINDFAC,RBOUND,PROFL2)
-----------------------------------------------------------------------------
MET1 Boundary conditions and water source
TEMAIR(MH) Air temperature at reference height above the canopy
VPAIR(MH) Vapor pressure at reference height above the canopy
PRECIP(MH) Precipitation (mm) for hourly increment
TEMSOL(MH) Soil temperature at the reference depth
WATSOL(MH) Soil volumetric water content at reference depth
(SKYIR,STRESS,PROFL2,INTERC,PENMAN)
-----------------------------------------------------------------------------
RESIS1 Leaf resistance related parameter
RHLEAF(20) Boundary layer resistance of leaves in a layer
RSLEAF(10,20) Stomatal resistance in angle class and layer
HPSI Water stress factor
(RADIAT,RBOUND,STOMA,LFEBAL,PROFL2)
-----------------------------------------------------------------------------
RESIS2 Input parameter for stomatal behaviour
RADN Light response characteristic of stomata
ANSTOM Light response characteristic of stomata
RCUT20 Light response characteristic of stomata
RSMIN Light response characteristic of stomata
TRSOPT Temperature response characteristic of stomata
TRSMIN Temperature response characteristic of stomata
TRSMAX Temperature response characteristic of stomata
RSM Temperature response characteristic of stomata
RSEXP Temperature response characteristic of stomata
PSI1 Leaf water potential parameter for stomata
PSI2 Leaf water potential parameter for stomata
(STOMA,PENMAN)
-----------------------------------------------------------------------------
LEAF1 Leaf transpiration parameter
DELT(10,20) Temp.difference betw. leaf and air (by layer and angle)
PSILF Leaf water potential (bar)
TRAN(10,20) Transpiration of a canopy layer and angle class
(RADIAT,LFEBAL,STRESS,DRIPS,ROOTEX)
-----------------------------------------------------------------------------
LEAF2 Fluxes across leaf in angle classes + layers
EVAP(10,20) Leaf E.T. of angle class and layer
GEVAP(10,20) Leaf E.T. for ground area for angle class + layer
HEAT(10,20) Leaf sensible heat for leaf area for " + "
GHEAT(10,20) Leaf sensible heat for ground area for " + "
ALAM Special conversion factor for wet bulb temperature
(RADIAT,LFEBAL,STRESS,PROFL2,DRIPS,ROOTEX)
-----------------------------------------------------------------------------
PROF1 Mostly profile related parameter through canopy ??
TAIR(20) Air temperature of layer in canopy
EAIR(20) Atmospheric vapor pressure of a layer in the canopy
PHIM =PHIMA, diabatic profile factor for momentum
PHIH =PHIMS, diabatic profile factor for heat
REFHTT Reference height for temperature (m)
REFHTE Reference height for vapor pressure (m)
RELH(20) Relative humidity in the canopy layer
NLABCY Number of layers above the canopy
NLBCPY Number of layers between the canopy and ground
(INPLNT,LFEBAL,HITE2,PROFL2)
-----------------------------------------------------------------------------
PROF2 Main results of profile calculations
TN(MH,50) Temperature profile through all layers
AKCPY(50) Eddy conductivity for canopy layers ???
CPCPY(50) Thermal capacity for canopy layers
UPROF(50) Upward diffuse radiation flux into canopy layer
Q(50) Sensible heat flux into canopy layers ???
ET(50) Evapotranspiration from canopy layers ???
EN(MH,50) Vapor pressure (soil moisture) profile through all layers
QCOND(50) Sensible heat flux through all layers
ECOND(50) Vertical water flux through all layers
TCHECK(50) ???
ESAT(50) Saturation vapor pressure for all layers ???
(PROFL2)
-----------------------------------------------------------------------------
PROF3 Some parameter for profile calculations
USTARA Friction velocity above the canopy
USTARS Friction velocity above the soil surface
PHIMA Diabatic profile factor for momentum
PHIMS Diabatic profile factor for heat
AKH(50) Therm.conductivity of layer div.by layer height incr.
CP(50) Thermal capacity of layer (MJm-2s-1K-1)
CP2(50) CP * deltz/dt
(PROFL2)
-----------------------------------------------------------------------------
PROF4 Coefficients for wind profile calculations
PSIMA Diabatic wind profile correction factor above the canopy
PSIMS Diabatic wind profile correction factor above the soil
ZDLA Height divided by the Monin-Obukov length above canopy
ZDLS Height divided by the Monin-Obukov length above the soil
(PROFL2)
-----------------------------------------------------------------------------
CPY1 Flux averages over layer
ETOTWT Evapotranspiration above canopy on weight basis (mgm-2s-1)
HTOT Sensible heat flux above the canopy
RPLNT Number of plants per row
EVAPG(20) Layer average for leaf evapotranspiration ???
HEATG(20) Layer avg for leaf sensible heat on leaf area basis
EAVE(20) Layer avg for leaf evapotransp. on ground area basis ???
HAVE(20) Layer avg for leaf sensible heat on ground area basis ???
(STRESS,PROFL2)
-----------------------------------------------------------------------------
CPY2 Fluxes and storage of soil and canopy
HSOIL Soil heat conduction flux
HCPYS Sensible heat flux from soil to canopy
EVTOT Evaporation flux from canopy on weight basis
CONTOT(20) Total condensation in a layer (gm-2)
SCOND(10,20) ????????
IHRWET(20) Accumulated hours of leaf wetness for a layer
ECPYS Evapotrans. from soil (calc. from canopy vapor flow)
CPHSTR Heat storage in the canopy
CPESTR Water vapor storage in the canopy
WCPYS Soil surface evaporation (calc. from soil water flux)
(LFEBAL,STRESS,PROFL2,PENMAN)
-----------------------------------------------------------------------------
SOIL1 Some soil layer properties
ZSOIL(25) Height of soil nodes (m)
WTI(50) Water content of soil layers
NDSOIL Number of soil layers
AKSOIL(50) Thermal conductivity of soil layer
CPSOIL(50) Thermal capacity of soil layer
TSN(50) Temperature of soil layer ???
(INSOIL,INITC,HITE2,PROFL2)
-----------------------------------------------------------------------------
SOIL2 Water flux related parameters for soil
AKSOL(50) ?????
AKW(50) Capillary conductivity of soil layer div. by depth incr.
CW(50) Capillary water capacity of soil layer div. by " "
WT(MH,50) Volumetric water content of soil layer
ESAVE(3,MH) ????
(PROFL2)
-----------------------------------------------------------------------------
SOIL3 Soil type input
CLODSZ Clod size for soil
(INSOIL,PROFL2)
-----------------------------------------------------------------------------
SOIL4 Soil property input data
PE ???
BX ???
BD Bulk density
AKS ???
AN ???
WS Volumetric water content
SANDFC Sand fraction
SILTFC Silt fraction
CLAYFC Clay fraction
QRTZFC Quartz fraction in fractions of total weight
(INSOIL,PROFL2)
----------------------------------------------------------------------------
WATER1
IPRECP ???
TPRECP Precipitation transmitted through the canopy
PN(MH,50) Soil layer water potential
WCOND(50) ???
WSTOR(50) ???
(PROFL2,DRIPS,ROOTEX)
-----------------------------------------------------------------------------
WATER2 Water flow through model boundaries
SW
DELD Convergence criterion for soil water flow equations
RHSLOP ???
DRHSFC Value to match soil water flow equations with canopy
DRAIN Drainage at lower boundary conditions
FILT Total infiltration
ETMM Canopy evapotranspiration
DRGRAV Drainage from soil by gravity
(PROFL2)
-----------------------------------------------------------------------------
WATER3
SWLAST ???
(WATER3)
-----------------------------------------------------------------------------
PSYN1 Summary of photosynthetic rates
PS1(10,20) Photosynthetic rate by angle class and layer
PSLAY(20) Photosynthetic rate of layer
PSCPYL Canopy photosynthesis based on leaf area
PSCPYG Canopy photosynthesis based on ground area
(RADIAT,PHOTO1)
-----------------------------------------------------------------------------
ITERAT Iteration counters
NOITER Loop counter for leaf temperature convergence
ITER2 Loop counter for leaf water potential convergence
ITER3 Loop counter for convergence of soil evaporation
ITERW Loop counter for conv. of water content with v.p. at sfc
ILOOP Iterations over soil water flow in NEWRAP
LOOPT NOT USED?
(PROFL2)
------------------------------------------------------------------------------
INTER1 Precipitation interception parameter
WTP(20) Weighting factor for canopy layer precip. intercept
FRWET(20) ???
FRWTMX Maximum fraction of leaf area wetted by precipitation
PINT(20) Precipitation intercepted by the canopy layer
PILAST(20) ???
(INPLNT,LFEBAL,DRIPS)
-----------------------------------------------------------------------------
INTER2 Intercepted precip. distribution
EVINT(20) ???
EVIMM(20) ???
PINTMX Maximum depth of water held on a wet leaf
FRSTEM Fraction of intercepted precip. that runs down the stem
DRIP(20) ???
STEM Intercepted precipitation running down the stem
(INPLNT,LFEBAL,DRIPS)
-----------------------------------------------------------------------------
ROOT1
FROOT(50) Fraction of root system extracting water in a layer
RESROT(50) Root resistance to water uptake?
ROOTSM ???
PSIXY Xylem water potential
PSITOP Leaf water potential
PROOT Percentage of root system extracting water in layer
ROOTUP(50) Water flow into root
CPYTR Canopy transpiration?
PSISUM Weighted effective soil water potential for root zone
(INPLNT,STOMA,PROFL2,ROOTEX)
-----------------------------------------------------------------------------
MITDAM Mite information
XMDCM(20)
XMDMIN
XMDMAX
RSDAM