CUPID Variables

Use Mosaic's FILE -> "Find In Current" option to locate a specific variable definition.

I/O tag indicates whether variable is input or output.


ABSPK  		I    	Number of layers above the sprinkler where
               		droplet evaporation is occurring.
               		
ADUM      		Ratio of upward to downward fluxes for diffuse radiation.

ADUM1       	O	Ratio of up VIS radiation flux in a layer           
                	divided by downward flux                  

ADUM2       	O	Ratio of up NIR radiation flux in a layer   
                	divided by downward flux                      

ADUM3       	O	Ratio of up THERM radiation flux in a layer    
                	divided by downward flux 

AKCPSO      	O	Thermal conductivity for soil layer   
                	(Jm-1s-1K-1), or canopy layer (m2s-1)

AKCPY(50)        	Eddy conductivity for canopy layers

AKCPY1      	O  	Eddy conductivity of layer 1 above canopy  
                	(Jm-1s-1C-1)

AKCPYS      	O	Eddy conductivity above the soil surface 
                	at the critical height (Jm-1s-1C-1)   
                	              	
AKH(50)         O	Therm.conductivity of layer div.by layer height incr.
			(ms-1)

AKSMMH      	O	Saturated hydraulic conductivity (mmh-1)
			
AKSOIL(50)       	Thermal conductivity of soil layer

AKW         	O	Capillary conductivity of a layer in soil 
                	divided by layer depth increment (kgsm-4)
                	
ALAM             	Special conversion factor for wet bulb temperature

ALEAF(3)         	Leaf absorptivity for VIS, NIR and thermal

AM          	O	Wind profiles in canopy factor from  
                	Landsberg, James & Jarvis

AMFULL          I       Canopy velocity profile coefficient (1.5)   
             	
ANIRR      	I	Polynomial coefficient for NIR reflectance
               		as a function of source incidence angle.

ANIRT      	I	Polynomial coefficient for NIR transmittance
               		as a function of source incidence angle.

ANSTOM     	I	Light response characteristics of stomata.               		  

APAR        	O	Absorbed PAR by canopy (Wm-2)  

APARD       	O	Absorbed PAR by canopy (Wm-2) on a daily
                	basis.

APRTOT

AROOT       	O	Rooting depth coefficient, exponent on root 
                	length density vs depth (usually 1-10)

AROTIN(IDAY)    	Rooting depth coefficient at day IDAY

ASOLR       	O	Absorbed solar radiation by canopy (Wm-2) 

ASOLRD      	O	Daily total absorbed solar radiation by the
                	canopy (MJm-2d-1)
                	
AVISR      	I	Polynomial coefficient for VIS reflectance as
			a function of source incidence angle.
			  
AVIST      	I	Polynomial coefficient for VIS transmittance
               		as a function of source incidence angle.
               		  
BD         	I	Bulk density

BDYANG     	I/O	Body angle of sprinkler jet (degrees)

BETA0      	I	Azimuth angle of minimum leaf area in azimuthal 
               		beta distribution.

BKV        	I	Coefficient for vapor pressure deficit effect on
               		stomata.

BMFLX1      	O	Flux of direct beam radiation above a    
                	layer (Wm-2)                               
     
BMFLX2      	O	Flux of direct beam radiation above a        
                	layer (Wm-2)                                
           
BMFLX3      	O	Flux of direct beam radiation above a   
                	layer (Wm-2) 

BN4         	O	Bidirectional reflectance factor for different
                	wavelengths.   
                	             	
BNIRR      	I	Polynomial coefficient for NIR reflectance as a 
               		function of source incidence angle.

BNIRT      	I	Polynomial coefficient for NIR transmittance as
               		a function of source incidence angle.

BSPK       	I	Number of layers below the sprinker where droplet
               		evaporation is occurring.

BVISR      	I	Polynomial coefficient for VIS reflectance as a 
               		function of source incidence angle.

BVIST      	I	Polynomial coefficient for VIS transmittance as
               		a function of source incidence angle.

BX          	O	Power in soil potential vs. water content   

C2          	O	Predicted daily soil heat storage for   
                	Nebraska modified Penman equation 
                	           		             
CAIR       	I	Ambient atmospheric CO2 concentration (ppm)
               		   
CETSUM     		..HETSUM parameters for droplet evaporation
                  	calculation (see parameter library)   

CIDCA       	O	Ratio of leaf internal CO2 concentration
                	over the CO2 concentration of the air.   

CIDCAI		I       If = -1, use CIDCA in code, otherwise use
                        CIDCA=CIDCAI
                                  
CILAY       	O	Leaf internal CO2 concentration (ppm)   
               	            		  
CILEAF(I,J)     O	Intercellular CO2 concentration of leaves in 
                 	angle class I and layer J (1E-6 mol/mol)  

CLAI        	O	Cumulative leaf area index (top down) 
                                                                        
CLAI1       	O	Cumulative leaf area index from the top     
                	of the canopy                          

CLAI2       	O	Cumulative leaf area index from the top  
                	of the canopy                         

CLAI3       	O	Cumulative leaf area index from the top  
                	of the canopy  
                	               	
CLAYFC     	I	Clay fraction                                                    
                                                                        
CLODSZ     	I   	Soil clod size index 
                 	             		  
CLUMP      	I	Clumping factor for canopy structure.  

CNIRR      	I	Polynomial coefficient for NIR refectance as a 
               		function of source incidence angle.

CNIRT      	I	Polynomial coefficient for NIR transmittance as
               		a function of source incidence angle.

CO2SOL		O	Soil surface CO2 flux (umol CO2 m-2 s-1) 

CONCPY      	O	Condensation on the canopy (gm-2)   
                                                                        
CONDDY      	O	Maximum canopy condensation during the day  
                	(mgm-2=mm)
                	
CONMIN     	I	Minimum stomatal conductance (moles m-2 s-1)  
           		 
CONTOT[JMAX]	O	Total condensation in a layer (gm-2)
			CONTOT < 0 means layer had condensation.
			
CP(50)          O 	Thermal capacity of layer (MJm-2s-1K-1)

CPCPSO      	O	Thermal capacity of soil layer (MJm-3K-1) 
                	or canopy layer (kJm-3K-1)        

