CUPID: SETTING UP AN INPUT FILE


                                                                    

	This documentation file is meant to assist someone in preparing
an input file for the Cupid model.  SECTION A is from file CUPINPUT.ASC
of the general documentation and provides a general over view. SECTION B
is from file IVRLIST.ASC and is a list of variables in the input. 
SECTION C is an example input file.




A. Overview Of A Cupid Input File

  1) General

  The Cupid input file consists of five types of information, necessary to 
  run and interpret the simulation of a plant-environment system.

     1.  The user  must decide on the information he or she wants from the 
         simulation. Cupid has the  capacity  to  provide more information 
         than could be useful for a single purpose. Hence, the Cupid input 
         requires a specification of the desired output. The first lines of 
         the input file represent code numbers that determine the Cupid 
         output.

     2.  The following lines of input are constant values specifying fixed
         properties of the plant-soil-atmosphere system. 
   
     3.  This is followed by a block of values that represent daily adjusted 
         properties of the system.

     4.  The next line specifies initial soil water status.

     5.  The rest is a block of environmental data as they are taken from
         weather stations. The input is typically hourly, but may also be
         in timesteps of varying length. Each line represents one time step. 


  2) The Code Numbers (Revised 930318)

  The code numbers are two digit integers. Each pair of numbers represents one 
  group of Cupid output variables. A separate paragraph has been written on the 
  code, so the code itself will not be discussed here. Also, the list of Cupid 
  output variables will be helpful to set up the code for a required output.
  The first integer of the input file is the number of code number pairs that 
  follow. Each pair of code numbers given here will result in the output of a 
  specific group of variables for every hour.

  The next line contains a flag ISPECL. If set to 1, the user has the option  
  to specify individual time steps for which additional output is wanted
  such as vertical profiles of quantities. The 
  following line specifies the number of code number pairs and the code 
  numbers for special hour output.

  The following line specifies the time steps for which the special output 
  applies. The time steps have to be entered individually and are led by the 
  number of special hours. These time step numbers are not identical with the 
  actual hour of the day. They are merely the indices of time step subscripted 
  variables. 



  3)  Fixed Header Data    

  The following parameters are best to be studied from the list of input 
  variables. They consist of various values to specify geographical location,
  soil properties, sprinkler set-up, plant properties, sensor locations
  and desired aspects of the structure and accuracy of the simulation.
  Refer to the list of input variables for more detailed discussion.


  4) Daily Parameters (revised 930318)

  The daily input variables are preceded by one line that indicates the number 
  of daily inputs that follow. It is, however, not necessary to specify values 
  for each day of the simulation. Daily parameters that are not given are 
  linearly interpolated between the preceding and the following value. The 
  daily parameter consists of the day number for which the values are valid, 
  the total plant LAI, the plant height, the rooting depth coefficient, the 
  relative height of the lowest leaf and the relative height of maximal
  leaf are density, the temperature and water content at the lower
  soil boundary.  The first day of this daily input must correspond with
  the first day of the input weather data.  The program does not check
  for this.

  5) Initial condition

  Cupid requires the initial condition of the soil water profile, which is 
  entered immediately after the daily input values. The numbers correspond to 
  the soil depths as specified above.


  6) Weather Station Data (Revised 930318)

  The time step data are preceded by the number of days for which the 
  simulation will run.  It may be practical to run fewer days than are actually 
  in the file. If all days are run it is necessary that the last line of  
  the time step data falls into the day following the last day specified. (For 
  a discussion of this refer to the paragraph on Cupid's time structure).
  The line preceding the time step block has two other values that belong in
  the category of Fixed Header Information: RATIO-DAY and RATIO-NITE, the
  ratios of potential to actual radiation, used only in case that the
  radiation flux data are not available.

  The time step data that follow consist of YEAR, DAY OF YEAR, LOCAL TIME,
  WIND SPEED (m/s), RADIATION FLUX (W/m2), AIR TEMPERATURE (C), AIR WATER
  VAPOR PRESSURE (mbar), PRECIPITATION (mm/time step) and IRRCHK, 
  an integer flag that determines the form of precipitation.

  Radiation input data is somewhat complicated because various
  combinations of PAR and NIR, direct and diffuse are possible.  The 5
  columns for radiation are TOTAL PAR, TOTAL NIR, PAR BEAM FRACTION
  NIR BEAM FRACTION, SKY THERMAL RADIATION.  Usually only total solar 
  radiation is available and then the input should be as follows:

           -1, Total Solar Radiation (W m-2), -1, -1, -1

  A "-1" indicates that those values must be calculated by the model.

  Precipitation codes: 0 = no precipitation
		       1 = overhead sprinkler irrigation
		       2 = furrow irrigation
		       3 = rainfall


B. List Of Input Variables (Revised930318)



       LIST OF CUPID INPUT VARIABLES IN THE ORDER OF HOW THEY ENTER THE
       PROGRAM (103 records):

       NAME       DEFINITION                                   (M,NN)  REF.

    Codes Controlling Output

       IPNT       Code numbers for output format (see key)

       JPNT       Code numbers for output format (see key)

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

       KPNT       Code number for output of vertical profiles
                  (see key)

       LPNT       Code number for output of vertical profiles
                  (see key)

       ISPN       Hours of special profile output

       MDAYIN     No. of day numbers input

       MDAY(LL)	  Day #'s for which profile output desired

    Site Factors

       XLAT       Latitude                                          1,01,5

       XLONG      Longitude                                         1,01,6

       STDLNG     Standard longitude (90 deg. central time zone)    1,01,7

    Soil Factors

       CLODSZ     Soil clod size index (1 m)                        

       NDSOIL     Number of depths in soil                          1,05,4
                  (number of soil layers)(13)

       REFDPT     Reference depth (2 m)                             1,03,6

       BD         Bulk density (1.4 Mg/m3)

       SANDFC     Sand fraction (0.7)

       SILTFC     Silt fraction (0.2)

       CLAYFC     Clay fraction (0.1)

       QRTZFC     Quartz fraction in fractions of total weight (0.2)

       DPE        Air entry potential (J/kg) (optional). If 0 is   
		  entered it is calculated by G. Cambbell's
		  equations.(-2.5J/kg)

       DBX        Exponent in the soil moisture release curve (optional).
		   (4.5)

       DAKS       Saturating water conductivity (optional)(0.7E-3kg-s/m3
		  .0007kg-s/m3 * 9.8m/s2 /1000kg/m3 = 6.86e-6m/s
		  =.00686mm/s=24.7mm/hour=1 inch/hour)

       ZSOIL      Depth of nodes in soil (m)                       

       DELD       Convergence criteria of diff. equation            7,13,4
                  solution (2.0E-5)

       DRHSFC     Value to match soil water flow equations          1,05,6
                  with canopy vapor flow equations (.01)

       RHSLOP     Value to match soil water flow equations
                  with canopy vapor flow equations (.05)

    Irrigation System Factors

       CTWA       ..GTWA parameters for irrigation water
                  temperature calculation (see parameter library)

       CETSUM     ..HETSUM parameters for droplet evaporation
                  calculation (see parameter library)

       CQSUM      ..GQSUM parameters for droplet heat exchange
                  calculation (see parameter library)

       CETFIX     .. FETFIX parameter for droplet evaporation
                  calculation (see parameter library)

       CQFIX      ..FQFIX parameter for droplet heat exchange
                  calculation (see parameter library)

       ZMAX       Maximum height of droplets over sprinkler (m)

       HTSPK      Height of sprinkler (m)                           1,09,1

       ABSPK      Number of sublayers above sprinkler in a
                  theoretical source profile

       BSPK       Number of sublayers below sprinkler in a
                  theoretical source profile

       SPKPRS     Sprinkler pressure (PSI)                          1,09,3

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

    Program Control Codes

       IBIDIR Code:       0:  no bidirectional radiation fluxes,
                          1:  bidirectional radiation fluxes


       ITASSL Code:       0:  assuming top canopy layer has no tassels
                          1:  assuming tassels present
    
    Canopy View Characteristics

       IPP       0: Calc view angles from viewzn and viewaz below
		 1: Calc view angles in principal plane of sun only

       NVZENP    No. of view zenith angles in principal plane

       VWZENP    View zenith angles in principal plane

       NOZENV    No. of view zenith angles if IPP=0

       VIEWZN    View zenith angles if IPP=0

       NOAZMV    No. of view azimuth angles if IPP=0

       VIEWAZ    View azimuth angles if IPP=0

    Canopy Factors

       ITOT       Number of leaf angle classes (9)                  1,01,1

       JMAX       Maximum number of canopy layers (19)              1,01,2

       JMIN       Minimum number of canopy layers (2)               1,01,3

       DFMIN      Thinnest leaf layer in terms of LAI units (0.2)   1,01,4

       CLUMP      Clumping factor for canopy structure.             5,01,6
                  Set to 1 for no clumping. (1.0)

       EMIS       Thermal emissivity of the leaves (0.95)           1,01,8

       EMISOL     Thermal emissivity of the soil (0.93)             1,01,9

       FACTIR     Factor to multiply total sky thermal flux by 
		  to get downward sky flux density in wavelength
		  interval sensed by IRT. If this is 0.0 then
		  the apparent emissivity is actual emissivity
		  of the canopy(emiscp). In 8-14 um band if this is 0.5
		  then a reasonable apparent sky infrared thermometer
		  temp. will be simulated for a clear sky.

       RSOIL      Soil Reflectance: (1): visible,                   1,02,1-3
                  (2): near infrared, (3): thermal

       RLEAF      Leaf reflectance: (1): visible,                   1,02,4-6
                  (2): near infrared, (3): thermal

       TLEAF      Leaf transmittance, (1): visible,                 1,02,7-9
                  (2) near infrared, (3) thermal

       AVISR      Non Lambertian coeff for leaf reflectance

       AVIST      NonLambertian coeff for leaf transmittance

       ISHPERE	  0: Use leaf andle dist input below
		  1: Spherical leaf angle dist, overrides and other
		     input.
       NALPHA	  Number of leaf angle classes input (9)

       GR(ALPHA)  Frac of leaf area in various leaf angle classes  1,04,1-9
		  if ISPHERE=0-leaf angle dist.

       IMUNU      0: Use GR(ALPHA) input above and calc. meu & neu

		  1: meu & neu input directly to calc leaf angle dist

		  2: mean and std. dev. of leaf angle dist. input

       XMEU       imunu=1: value of meu
		  imunu=2: mean angle of leaf angle dist(0=hor leaves)

       XNEU       imunu=1: value of neu
		  imunu=2: std. dev. of leaf angle dist.

                  Canopy Type	meu	neu (See Subroutine LAD)
		  planophile	2.770	1.172
		  erectophile	1.172	2.770
		  plagiophile	3.326	3.326
		  extremophile	0.433	0.433
		  uniform	1.000	1.000
		  spherical	1.101	1.930

       IMUNUA     0: Leaves symmetric about azimuth (meu=neu=1.0)
		  1: Mean and std. dev. of azimuth angle dist

       BETA0      Angle of minimum in leaf area if azimuth dist is
		  asymmetric (37.62)

       NBETA      No. of azimuth angle classes (50)

       ROWSPC     Row spacing (0.76m)                               1,05,1

       PLSPC      Plant spacing (0.15m)                             1,05,2

       DMAX       Maximum width of the leaf to calc leaf B.L.R.(0.1m)

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

       PLTPM2     Number of plants per square meter

       REFHTW     Reference height for wind (= temp and vp)         1,03,7

       ZDHCR      =.2, critical layer height/canopy height

       AMFULL	  Canopy velocity profile coefficient (1.5)

       NLABCY     Number of layers above the canopy

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


       NLBCPY     Number of layers below the canopy

    Parameters for describing stom resis--no longer used

       ZBCPY      Height of layers below the canopy (m)

       RCUT20     Light response characteristics of stomata

       RSMIN      Light response characteristics of stomata

       ANSTOM     Light response characteristics of stomata

       RADN       Light response characteristics of stomata

       TRSOPT     Temp. response characteristics of stomata

       TRSMAX     Temp. response characteristics of stomata

       TRSMIN     Temp. response characteristics of stomata

       RSEXP      Temp. response characteristics of stomata

       RSM        Temp. response characteristics of stomata

    End obsolete stomatal resistance parameters

       PSI1       Leaf water pot. at which stomata begin to close   2,04,4
		  (-15 bar)

       PSI2       Leaf water pot. at which stomata closed (-20 bar) 2,04,5

       IC3C4      Code for C3 plant (0), C4 plant (1)

		   0: C3 Farquhar
		   1: C4 Corn		Berry model w/ Index for stom
		   2: C4 Big Blue Stem  Berry model w/ Index for stom
		   3: C3 Soybean        Berry model w/ Index for stom
		  10: C3 Farquhar
		  11: C4 Corn           Farquhar model stom use D1 & D2
		  12: C4 Big Blue Stem  Farquhar model stom use D1 & D2
		  13: C3 Soybean	Farquhar model stom use D1 & D2

       CAIR       Ambient atmospheric CO2 concentration (330 vpm)

       FS         Ratio of stomatal resistance on top to
                  bottom of leaf (1.0)

       CUCOND     Cuticular conductivity to water vapor
                  (.0137 mol/m2/sec)

       D1         VPD for start of stomatal closure (0.0 mb)

       bkv        Slope of stom conduct. vs vpd func (4.0)

       CONMIN     Used with D1 and vpd but not used now
                                       
       RXCHAM     Boundary layer resistance of one leaf side in
                  chamber for measuring P.S. (30. s/m)

       NOTEMP     0: Normal temperature response for photosyn       1,03,8
		  1: Photosyn. independent of temp. and set to
		     value at temp TFIX. Dark Resp. still temp. dep.
       TFIX       Photosyn. evaluated at this temp for all leaves   1,03,9
		  (35C)

       RROOT      Hydraulic root resistance (3.E6 m4kg-1s-1)        1,03,7

       FRSTEM     Fraction of intercepted precip. that runs         1,05,7
                  down the stem (.5)

       FRWTMX     Max. fraction of leaf area wetted by precip.      1,05,8
		  (0.2)

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

       INDAYS     Number of days for which daily input values
                  are listed

       LDAYIN     Day of year number of daily values to follow

       TLAI       Total leaf area index at day LDAYIN (3.0)         2,01,5

       HT         Height of canopy at day LDAYIN (2.5 m)            2,01,4

       AROTIN     Rooting depth coefficient at day LDAYIN (2.4)     2,01,8

       ZLDH       Relative height of the lowest green leaf at
                  day LDAYIN (0.2)

       ZMDH       Relative height of the maximum leaf area density
                  at day LDAYIN (0.55)

       TEMSOL     Soil temperature at soil boundary condition
                  at day LDAYIN (24.0 C)

       WATSOL     Volumetric water content at soil boundary condition
                  at day LDAYIN (0.13)

       WATER      Initial water content profile (fraction by volume)

       NODAYS     Number of days for the simulation

       RATIOD     Ratio of actual/potential solar radiation
                  during the day

       RATION     Ratio of actual/potential solar radiation
                  used at night to estimate cloud cover (1.0)

    Weather Data

       IYEAR      Year                                              2,01,1

       ICUMDY     Cumulative day                                    2,01,2

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

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

       RADTOP     Solar radiation above canopy, 3 values for VIS.   4,01,5-7
                  NIR and thermal (RADTP1,2,3 in output)


       FBEAM      Fraction of total radiation above canopy,         4,01,8-9
                  2 values for VIS. and NIR (FBEAM1,2 in output)

       TEMAIR     Air temperature  (C)                              4,02,1

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

       PRECIP     Precipitation (mm)                                4,02,4

       IRRCHK     Flag to specify type of precipitation.
                  0: no precipitation
                  1: sprinkler 
                  2: irrigation below canopy (drop tubes)
                  3: rain

       Note: Last line of file must have day number on it that is one
	     larger than the last day of the simulation.

C. Example Input File called CUPID6.KSUIN81


file cupid6.ksuin81  day 201-208 1981  
Data from Kansas State University 
48 1 1  1 2  1 3  1 4  1 5  1 6  1 7  1 8  1 9  1 10                            
2 1  2 2  2 3  2 4  2 5  2 6  2 7  2 8  2 9  2 10  2 11  2 12  2 13  2 14
4 1  4 2  4 3
7 1  7 2  7 3  7 4  7 5  7 7  7 8  7 9  7 10  7 11  7 12  7 13  7 14  7 15
7 16  7 17  7 18  7 19  7 20  7 21  7 22
1      (ISPECL)                                                                 
10  5 1  5 2  5 3  5 4  5 5  5 6  5 7  8 1  8 2  8 3                            
6 6 7 11 13 15 20 (NN2PNT,(ISPN(L),L=1,NN2PNT)                                  
1  1                              (MDAYIN,MDAY(LL),LL=1,MDAYIN)                 
43 103 90   (LAT  LONG  STD LONG)          
1. 13  2.0       ( CLODSZ NDSOIL REFDPT )                                       
1.4 .7 .2 .1 .2 -2.5    4.5 .7E-3  (BD,SAND,SILT,CLAY,QRTZ,PE,BX,DAKS)          
0. .02 .06 .21 .37 .52 .67 .82 .97 1.12 1.27 1.42 2.0   (ZSOIL)                 
2.E-5 0.01 0.05       (DELD DRHSFC RHSLOP)                                      
4.019876  0.257870 0.028933  0.001238 0.000077             (TIRRIG)             
10882.7266 287.241455 -3.0882 -75.019791 -2.499599 -0.1965 (ETSUM)              
5619.871090 -1266.63257 47.500488  3.507599  0.083400      (QSUM)               
173.876083  3.428699 -0.770100 -0.105900                   (ETFIX)              
-68.959198 -18.997482  0.410400  0.192600                  (QFIX)               
2.8 4.5 28. 24. 53. 27. 13.5 (ZMAX,HTSPK,ABSPK,BSPK,SPKPRS,BDYANG,TDROPI        
0 0 (ibidir itassl)
1  ipp, =1 only reflectance in principle plane is calculated
4  nvzenp: number of view zenith angles in principle plane
0.  20.  35.  50.  vwzenp(i),i=1,nvzenp: view zenith angles
6  0.5 15. 30. 45. 60. 75.    (NOZENV,VIEWZN(L),L=1,NOZENV)
8  0. 45. 90. 135. 180. 225. 270. 315.   (NOAZMV,VIEWAZ(L),L=1,NOAZMV)
9 19  2  .2  1.0  3 (ITOT JMAX JMIN DFMIN CLUMP kmax)  
.95 .93  .2     (EMIS EMISOL factir)   
.1 .3 .07 .09 .39  .05  .04  .45  .0   (RSOIL RLEAF  TLEAF CORN)   
0.        0.         8.16 0.       0.      41.18  (AVISR,BVISR,CVISR,ANIRR,BNIRR
0.        0.      3.22  0.       0.       57.23   (AVIST,BVIST,CVIST,ANIRT,BNIRT
1 9 .07 .16 .29 .30 .18 .18 .18 .18 .18     (ISPHER NALPHA GR(ALPHA) )
1 1.172  2.770 (IMUNU  XMEU or angmn   XNEU or angsd     IMUNU=1 or 2)
0  180.   15.  37.62 50  (imunua xmeuaz/angamn xneuaz/angasd beta0 nbeta)
.75  .15 0.10   (ROWSPC  PLSPC DMAX)                                    
3 7    (ILAIUP PLTPM2)                                                          
4. .20 1.5  (REFHTW ZDHCR amfull)                                  
5  4. 3.59 3.18 2.77 2.36            (NLABCY,(Z(L),L=1,NLABCY))                 
3  .28 .195 .099                  (NLBCPY,(Z(L),L=1,NLBCPY))                    
3000. 100. 3. 100.  (RS VS LIGHT   RCUT20 RSMIN ANSTOM RADN)                    
30. 45. -2. 6. 5.     (RS VS TEMP   TRSOPT TRSMAX TRSMIN RSEXP RSM)             
-15.0 -20.0  .5  (PSI1 PSI2  cidcai, -1 to use cidca in code)  
3  330.  1.  .0137   0.  4.0 .08 30. (IC3C4 CAIR FS CUCOND D1 bkv conmin RXCHAM
0 35.            (NOTEMP=1 IF NO TEMP EFFECT IS PS AND KS ELSE=0, TFIX)         
3.E6      (RROOT)                                                               
0.5  0.2 0.15      (FRSTEM FRWTMX PINTMX)                                       
8                  (INDAYS)
201 3.0  2.5 2.4 .2  .55 24. .13   (LDAYIN,XTLAI,XHT,XAROT,XZLDH,XZMDH,         
202 3.0  2.5 2.4 .20 .55 24. .13    TSOLBC,WSOLBC)                              
203 3.0  2.5 2.4 .20 .55 24. .13                                                
204 3.0  2.5 2.4 .20 .55 24. .13                                                
205 3.0  2.5 2.4 .20 .55 24. .13                                                
206 3.0  2.5 2.4 .20 .55 24. .13                                                
207 3.0  2.5 2.4 .20 .55 24. .13                                                
208 3.0  2.5 2.4 .20 .55 24. .13                                                
.03 .03 .03 .24 .26 .22 .2 .2 .2 .19 .17 .15 .12     (WATER)                    
  1  1.0 1.0      ( NODAYS  RATIO-DAY RATIO-NITE)                               
 1981 201  0.9 2.70   -1.0    0.0   -1.0   -1.0   -1.0  22.9  20.8   0.0 0      
 1981 201  1.9 2.00   -1.0    0.0   -1.0   -1.0   -1.0  23.3  20.4   0.0 0      
 1981 201  2.9 3.90   -1.0    0.0   -1.0   -1.0   -1.0  23.1  20.0   0.0 0      
 1981 201  3.9 4.20   -1.0    0.0   -1.0   -1.0   -1.0  22.1  18.8   0.0 0      
 1981 201  4.9 1.10   -1.0    0.0   -1.0   -1.0   -1.0  19.1  18.6   0.0 0      
 1981 201  5.9 0.90   -1.0    2.3   -1.0   -1.0   -1.0  18.2  19.2   0.0 0      
 1981 201  6.9 1.20   -1.0   88.3   -1.0   -1.0   -1.0  19.8  19.9   0.0 0      
 1981 201  7.9 2.10   -1.0  264.9   -1.0   -1.0   -1.0  23.6  20.2   0.0 0      
 1981 201  8.9 3.60   -1.0  464.7   -1.0   -1.0   -1.0  26.7  21.5   0.0 0      
 1981 201  9.9 3.10   -1.0  651.7   -1.0   -1.0   -1.0  29.0  22.5   0.0 0      
 1981 201 10.9 1.40   -1.0  809.7   -1.0   -1.0   -1.0  30.6  23.5   0.0 0      
 1981 201 11.9 1.00   -1.0  907.3   -1.0   -1.0   -1.0  31.8  24.5   0.0 0      
 1981 201 12.9 1.20   -1.0  973.5   -1.0   -1.0   -1.0  33.5  24.5   0.0 0      
 1981 201 13.9 1.10   -1.0  971.2   -1.0   -1.0   -1.0  34.2  24.5   0.0 0      
 1981 201 14.9 1.40   -1.0  903.8   -1.0   -1.0   -1.0  34.6  24.5   0.0 0      
 1981 201 15.9 1.30   -1.0  789.9   -1.0   -1.0   -1.0  35.0  24.5   0.0 0      
 1981 201 16.9 2.20   -1.0  630.8   -1.0   -1.0   -1.0  34.5  24.5   0.0 0      
 1981 201 17.9 2.30   -1.0  472.8   -1.0   -1.0   -1.0  34.1  24.0   0.0 0      
 1981 201 18.9 1.90   -1.0  245.1   -1.0   -1.0   -1.0  33.0  23.5   0.0 0      
 1981 201 19.9 1.50   -1.0   67.4   -1.0   -1.0   -1.0  29.5  22.5   0.0 0      
 1981 201 20.9 1.70   -1.0    1.2   -1.0   -1.0   -1.0  26.6  23.0   0.0 0      
 1981 201 21.9 2.10   -1.0    0.0   -1.0   -1.0   -1.0  25.2  24.0   0.0 0      
 1981 201 22.9 2.40   -1.0    0.0   -1.0   -1.0   -1.0  25.1  25.0   0.0 0      
 1981 201 23.9 2.10   -1.0    0.0   -1.0   -1.0   -1.0  23.8  25.0   0.0 0      
 1981 202  0.9 2.40   -1.0    0.0   -1.0   -1.0   -1.0  23.3  24.0   0.0 0      
 1981 202  1.9 3.00   -1.0    0.0   -1.0   -1.0   -1.0  23.0  23.0   0.0 0      
 1981 202  2.9 3.40   -1.0    0.0   -1.0   -1.0   -1.0  23.1  22.0   0.0 0      
 1981 202  3.9 3.20   -1.0    0.0   -1.0   -1.0   -1.0  22.5  21.0   0.0 0      
 1981 202  4.9 3.20   -1.0    0.0   -1.0   -1.0   -1.0  22.3  20.5   0.0 0      
 1981 202  5.9 2.70   -1.0    1.2   -1.0   -1.0   -1.0  21.6  20.1   0.0 0      
 1981 202  6.9 2.20   -1.0   59.2   -1.0   -1.0   -1.0  21.2  20.5   0.0 0      
 1981 202  7.9 3.70   -1.0  232.3   -1.0   -1.0   -1.0  22.8  21.1   0.0 0      
 1981 202  8.9 4.30   -1.0  461.2   -1.0   -1.0   -1.0  25.8  22.5   0.0 0      
 1981 202  9.9 5.30   -1.0  655.2   -1.0   -1.0   -1.0  29.1  23.5   0.0 0      
 1981 202 10.9 3.80   -1.0  819.0   -1.0   -1.0   -1.0  31.6  24.5   0.0 0      
 1981 202 11.9 3.80   -1.0  937.5   -1.0   -1.0   -1.0  34.6  24.5   0.0 0      
 1981 202 12.9 3.10   -1.0  983.9   -1.0   -1.0   -1.0  35.8  25.0   0.0 0      
 1981 202 13.9 1.70   -1.0  983.9   -1.0   -1.0   -1.0  37.5  25.0   0.0 0      
 1981 202 14.9 1.60   -1.0  909.6   -1.0   -1.0   -1.0  37.5  25.0   0.0 0      
 1981 202 15.9 1.60   -1.0  803.9   -1.0   -1.0   -1.0  37.0  24.5   0.0 0      
 1981 202 16.9 2.00   -1.0  626.1   -1.0   -1.0   -1.0  36.0  24.0   0.0 0      
 1981 202 17.9 5.00   -1.0  435.6   -1.0   -1.0   -1.0  32.8  24.0   0.0 0      
 1981 202 18.9 4.40   -1.0  242.8   -1.0   -1.0   -1.0  31.0  24.0   0.0 0      
 1981 202 19.9 3.70   -1.0   68.5   -1.0   -1.0   -1.0  28.8  23.5   0.0 0      
 1981 202 20.9 3.00   -1.0    1.2   -1.0   -1.0   -1.0  25.6  23.5   0.0 0      
 1981 202 21.9 2.90   -1.0    0.0   -1.0   -1.0   -1.0  24.2  23.5   0.0 0      
 1981 202 22.9 2.80   -1.0    0.0   -1.0   -1.0   -1.0  22.9  23.5   0.0 0      
 1981 202 23.9 2.80   -1.0    0.0   -1.0   -1.0   -1.0  21.7  23.2   0.0 0      
 1981 203  0.9 2.90   -1.0    0.0   -1.0   -1.0   -1.0  20.1  23.2   0.0 0      
 1981 203  1.9 2.70   -1.0    0.0   -1.0   -1.0   -1.0  20.5  23.0   0.0 0      
 1981 203  2.9 2.80   -1.0    0.0   -1.0   -1.0   -1.0  19.7  22.8   0.0 0      
 1981 203  3.9 3.10   -1.0    0.0   -1.0   -1.0   -1.0  19.0  21.7   0.0 0      
 1981 203  4.9 3.10   -1.0    0.0   -1.0   -1.0   -1.0  18.8  21.3   0.0 0      
 1981 203  5.9 3.40   -1.0    2.3   -1.0   -1.0   -1.0  18.7  21.4   0.0 0      
 1981 203  6.9 3.50   -1.0   39.5   -1.0   -1.0   -1.0  18.4  20.9   0.0 0      
 1981 203  7.9 2.30   -1.0  124.3   -1.0   -1.0   -1.0  18.9  20.9  12.0 1      
 1981 203  8.9 1.50   -1.0  333.4   -1.0   -1.0   -1.0  20.5  21.1  12.1 1      
 1981 203  9.9 2.70   -1.0  375.2   -1.0   -1.0   -1.0  22.6  21.8  12.0 1      
 1981 203 10.9 3.40   -1.0  474.0   -1.0   -1.0   -1.0  22.8  22.5   0.0 0      
 1981 203 11.9 3.40   -1.0  225.4   -1.0   -1.0   -1.0  21.8  21.5   0.0 0      
 1981 203 12.9 3.20   -1.0  196.3   -1.0   -1.0   -1.0  21.8  22.0   0.0 0      
 1981 203 13.9 2.50   -1.0  262.5   -1.0   -1.0   -1.0  22.6  22.0   0.0 0      
 1981 203 14.9 2.20   -1.0  367.1   -1.0   -1.0   -1.0  23.9  22.5   0.0 0      
 1981 203 15.9 2.00   -1.0  537.9   -1.0   -1.0   -1.0  25.3  23.2   0.0 0      
 1981 203 16.9 2.50   -1.0  535.5   -1.0   -1.0   -1.0  25.8  23.5   0.0 0      
 1981 203 17.9 2.50   -1.0  357.8   -1.0   -1.0   -1.0  26.1  23.8   0.0 0      
 1981 203 18.9 1.70   -1.0   82.5   -1.0   -1.0   -1.0  25.3  23.5   0.0 0      
 1981 203 19.9 2.30   -1.0   13.9   -1.0   -1.0   -1.0  24.8  23.2   0.0 0      
 1981 203 20.9 2.80   -1.0    0.0   -1.0   -1.0   -1.0  24.2  23.2   0.0 0      
 1981 203 21.9 2.60   -1.0    0.0   -1.0   -1.0   -1.0  23.9  23.0   0.0 0      
 1981 203 22.9 1.80   -1.0    0.0   -1.0   -1.0   -1.0  22.3  22.9   0.0 0      
 1981 203 23.9 1.90   -1.0    0.0   -1.0   -1.0   -1.0  21.6  22.4   0.0 0      
 1981 204  0.9 1.00   -1.0    0.0   -1.0   -1.0   -1.0  21.4  22.5   0.0 0      
 1981 204  1.9 1.30   -1.0    0.0   -1.0   -1.0   -1.0  21.0  22.3   0.0 0      
 1981 204  2.9 1.30   -1.0    0.0   -1.0   -1.0   -1.0  20.8  22.2   0.0 0      
 1981 204  3.9 1.70   -1.0    0.0   -1.0   -1.0   -1.0  20.4  22.2   0.0 0      
 1981 204  4.9 1.50   -1.0    0.0   -1.0   -1.0   -1.0  19.9  21.6   0.0 0      
 1981 204  5.9 1.90   -1.0    2.3   -1.0   -1.0   -1.0  21.3  22.0   0.0 0      
 1981 204  6.9 2.40   -1.0   47.6   -1.0   -1.0   -1.0  21.4  22.2   0.0 0      
 1981 204  7.9 2.40   -1.0  163.8   -1.0   -1.0   -1.0  22.5  23.0   0.0 0      
 1981 204  8.9 2.20   -1.0  345.0   -1.0   -1.0   -1.0  24.9  24.0   0.0 0      
 1981 204  9.9 2.00   -1.0  534.4   -1.0   -1.0   -1.0  27.7  25.0   0.0 0      
 1981 204 10.9 2.70   -1.0  792.3   -1.0   -1.0   -1.0  29.7  25.5   0.0 0      
 1981 204 11.9 2.40   -1.0  844.5   -1.0   -1.0   -1.0  30.6  25.5   0.0 0      
 1981 204 12.9 2.00   -1.0  921.2   -1.0   -1.0   -1.0  31.6  25.5   0.0 0      
 1981 204 13.9 2.00   -1.0  781.8   -1.0   -1.0   -1.0  31.8  25.0   0.0 0      
 1981 204 14.9 2.00   -1.0  569.2   -1.0   -1.0   -1.0  32.5  24.5   0.0 0      
 1981 204 15.9 2.40   -1.0  464.7   -1.0   -1.0   -1.0  32.5  24.0   0.0 0      
 1981 204 16.9 2.10   -1.0  647.0   -1.0   -1.0   -1.0  33.4  24.0   0.0 0      
 1981 204 17.9 2.20   -1.0  252.1   -1.0   -1.0   -1.0  31.5  23.5   0.0 0      
 1981 204 18.9 2.40   -1.0  124.3   -1.0   -1.0   -1.0  30.2  23.5   0.0 0      
 1981 204 19.9 2.20   -1.0   17.4   -1.0   -1.0   -1.0  28.7  23.0   0.0 0      
 1981 204 20.9 1.50   -1.0    1.2   -1.0   -1.0   -1.0  27.2  22.0   0.0 0      
 1981 204 21.9 4.50   -1.0    0.0   -1.0   -1.0   -1.0  25.9  21.0   0.0 0      
 1981 204 22.9 1.40   -1.0    0.0   -1.0   -1.0   -1.0  26.0  21.0   0.0 0      
 1981 204 23.9 1.60   -1.0    0.0   -1.0   -1.0   -1.0  24.5  21.0   0.0 0      
 1981 205  0.9 1.60   -1.0    0.0   -1.0   -1.0   -1.0  23.7  22.0   0.0 0      
 1981 205  1.9 1.30   -1.0    0.0   -1.0   -1.0   -1.0  23.5  22.0   0.0 0      
 1981 205  2.9 1.70   -1.0    0.0   -1.0   -1.0   -1.0  24.7  22.0   0.0 0      
 1981 205  3.9 1.80   -1.0    0.0   -1.0   -1.0   -1.0  24.4  22.0   0.0 0      
 1981 205  4.9 0.90   -1.0    0.0   -1.0   -1.0   -1.0  23.1  22.0   0.0 0      
 1981 205  5.9 0.90   -1.0    1.2   -1.0   -1.0   -1.0  22.3  22.5   0.0 0      
 1981 205  6.9 1.20   -1.0   34.8   -1.0   -1.0   -1.0  23.0  22.5   0.0 0      
 1981 205  7.9 1.30   -1.0  127.8   -1.0   -1.0   -1.0  24.1  22.7   0.0 0      
 1981 205  8.9 2.00   -1.0  237.0   -1.0   -1.0   -1.0  26.0  22.8   0.0 0      
 1981 205  9.9 2.80   -1.0  427.5   -1.0   -1.0   -1.0  28.1  22.5   0.0 0      
 1981 205 10.9 3.10   -1.0  530.9   -1.0   -1.0   -1.0  29.4  23.0   0.0 0      
 1981 205 11.9 2.90   -1.0  778.3   -1.0   -1.0   -1.0  32.1  23.5   0.0 0      
 1981 205 12.9 4.10   -1.0  812.0   -1.0   -1.0   -1.0  33.8  24.0   0.0 0      
 1981 205 13.9 4.80   -1.0  647.0   -1.0   -1.0   -1.0  34.3  24.0   0.0 0      
 1981 205 14.9 4.60   -1.0  670.3   -1.0   -1.0   -1.0  33.7  24.0   0.0 0      
 1981 205 15.9 4.60   -1.0  406.6   -1.0   -1.0   -1.0  33.7  23.5   0.0 0      
 1981 205 16.9 4.80   -1.0  460.0   -1.0   -1.0   -1.0  34.5  23.0   0.0 0      
 1981 205 17.9 3.40   -1.0  174.2   -1.0   -1.0   -1.0  33.5  22.5   0.0 0      
 1981 205 18.9 2.20   -1.0  110.4   -1.0   -1.0   -1.0  32.1  22.0   0.0 0      
 1981 205 19.9 2.10   -1.0    8.1   -1.0   -1.0   -1.0  28.0  21.5   0.0 0      
 1981 205 20.9 3.40   -1.0    0.0   -1.0   -1.0   -1.0  25.0  21.0   0.0 0      
 1981 205 21.9 1.20   -1.0    0.0   -1.0   -1.0   -1.0  23.7  21.0   0.0 0      
 1981 205 22.9 1.50   -1.0    0.0   -1.0   -1.0   -1.0  23.5  21.0   0.0 0      
 1981 205 23.9 2.30   -1.0    0.0   -1.0   -1.0   -1.0  22.5  21.5   0.0 0      
 1981 206  0.9 1.30   -1.0    0.0   -1.0   -1.0   -1.0  22.3  22.0   0.0 0      
 1981 206  1.9 1.30   -1.0    0.0   -1.0   -1.0   -1.0  21.5  22.0   0.0 0      
 1981 206  2.9 2.30   -1.0    0.0   -1.0   -1.0   -1.0  21.3  22.0   0.0 0      
 1981 206  3.9 2.80   -1.0    0.0   -1.0   -1.0   -1.0  21.0  21.5   0.0 0      
 1981 206  4.9 3.30   -1.0    0.0   -1.0   -1.0   -1.0  19.6  21.0   0.0 0      
 1981 206  5.9 3.40   -1.0    1.2   -1.0   -1.0   -1.0  17.8  20.5   0.0 0      
 1981 206  6.9 2.40   -1.0   65.1   -1.0   -1.0   -1.0  17.9  20.5   0.0 0      
 1981 206  7.9 2.80   -1.0  256.7   -1.0   -1.0   -1.0  19.5  21.0   0.0 0      
 1981 206  8.9 2.70   -1.0  369.4   -1.0   -1.0   -1.0  20.1  21.5   0.0 0      
 1981 206  9.9 1.90   -1.0  628.5   -1.0   -1.0   -1.0  22.6  22.0   0.0 0      
 1981 206 10.9 1.20   -1.0  626.1   -1.0   -1.0   -1.0  25.2  22.5   0.0 0      
 1981 206 11.9 1.20   -1.0  709.8   -1.0   -1.0   -1.0  26.3  23.0   0.0 0      
 1981 206 12.9 2.10   -1.0  957.2   -1.0   -1.0   -1.0  27.8  23.5   0.0 0      
 1981 206 13.9 1.40   -1.0  939.8   -1.0   -1.0   -1.0  29.9  24.0   0.0 0      
 1981 206 14.9 1.20   -1.0  470.5   -1.0   -1.0   -1.0  29.2  24.0   0.0 0      
 1981 206 15.9 1.70   -1.0  297.4   -1.0   -1.0   -1.0  28.7  24.0   6.4 3      
 1981 206 16.9 1.00   -1.0  139.4   -1.0   -1.0   -1.0  27.8  24.0   0.0 0      
 1981 206 17.9 1.20   -1.0  151.0   -1.0   -1.0   -1.0  28.0  24.0   0.0 0      
 1981 206 18.9 1.90   -1.0  209.1    1.0   -1.0   -1.0  28.2  24.0   0.0 0      
 1981 206 19.9 1.00   -1.0   37.2   -1.0   -1.0   -1.0  24.7  24.0   0.0 0      
 1981 206 20.9 0.90   -1.0    1.2   -1.0   -1.0   -1.0  22.7  23.7   0.0 0      
 1981 206 21.9 0.80   -1.0    0.0   -1.0   -1.0   -1.0  22.4  23.5   0.0 0      
 1981 206 22.9 0.90   -1.0    0.0   -1.0   -1.0   -1.0  22.5  23.2   0.0 0      
 1981 206 23.9 0.90   -1.0    0.0   -1.0   -1.0   -1.0  22.3  23.2   0.0 0      
 1981 207  0.9 1.80   -1.0    0.0   -1.0   -1.0   -1.0  21.4  23.2   0.0 0      
 1981 207  1.9 1.60   -1.0    0.0   -1.0   -1.0   -1.0  20.5  23.2   0.0 0      
 1981 207  2.9 0.90   -1.0    0.0   -1.0   -1.0   -1.0  20.0  22.9   0.0 0      
 1981 207  3.9 1.00   -1.0    0.0   -1.0   -1.0   -1.0  20.5  23.0   0.0 0      
 1981 207  4.9 1.00   -1.0    0.0   -1.0   -1.0   -1.0  20.5  23.2   0.0 0      
 1981 207  5.9 1.40   -1.0    1.2   -1.0   -1.0   -1.0  20.4  23.2   0.0 0      
 1981 207  6.9 1.30   -1.0    7.0   -1.0   -1.0   -1.0  20.4  23.5   0.0 0      
 1981 207  7.9 2.20   -1.0   82.5   -1.0   -1.0   -1.0  20.7  23.5   0.0 0      
 1981 207  8.9 2.80   -1.0  248.6   -1.0   -1.0   -1.0  21.2  24.0   0.0 0      
 1981 207  9.9 2.30   -1.0  335.7   -1.0   -1.0   -1.0  21.9  24.0   0.0 0      
 1981 207 10.9 2.10   -1.0  329.9   -1.0   -1.0   -1.0  22.3  24.5   0.0 0      
 1981 207 11.9 2.10   -1.0  353.1   -1.0   -1.0   -1.0  22.7  24.5   0.0 0      
 1981 207 12.9 2.00   -1.0  770.2   -1.0   -1.0   -1.0  24.7  25.0   0.0 0      
 1981 207 13.9 2.60   -1.0  604.1   -1.0   -1.0   -1.0  25.3  25.1   0.0 0      
 1981 207 14.9 2.40   -1.0  374.1   -1.0   -1.0   -1.0  24.6  25.5   0.0 0      
 1981 207 15.9 2.80   -1.0  427.5   -1.0   -1.0   -1.0  24.9  25.5   0.0 0      
 1981 207 16.9 3.00   -1.0  491.4   -1.0   -1.0   -1.0  25.5  25.0   6.1 3      
 1981 207 17.9 3.60   -1.0  167.3   -1.0   -1.0   -1.0  23.8  24.5   6.1 3      
 1981 207 18.9 3.80   -1.0   55.8   -1.0   -1.0   -1.0  22.8  24.0   6.1 3      
 1981 207 19.9 3.00   -1.0    5.8   -1.0   -1.0   -1.0  21.8  23.5   6.1 3      
 1981 207 20.9 2.50   -1.0    0.0   -1.0   -1.0   -1.0  21.6  23.5   6.1 3      
 1981 207 21.9 2.30   -1.0    0.0   -1.0   -1.0   -1.0  20.0  23.0   6.1 3      
 1981 207 22.9 2.50   -1.0    0.0   -1.0   -1.0   -1.0  19.8  23.0   0.0 0      
 1981 207 23.9 3.20   -1.0    0.0   -1.0   -1.0   -1.0  19.2  22.5   0.0 0      
 1981 208  0.9 2.50   -1.0    0.0   -1.0   -1.0   -1.0  19.1  22.5   0.0 0      
 1981 208  1.9 2.40   -1.0    0.0   -1.0   -1.0   -1.0  18.3  22.0   0.0 0      
 1981 208  2.9 1.60   -1.0    0.0   -1.0   -1.0   -1.0  18.2  21.5   0.0 0      
 1981 208  3.9 1.10   -1.0    0.0   -1.0   -1.0   -1.0  18.1  21.0   0.0 0      
 1981 208  4.9 1.40   -1.0    0.0   -1.0   -1.0   -1.0  18.1  20.5   0.0 0      
 1981 208  5.9 1.80   -1.0    1.2   -1.0   -1.0   -1.0  17.8  20.5   0.0 0      
 1981 208  6.9 2.10   -1.0   40.7   -1.0   -1.0   -1.0  18.1  20.5   0.0 0      
 1981 208  7.9 2.00   -1.0   97.6   -1.0   -1.0   -1.0  18.7  21.0   0.0 0      
 1981 208  8.9 3.60   -1.0  160.3   -1.0   -1.0   -1.0  19.5  21.0   0.0 0      
 1981 208  9.9 2.70   -1.0  163.8   -1.0   -1.0   -1.0  20.6  21.5   0.0 0      
 1981 208 10.9 2.40   -1.0  278.8   -1.0   -1.0   -1.0  22.0  22.0   0.0 0      
 1981 208 11.9 3.10   -1.0  535.5   -1.0   -1.0   -1.0  24.4  23.0   0.0 0      
 1981 208 12.9 2.70   -1.0  447.2   -1.0   -1.0   -1.0  24.6  23.5   0.0 0      
 1981 208 13.9 2.10   -1.0  713.3   -1.0   -1.0   -1.0  26.0  24.0   0.0 0      
 1981 208 14.9 1.10   -1.0  564.6   -1.0   -1.0   -1.0  26.4  24.5   0.0 0      
 1981 208 15.9 0.90   -1.0  658.7   -1.0   -1.0   -1.0  27.8  24.5   0.0 0      
 1981 208 16.9 1.50   -1.0  520.4   -1.0   -1.0   -1.0  27.7  24.0   0.0 0      
 1981 208 17.9 1.90   -1.0  451.9   -1.0   -1.0   -1.0  27.8  23.5   0.0 0      
 1981 208 18.9 2.60   -1.0   66.2   -1.0   -1.0   -1.0  25.2  22.5   0.0 0      
 1981 208 19.9 3.10   -1.0    5.8   -1.0   -1.0   -1.0  23.0  21.5   0.0 0      
 1981 208 20.9 2.50   -1.0    2.3   -1.0   -1.0   -1.0  21.3  21.5   0.0 0      
 1981 208 21.9 2.80   -1.0    0.0   -1.0   -1.0   -1.0  21.6  21.5   0.0 0      
 1981 208 22.9 2.60   -1.0    0.0   -1.0   -1.0   -1.0  19.9  21.5   0.0 0      
 1981 208 23.9 2.60   -1.0    0.0   -1.0   -1.0   -1.0  19.9  21.5   0.0 0      
 1981 209  0.9 2.60   -1.0    0.0   -1.0   -1.0   -1.0  19.9  21.5   0.0 0      
 year,day,std time,wind,PAR,NIR,Fbeam PAR,Fbeam NIR, sky
 thermal,Tair,eair,precip,code=3 for rain.