Cupid Output Variables

Sorted by code number


    

    LIST OF CUPID OUTPUT VARIABLES IN THE ORDER OF HOW THEY FOLLOW IN THE
    OUTPUT FILE 

    NAME         DEFINITION                                  (M,NN,C)   PROGREF
    ____         __________                                  ________   ______
    
    
    
			***** FIXED HEADER DATA ****
			
    
    ITOT         Number of leaf angle classes                  1,01,1   
    
    JMAX	 Maximum allowed number of canopy layers       1,01,2
    
    JMIN	 Minimum allowed number of canopy layers       1,01,3 
    
    DFMIN        Minimum increment in LAI units                1,01,4 
    
    XLAT         Site latitude                                 1,01,5   
                                                                        
    XLONG        Site longitude                                1,01,6   
    
    STDLNG       Standard longitude for time zone              1,01,7  
    
    EMIS         Thermal emissivity of leaf                    1,01,8 
    
    EMISOL       Thermal emissivity of soil                    1,01,9  
    
    RSOIL1       Soil reflectance, visible radiation           1,02,1    
                                                                        
    RSOIL2       Soil reflectance, near infrared radiation     1,02,2
                                                                        
    RSOIL3       Soil reflectance, thermal radiation           1,02,3   
    
    RLIVE1       Live leaf reflectance, visible radiation      1,02,4  
                                                                        
    RLIVE2       Live leaf reflectance, near infrared          1,02,5    
                                                                        
    RLIVE3       Live leaf reflectance, thermal                1,02,6 
    
    TLIVE1       Live leaf transmittance, visible              1,02,7
                                                                        
    TLIVE2       Live leaf transmittance, near infrared        1,02,8    
                                                                        
    TLIVE3       Live leaf transmittance, thermal              1,02,9                       

    ILAIUP       Code for updating canopy structure data       1,03,1
                 (=3 for reading in)

    NODAY        Number of days of data input                  1,03,2

    REFHTT       Reference height for temperature (m)          1,03,3
    
    REFHT						       1,03,4    

    REFHTW       Reference height for wind speed               1,03,5   
                 measurement (m)

    REFDPT       Reference depth for soil temperature          1,03,6 
                 measurement at lower boundary condition (m)

    RROOT        Hydraulic plant root resistance  (m4kg-1s-1)  1,03,7

    NOTEMP       =1 if not temp effect on photosynthesis       1,03,8
		 =0 if normal temp effect on photosynthesis 

    TFIX         if notemp=1 then this is the temp 	       1,03,9
		 photosynthesis is evaluated at.			      

    FR(I)        Fraction of leaf area in the I'th leaf        1,04,1-9   
                 angle class, also frac in each of 50 azm 
                 classes

    ROWSPC       Row spacing (m)                               1,05,1

    PLSPC        Plant spacing (m)                             1,05,2

    Z0SOIL       Soil roughness parameter                      1,05,3

    NDSOIL       Number of soil nodes including the surface    1,05,4

    DELD         Convergence criteria for soil water flow      1,05,5
                 equations  (kgm-2s-1)

    DRHSFC       Value to match soil water flow equations      1,05,6
                 with canopy vapor flow equations

    FRSTEM       Fraction of intercepted precipitation         1,05,7
                 that runs down the stem

    FRWTMX       Maximum fraction of leaf area wetted by       1,05,8
                 precipitation

    PINTMX       Maximum depth of water (mm) that can be       1,05,9
                 held on a wet leaf in a uniform film

    IC3C4        Index to indicate mode of plant               1,06,1
                 photosynthesis  (1:C4, 0:C3)
		 0   Farquhar C3
		 1   Berry corn
		 2   Berry prairie grass
		 3   Berry soybean
		 10  Farquhar C3
		 11  Farquhar corn
		 12  Farquhar soybean

    Leaf photosynthesis coefficients from Farquahar,1980 
    (1,06-07):
    
    T0		 Reference temperature for respiration         1,06,2
	  	 and k's (C)

    XK1	    	 Mult coefficient for kc (microbars)           1,06,3
    
    XK2	 	 Exp coefficient for kc	(J/micromol)           1,06,4
    
    XK3		 Mult coefficient for ko (millibars)	       1,06,5
    
    XK4	 	 Exp coefficient for ko (J/micromol)           1,06,6
    
    XK5	 	 Exp coefficient for rdark (J/micromol)        1,06,7
    
    XK6	 	 Exp coefficient for vm (j/micromol)           1,06,8
    
    RO		 Dark respiration rate at t0 		       1,06,9
                 (micromol/m2/s)					 
    
    VM           Maximum velocity of carboxylation	       1,07,1
                 t0 (micromol/m2/s) 

    XJ2		 Maximum rate electron transfer 	       1,07,2
                 (micro eq/m2/s)
    
    DELHA        Coefficient in the temperature dependent      1,07,3
                 portion of electron transfer in 
                 photosynthesis.   
                  	 
    DELHL        Same as DELHA.				       1,07,4

    DELHH        Same as DELHA.				       1,07,5
    
    THALFL       Coefficient in temperature dependence of      1,07,6
                 electron transfer in photosynthesis.    
                 
    THALFH       Same as THALF.				       1,07,7  
    
    CIDCA        Ratio of leaf internal CO2 concentration      1,07,8
                 to the CO2 concentration of the air.  
                 
    GSIN         Input CO2 concentration point in the          1,07,9
                 absence of dark respiration 
                 (micromoles per mole).                                         

    PWP		 Soil permanent wilting percentage             1,08,1

    WT3.0	 Soil Volumetric water content at -3 bars      1,08,2

    WT0.3	 Soil Volumetric water content at -.3 bars     1,08,3

    WT0.1	 Soil volumetric water content at -.1 bars     1,08,4

    WT.05	 Soil volumetric water content at -.05 bars    1,08,5

    WT.01	 Soil volumetric water content at -.01 bars    1,08,6
    
    HTSPK        Height of sprinkler (m)                       1,09,1

    HTDROP       Maximum height of droplets (m)                1,09,2

    SPKPRS       Sprinkler pressure (PSI)                      1,09,3

    BDYANG       Sprinkler body angle (degrees)                1,09,4

    TDROPI       Irrigation water temperature		       1,09,5
    
    OX		 Oxygen concentration (millibars)	       1,10,1
    
    CAIR         Ambient atmospheric CO2 concentration (ppm)   1,10,2
    
    RASTOM       Ratio of stomatal resistance on top to	       1,10,3
                 bottom of leaf

    CUCOND       Cuticular conductivity for water vapor        1,10,4
                 (micromol/m2/sec)                 
                 
    RXCHAM       Boundary layer resistance in the chamber      1,10,5
                 used to measure leaf photosynthesis (s/m)     

    FACJ         Coefficient in light dependent portion of     1,10,6
                 photosynthesis. 
                 
    ZNON         Non-rectangular hyperbola coefficient         1,10,7
                 from Farquhar and Wong, 1984  
                 
    D1		 Vapor pressure deficit for start of           1,10,8
                 stomatal closure (mbar)
    
    BKV          Parameter for vapor pressure deficit          1,10,9
                 response of stomata                    
                                   
    PE(JZ1)	 Air entry potential for soil (J/kg)	       1,11,1

    BX(JZ1)	 Power in moisture release curve eqn.	       1,11,2

    AKSMMS(JZ1)  Saturated water conductivity in mm/second     1,11,3

    WS(JZ1)	 Saturated soil water content(m3/m3)	       1,11,4 
    
    SLOPE        Slope of site				       1,12,1
    
    ASPECT	 Aspect angle of site			       1,12,2   

    RDEAD1       Dead leaf reflectance, visible radiation      1,13,1  
                                                                        
    RDEAD2       Dead leaf reflectance, near infrared          1,13,2    
                                                                        
    RDEAD3       Dead leaf reflectance, thermal                1,13,3 
    
    TDEAD1       Dead leaf transmittance, visible              1,13,4
                                                                        
    TDEAD2       Dead leaf transmittance, near infrared        1,13,5    
                                                                        
    TDEAD3       Dead leaf transmittance, thermal              1,13,6 
    
    
    
    
    			***** DAILY HEADER DATA *****  
    			                   

    IYEAR        Year                                          2,01,1          

    ICUMDY       Cumulative day number (JAN 1= 1)              2,01,2      

    MONTH        Month                                         2,01,3     

    H            Height of canopy (m)                          2,01,4      

    TOTLAI       Total leaf area index                         2,01,5      

    IDECL        Solar declination                             2,01,6

    IEQTM        Equation of time for the sun                  2,01,7

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

    HFDAY        Half day length  (hours)                      2,02,1

    STRT         Sunrise (hours.tenths)                        2,02,2

    END          Sunset (hours.tenths)                         2,02,3

    IDAY         Daily loop variable (from 1 to NODAYS)        2,02,4

    FACTIR       Fraction that multiplies thermal IR sky       2,02,5
		 radiation from 4-100 um to simulate the
		 sky thermal radiation in a waveband of an
		 infrared thermometer that might be used to
		 determine canopy temperature.  For example
		 8-14 um.  If this is zero than emiscp is the
		 true canopy emissivity.

    SIZELF       Maximum length of leaf (m)                    2,03,1

    Z0           Roughness length of canopy (m)                2,03,2

    DISP         Canopy displacement heigth (m)                2,03,3

    AM           Wind profiles in canopy factor from           2,03,4
                 Landsberg, James & Jarvis

    ZLDH         Fraction of canopy height below which         2,03,5
                 there are no leaves

    ZMDH         Fraction of canopy height with maximum leaf   2,03,6
                 density

    JTOT         Number of canopy layers +1                    2,03,7
                 (to include soil surface)

    NDSOIL       Number of depths in soil including the        2,03,8
                 surface and lower boundary condition

    JZCPY        Value of the index number of JZ of the        2,04,1
                 lowest layer containing leaves

    JZFSC        Value of the index number of JZ of the node   2,04,2
                 immediately above the soil surface

    JZBOT        Value of the index number of JZ of the lower  2,04,3
                 soil boundary condition
                 
    NLABCY	 Number of layers above the canopy             2,04,4	
    
    NLBCPY      Number of layers below the canopy              2,04,5          

    PSI1         Leaf water potential where stomata            2,04,6
                 begin to close (bars)

    PSI2         Leaf water potential where stomata            2,04,7
                 completely close (bars)

    TAIRD        Daily mean air temperature at reference       2,05,1
                 height (C)

    TAIRMX       Maximum temperature of air at reference       2,05,2
                 height for day (C)

    TAIRMN       Minimum temperature of air at reference       2,05,3
                 height for day (C)

    VPAIRD       Daily mean air vapor pressure at reference    2,05,4
                 height (mb)

    WINDDY       Daily mean wind speed at reference height     2,05,5
                 (m/s)

    SOLARD       Daily accumulated solar radiation (MJm-2d-1)  2,05,6

    PARD         Daily accumulated PAR (MJm-2d-1)              2,05,7

    PRECPD       Daily accumulated precipitation (mm)          2,05,8

    TCPYD        Daily mean canopy temperature (C)             2,06,1

    TSFCSD       Daily mean soil surface temperature (C)       2,06,2

    RHMAX        Daily maximum relative humidity at            2,06,3
                 the reference height

    RHMIN        Daily minimum relative humidity at            2,06,4
                 the reference height

    RNCPYD       Daily mean net radiation above canopy         2,07,1
                 (Wm-2) (+ for radiation input)

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

    ECPYD        Daily mean evapotranspiration above canopy    2,07,3
                 (Wm-2) (+ for evap. - for cond.)

    HSOILD       Daily mean soil thermal conduction            2,07,4
                 into the soil (Wm-2)

    RNSOLD       Daily mean net radiation above the soil       2,07,5
                 (Wm-2)

    HCPYSD       Daily mean sensible heat above soil (Wm-2)    2,07,6

    WCPYSD       Daily mean soil surface evaporation (Wm-2)    2,07,7

    WSMMD        Daily mean soil surface evaporation (mm/day)  2,07,8

    CONDDY       Maximum canopy condensation during the day    2,07,9
                 (mgm-2=mm)

    SWDAY        Daily mean soil water storage from the        2,08,1
                 previous day (mm)

    DELWD        Change in soil water storage from the         2,08,2
                 previous day (mm)

    DRAIND       Daily total drainage at the lower boundary    2,08,3
                 condition (mm)

    FILTD        Daily total infiltration (mm)                 2,08,4

    ETMMD        Daily ET (mm)                                 2,08,5

    PREEVD       Daily intercepted precipitation which has     2,08,6
                 evaporated (mm)

    TRMMD        Daily total transpiration (mm)                2,08,7

    PARID        Daily intercepted PAR (Wm-2)                  2,08,8

    SOLRID       Daily intercepted solar radiation (Wm-2)      2,08,9

    ETPEND       Daily ET from Penman equation (mm)            2,09,1

    ETPTD        Daily ET from Priestly-Taylor eq. (mm)        2,09,2

    ETJHD        Daily ET from Jensen-Haise eq. (mm)           2,09,3

    ETBCD        Daily ET from Blaney-Criddle eq. (mm)         2,09,4

    KMPDAY       Daily windrun (kmday-1)                       2,09,5

    RNETMM       Daily net radiation equivalent (mmday-1)      2,09,6

    LYPDAY       Daily total solar radiation (Langleys/day)    2,09,7

    ETPMD        Daily ET from Penman-Monteith eq. (mm)        2,09,8

    ETPNEB       Daily ET from Nebraska modified Penman        2,09,9
                 equation (mm)

    HRWET        Daily cumulative hours of leaf wetness        2,10,1

    WETAM        Hours of leaf wetness between midnight and    2,10,2
                 noon

    WETPM        Hours of leaf wetness between noon and        2,10,3
                 midnight

    TSOL24       Soil temperature profiles at midnight         2,11,1

    WSOL24       Soil volumetric water content at midnight     2,11,2

    RN2          Predicted net daily radiation for Nebraska    2,12,1
                 modified Penman equation (ly/day)

    G2           Predicted daily soil heat storage for         2,12,2
                 Nebraska modified Penman equation (ly/day)

    CRCOEF       Crop coefficient for simple ET calculations   2,12,3

    PARID        Daily intercepted PAR (Wm-2)                  2,13,1

    SOLRID       Daily intercepted solar radiation (Wm-2)      2,13,2

    APARD        Daily absorbed PAR  radiation (Wm-2)          2,13,3

    ASOLRD       Daily absorbed    solar radiation (Wm-2)      2,13,4

    PSYNGD       Daily canopy photosynthetic rate              2,13,5

    ELUPAD       Daily PAR absorbed light use efficiency       2,13,6
                 (g/MJ)
    
    ELUSAD       Daily SOL absorbed light use efficiency       2,13,7
 		 (g/MJ)
    
    ELUPID       Daily PAR intcpted light use efficiency       2,13,8
  		 (g/MJ)
    
    ELUSID       Daily SOL absorbed light use efficiency       2,13,9
		 (g/MJ)
		
    RMSTMD	 Daily total of maintenance respiration        2,14,1
  		 of stem (micromole m-2 ground area s-1) 					       
    
    RMLFD	 Daily total of maintenance respiration        2,14,2
  		 of leaf (micromole m-2 ground area s-1) 
  		    
    RMROTD	 Daily total of maintenance respiration        2,14,3
    		 of root (micromole m-2 ground area s-1)					       
    
    RMGRD	 Daily total of maintenance respiration        2,14,4
   		 of grain (micromole m-2 ground area s-1)   
   		  
    RMTOTD	 Daily total of maintenance respiration        2,14,5
    		 of total plant 
    		 (micromole m-2 ground area s-1)
    		 
    RGTOTD	 Daily total of total growth respiration       2,14,6
    
    PSMRMD	 Daily average of photosynthesis minus         2,14,7
                 maintenance respiration
    
    DMINCD	 Daily dry matter increment		       2,14,8	       
                 (mg of dry weight m-2d-1)	    
    
    RDDAY						       2,15,1
    
    RDDAY2						       2,15,2 
       
    PSS1AD						       2,15,3 
      
    PSS1BD						       2,15,4 
     
    PSS1CD						       2,15,5 
      
    PSS1DD						       2,15,6 
    
    PSS2AD						       2,15,7 
      
    PSS3AD						       2,15,8 
      
    PSS3BD						       2,15,9 
    
    PSYNGD       Canopy photosynthetic rate based on ground    2,16,1
                 area averaged over a day.	 
    
    PSS4D						       2,16,2 
       
    PSS2BD						       2,16,3  
    
    RNLFD	 					       2,17,1 
     
    EVLFD						       2,17,2 
      
    QLFD						       2,17,3  
    
    RNCPYD       Daily mean net radiation above canopy         2,18,1
    		 (mj m-2 day-1)
                 (+ for radiation input)    
    
    HCPYD        Daily mean sensible heat above canopy         2,18,2
                 (mj m-2 day-1)
                 (+ for canopy warmer than air)

    ECPYD        Daily mean evapotranspiration above canopy    2,18,3
    	  	 (mj m-2 day-1)
                 (+ for evap. - for cond.)

    HSOILD       Daily mean soil thermal conduction            2,18,4
                 into the soil (mj m-2 day-1)

    RNSOLD       Daily mean net radiation above the soil       2,18,5
                 (mj m-2 day-1)

    HCPYSD       Daily mean sensible heat above soil           2,18,6
    		 (mj m-2 day-1)

    WCPYSD       Daily mean soil surface evaporation           2,18,7 
    		 (mj m-2 day-1)
    		 
    DRAN5D	 Daily mean drainage from soil node 5 up       2,20,1
                 from jzbot		               
                     
    HCPYDL       HCPY averaged over the daylight hours         2,21,1
    
    
    
    
    			***** HOURLY HEADER DATA *****
    			

    TIMLOC       Local time (hours.tenths)                     4,01,1

    TIMSUN       Solar time (hours.tenths)                     4,01,2

    ZENANG       Solar zenith angle (degree)                   4,01,3

    SUNAZM       Solar azimuth angle (0 degrees = south)       4,01,4

    RADTP1       Visible radiation above the canopy on a       4,01,5
                 horizontal plane (Wm-2)

    RADTP2       Near infrared radiation above the canopy      4,01,6
                 on a horizontal plane (Wm-2)

    RADTP3       Thermal radiation above the canopy on a       4,01,7
                 horizontal plane (Wm-2)

    FBEAM1       Fraction of total visible radiation above     4,01,8
                 canopy (beam as sensed by horizontal sensor)

    FBEAM2       Fraction of total NIR radiation above the     4,01,9
                 canopy

    TEMAIR       Air temperature (C)                           4,02,1

    VPAIR        Atmospheric vapor pressure (mb)               4,02,2

    WIND         Wind speed (m/s)                              4,02,3

    PRECIP       Precipitation (mm)                            4,02,4
                 (Set to 0 for surface irrigation)

    TEMSOL       Soil temperature at the lower boundary (C)    4,02,5

    WATSOL       Soil water content at the lower boundary      4,02,6
                 (volume fraction)

    SOLAR        Incident solar radiation on a horizontal      4,02,7
                 plane above the canopy (Wm-2)

    RATIO        Ratio of actual solar radiation to potential  4,02,8
                 clear-day solar radiation

    CUMHR        Plottting variable with cumulative day and    4,02,9
                 decimal fractions of hours

    VPDIN        Water vapor pressure deficit at upper         4,03,1
                 boundary (mb)

    RELHUM       Relative humidity at upper boundary	       4,03,2

    IRRCHK       Irrigation code                               4,03,3
		  0   no irrigation or precip
		  1   above canopy sprinkler
		  2   furrow irrigation
		  3   rain

    FACTV        Factor that wind value in input file is       4,03,4
		 multiplied by to get wind value used by 
		 model
		
    PATH         path*h = path length of sunbeam through       4,03,5
                 canopy	
                 
    ZENSLP       Angle of sunbeam perpendicular to slope       4,04,1
    
    RADAB1						       4,04,2
    
    RADAB2 						       4,04,3
    
    RADAB3						       4,04,4
    
    FBEAM1						       4,04,5
    
    FBEAM2 						       4,04,6  
    
    
    
    
    			***** HOURLY SUMMARY BY LAYER *****   
    			      

                 5,1-3 are two lines of output each, hence the
                 duplicitous listings.  Subscripts 1,2, and 3
                 refer to VIS, NIR and THERM, respectfully.
                  
    ITER1        Number of iterations required for convergence 5,01,1
                 of scattered VIS radiation equations

    EXDIR1       Fraction of direct VIS beam that passes       5,01,2
                 unintercepted through a single canopy layer

    EXDIF1       Fraction of diffuse VIS radiation that passes 5,01,3
                 unintercepted through a single canopy layer

    CLUMP1       Canopy clumping factor after Nilson(1971)     5,01,4

    CLAI1(J)     Cumulative leaf area index from the top       5,01,1
                 of the canopy 
                    
    RLAYER1(J)   Reflectance of diffuse VIS radiation from     5,01,2  
                 a canopy layer    

    TLAYER1(J)   Transmittance of diffuse VIS radiation for    5,01,3 
                 a canopy layer

    D1(J)        Downward VIS radiation flux above a layer     5,01,4
                 (Wm-2)

    U1(J)        Upward VIS radiation flux above a layer       5,01,5
                 (Wm-2)

    TBEAM1(J)    Transmittance of direct beam radiation to     5,01,6
                 above a layer (above canopy: TBEAM=FBEAM)

    BMFLX1(J)    Flux of direct beam VIS radiation above a     5,01,7
                 layer (Wm-2)

    RN1231(J)    Net VIS radiation flux into a layer           5,01,8
                 (D1-U1)(Wm-2)

    RNET1(J)     Net all-wavelength flux above a layer         5,01,9
                 (Wm-2)

    ITER2        Number of iterations required for convergence 5,02,1
                 of scattered NIR radiation equations

    EXDIR2       Fraction of direct NIR beam that passes       5,02,2
                 unintercepted through a single canopy layer

    EXDIF2       Fraction of diffuse NIR radiation that passes 5,02,3
                 unintercepted through a single canopy layer

    CLUMP2       Canopy clumping factor after Nilson(1971)     5,02,4

    CLAI2(J)     Cumulative leaf area index from the top       5,02,1
                 of the canopy 
                    
    RLAYER2(J)   Reflectance of diffuse NIR radiation from     5,02,2  
                 a canopy layer    

    TLAYER2(J)   Transmittance of diffuse NIR radiation for    5,02,3 
                 a canopy layer

    D2(J)        Downward NIR radiation flux above a layer     5,02,4
                 (Wm-2)

    U2(J)        Upward NIR radiation flux above a layer       5,02,5
                 (Wm-2)

    TBEAM2(J)    Transmittance of direct beam radiation to     5,02,6
                 above a layer (above canopy: TBEAM=FBEAM)

    BMFLX2(J)    Flux of direct beam NIR radiation above a     5,02,7
                 layer (Wm-2)

    RN1232(J)    Net NIR radiation flux into a layer           5,02,8
                 (D2-U2)(Wm-2)

    RNET2(J)     Net all-wavelength flux above a layer         5,02,9
                 (Wm-2)

    ITER3        Number of iterations required for convergence 5,03,1
                 of scattered THERM radiation equations

    EXDIR3       Fraction of direct THERM beam that passes     5,03,2
                 unintercepted through a single canopy layer

    EXDIF3       Fraction of diffuse THERM radiation that      5,03,3
                 passes unintercepted through a single canopy 
                 layer

    CLUMP3       Canopy clumping factor after Nilson(1971)     5,03,4

    CLAI3(J)     Cumulative leaf area index from the top       5,03,1
                 of the canopy 
                    
    RLAYER3(J)   Reflectance of diffuse THERM radiation from   5,03,2  
                 a canopy layer    

    TLAYER3(J)   Transmittance of diffuse THERM radiation for  5,03,3 
                 a canopy layer

    D3(J)        Downward THERM radiation flux above a layer   5,03,4
                 (Wm-2)

    U3(J)        Upward THERM radiation flux above a layer     5,03,5
                 (Wm-2)

    TBEAM3(J)    Transmittance of direct beam radiation to     5,03,6
                 above a layer (above canopy: TBEAM=FBEAM)

    BMFLX3(J)    Flux of direct beam THERM radiation above a   5,03,7
                 layer (Wm-2)

    RN1233(J)    Net THERM radiation flux into a layer         5,03,8
                 (D3-U3)(Wm-2)

    RNET3(J)     Net all-wavelength flux above a layer         5,03,9
                 (Wm-2)

    HEIGHT       Height of the top of the layer (m)            5,04,1 

    ZDH          Height of top of the layer divided by         5,04,2 
                 the canopy height (canopy layers only)

    CONTOT       Total condensation for a layer                5,04,3 

    ETOTW        Evapotranspiration in a canopy layer          5,04,4         
                 (mgm-2s-1)

    VPD          Vapor pressure deficit (mb)                   5,04,5       

    QCOND        Sensible heat flux by the soil and            5,04,6        
                 canopy layer (Wm-2)

    TLFAVE       Average leaf temperature in canopy (C)        5,04,7

    ECOND        Vertical water flux in each layer (Wm-2)      5,04,8

    PSLAY        Photosynthetic rate in each layer             5,04,9       
                 (1E-6 mol CO2 m-2 s-1)

    RNDIV        Net radiation divergence in layer (Wm-2)      5,05,1        

    PINT         Precipitation intercepted by a canopy layer   5,05,2
                 (mm)

    TRAN         Transpiration of a canopy layer (Wm-2)        5,05,3         

    EVAPG        Evaporation/ground area by layer (Wm-2)       5,05,4         

    HEATG        Sensible heat/ground area (Wm-2)              5,05,5       

    TAIR         Air temperature in canopy (C)                 5,05,6

    U            Wind speed (ms-1)                             5,05,7

    EAIR         Air vapor pressure in canopy (mb)             5,05,8

    RH           Relative humidity of air (percent)            5,05,9

    ZMID         Height of midpoint of each layer (m)          5,06,1

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

    CPCPSO       Thermal capacity of soil layer (MJm-3K-1)     5,06,3
                 or the canopy layer (kJm-3K-1)

    AKH          Thermal conductivity of soil or atmosphere    5,06,4
                 in canopy div. by layer height incr. (ms-1)

    CP           Thermal capacity of soil or atmosphere in     5,06,5
                 the canopy (Jm-2s-1K-1)

    RSAVE        Average stomatal resistance of a layer of     5,06,6
                 leaves (sm-1)

    RHLEAF       Boundary layer resistance of leaves in        5,06,7
                 a layer (s/m)

    AKW          Capillary conductivity of a layer in soil     5,06,8
                 divided by a layer depth increment (kgsm-4)

    CW           Capillary water capacity of soil divided      5,06,9
                 by a layer depth increment (kgsm-4)

    CLAI         Cumulative leaf area index (from top down)    5,07,1

    CSLAY        Leaf conductivity (stomatal and bound.layer)  5,07,2
                 (mol-1 H2O m2s)

    RSLAY        Stomatal resistance (s/m)                     5,07,3

    PSLAY        Leaf photosynthetic rate (E-6 mol m-2s-1)     5,07,4

    CILAY        Leaf internal CO2 concentration (ppm)         5,07,5

    FEHIST       Light history factor for photosynthesis       5,07,6
                 model

    RDLAY        Coefficient to adjust dark respiration for    5,07,7
                 light history

    ECOMPL       Light compensation point                      5,07,8
    
    
    
    
    			***** HOURLY SUMMARY DATA *****
    			

    TSFCSO       Soil surface temperature (C)                  7,01,1

    PSILF        Leaf water potential (bar)                    7,01,2

    TCPY         Average canopy temperature (C)                7,01,3

    ETOTWT       Evaporation flux above the canopy on          7,01,4
                 a weight basis (mgm-2s-1)

    EVTOT        Evaporation flux from the canopy              7,01,5
                 (Wm-2)

    HTOT         Sensible heat flux above canopy (Wm-2)        7,01,6

    HSOIL        Soil heat conduction flux (Wm-2)              7,01,7
                 (+ for into the soil)

    CPHSTR       Heat storage in canopy (Wm-2)                 7,01,8

    HPSI         Water stress factor (Obsolete)                7,01,9           

    USTARA       Friction velocity above the canopy (m/s)      7,02,1

    USTARS       Friction velocity above the soil surface      7,02,2
                 (m/s)

    PHIMA        PHIm, diabatic profile factor for momentum    7,02,3
                 above the canopy

    PHIMS        PHIm, diabatic profile factor for momentum    7,02,4        
                 above the soil surface

    AKCPY1       Eddy conductivity of layer 1 above            7,02,5
                 the canopy (Jm-1s-1C-1)

    AKCPYS       Eddy conductivity above the soil surface      7,02,6
                 at the critical height (Jm-1s-1C-1)

    HCPYS        Sensible heat flux from the soil to canopy    7,02,7
                 (Wm-2)

    ECPYS        Evaporation flux from soil calculated from    7,02,8
                 the canopy vapor flow equations (Wm-2)

    RNSOIL       Net radiation above the soil (Wm-2)           7,02,9

    JZCRT        Critical layer in the canopy where influence  7,03,1
                 of soil is significant on wind speed profile
                 within the canopy

    ZCRIT        Height of the layer JZCRIT  (m)               7,03,2

    PSIMA        Diabatic wind profile correction factor for   7,03,3
                 profiles above the canopy

    PSIMS        Diabatic wind profile correction factor for   7,03,4
                 profiles above the soil surface 

    ZDLA         Height divided by Monin-Obukov length above   7,03,5
                 the canopy

    ZDLS         Height divided by Monin-Obukov length above   7,03,6
                 the soil

    CPESTR       Canopy storage of water vapor                 7,03,7

    DRGRAV       Drainage from the soil by gravity             7,03,8

    WCPYS        Soil surface evaporation calculated from      7,03,9
                 soil water flux equations (Wm-2)

    RNCPY        Net radiation above the canopy (Wm-2)         7,04,1

    HCPY         Net sensible heat from the canopy and soil    7,04,2
                 to the atmosphere (Wm-2)

    ECPY         Evapotranspiration from the canopy (Wm-2)     7,04,3

    CONCPY       Condensation on the canopy (gm-2)             7,04,4

    MMH20        Total water stored in root zone (mm)          7,04,5

    DELH20       Change in soil water storage from the last    7,04,6
                 time step (mm)

    DRAIN        Drainage out of the bottom of the root zone   7,04,7
                 during this timestep (mm)

    FILT         Infiltration into the soil during this        7,04,8
                 timestep (mm)

    ETMIC        Canopy evapotranspiration during this         7,04,9
                 timestep (microns)

    TPRECP       Precip. transmitted through the canopy (mm/s) 7,05,1

    PREINT       Precipitation intercepted by canopy (mm)      7,05,2

    PEVMIC       Intercepted precipitation that evaporated     7,05,3
                 during this time step (microns)

    STEM         Stem flow of intercepted precipitation (mm)   7,05,4

    TRMIC        Canopy transpiration during this timestep     7,05,5
                 (microns)
    
    EVSMIC       Soil evaporation during this timestep         7,05,6
                 (microns)
    
    DPEVP        Droplet evaporation from sprinkler            7,05,7
                 irrigation.        

    PSCPYL       Photosynthetic rate of canopy on leaf area    7,05,8
                 basis (umol CO2 m-2 s-1)

    The interpretation of 7,7-9 depends on the value of FACTIR:
                 
    FACTIR is the factor that the thermal sky flux calc'd from subroutine
    skyir is multiplied by to get sky thermal in band of interest.
    
    FACTIR=0.99 means calculate the average temps (as would be measured 
    by a broad band hemispherical radiometer). 

    FACTIR=0.50 means calculate the directional temps (as would be 
    measured by a directional narrow band IR thermometer). 

    If FACTIR=0.99:
                 
    T1		 Average vegetation kinetic temperature	       7,07,1		 
    
    T4		 Hemispherical infrared temperature            7,07,2  
    
    If FACTIR=0.50:                   

    VIEWAZ       Azimuth view angle usually wrt north          7,07-9,1
                     
    TB4 NN       Average directional kinetic temperature.
                 Directional temp. of cpy if all elements      7,07,2-9
		 were black bodies at view zenith angle=NN   
		 deg. Weighting is T^4. (Also known at T2.)

    EM4 NN       Apparent directional infrared temperature.
    		 Directional apparent emissivity of cpy based  7,08,2-9 
		 on apparent IR temp simulated for infrared 
		 thermometer and black body canopy temp.  
		 This vbl is emiscp and thus depends on sky 
		 irradiation.  To get true cpy emissivity  
		 factir must be 0.  if factir=1 then reflected 
		 sky rad makes apparent emissivity look close 
		 to 1.0. 

    TA4 NN       Directional brightness temperature of canopy  7,09,2-9
		 simulated for an infrared thermometer 
		 sensitive to some waveband like 8-14 um.  
		 This includes emission from the canopy and 
		 reflected sky irradiation so typically 
		 factir=.5 for a clear sky.  If factir is set 
		 to zero to get true emissivity then this temp 
		 will be below what an infrared thermometer 
		 would measure because the reflected sky 
		 irradiation will be zero. (Also known as T5.)

    VIS I        Visible radiation reflectance factors for     7,11,1-9
                 calculation of bidirectional fluxes

    NIR I        NIR reflectance factors for calculation of    7,12,1-9
                 bidirectional fluxes

    PSIXY        Root xylem water potential (bars)             7,13,1

    PSITOP       Leaf water potential (bars)                   7,13,2

    PSISUM       Root length weighted soil water potential     7,13,3
                 (bars)

    DELD         Convergence factor for soil water equations   7,13,4

    WSMM         Soil surface evaporation (mm)                 7,13,5

    RSCPY        Effective canopy stomatal resistance          7,13,6
                 (sm-1)

    RVCPY        Canopy resistance to water vapor              7,13,7
                 (sm-1)

    RHCPY        Canopy resistance to heat flow                7,13,8
                 (sm-1)

    QNOITR       Number of iterations in leaf energy budget    7,14,1

    QITER2       Iterations over canopy temperature and V.P.   7,14,2
                 profiles to get convergence in main program

    QITER3       Iterations over soil evap.to merge soil water 7,14,3
                 flow eq. & canopy vapor eq. in PROFL2 
                 subroutine

    QITERW       Iterations to match soil sfc W.C. with V.P.   7,14,4
                 at soil surface through log relationship

    QLOOPE       Iterations over canopy water vapor flow       7,14,5
                 equations alone in Newton-Raphson

    QITERH       Iterations over water stress in canopy        7,14,6

    QLOOPT       Iterations over soil heat flow equations      7,14.7  
                 alone in Newton-Raphson

    QLOOPW       Iterations over soil water flow equations     7,14,8
                 by Newton-Raphson in PROFL2

    ETPM2        Canopy evapotranspiration by Penman-Monteith  7,15,1
                 equations neglecting soil heat flux

    RNETDV       Net radiation divergence in canopy   	       7,15,2
                 

    ETPM         Canopy ET by Penman-Monteith equ. using soil  7,15,3
                 heat flux from CUPID

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

    RHEAT        Aerodynamic resistance to heat and water      7,15,5
                 flow from diabatic wind profile equation

    SENS         Canopy sensible heat flux by Penman-Monteith  7,15,6
                 equations using soil heat flux from CUPID

    USTR         Friction velocity used in Penman-Motheith     7,15,7
                 eq. for the calculation of aerodynamic 
                 resistance to heat and water flow

    TCROP        Crop temperature from P.-M. equations         7,15,8

    ZDL2         Height divided by Monin-Obukov length used    7,15,9
                 for diabatic profile corrections in 
                 Penman-Monteith equations

    PARI         Intercepted PAR by the canopy (Wm-2)          7,16,1

    SOLRI        Intercepted solar radiation by canopy (Wm-2)  7,16,2

    PSCPYG       Photosynthetic rate of canopy on ground area  7,16,3
                 basis (umol CO2 m-2s-1)

    APAR         Absorbed PAR by the canopy (Wm-2)             7,16,4

    ASOLR        Absorbed solar radiation by canopy (Wm-2)     7,16,5

    ELUAP        Efficiency of light use based on APAR         7,16,6
                 (g CO2 MJ-1)

    ELUAS        Efficiency of light use based on ASOLR        7,16,7
                 (g CO2 MJ-1)

    ELUIP        Efficiency of light use based on PARI         7,16,8
                 (g CO2 MJ-1)

    ELUIS        Efficiency of light use based on SOLRI        7,16,9
                 (g CO2 MJ-1)  

    TIRRIG       Irrigation water temperature at canopy top    7,17,1

    ETTOT        Total droplet evaporation (+) or condensa-    7,17,2
                 tion (-) above the canopy (Wm-2)

    QTOT         Total droplet sensible heat exchange above    7,17,3
                 the canopy (+ droplet cooling) (Wm-2)
                 
    DPSTOR	 Heat storage in droplets (proportional to     7,17,4
                 TDROPI-TIRRIG)
    
    TIRETQ						       7,17,5

    RDCPY        Canopy dark respiration (umol co2 m-2 s-1)    7,18,1

    WTSFC        Volumetric water content of top soil layer    7,18,2

    CO2SOL       Soil surface CO2 flux (umol CO2 m-2 s-1)      7,18,3
    
    RMSTEM       Maintenance respiration of stem.	       7,19,1
                 (micromole m-2 ground area s-1)    

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

    RMROOT       Maintenance respiration of root.	       7,19,3
                 (micromole m-2 ground area s-1)     

    RMGR         Maintenance respiration of grain.  	       7,19,4
                 (micromole m-2 ground area s-1)   

    RMTOT        Maintenance respiration of total plant.       7,19,5
                 (micromole m-2 ground area s-1)

    RGTOT        Growth respiration of total.  		       7,19,6  

    PSMRM        Photosyntheis minus maintenance respiration   7,19,7  
    
    DMINC        Dry matter increment 			       7,19,8
                 (mg of dry weight m-2s-1)   

    TROOTM       Root-weighted mean temperature of soil (C).   7,19,9  
                                                                                
    DPVMMH       Droplet evaporation in mm/hr                  7,20,1

    PEVMMH       Interception evaporation in mm.hr             7,20,2

    TRMMH        Canopy transpiration in mm/hr		       7,20,3

    EVSMMH       Soil evap in mm/hr                            7,20,4

    ETMMH        Canopy + soil evapotranspiration in mm/hr     7,20,5

    ECPYDV       Canopy only evapotranspiration in mm/hr       7,20,6
    
    WTIRSL						       7,21,1
    
    TSFC						       7,21,2
    
    TLFMAX	 Maximum allowable leaf temperature.	       7,21,3
    
    TLFMIN	 Minimum allowable leaf temperature.	       7,21,4
    
    TIRCPY	 Average temperature of canopy for sunlit      7,21,5
                 leaves.	       
    
    TIRHOT						       7,21,6
    
    TAPHOT						       7,21,7
    
    EMIHOT						       7,21,8
    
    TAREOM	 Canopy aerodynamic temperature using 	       7,22,1
                 the zero-plane displacement height for
                 momentum (C).

    TAEROH       Canopy aerodynamic temperature using	       7,22,2
                 Z  for heat (C).
                  h    
                  
    SAUER						       7,22,3 
    
    TANNER						       7,22,4
    
    PSS1A						       7,23,1  
    
    PSS1B						       7,23,2 
     
    PSS2A						       7,23,3 
     
    PSS2B						       7,23,4 
     
    EVLF						       7,23,5 
    
    QLF 	 Sensible heat flux from leaf surface??	       7,23,6 
     
    RNLF						       7,23,7 
     
    LOOP						       7,23,8  
     
    RNLF2						       7,23,9 
    
    PSS3A						       7,24,1  
    
    RDLF						       7,24,2  
    
    PSS4						       7,24,3 
     
    ZDL 						       7,24,4 
     
    PSIH						       7,24,5 
       
    PCE 	 					       7,24,6 
     
    PCEA						       7,24,7 
      
    BLRLF						       7,24,8 
     
    RSLF						       7,24,9 
    
    TEMAIR       Air temperature (C)                           7,25,1   
    
    TDISP	 Temperature at the displacement height??      7,25,2
    
    TLF2						       7,25,3
    
    VPAIR        Air vapor pressure (mb)                       7,25,4
    
    EDISP	 Vapor pressure at displacement height??       7,25,5
    
    ELF		 Vapor pressure at leaf surface		       7,25,6
       
    CAIR         Ambient atmospheric CO2 concentration (ppm)   7,25,7
    
    CDISP	 C02 concentration at displacement height      7,25,8
    
    CLF							       7,25,9
    
    DRAIN5	 Drainage from 5th soil node up from jzbot     7,26,1
    
    ALBEDO						       7,27,1
    
    VISRFL						       7,27,2
    
    NIRRFL						       7,27,3
    
    ALB2						       7,27,4
    
    VISR2						       7,27,5
    
    NIRR2						       7,27,6
    
    
    
    
    
    		***** HOURLY SUMMARY BY LAYER AND LEAF ANGLE *****
    		
                                                        
    DSTNET       Net all-wave radiation in angle class I and   8,01,1
                 layer J (Wm-2)

    DSTVIS       Incident visible radiation in angle class I   8,01,2
                 and layer J (Wm-2)

    DSTNIR       Incident NIR in angle class I and layer J     8,01,3
                 (Wm-2)

    DSTTHR       Incident thermal radiation in angle class I   8,01,4
                 and layer J (Wm-2)

    PRAREA       Percent of leaf area in angle class I and     8,01,5
                 layer J

    PS1          Net photosynthetic rate of a leaf in layer I  8,01,6
                 and angle class I  (1E-6 mol CO2 m-2 s-1)

    TEMPLF       Leaf temperature in angle class I and         8,02,1
                 layer J (C)

    DELT         Temperature difference between leaf in angle  8,02,2
                 class I and layer J and air in layer J (C)

    EVAP         Leaf ET in angle class I and layer J (mm)     8,02,3

    HEAT         Leaf sensible heat based on leaf area in      8,02,4
                 angle class I and layer J (Wm-2)

    GEVAP        Leaf ET based on ground area for angle class  8,02,5
                 I and layer J

    GHEAT        Leaf sensible heat based on ground area for   8,02,6
                 angle class I and layer J

    RSLEAF       Stomatal resistance for angle class I and     8,02,7
                 layer J (s/m)
                 
    RSNOVP       Stomatal resistance with no direct vapor      8,02,8
                 pressure deficit effect (sm-1).
                 
    CSLEAF       Leaf stomatal and boundary layer conductance  8,03,1
                 for angle class I and layer J (mol-1 H2O m2s)

    RSLEAF       Leaf stomatal resistance for angle class I    8,03,2
                 and layer j (mol H2O m-2 s-1 )

    PSLEAF       Leaf photosynthetic rate for angle class I    8,03,3
                 and layer J (1E-6 mol m-2s-1)

    CILEAF       Intercellular CO2 concentration of leaves in  8,03,4
                 angle class I and layer J (1E-6 mol/mol)

    RDRK         Dark respiration (umol CO2 m-2 s-1)           8,03,5

    RHLEAF       Leaf boundary layer resistance (s m-1)        8,03,6

    ABSPAR       Absorbed PAR by leaf(umol photon m-2 s-1)     8,03,7

    CSFCLF       CO2 conc at surface of leaf (umol/mol)        8,03,8

    HSLEAF       Rel Hum at the surface of the leaf            8,03,9
    
    
    
    
    		***** HOURLY BOUNDARY LAYER SUMMARY *****
    			 (Cupid-PBL version only)
    		

    PRESS(K)     Atmospheric pressure at level k (mbar)	       9,01,1
    
    TEMP(K)      Air temperature at level k (C)                9,01,2
    
    POTMP(K)     Potential temperature at level k (C)          9,01,3
    
    MIXRAT(K)    Mixing ratio at level k (g/kg)		       9,01,4
    
    WINDU(K)	 U Wind speed component at level k (m/s)       9,01,5
    
    WINDV(K)	 V Wind speed component at level k (m/s)       9,01,6
    
    WIND(K)	 Wind speed at level k (m/s)		       9,01,7
    
    GPHITE(K)    Geopotential height at level k (m)	       9,02,1
    
    EKMH(K)	 Exchange coefficient for momentum	       9,02,2
    
    EKHH(K)      Exchange coefficient for heat                 9,02,3
    
    STABRAT(K)   Stability criterion: (dz^2/dt)/max(ekmh,ekhh) 9,02,4
    
    RI(K)	 Richardson's number