CPESTR      	O	Canopy storage of water vapor          
                                                                        
CPHSTR      	O	Heat storage in the canopy (Wm-2)

CPSOIL(50)       	Thermal capacity of soil layer

CPYTR            	Canopy transpiration?

CRCOEF      	O	Crop coefficient for simple ET calculations

CSFCLF       	O	CO2 conc at surface of leaf (umol/mol)	

CSLAY       	O	Leaf conductivity (stomatal and bound.lay.)     
                	(mol-1 H2O m2s)
                	
CSLEAF       		Leaf stomatal and boundary layer conductance 
                 	for angle class I and layer J (mol-1 H2O m2s)
                 	
CT(i)      		Fraction of beam flux density into angle class i in 
			beam flux as would be measured perpendicular to the 
			incidence angle.
			
CTWA       		..GTWA parameters for irrigation water
                  	temperature calculation (see parameter library)
                  	
CUCOND		I/O	Cuticilar conductivity for water vapor  (mol/m2/sec) 

CUMHR       	O	Plottting variable with cumulative day and    
                	decimal fractions of hours
                	
CVISR      	I	Polynomial coefficient for VIS reflectance as a 
               		function of source incidence angle.

CVIST      	I	Polynomial coefficient for VIS transmittance as 
               		a function of source incidence angle.

CW          	O	Capillary water capacity of soil divided      
                	by a layer depth increment (kgsm-4)   
                	            		
D(K,J)    		Downward diffuse radiation above sublayer J divided 
			by the total radiation above the canopy, 
			subscript wavelength band K.
			
D1         	I	VPD for start of stomatal closure (mb)                         
                                                                        
D2         	I	VPD for end of stomatal closure (mb)                           

D1          	O	Downward VIS radiation flux above a layer   
                	(Wm-2)                                     
         
D2          	O	Downward NIR radiation flux above a layer     
                	(Wm-2)                                           

D3          	O	Downward THERM radiation flux above a layer 
                	(Wm-2)   
                	                                                                     
DAKS       	I	Saturating water conductivity (optional).                      
                                                                        
DBX        	I	Exponent in the soil moisture release curve (optional).        
                                                                        
DELD       	I	Convergence criteria of diff. equation          
               		solution   
DELD        	O	Convergence criteria for soil water flow 
                	equations  (kgm-2s-1)                

DELD        	O	Convergence factor for soil water equations   
                                                                        
DELH20      	O	Change in soil water storage from last      
                	time step (mm) 

DELHA       	O	Coefficient in the temperature dependant
                	portion of electron transfer in photosynthesis.

DELHL       	O	Same as DELHA.

DELHH       	O	Same as DELHA.
    
DELT        	O	Temperature difference between leaf in angle 
                	class I and layer J and air in layer J (C)
			
DELTAE

DELTAT

DELWD       	O	Change in soil water storage from previous 
                	day (mm) 

DF         		Leaf area index increment of one canopy layer.

DFMIN      	I/O	Thinnest leaf layer in terms of LAI units

DIANOZ     	I	Sprinkler nozzle diameter (mm)

DISP        	O	Canopy displacement height (m) 

DISPDH		I	DISP/H

DISTLS(angle,hour)	Distribution of leaf-normal to sun angle for a given
			hour (sun direction).

DMAX       	I	Maximum size of leaf	

DMINC       	O	Dry matter increment (mg of dry weight m-2s-1)

DMINCD      	O	Daily dry matter increment
                	(mg of dry weight m-2d-1)
                	
DPDMIN     	I	For sprinkler irrigation: minimum droplet                        
               		diameter (mm) 	
               			  
DPE        	I	Air entry potential (J/kg) (optional). If 0 is
                  	entered it is calculated by G. Cambbell's
                  	equations.(-2.5J/kg)	

DPEVP       	O	Droplet evaporation from sprinkler irrigation.    

DPSTOR		O 	Heat storage in droplets (prop to TDROPI-TIRRIG) 

DPVMMH		O	Droplet evaporation in mm/hr	

DRAIND      	O	Daily total drainage at lower boundary  
                	condition (mm)                           

DRAIN       	O	Drainage out bottom of root zone  
                	(kgm-2s-1=mms-1)
                                         
DRGRAV      	O	Drainage from soil by gravity    
          	
DRHSFC     	I/O	Value to match soil water flow equations          
               		with canopy vapor flow equations

DRIP(20)         	???  
                			  
DSTNET(i,j)  	O	Net total radiation of all wavelengths into layer j
                  	and angle class i.
                  				  
DSTRAD(k,i,j)       	Total radiation incident in layer j and angle class 
                        j of wavelength k (diffuse+direct beam).

DSTNET      	O	Net all-wave radiation in angle class I and 
                	layer J (Wm-2)     

DSTNG			???

DSTNIR      	O	Incident NIR in angle class I and layer J   
                	(Wm-2)                   

DSTTHR      	O	Incident thermal radiation in angle class I  
                	and layer J (Wm-2)                      

DSTVIS      	O	Incident visible radiation in angle class I  
                	and layer J (Wm-2) 
                	                       
DT               	Length of time step for a single computation    

DTINT      	I	Length of the time step (sec) for ballistic                      
               		droplet calculations                                             

DTNEW      	I	Time at which to go to longer intervals
                      
EAIR         	O	Air vapor pressure in canopy (mb)

ECOMPL      	O	Light compensation point (irradiation level where
			C02 fixation through photosynthesis equals release 
			by respiration).

ECOND       	O	Vertical water flux in each layer (Wm-2)  
                                                                        
ECPY        	O	Evapotranspiration from canopy (Wm-2)      
                                                                        
ECPYD       	O	Daily mean evapotranspiration above canopy   
                	(Wm-2) (+ for evap. and - for cond.)

ECPYDV 		O       Canopy only evapotranspiration in mm/hr
                				
ECPYS       		Evaporation flux from soil calculated from 
                	canopy vapor flow equations (Wm-2)

ELUAP       	O	Efficiency of light use based on APAR  
                	(g CO2 MJ-1) 
                                                                        
ELUAS       	O	Efficiency of light use based on ASOLR     
                	(g CO2 MJ-1)
                                                                        
ELUIP       	O	Efficiency of light use based on PARI   
                	(g CO2 MJ-1)
                                                                        
ELUIS       	O	Efficiency of light use based on SOLRI      
                	(g CO2 MJ-1) 

ELUPAD      	O	Same as ELUAP except it is a daily total.

ELUPID      	O	Same as ELUIP except it is a daily total.

ELUSAD      	O	Same as ELUAS except it is a daily total.

ELUSID      	O	Same as ELUIS except it is a daily total.   
             	               				  
EMIS        	I/O	Thermal emissivity of leaf

EMISOL     	I/O	Thermal emissivity of the soil 	
		  
EN(MH,50)        	Vapor pressure (soil moisture) profile through
			all layers
			
END         	O	Sunset (hours.tenths)	

ESOIL      		Soil thermal radiation flux density

ET(50)           	Evapotranspiration from canopy layers      ???

ETBCD       	O	Daily ET from Blaney Criddle eq. (mm)

ETFIX      	I	Polynomial coefficients for droplet evaporation
               		rate for layers below the sprinkler height.

ETJHD       	O	Daily ET from Jensen-Haise eq. (mm)     
                                                                        
ETMIC       	O	Canopy evapotranspiration (1E-6 m)          
                                                                        
ETMMD       	O	Daily ET (mm)

ETMMH    	O	Canopy + soil evapotranspiration in mm/hr

ETOTW       	O	Evapotranspiration in a canopy layer 
                	(mgm-2s-1) (are these units right??? how is this
                	diff from ETOTWT???)

ETOTWT      	O	Evaporation flux above the canopy on a 
                	weight basis (mgm-2s-1)               		

ETPEND      	O	Daily ET from Penman equation (mm)

ETPM         	O	Canopy ET by Penman-Monteith equ. using soil
                 	heat flux from CUPID

ETPM2       	O	Canopy evapotranspiration by Penman-Monteith  
                	equations neglecting soil heat flux   
                	              	
ETPMD        	O	Daily ET from Penman-Monteith eq. (mm)

ETPNEB      	O	Daily ET from Nebraska modified Penman      
                	equation (mm)                    

ETPTD       	O	Daily ET from Priestly-Taylor eq. (mm)

ETSUM      	I	Polynomial coefficients for droplet evaporation
               		flux summed over all layers above the sprinkler
               		height.
               		
ETTOT       	O	Total droplet evaporation (+) or condensa-
                	tion (-) above the canopy (Wm-2)

EVAP        	O	Leaf ET in angle class I and layer J (mm)

EVAPG(J)        O	Layer average for leaf evapotranspiration (Wm-2)
			(= GEVAP(I,J) summed over angle class)
              	
EVINT(J)		???	

EVSMIC      	O	Soil evaporation (microns).

EVSMMH		O	Soil evap in mm/hr 

EVTOT       	O	Evaporation flux from the canopy on a    
                	weight basis (mgm-2s-1)	  
                	
EXPDIF    	O	Fraction of diffuse radiation that passes unintercepted 
               		through one layer.
               		
EXPDIR		O	The direct beam non-interception factor  

FACJ        	O	Coefficient in light dependent portion of 
                	photosynthesis. 
                	            		  
FACTE

FACTIR     	I	Factor that multiplies sky emissivity calculated
               		in subroutine SKYIR from Brusaert's equation.
               		See further notes on this topic.
               		
FACTV			Factor that wind value in input file is 
                  	multiplied by to get wind value used by model

FBEAM1      	O	Fraction of total visible radiation above   
                	canopy (beam as sensed by horizontal sensor)    

FBEAM2      	O	Fraction of total NIR radiation above     
                	the canopy 
                	
FBEAM(k)   	I	Fraction of direct beam radiation of wavelengths k in 
			total radiation of those wavelengths above canopy 
			(horiz.sensor).
			
FBEAM1		I	???			

FC          	O	Field capacity of -1/3 bar (% by volume).	
		
FEHIST      	O	Light history factor for photosynthesis
               	 	modeL	

FILT        	O	Infiltration into soil (mm)          
                                                                        
FILTD       	O	Daily total infiltration (mm)   
             	 			  
FR(i)      	O	Fraction of leaf area in each I'th leaf       
                 	inclination angle class including all azimuth angles.
                   	A spherical distribution is assumed.
                   	
FRAREA(i,j):         	Fraction of sunlit leaf area in layer j that is 
 			contained in angle class i (FRAREA of one layer
                        adds up to SUNLIT). 
                                          	  
FRAZ(I)     	O	Fraction of leaf area in each I'th leaf
                	azimuth angle class. 
                	
FREQR                	  

FRLIVE(IDAY)		Fraction of LAI which is live vegetation on day IDAY

FROOT(50)        	Fraction of root system extracting water in a layer 

FRSTEM     	I/O	Fraction of intercepted precip. that runs           
               		down the stem
               		
FRWET(20)        	???

FRWTMX          I/O	Maximum fraction of leaf area wetted by precipitation

FS         	I	Ratio of stomatal conductance on the top and
               		bottom of the leaf.

G2          	O	???

GEVAP(I,J)      O 	Evapotranspiration per ground area in angle class i
                	and layer j.  (Wm-2)
 
GHEAT(I,J)      O	Sensible heat flux per ground area in angle
                	class i and layer j.  (Wm-2)                	

GHWATR(I,J)		???

GR(ALPHA)  	I	Fraction of leaf area input at each angle alpha.  

GSIN        	O	Input CO2 concentration point in the absense of 
                	dark respirtaion (micromoles per mole).

H           	O	Height of canopy (m)                         
                                                                        
HCPY        	O	Sensible heat from soil and canopy to atmosphere         
                	(Wm-2)  (= QCOND(2))                      

HCPYD       	O	Daily mean sensible heat above canopy (Wm-2)  
                	(+ for canopy warmer than air)            

HCPYDL       	O	Daily mean sensible heat above canopy   
                	averaged over daylight hours (Wm-2)    
                	       
HCPYSD      	O	Daily mean sensible heat above soil (Wm-2)     
                                                                        
HCPYS       	O	Sensible heat flux from soil to canopy        
                	(Wm-2)                                    

HDSOIL      	O	Daily mean soil thermal conduction into        
                	soil (Wm-2) 
                	                	
HEATG(J)	O	Sensible heat flux per ground area in layer j.  (Wm-2)
			(=GHEAT(I,J) summed over angle class)

HEIGHT      	O	Height of the top of layer (m)            
                                                                        
HFDAY       	O	Half day length  (hours) 
                   	
HH 			Total conductance for evaporating water vapor from 
         		the leaves.

HPSI        	O	Water stress factor  

HSOIL       	O	Soil heat conduction flux (Wm-2)         
                	(+ for into soil)                    

HSOILD      	O	Same as HSOIL except it is a daily total.  
      		             		
HSLEAF       		Rel Hum at the surface of the leaf    
           	  
HT(IDAY)        	Height of canopy at day	IDAY	

HTDROP      	O	Maximum height of droplets (m)          
                                                                        
HTOT        	O	Sensible heat flux above canopy (Wm-2) 
			( = HEATG(J) summed over all canopy layers )

HTSPK      	I/O	Height of sprinkler nozzle (m)

IBIDIR Code:    I   	0:  no bidirectional radiation fluxes                   
                       	1:  bidirectional radiation fluxes
                       	
IC3C4       	I/O	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

ICUMDY      	I/O	Cumulative day                       	

IDAY        	O	Daily loop variable (from 1 to NODAYS)  
                                                                        
IDECL       	O	Solar declination 

IDOROC		I

IEQTM       	O	Equation of time for the sun

IHRWET[JMAX]		Accumulated hours of leaf wetness for a layer	

ILAIUP		I/O	Code for updating canopy structure data
			(=3 for reading in).

IMUNU       	I	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

IMUNUA      	I	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                 

INDAT       	O	Code for the frequency of data           
                	(=1 for hourly) 
                	
INOINT  	I                                                                
                                                                        
INTSKY			Flag indicating RADTOP(3) was not specified in input
            		so we should calculate it.
            		
IPOT     		Index indicating method of calculation for RATIO(IHR),
			=1: fixed

IPNT        	I	Code numbers for output format (see key)
    			
			
IPRTIN		I	
		
IROCLY	
		  
IRRCHK     	O	Flag to specify type of precipitation:
			  0 = no precipitation
                       	  1 = overhead sprinkler irrigation
                       	  2 = furrow irrigation
                       	  3 = rainfall

ISPEC Code  	I	(0: no vertical profiles printed,                           
                 	 1:prints vertical profiles (see key))  

ISPHER      	I	1: Leaf angle distribution is spherical
                	0: Other inputs are used  

ISPKTY      	I	Sprinkler type: -1: Impact, 0: Spray, 1: CDS     
                                                                        
ISPN        	I	Hours of special profile output          
                                                                        
ITASSL Code     I  	0: assuming top canopy layer has no tassels      
                      	1: assuming tassels are present 

ITER1       	O	Number of iterations required for convergence 
                	of scattered VIS radiation equations    

ITER2       	O	Number of iterations required for convergence 
                	of scattered NIR radiation equations    

ITER3       	O	Number of iterations required for convergence 
                	of scattered THERM radiation equations  
                	                    	
ITOT		I	# leaf (inclination?) angle classes

ITOTP1			ITOT + 1

IYEAR       	I/O	Year

JMAX       	I/O	Maximum number of canopy layers (pres. 19)        
                                                                        
JMIN        	I/O	Minimum number of canopy layers (pres. 2)         
                                                                        
JPNT        	I	Code numbers for output format (see key)

JTOT            O 	Number of canopy layers

KMAX        	I	Maximum number of wavelengths.

KMPDAY      	O	Daily windrun (kmday-1)

KPNT        	I	Code number for output of vertical profiles           
                	(see key)
                	
KSTRT      		Index indicating if IHR before or after sunrise/sunset:
                	(=3 for nighttime and =1 for daytime).
                	
LDAYIN      	I	Julian Day number of daily values to follow     
                                                                        
LPNT        	I	Code number for output of vertical profiles   
                	(see key)    

LYPDAY      	O	Daily total solar radiation (Langleys/day)

MDAY(LL)    	I	Day number associated with MDAYIN, where MDAY=1
                        signifies the first day of the simulation.

MDAYIN      	I	Number of days that profile output (or previous
                	line of input file) are printed. (MDAYIN=1 ->
                	output profile for all days)  
                	              	
MH			Related to hours?

MMH20       	O	Total water stored in root zone (mm)     
                                                                        
MONTH       	O	Month

NALPHA      	I	Number of leaf angle classes input in 
                	GR(ALPHA)
                	
NBETA      	I	No. of azimuth angle classes (50)

NDSOIL      	I/O	Number of depths in the soil
                	(number of soil layers)

NIR I       	O	NIR reflectance factors for calculation of  
                	bidirectional fluxes 
                	
NLABCY      	I	Number of layers above the canopy                  
                                                                        
NLBCPY      	I	Number of layers below the canopy.

NN2PNT      	I	Number of hours during the day when detailed
                	profile output is printed if ISPECL=1.         

NOAZMV      	I	Number of azimuth angles of view for BRDF.

NODAY       	O	Number of days of data input   
                                                                     
NODAYS      	I	Number of days for the simulation                
                                                                        
NODYIN     	I	Number of days for which daily input values   
                	are listed                                

NOHRS       	O	Number of hours of data per day     

NOTEMP      	O	Variable to override the temperature
                	dependence of photosynthesis
                	(=0, normal temperature dependence)
                	(=1, temperature = TFIX)
                	
NOZENV      	I	Number of zenith angles of view for BRDF.   
              	
NXINTZ			Dimension of XINTZ (90)

PA			Stomatal resistance water potential

PARD        	O	Daily accumulated PAR (MJm-2d-1)           
                                                                        
PARI        	O	Intercepted PAR by the canopy (Wm-2)     
                                                                        
PARID       	O	Daily intercepted PAR (Wm-2)

PATH		O	path*h = path length of sunbeam through canopy

PE          	O	???

PEVMIC      	O	Precipitation interception evaporation (microns)

PEVMMH		O	Interception evaporation in mm/hr

PHIMA           O 	Diabatic profile factor for momentum above the
			canopy

PHIMS           O 	Diabatic profile factor for momentum above the
			surface

PILAST(20)       	Value of PINT on last iteration.

PINT(20)        O 	Precipitation intercepted by the canopy layer

PINT1

PINTMX      	I/O	Maximum layer of water (mm) that can be held      
                	on a wet leaf in a uniform film 

PLSPC       	I/O	Plant spacing (m)   
             	
PLTPM2     	I	Number of plants per square meter

PN(MH,50)        	Soil layer water potential

POTBM1			Potential direct beam visible radiation at the earth's
 			surface.
 			
POTDIF  		Potential diffuse visible radiation at the earth's surface.

POTVIS	  		Potential total visible radiation at the earth's surface.			  

PRAREA      	O	Percent of leaf area in angle class I and 
                	layer J
                	
PRECIP(MH)      O 	Precipitation (mm) for hourly increment

PRECPD      	O	Daily accumulated precipitation (mm)      
                                                                        
PREEVD      	O	Daily intercepted precipitation evaporated    
                	(mm)                                       
           
PREEV       	O	Intercepted precipitation that evaporated     
                	this time step (mm)                        
           
PREINT      	O	Precipitation intercepted by canopy (mm)     
                                                                        
PROOT       	O	Hydraulic plant root resistance  (m4kg-1s-1) 
                                                                        
PS1         	O	Net photosynthetic rate of a leaf in layer I  
                	and angle class I  (1E-6 mol CO2 m-2 s-1)             

PSCPYL      	O	Photosynthetic rate of canopy on leaf area    
                	basis (umol CO2 m-2 s-1)       

PSCPYG      	O	Photosynthetic rate of canopy on ground area  
                	basis (umol CO2 m-2s-1)
                	
PSI1        	I/O	Leaf water potential where stomata 
                	begin to close (bars)
                	
PSI2        	I/O	Leaf water potential where stomata 
                	completely closed (bars) 

PSILF       	O	Leaf water potential (bars) 

PSIMA           O	Diabatic wind profile correction factor above the canopy                	          

PSIMS           O 	Diabatic wind profile correction factor above the soil

PSISUM      	O	Root length weighted soil water potential 
                	(bars)
                	
PSITOP      	O	Leaf water potential (bars) (at the top of the 
			canopy???)

PSITP1			PSITOP from previous iteration.

PSIXY        	O	Root xylem water potential (bars)

PSLAY       	O	Photosynthetic rate in each layer        
                	(1E-6 mol CO2 m-2 s-1)
                	
PSLEAF       	O	Leaf photosynthetic rate for angle class I
                 	and layer J (1E-6 mol m-2s-1)

PSMRM       	O	Photosyntheis minus maintenance respiration 

PSMRMD      	O	Same as PSMRM except it is a daily average.

PSYNGD      	O	Canopy photosynthetic rate based on ground 
                	area averaged over a day.

PWP         	O	Permanent wilting percentage at -15 bars
                	(% by volume)
                	
Q(50)            	Sensible heat flux into canopy layers      ???

QCOND       	O	Sensible heat flux by soil and canopy      
                	layer (Wm-2) 

QFIX        	I	Polynomial coefficients for sensible heat flux
                	for layers below the sprinkler height.        

QHRWET      	O	Daily cumulative hours of leaf wetness      
               
QITER2      	O	Iterations over canopy temperature and V.P.  
                	profiles to get convergence in main program
         
QITER3      	O	Iterations over soil evap.to merge soil water 
                	flow eq. & canopy vapor eq. in PROFL2 subr.
        
QITERW      	O	Iterations to match soil sfc W.C. with V.P.   
                	at the soil surface through log relationship   
        
QITERH      	O	Iterations over water stress in canopy        

QITOT       	O	Number of leaf angle classes                
  
QJMAX       	O	???

QJMIN       	O	???

QJZBOT      	O	Value of the index number of JZ of the lower  
                	soil boundary condition     

QJZCPY      	O	Value of the index number of JZ of the      
                	lowest layer containing leaves       

QJZCRT      	O	Critical layer in canopy where influence     
                	of the soil is significant on the wind 
                	speed profile in the canopy

QJZFSC      	O	Value of the index number of JZ of the node   
                	immediately above the soil surface      

QLAIUP      	O	???
                                                                        
QLOOPE      	O	Iterations over canopy water vapor flow      
                	equations alone in Newton-Raphson    

QLOOPT      	O	Iterations over soil heat flow equations      
                	alone in Newton-Raphson       

QLOOPW      	O	Iterations over soil water flow equations     
                	by Newton-Raphson in PROFL2    

QNODAY      	O	???

QNOITR      	O	Number of iterations in leaf energy budget                	
                                                                        
QRTZFC      	I	Quartz fraction in fractions of total weight

QSUM        	I	Polynomial coefficients for sensible heat flux
                	summed over all layers above the sprinkler height.

QTOT        	O	Total droplet sensible heat exchange above 
                	canopy (+ droplet cooling) (Wm-2)

QWETAM      	O	Hours of leaf wetness between midnight and
                	noon                      

QWETPM      	O	Hours of leaf wetness between noon and     
                	midnight
                	
RADSOL

RADABV(1,IHR)  		Ratio of the actual total visible radiation over
                     	the potential total visible radiation.
                     	
RADABV(2,IHR)  		Ratio of the actual total NIR radiation over the
                      	potential total NIR radiation.

RADN        	I	Light response characteristics of stomata ?????
			I THINK THIS MAY BE NET RADIATION.  
                  	
RADTOP(3,MH)    I 	Radiation onto a horizontal plane above the canopy
			1st index = wavelength: VIS, NIR, Thermal
			2nd index = hour.

RADTP1      	O	Visible radiation above canopy on horizontal  
                	(Wm-2)            

RADTP2      	O	Near infrared radiation above the canopy on   
                	a horizontal plane (Wm-2)                

RADTP3      	O	Thermal radiation above the canopy on a      
                	horizontal plane (Wm-2)
                	
RASTOM      	I/O	Ratio of stomatal resistance on top to   
                	bottom of leaf	
                			
RATIO   	O	Ratio of actual total NIR radiation over the potential
              		total visible and NIR radiation.

RATIOD      	I	Ratio of actual/potential solar radiation     
                	during the day 

RATION      	I	Ratio of actual/potential solar radiation   
                	at night                                 

RCUT20      	I 	Light response characteristics of stomata   

RDCPY       	O	Dark respiration of canopy.               
            
RDLAY       	O	Coefficient to adjust dark respiration for  
                	light history     
                	        	             		
RDRK         	O	dark respiration (umol CO2 m-2 s-1) 

REFDPT      	I/O	Reference depth    

REFHTT      	O	Reference height for temperature (m)  
       		  			
REFHTW      	I/O	Reference height for wind (= temp and vp)

RELHUM      	O	Relative humidity of the atmosphere

RGTOT       	O	Growth respiration of total.
                                                                        
RHCPY       	O	Canopy resistance to heat flow     
                	(sm-1)                                     
               
RHEAT       	O	Aerodynamic resistance to heat and water  
                	flow from diabatic wind profile equation 
                	
RESROT(50)       	Root resistance to water uptake?

RHLEAF      	O	Boundary layer resistance of leaves in   
                	a layer (s/m) 

RHMAX       	O	Daily maximum relative humidity at the  
                	reference height                           
             
RHMIN       	O	Daily minimum relative humidity at the     
                	reference height   
                	             	  
RHOCPY

RHSLOP      	I	Value to match soil water flow equations     
                	with canopy vapor flow equations

RLAY1       	O	Reflectance of diffuse VIS radiation from     
                	a canopy layer                             
             
RLAY2       	O	Reflectance of diffuse NIR radiation from     
                	a canopy layer                             
           
RLAY3       	O	Reflectance of diffuse THERM. radiation      
                	from a canopy layer     
                	           	
RLAYR(K)    		Reflectance factor for diffuse radiation of a layer. 
               		Subscripted by K for wavelength band.
               		
RLEAF1      	I/O	Leaf reflectance, visible radiation 

RLEAF2      	I/O	Leaf reflectance, near infrared  
             
RLEAF3      	I/O	Leaf reflectance, thermal

RLFDIF

RLFDIR			Calculated in DSTLIT (?)  Weighted mean sunlit leaf
			reflectance.  Perhaps this is for direct light, and
			above is for diffuse light.
			
RLFHEM			Leaf hemispherical reflectance

RMGR        	O	Maintenance respiration of grain.
                	(micromole m-2 ground area s-1)

RMGRD       	O	Same as RMGR except for a daily total.

RMLEAF      	O	Maintenance respiration of leaf.
                	(micromole m-2 ground area s-1)

RMLFD       	O	Same as RMLEAF except for a daily total.

RMROOT      	O	Maintenance respiration of root.
                	(micromole m-2 ground area s-1)

RMROTD      	O	Same as RMROOT except for a daily total.

RMSTEM      	O	Maintenance respiration of stem.
                	(micromole m-2 ground area s-1)

RMSTMD      	O	Same as RMSTEM except for a daily total.

RMTOT       	O	Maintenance respiration of total plant.
                	(micromole m-2 ground area s-1)

RMTOTD      	O	Same as RMTOT except for a daily total.
                                                                        
RN1         	O	Net VIS radiation flux above a layer 
                	(D1-U1)(Wm-2)                         

RN2         	O	Predicted net daily radiation for Nebraska   
                	modified Penman equation   

RN2         	O	Net NIR radiation flux above a layer         
                	(D1-U1)(Wm-2)                         

RN3         	O	Net THERM radiation flux above a layer        
                	(D1-U1)(Wm-2)                           

RNCPYD      	O	Daily mean net radiation above canopy       
                	(Wm-2) (+ for radiation input)        

RNCPY       	O	Net radiation above canopy (Wm-2)          
                                                                        
RNDIV       	O	Net radiation divergence in layer (Wm-2)    
                                                                        
RNET1       	O	Net all-wavelength flux above a layer        
                	(Wm-2)                                  

RNET2       	O	Net all-wavelength flux above a layer     
                	(Wm-2)                                         

RNET3       	O	Net all-wavelength flux above a layer      
                	(Wm-2)                                     

RNETMM      	O	Daily net radiation equivalent (mmday-1) 

RNDIV(j):   		Net total radiation flux density of all wavelengths
                	lost in layer j.
                	
RNLAM(k,j) 		Net radiation flux density into layer j of wavelength k
                	diffuse and direct beam.  (+) downward.  
                	                                                                      
RNSOIL      	O	Net radiation above soil (Wm-2)           
                                                                        
RNSOLD      	O	Daily mean net radiation above soil         
                	(Wm-2) 
                
ROOTSM			Sum of reciprocals of root resistances

ROOTUP(50)       	Water flow into root

ROWSPC      	I/O	Row spacing (m)

RROOT      	I	Hydraulic root resistance (3.E6 m4kg-1s-1)

RSAVE       	O	Average stomatal resistance of a layer of   
                	leaves (sm-1)

RSCPY       	O	Effective canopy stomatal resistance   
                	(sm-1)      
                	           	
RSEXP       	I	Temp. response characteristics of stomata

RSFAC		   	Factor to adjust leaf boundary layer resistance
			for stomatal conductance or resistance ratio
			
RSLAY       	O	Stomatal resistance (s/m)

RSLEAF      	O	Leaf stomatal resistance in leaf angle 
                	class i and layer j (sm-1).
                	
RSM         	I	Temp. response characteristics of stomata  
                                                                        
RSMIN       	I	Light response characteristics of stomata

RSNOVP      	O	Stomatal resistance with no direct vapor
                	pressure deficit effect (sm-1).

RSOIL       	I/O	Soil Reflectance: (1): visible,                  
                	(2): near infrared, (3): thermal   
                	             	
RSUBL     		Reflectance factor for diffuse radiation of a sublayer. 
               		Subscripted by K for wavelength band.

RSTOT        	O	Canopy resistance to water vapor flow 
                 	calculated from CUPID. Set=0 if precipitation
                 	is more than 0.5 mm  

RVCPY       	O	Canopy resistance to water vapor        
                	(sm-1)             
                	    	 
RXCHAM      	I/O	Boundary layer resitance of one leaf side in   
                	chamber for measuring P.S. (s/m)          

SANDFC      	I	Sand fraction   
              	              		  
SCOND[ITOTPL,JMAX]	Stomatal conductivity for water vapor?

SDN(j):     		Fraction of downward diffuse radiation in total 
			radiation above the canopy that originated in that 
			layer from direct beam interception.
			
SDNSUB    		Fraction of intercepted direct beam that is scattered 
               		downward in one sublayer.

SENS        	O	Canopy sensible heat flux by Penman-Monteith 
                	equations using soil heat flux from CUPID  
          
SENS2       	O	Canopy sensible heat by Penman-Monteith eq.  
                	neglecting soil heat flux
                	
SILTFC      	I	Silt fraction

SIZELF      	O	Maximum length of leaf (m)                   
                                                                        
SOLARD      	O	Daily accumulated solar radiation (MJm-2d-1)  
                                                                        
SOLAR       	O	Incident solar radiation on a horizontal      
                	plane above the canopy (Wm-2)  

SOLRID      	O	Daily intercepted solar radiation (Wm-2)    
                                                                        
SOLRI       	O	Intercepted solar radiation by canopy (Wm-2) 
               	              		
SOURDN(K,IHR)    	Source term for downward diffuse radiation for the 
               		wavelength band K and IHR for one layer.
               		
SOURUP(K,IHR)    	Source term for upward diffuse radiation for the 
               		wavelength band K and IHR for one layer.

SPKPRS      	I/O	Sprinkler pressure (kPa)
               		
STDLNG      	I/O	Standard longitude (90 deg. central time zone)

STEM        	O	Stem flow of intercepted precipitation (mm)

STRT         	O	Sunrise (hours.tenths) 

SUNAZM      	O	Solar azimuth (0 degrees = south)

SUNLIT     		LAI of a layer that is sunlit. 

SUP(j):     		Fraction of upward diffuse radiation in total radiation 
			above the canopy that originated in that layer from 
			direct beam interception.             		  

SUPSUB    		Fraction of intercepted direct beam that is reflected 
               		upward in one sublayer.  

SWDAY       	O	Daily mean soil water storage from      
                	the previous day (mm)

T1		O	Average vegetation kinetic temperature

T2(THETV,PHIV)	O	Average directional kinetic temperature at view
			zenith angle THETV, view az angle PHIV
			
T4		O	Hemispherical infrared temperature

T5(THETV,PHIV)	O	Apparent directional infrared temperature at view
			zenith angle THETV, view az angle PHIV
				                	
TAEROH      	O	Canopy aerodynamic temperature using
                	Z  for heat (C).
                          h
                          
TAEROM      	O	Canopy aerodynamic temperature using the zero-plane 
			displacement height for momentum (C).

TAIR        	O	Air temperature in the canopy (C)   
                                                                        
TAIRD       	O	Daily mean air temperature at reference    
                	height (C)                    

TAIRMX     	O	 Maximum temperature of air at reference   
                	height for day (C)              

TAIRMN      	O	Minimum temperature of air at reference   
                	height for day (C) 
                	
TBEAM(j)   	O	Fraction of direct beam in direct beam above the canopy
                	passing through layer j.  It is decremented in each
                	layer due to physical interception by the leaves in
                	that layer - it has nothing to do with the radiative
                	properties of the leaves themselves.              		               		

TCPY        	O	Average canopy temperature (C)  
                                                                        
TCPYD       	O	Daily mean canopy temperature (C)          
                                                                        
TCROP       	O	Crop temperature from P.-M. equations

TDROPI      	I	Initial droplet temperature at sprinkler (C)

TEMAIR     	I/O	Air temperature  (C)

TEMPLF(i,j) 	O	Leaf temperature in angle class i and layer j.

TEMPLF(ITOTP1,j)  	Average leaf temperature over all angle classes
                       	in layer j.

TEMSOL      	O	Soil temperature at the lower boundary (C)

TFIX        	O	Fixed temperature of photosynthesis if the
                	temperature dependence is overriden by NOTEMP.

THALFL      	O	Coefficient in temperature dependence of
                	electron transfer in photosynthesis.

THALFH      	O	Same as THALF.

TIMLOC      	I/O	Local time (hours.tenth)

TIMSUN       	O	Solar time (hours.tenths)

TIR #       	O	Apparent temperature of canopy (IR-thermom.) 
                	at # deg. view angle from zenith.
                	
TIRETQ

TIRRIG      	I/O	Average temperature of water droplets when they
                	reach the top of the canopy (C).

TLAY1       	O	Transmittance of diffuse VIS radiation for   
                	a canopy layer       

TLAY2       	O	Transmittance of diffuse NIR radiation for   
                	a canopy layer    

TLAY3       	O	Transmittance of diffuse THERM. radiation   
                	for a canopy layer 
                	
TLAYR(K)     		Transmittance factors for diffuse radiation of a layer. 
               		Subscripted by K for wavelength band. 

TLEAF1      	I/O	Leaf transmittance, visible                

TLEAF2      	I/O	Leaf transmittance, near infrared            

TLEAF3      	I/O	Leaf transmittance, thermal

TLFAVE      	O	Average leaf temperature in canopy (C)

TLFDIF

TLFDIR			Calculated in DSTLIT (?) Weighted mean sunlit leaf
			transmittance.  Perhaps this is for direct light,
			and above is for diffuse light.

TLFHEM			Leaf hemispherical transmittance

TOTLAI		O	Total leaf area index

TLAI(IDAY)      	Total leaf area index at day IDAY

TN(MH,50)        	Temperature profile through all layers

TPRECP      	O	Precip. transmitted through the canopy

TRAN(10,20)     O 	Transpiration of a canopy layer and angle class

TRMIC       	O	Canopy transpiration (1E-6 m)   
                                                                        
TRMMD       	O	Daily total transpiration (mm)   

TRMMH		O	Canopy transpiration in mm/hr   

TROOTM      	O	Root-weighted mean temperature of soil (C).

TRSMAX      	I	Temp. response characteristics of stomata 
                                                                        
TRSMIN      	I	Temp. response characteristics of stomata              
                                                                        
TRSOPT      	I	Temp. response characteristics of stomata

TSBC(IDAY)    		Soil temperature at soil boundary condition
              		at day IDAY 
              		 
TSFC       		Soil surface temperature.

TSFCSD      	O	Daily mean soil surface temperature (C)   
                                                                        
TSFCSO      	O	Soil surface temperature (C)

TSN(50)          	Temperature of soil layer   ???

TSOL24      	O	Soil temperature profiles at midnight

TSOLBC      	I	Soil temperature at soil boundary cond. at day         
                	LDAYIN
                	
TSUBL     		Transmittance factors for diffuse radiation of a sublayer. 
               		Subscripted by K for wavelength band.

TWB1

U

U(K,J)   	O	Upward diffuse radiation above sublayer J divided by 
			the total radiation above the canopy, subscript 
			wavelength band K.
			
USTARA       	O	Friction velocity above the canopy (m/s)

USTARS          O 	Friction velocity above the soil surface			  

USTR        	O	Friction velocity used in P.-M. eq. for   
                	calculation of aerodynamic resistance
                	to heat and water flow 
                	
VIEWAZ      	I/O	View azimuth angle

VIEWZN      	I	View zenith angle 

VIS I       	O	Visible radiation reflectance factors for   
                	calculation of bidirectional fluxes.

VM          	O	Maximum velocity of carboxylation. 

VPAIR(MH)       O 	Air vapor pressure at reference height above the canopy

VPAIRD      	O	Daily mean air vapor pressure at reference
                	height (mb) 
 
VPD         	O	Vapor pressure deficit (mb) 
               	
VPDIN        	O	Water vapor pressure deficit at upper boundary (mb)

WATER       	I	Initial water content profile               

WATSOL      	O	Soil water content at the lower boundary
                	(volume fraction)  

WCPYSD      	O	Daily mean soil surface evaporation (Wm-2)   
                                                                        
WCPYS       	O	Soil surface evaporation calculated from     
                	soil water flux equations (Wm-2)  
                	              	                                                                      
WIND        	I/O	Wind speed  (m/s)

WINDDY      	O	Daily mean wind speed at reference height      
                	(m/s)               

WS          	O	Saturation water content (% by volume)

WSBC(IDAY)    		Volumetric water content at soil boundary condition
              		at day IDAY
              		
WSMM        	O	Soil surface evaporation (mm)              
                                                                        
WSMMD       	O	Daily mean soil surface evaporation (mm/day)

WSOL24      	O	Soil volumetric water content at midnight 

WSOLBC      	I	Volumetric water content at soil boundary cond.
                	at day LDAYIN

WT(MH,50)        	Volumetric water content of soil layer

WT3BAR      	O	3-bar water content (% by volume)

WTI(50)          	Water content of soil layers

WTIR(J,IANGLE)		Weighting factor for thermal IR interception by
                        each canopy layer.
                        
WTP(20)          	Weighting factor for canopy layer precip. intercept

WTSFC		O	Volumetric water content of top soil layer

XAROT       	I	Rooting depth coefficient at day LDAYIN

XFRVL           I	Fraction of green vegetation at day LDAYIN

XHT         	I	Height of canopy at day LDAYIN

XINT(ihr)		Weighting factor for canopy radiation interception
			calculations.  It involves the sun->leaf-normal 
			angle.

XINTV(THETV,PHIV)       Extinc. coef. at different viewing angles			

XINTZ(NXINTZ)		extin. coef. at different zenith averaged over azimuth

XLAT        	I/O	Latitude

XLONG       	I/O	Longitude

XMEU        	I	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)

XMEUAZ      	I	One parameter of the beta distribution for the 
                	leaf azimuthal angle distribution.

XNEU        	I	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

XNEUAZ      	I	Second parameter of the beta distribution for the
                	leaf azimuthal angle distribution.
                                                                        
XTLAI       	I	Total leaf area index at day LDAYIN               
                                                                        
XZLDH       	I	Relative height of the lowest green leaf at
                	day LDAYIN                      

XZMDH       	I	Relative height of the maximum leaf area density    
                	at day LDAYIN

Z0          	O	Roughness length of canopy (m)  
              	
Z0DH            I       Z0/H 

ZABCPY      	I	Height of layers above canopy. Top one = upper 
                	boundary condition.(m)          

ZBCPY       	I	Height of layers below the canopy (m)

ZCRIT       	O	Height of the layer JZCRIT  (m)

ZDH         	O	Height of top of layer divided by the   
                	canopy height
                	
ZDHCR       	I	=0.2, critical layer height/canopy height

ZDL2        	O	Height divided by Monin-Obukov length used  
                	for diabatic profile corrections in P.-M. calc. 

ZDLA        	O	Height divided by Monin-Obukov length above          
                	the canopy                                 

ZDLS        	O	Height divided by Monin-Obukov length above  
                	the soil                                                       

ZENANG      	O	Solar zenith angle (degree)                     
                                                                        
ZLDH        	O	Fraction of canopy height below which there  
                	are no leaves
                	
ZLDH1(IDAY)		Relative height of the lowest green leaf at
                        day IDAY

ZMAX        	I	Maximum height of droplet trajectory (m)

ZMDH        	O	Fraction of canopy height with maximum leaf
                	density
                	
ZMDH1(IDAY)    		Relative height of the highest green leaf at
              		day IDAY

ZMID        	O	Height of midpoint of each layer (m)     

ZNON        	O	Non-rectangular hyperbola coefficient

ZSOIL       	I	Depth of nodes in soil in meters (ZOSOIL in output)

Further notes on selected variables

FACTIR

FACTIR is the factor that the thermal sky flux, as calulated in the subroutine skyir, is multiplied by to get the sky thermal flux in the band of interest. It is used in the subroutine bdrtm as such:

d(kt,jtot)=radabv(3)*factir

Here are some values it can take:

Return to definition