WPC 2BSc Z7#|d2cpi (M)h8) (J)Prestige Elite 12cpi (M)Prestige Elite 12cpi Bold (J)Prestige Elite 16.67cpi (J)HP LaserJet Series II (Additional)HLSEIIAD.PRSd4 @,t0kG@Prestige Elite 12cpi (M)Prestige Elite 12cpi Bold (J):0*((@@  uT . REFERENCES ă X Acock, B., V. R. Reddy, F. D. Whisler, D. N. Baker, J. M. McKinion, H. F. Hodges, and K. J. Boote. 1982. The soybean crop simulator GLYCIM: Model documentation. USDAARS Crop Simulation Research Unit and Department of Agronomy, Mississippi State University, Mississippi State, MS. 316 pp.(# X Baker, D. N., J. R. Lambert, and J. M. McKinion. 1983. GOSSYM: A Simulator of Cotton Crop Growth and Yield. Tech. Bull. 1089, S. Carolina Agric. Expt. Sta., Clemson, SC. 134 pp.(# X Bouwer, H., and R. C. Rice. 1965. Field measurement of critical pressure head of soils. 1965 Annual Report, U. S. Water Conservation Laboratory, Phoenix, AZ. 411 to 4119.(# X Brooks, R. H., and A. T. Corey. 1964. Hydraulic properties of porous media. Hydrology Papers, Colorado State University, Ft. Collins, CO. 3: 127.(# X Brust, K. J., C. H. M. van Bavel, and G. B. Stirk. 1968. Hydraulic properties of a clay loam soil and the field measurement of water uptake by roots: III. Comparison of field and laboratory data on retention and of measured and calculated conductivities. Soil Sci. Soc. Amer. Proc. 32: 322326.(# X Cataldo, D. A., M. Haroon, L. E. Schrader, and V. L. Youngs. 1975. Rapid colorimetric determination of nitrate in plant tissue by nitration of salicylic acid. Commun. Soil Sci. Plant Anal. 6:7180.(# X Doerge, T. A., C. Farr, and J. Watson. 1986. Survey of residual soil nitratenitrogen in three cottongrowing areas of Maricopa County, Arizona. Bull. 8657 Cooperative Extension Service, University of Arizona, Tucson, AZ.(# X Erie, L. J., O. F. French, D. A. Bucks, and K. Harris. 1982. Consumptive use of water by major crops in the southwestern United States. Conservation Res. Rept. No. 29, U. S. Dept. of Agric., Agric. Res. Ser., Wash. DC. 40 pp.(# X Gardner, W. R., and M. S. Mayhugh. 1958. Solutions and tests of the diffusion equation for the movement of water in soil. Soil Sci. Soc. Amer. Proc. 22: 197201.(# X Fink, D. H., and R. D. Jackson. 1973. An equation for describing water vapor adsorption isotherms of soils. Soil Sci. 116: 256261.(# X IBSNAT, the International Benchmark Sites Network for Agrotechnology Transfer. 1986. IBSNAT Technical Report 5, Decision Support System for Agrotechnology Transfer (DSSAT), Documentation for IBSNAT Crop Model Input and Output Files, Version 1.0. Department of Agronomy and Soil Science, University of Hawaii, Honolulu, HA. 53 pp.(# X IBSNAT, the International Benchmark Sites Network for Agrotechnology Transfer. 1988. IBSNAT Technical Report 1, Experimental Design and Data Collection Procedures for IBSNAT, Third Edition, Department of Agronomy and Soil'?0*((@@ Science, University of Hawaii, Honolulu, HA. 73 pp.(# X Idso, S. B., R. D. Jackson, R. J. Reginato, B. A. Kimball, and F. S. Nakayama. 1975. the dependence of bare soil albedo on soil water content. J. Appl. Meteorol. 14: 109113.(# X Jackson, R. D. 1972. On the calculation of hydraulic conductivity. Soil Sci. Soc. Amer. Proc. 36:380382.(# X Jackson, R. D. 1973. Diurnal changes in soil water content during drying. In R. R. Bruce et al. (ed.), Field soil water regime. Special Pub. 5, Soil Sci. Soc. Amer., Madison, WI. 3755.(# X Jackson, R. D., R. J. Reginato, B. A. Kimball, and F. S. Nakayama. 1974. Diurnal soil water evaporation: Comparison of measured and calculated soilwater fluxes. Soil Sci. Soc. Amer. Proc. 38:861866.(# X Jones, C. A., and J. R. Kiniry (eds.). 1986. CERESMaize: A simulation model of maize growth and development. Texas A & M University Press, College Station, TX. 194 pp.(# X Keeney, D. R., and D. W. Nelson. 1982. Nitrogen inorganic forms. In A. L. Page et al. (eds.), Methods of soil analysis, Part 2, Chemical and microbiological properties, Amer. Soc. Agronomy and Soil Sci. Soc. Amer., Madison, WI. 595624.(# X Kimball, B. A., S. T. Mitchell, and G. Brooks. 1981. Annual Report, U. S. Water Conservation Laboratory, Phoenix, AZ. (# X Kimball, B. A., and R. D. Jackson. 1971. Seasonal effects on soil drying after irrigation. Hydrology and Water Resources in Arizona and the Southwest, Vol. I, Univ. of Arizona Press, Tucson, AZ. 8598.(# X Kimball, B. A., J. R. Mauney, G. Guinn, F. S. Nakayama, P. J. Pinter, Jr., K. L. Clawson, R. J. Reginato, and S. B. Idso. 1983. Effects of Increasing At tT  mospheric CO2 on the Yield and Water Use of Crops. No. 021, U. S. Dept. of Energy Series, Response of Vegetation to Carbon Dioxide. Agricultural Research Service, U. S. Dept. of Agriculture, Washington, DC. 37 pp.(# X Kimball, B. A., J. R. Mauney, G. Guinn, F. S. Nakayama, P. J. Pinter, Jr., K. L. Clawson, S. B. Idso, G. D. Butler, Jr., and J. W. Radin. 1984. Effects  tT of Increasing Atmospheric CO2 on the Yield and Water Use of Crops. ; No. 023, Response of Vegetation to Carbon Dioxide, U. S. Dept. of Energy, Carbon Dioxide Research Division, and the U. S. Dept. of Agriculture, Agricultural Research Service, Washington, DC. 60 pp.(# X Kimball, B. A., J. R. Mauney, G. Guinn, F. S. Nakayama, S. B. Idso, J. W. Radin, D. L. Hendrix, G. D. Butler, Jr., T. I. Zarembinski, and P. E. Nixon, III.  tTH& 1985. Effects of Increasing Atmospheric CO2 on the Yield and Water Use of Crops. No. 027, Response of Vegetation to Carbon Dioxide, U. S. Dept. of Energy, Carbon Dioxide Research Division, and the U. S. Dept. of Agricul'@0*((@@Ԯture, Agricultural Research Service, Washington, DC. 75 pp.(# X Kimball, B. A., J. R. Mauney, J. W. Radin, F. S. Nakayama, S. B. Idso, D. L. Hendrix, D. H. Akey, S. G. Allen, M. G. Anderson, and W. Hartung. 1986.  tT  Effects of Increasing Atmospheric CO2 on the Growth, Water Relations, and Physiology of Plants Grown under Optimal and Limiting Levels of Water and Nitrogen. No. 039, Response of Vegetation to Carbon Dioxide, U. S. Dept. of Energy, Carbon Dioxide Research Division, and the U. S. Dept. of Agriculture, Agricultural Research Service, Washington, DC. 125 pp.(# X Kimball, B. A., J. R. Mauney, D. H. Akey, D. L. Hendrix, S. G. Allen, S. B. Idso, J. W. Radin, and E. A. Lakatos. 1987. Effects of Increasing At tT` mospheric CO2 on the Growth, Water Relations, and Physiology of Plants Grown under Optimal and Limiting Levels of Water and Nitrogen. No. 049, Response of Vegetation to Carbon Dioxide, U. S. Dept. of Energy, Carbon Dioxide Research Division, and the U. S. Dept. of Agriculture, Agricultural Research Service, Washington, DC. 124 pp.(# X Mualem, Y. 1976. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resources Res. 12: 513522.(# X Nakayama, F. S. 1965. Dispersion and flocculation of soil and clay minerals as related to the Na and Ca status of the ambient solution. 1965 Annual Report, U. S. Water Conservation Laboratory, Phoenix, AZ. 401 to 409.(# X Nielsen, D. R., R. D. Jackson, J. W. Cary, D. D. Evans et al. 1972. Soil water. Amer. Soc. of Agronomy and Soil Sci. Soc. Amer., Madison, WI.(# X Ritchie, J. T. 1972. Model for predicting evaporation from a row crop with incomplete cover. Water Resources Research 8: 12041213.(# X Ritchie, J. T. 1985. A userorientated model of the soil water balance. In W. Day and R. K. Arkin (eds.), Wheat growth and modeling, Plenum Publishing Corp. 293305.(# X Ritchie, J. T., and S. Otter. 1985. Description and performance of CERESWheat: A useroriented wheat yield model. USDAARS. ARS38. 159175.(# X Ritchie, J. T., J. R. Kiniry, C. A. Jones, and P. T. Dyke. 1986. Model inputs. In C. A. Jones and J. R. Kiniry, CERESMaize: A simulation model of maize growth and development, Texas A&M University Press, College Station, TX. 3748.(# X Rosenberg, N. J., B. L. Blad, and S. B. Verma. 1983. Microclimate: The biological environment. John Wiley and Sons, New York. 495 pp.(# X Schmidli, R. J. 1986. Climate of Phoenix, Arizona. NOAA Technical Memorandum NWS WR177, National Weather Service, National Oceanic and Atmospheric Administration, U. S. Dept. of Commerce, Salt Lake City, UT. 126 pp.(# X Soil Improvement Committee, California Fertilizer Association. 1985. Western fertilizer handbook. Interstate Printers and Publishers, Inc., Danville,'A0*((@@ IL.(# X Schnell, S. M. 1983. Field psychrometer construction and evaluation. U. S. Water Conservation Laboratory, Phoenix, AZ. 10 pp.(# X Thompson, L. M. 1957. Soils and soil fertility. McGrawHill Book Co., Inc., New York, NY.(# X Van Bavel, C. H. M., G. B. Stirk, and K. J. Brust. 1968a. Hydraulic properties of a clay loam soil and the field measurement of water uptake by roots: I. Interpretation of water content and pressure profiles. Soil Sci. Soc. Amer. Proc. 32: 310317.(# X Van Bavel, C. H. M., K. J. Brust, and G. B. Stirk. 1968b. Hydraulic properties of a clay loam soil and the field measurement of water uptake by roots: II. The water balance of the root zone. Soil Sci. Soc. Amer. Proc. 32: 317321.(# X Van Genuchten, M. Th. 1980. A closedform equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Amer. J. 44: 892898.(# X Van Genuchten, M. Th., and D. R. Nielsen. 1985. On describing and predicting the hydraulic properties of unsaturated soils. Annales Geophysicae 3: 615628.(# X Van Genuchten, M. Th., F. J. Leij, and S. R. Yates. 1992. The RETC code for quantifying the hydraulic functions of unsaturated soils. EPA Report ____, Robert S. Kerr Environmental Res. Lab., U. S. Environmental Protection Agency, Ada, OK. (In press).(# X Whisler, F. D. 1976. Calculating the unsaturated hydraulic conductivity and diffusivity. Soil Sci. Soc. Amer. J. 40: 150151.(# X Whisler, F. D. 1982. Soil data preparation for crop growth models. Agronomic Series No. 1, Mississippi State University, Mississippi State, MS. 9 pp.(# X Wilkerson, G. G., J. W. Jones, K. J. Boote, K. T. Ingram, and J. W. Mishoe. 1983. Modeling soybean growth for management. Trans. of the ASAE 26:6373.(# X Wosten, J. H. M., and M. Th. van Genuchten. 1988. Using texture and other soil properties to predict unsaturated hydraulic functions. Soil Sci. Soc. Am. J. 52: 17621770.(#  B0*((@@   Ђ\/GLOSSARY ALPHA(L)` `  parameter in the van Genuchten soil water retention curve ` `  (van Genuchten, 1980; van Genuchten and Nielsen, 1985;  tT ` `  van Genuchten et al., 1992) (mé1). AFERT` `  amount of nitrogen fertilizer added on JFDAY(J) (kg N/ha) ` `  (same as FERTN). AIRDR` `  volumetric water content of air dry soil in soil layer L  tT ` `  (m3/m3). AMP` `  annual amplitude (maximum minus minimum) in mean monthly  tT ` `  temperature ($C). AMTIRR` `  amount of irrigation added on JDIRR(J) (mm). A00` `  character representation of day of year.  tT BD(L)` `  moist bulk density of the soil in soil layer L (Mg/m3). BDRATO` `  ratio of initial to final volumetric water content of the ` `  bottom soil layer. BDSLOP` `  rate of change of water content in bottom soil layer after day  tTP ` `  of first bloom ((m3/m3)/day). BEGDATE` `  beginning date in weather file (8 characters). BETA(L)` `  slope of graph of log (hydraulic diffusivity) vs. volumetric ` `  water content for soil layer L (Gardner and Mayhugh 1958). CATEXC(L)` `  cation exchange capacity in soil layer L, (meq/kg). CLAY(L)` `  percentage of clay in soil layer L (%). CN2` `  SCS curve number used to calculate daily runoff.  tTx COND` `  soil hydraulic conductivity (m/s).  tT CO2` `  mean daily carbon dioxide concentration (mol/mol). CO2DAT` `  switch to indicate if daily carbon dioxide concentration data ` `  are available (Yes=1,No=0).  tT# CO2YR` `  mean seasonal CO2 concentration (mol/mol). CTEXPyy.DIR` `  file name of directory of cotton experiments where yy is year. DATAID` `  soil classification. 'C0*((@@ԌDEWDAT` `  switch to indicate if dew point temperature data are available ` `  (Yes=1,No=0).  tTX DEWPT` `  dew point temperature ($C). DFERT` `  depth of incorporation of fertilizer application (cm).  tTx DIFF` `  soil water diffusivity (m2/s). DIFF0` `  diffusivity of soil layer L at a low reference water content  tT ` `  for Gardner and Mayhugh (1958) equation (m2/s). DIFRN` `  path to Experiment Directory File DLAYR(L)` `  thickness of soil layer L (cm). DMOD` `  zerotounity factor which reduces the rate constant for ` `  mineralization of humus pool of soil for which organic ` `  matter is chemically or physically protected (default =1). DNH4(L)` `  default soil ammonium in soil layer L (mg of N/kg of soil). DNO3(L)` `  default soil nitrate in soil layer L (mg of N/kg of soil). DPH(L)` `  default pH of soil layer L. DSOIL` `  irrigation management depth (cm).  tT DSW(L)` `  default water content of soil layer L (m3/m3).  tTp DUL(L)` `  drained upper limit soil water content for soil layer L (m3/m3). EFFIRR` `  irrigation system efficiency (fraction). ENDDATE` `  ending date in weather file (8 characters). ETA(L)` `  soil characteristic parameter (Brooks & Corey, 1964) relating (#(#Z` `  volumetric water content to soil water potential (matric) ` `  for soil inhh*layer L. EXPDES` `  experiment description (40 characters). EXPID` `  experiment identifier (8 characters) with in institute code, ` `  site code, year of the experiment, and experiment number. EXPTNO` `  experiment number. FC(L)` `  "field capacity" defined as the volumetric water content at a ` `  soil water potential (matric) of 33 kPa or other  tTH& ` `  appropriate value, as specified in PSISFC (m3/m3). FERCOD` `  code for placement of fertilizer application, as defined by'D0*((@@Ԍ` `  IBSNAT (1988) where: ` `  01 broadcast, not incorporated ` `  02 broadcast, incorporated ` `  03 banded on surface ` `  04 banded beneath surface ` `  05 applied in irrigation water ` `  06 foliar spray ` `  07 bottom of hole ` `  08 on the seed ` `  99 other FERTIN` `  effective amount of inoculants and amendments (kg/ha). FERTK` `  effective amount of potassium in fertilizer added (kg/ha). FERTN` `  effective amount of nitrogen in fertilizer added (kg/ha). FERTP` `  effective amount of phosphorus in fertilizer added (kg/ha). FILE1 to` `  names of the data files as defined in the 'File Structures FILE9` `  and Data Formats' section. FILEA FILEB  tT GBD(L)` `  bulk density of soil layer L (g/cm3). GH2OC( )` `  gravimetric water content of a particular bulk densityroot ` `  impedance curve (kg/kg). GLAYR(L)` `  thickness of soil layer L (cm).  tT8 GRAV` `  acceleration of gravity (m/s2). H` `  soil water pressure head (m). HISTRY` `  brief description of previous cropping history for site. ID` `  is a combination of codes for institute, site, year (last ` `  two digits), and experiment number. IDUMSL` `  number assigned to a soil type. IFTYPE` `  code number for the type of fertilizer added on JDFERT, as ` `  defined by IBSNAT (1988). IITYPE` `  code number for the type of inoculants and amendments added on ` `  JDFERT, as defined by IBSNAT (1988). IIRR` `  switch describing irrigation (default = 1). ` `  1: no irrigation ` `  2: irrigation applied using field schedule ` `  3: automatically irrigated at threshold soil water'E0*((@@Ԍ` `  4: assume no water stress, water balance not used IKTYPE` `  code number for the type of potassium fertilizer, as defined ` `  by IBSNAT (1988). IMERGE` `  emergence date, Julian day of the year. INRIM( )` `  number of pairs of bulk density / root impedance values defining ` `  a curve at a particular water content. INSTE` `  code for institute ID in files FILE6, FILE7, FILE8, FILEA, ` `  FILEB, and weather files. INSTS` `  code for institute ID in soil files FILE4 and FILE5. INSTW` `  code for institute ID in weather files FILE1. INTYPE` `  code number for the type of nitrogen fertilizer, as defined by ` `  IBSNAT (1988), where: ` `  05 urea ` `  08 calcium nitrate ` `  10 urea ammonium nitrate solution. IPTYPE` `  code number for the type of potassium fertilizer, as defined by ` `  IBSNAT (1988). IRRCOD` `  type of irrigation system, as defined by IBSNAT (1988) where: ` `  01 furrow ` `  02 alternating furrows ` `  03 flood ` `  04 sprinkler ` `  05 drip or trickle ` `  99 other ISIM` `  Julian day of year simulation begins. ISOILT` `  soil number for this treatment. ISOW` `  sowing date, Julian day of the year. ISWDIS` `  switch to indicate if disease data are available (Yes=1,No=0). ISWEED` `  switch to indicate if weed data are available (Yes=1,No=0). ISWINS` `  switch to indicate if insect data are available (Yes=1,No=0). ISWNEM` `  switch to indicate if nematode data are available (Yes=1,No=0). ISWNIT` `  switch to indicate if nitrogen routines are used (default=0) ` `  0: nitrogen subroutines are not used, assumes adequate (#(#Z` `   nitrogen. ` `  1: nitrogen subroutines are used.'F0*((@@Ԍ IVARTY` `  cultivar number for this treatment. IYR` `  year for which weather data is being read.  tT J` `  soil water flux (m3 of water mé2 sé1). JDFERT` `  Julian day of the year of fertilizer application ` `  (same as JFDAY). JDIRR` `  Julian day of the year of irrigation event J. JDOY` `  Julian day of year. JEMRGD` `  day of emergence, which is the day of year that 50% of the ` `  plants emerge from the soil (day of year). JFDAY` `  Julian day of the year of fertilizer application. JFLRJD` `  day of first flower, which is the day of year that 50% of the ` `  plants display their first flower (day of year).  tT0 JNABSM` `  measured number of abscised sites per m2.  tT JNFLWM` `  measured number of flowers per m2.  tTP JNGBLM` `  measured number of green bolls per m2.  tT JNMBLM` `  measured number of mature bolls per m2.  tTp JNNODM` `  measured number of nodes per m2.  tT JNSQRM` `  measured number of squares per m2. JSQRJD` `  day of first square, which is the day of year that 50% of the ` `  plants display their first square (day of year.) JUL` `  Julian date of weather record in data file. L` `  soil layer index. LISTID` `  Subroutine called by main program of RETRVE to identify which ` `  files are to be retrieved (read) and sets the path to these ` `  files. LL(L)` `  lower limit of plantextractable soil water for soil layer L, (#(#Z tT# ` `  (m3/m3). LYRSOL` `  number of soil layers. NCURVE` `  number of bulk density / root impedance curves, ` `  one for each water content.'G0*((@@ԌNLAYR` `  number of soil layers. NH4(L)` `  soil ammonium in soil layer L (mg N/kg soil). NOVAR` `  number of "state" variables for which there are intermediate ` `  growth data in FILEB. NO3(L)` `  soil nitrate in soil layer L (mg N/kg soil). NPATH` `  number of characters in path given by PTODAT. NV( )` `  array of pointers to indicate which "state" variables for which ` `  there are data in FILEB. OC(L)` `  organic carbon concentration in soil layer L (%). OUT1 to` `  output file names. OUT4 PARDAT` `  switch to indicate if PAR data are available (Yes=1,No=0).  tTh PARFAC` `  factor to convert MJ/m2 to mol/m2 (2.30 if unknown). PATHINP.INC` `  four line file of pathnames for input to Program RETRVE. PEDON` `  SCS pedon number. PH(L)` `  pH of soil in soil layer L in a 1:1 soil water ` `  slurry. PHFAC3` `  variable to reduce apparent photosynthesis attributed to ` `  soil fertility (for grain legume models, default = 1.00). PHOSPH` `  soil phosphorus in soil layer L (mg P/kg soil).  tTX PLANTS` `  plant population (plants/m2) POTASS` `  soil phosphorus in soil layer L (mg K/kg soil). PROCK` `  procedure used to determine soil K. PROCN` `  procedure used to determine soil N. PROCP` `  procedure used to determine soil P. PSI` `  soil water potential expressed as pressure head (m). PSIBUB(L)` `  soil water potential (matric) near saturation taken to be ` `  bubbling pressure at which air starts to enter the soil ` `  as it$dries (kPa, positive number). PSISFC` `  soil water potential (matric) used to define field capacity'H0*((@@Ԍ` `  (kPa,$positive number). PSIR` `  reference low (dry) soil water potential (matric) in Marani ` `  equation corresponding to THETAR (kPa). PTODAT` `  path to directory containing the data files. READDA` `  Subroutine in Program RETRVE that reads data files. RETRVE.FOR` `  FORTRAN program to read data files. RETRVE.OUT` `  output file from Program RETRVE which is an echo of the ` `  data files read by RETRVE.  tT RHOS` `  density of solid soil particles (Mg/m3).  tT RHOW` `  density of water (Mg/m3). ROCK(L)` `  coarse rock fragments in soil layer L (%). ROOT` `  dry weight of root residue of previous crop (kg/ha, ` `  default=500). ROWSPC` `  row spacing (m). RSAT` `  relative saturation as defined by Eq. 14. RWUMX` `  maximum daily root water uptake per unit root length  tT ` `  (cm3/cm rootday, default = 0.03). SALB` `  bare soil albedo. SAND(L)` `  sand in soil layer L (%).  tT SAT(L)` `  saturated water content for soil layer L (m3/m3). SCN` `  C:N ratio of surface residue of previous crop ` `  (kg C/kg N, default = 75.).  tTx SCOND(L)` `  saturated hydraulic conductivity (m/s, or m/s in the data ` `  files). SDEP` `  depth of surface residue incorporation (cm). SDEPTH` `  sowing depth L (cm). SILT(L)` `  silt in soil layer L (%). SITEE` `  code for site ID in files FILE6, FILE7, FILE8, FILEA, and FILEB. SITES` `  code for site ID in soil files FILE4 and FILE5. 'I0*((@@ԌSNAME` `  description of soil for root impedance data.  tT SOLRAD` `  daily total solar radiation (MJ/m2). STATW` `  code for weather station ID (2 characters). STDAT` `  switch to indicate if soil temperature data are available ` `  (Yes=1,No=0). STMAX` `  daily value of maximum temperature at 10 cm depth in bare soil  tT ` `  ($C). STMIN` `  daily value of minimum temperature at 10 cm depth in bare soil  tT( ` `  ($C). STRAW` `  weight of organic residue of previous crop and/or added green ` `  manure (kg/ha).  tT SW(L)` `  soil water content for soil layer L (m3/m3). SWCON` `  soil water drainage constant, (fraction drained/day). SWCON1` `  coefficient in the steadystate solution of the radialflow  tT ` `  root uptake equation (cm3/cm rootday, default=0.00267). SWCON2` `  coefficient in the steadystate solution of the radialflow  tTP ` `  root uptake equation (cm3/cm rootday, default=58.0). SWCON3` `  coefficient in the steadystate solution of the radialflow  tT ` `  root uptake equation (cm3/cm rootday, default=6.68).  tT8 TAV` `  annual average ambient air temperature ($C). TAXON` `  soil classification. TD` `  number of days after first bloom that water content of the ` `  bottom soil layer begins to change (days).  tT THETA` `  volumetric soil water content (m3/m3). THETAC` `  percent available water which triggers irrigation (%). THETAI` `  initial water content of bottom layer in soil profile for some  tT! ` `  sites and data sets where this is known (m3/m3). THETAR(L)` `  "residual" soil water content of Brooks and Corey (1964) or ` `  "reference" soil water content near the wilting point for ` `  the Marani equation, taken as THETA0 slightly reduced  tT% ` `  (typically by 0.001) (m3/m3). THETAS(L)` `  volumetric water content of soil in layer L at "saturation". ` `  Assumed to be the same as the water content at which air'J0*((@@Ԍ` `  starts to enter the soil as it dries (ie. at the bubbling  tT ` `  pressure) (m3/m3). THETA0(L)` `  reference volumetric water content of soil in layer L for ` `  Gardner and Mayhugh (1958) equation, generally at a low  tT ` `  (wilting point) soil water potential of 1500 kPa (m3/m3). TITLET` `  title of treatment. TSTBD( , )` `  bulk density at a point on bulk density / root impedance curve  tT ` `  (Mg/m3). TRTNO` `  treatment number. TSTIMP( , )` `  root impedance at a point on the bulk density / root impedance  tT ` `  curve (kg/cm2). U` `  upper limit of stage 1 soil evaporation (mm). VALFAB` `  subroutine in Program RETRVE for reading validation data from ` `  FILEA and FILEB. VGTHR(L)` `  "residual" soil water content fitted to the extrapolated dry ` `  end of the retention curve for the van Genuchten equation. ` `  (van Genuchten, 1980; van Genuchten and Nielsen, 1985;  tT ` `  van Genuchten et al., 1992) (m3/m3). VGTHS(L)` `  "saturated" soil water content fitted to the extrapolated wet ` `  end of the retention curve for the van Genuchten equation. ` `  (van Genuchten, 1980; van Genuchten and Nielsen, 1985;  tTp ` `  van Genuchten et al., 1992) (m3/m3). WIND` `  daily total wind run (km). WINDYR` `  longterm average daily total wind run for the experimental ` `  site (km). WNDDAT` `  switch to indicate if daily wind data are available ` `  (Yes=1, No=0). WR(L)` `  weighting factor for soil depth L to determine new root growth ` `  distribution. WTHES` `  weather station description (40 characters). WTHID` `  weather station identifier (4 characters). WTHyy.DIR` `  file name of directory of weather files where yy is year. XAPTNP` `  measured nitrogen content of stems plus burrs plus lint ` `  at maturity (kg/ha). 'K0*((@@ԌXBIOM` `  fieldmeasured aboveground biomass at maturity (kg/ha).  tT XBLSM` `  fieldmeasured mature boll number (bolls/m2).  tT  XLAI` `  measured leaf area index (m2/m2).  tT XLAIMX` `  maximum leaf area index during season (m2/m2). XLAT` `  latitude of station (degrees). XLONG` `  longitude of station (degrees). XLTYLD` `  final dry lint yield attainable by hand harvesting (kg/ha). XM(L)` `  parameter in the van Genuchten soil water retention curve. ` `  (van Genuchten, 1980; van Genuchten and Nielsen, 1985; ` `  van Genuchten et al., 1992) XN(L)` `  parameter in the van Genuchten soil water retention curve. ` `  (van Genuchten, 1980; van Genuchten and Nielsen, 1985; ` `  van Genuchten et al., 1992)  tT0 XPAR` `  daily total PAR (mol/m2). XPLTHT` `  measured plant height (cm). XRAIN` `  daily total precipitation (mm/day). XSDN` `  measured nitrogen content of seed at maturity (kg/ha). XSDTN` `  measured nitrogen concentration in seed at maturity (%). XSDWT` `  measured seed dry weight (g/seed). XSDYLD` `  final dry seed yield (kg/ha). XSPB` `  measured number of seeds per boll (seeds/boll). XSTBR` `  fieldmeasured stem plus burr weight at maturity (kg/ha).  tT@ XTMAX` `  daily value of maximum air temperature ($C).  tT XTMIN` `  daily value of minimum air temperature ($C). XTOTNP` `  measured total crop nitrogen content at maturity (kg/ha). XWBURR` `  measured dry weight of burrs (kg/ha). XWGBLH` `  measured dry weight of green bolls (kg/ha). XWLEFH` `  measured dry weight of leaves (kg/ha). 'L0*((@@ԌXWLINT` `  measured dry weight of lint (kg/ha). XWMBLH` `  measured dry weight of mature bolls (kg/ha). XWROTH` `  measured dry weight of roots (kg/ha). XWSEED` `  measured dry weight of seed (kg/ha). XWSTMH` `  measured dry weight of stems (kg/ha). X1FERT` `  distance from plant row to closest edge of fertilizer ` `  application area (cm). X1IRR` `  distance from plant row to closest edge of water application ` `  area (cm). X2FERT` `  distance from plant row to furthest edge of fertilizer ` `  application area (cm). X2IRR` `  distance from plant row to furthest edge of water application ` `  area (cm). YEAR` `  year, last 2 digits Z` `  depth or length variable (m). Z1FERT` `  topmost depth in soil where fertilizer is applied (cm). Z1IRR` `  topmost depth in soil where water is applied (cm). For flood, ` `  sprinkler, or furrow irrigation, the appropriate entry would ` `  be zero. Z2FERT` `  lowest depth in soil where fertilizer is applied (cm). Z2IRR` `  lowest depth in soil where water is applied (cm). XM0*((@@  uS   #FILE STRUCTURES AND DATA FORMATS ă  uT   # Modifications to IBSNAT Standard ă In 1986 a group of plant growth modelers collaborated in an effort to standardize crop model input and output file structures. Operating under the auspices of the International Benchmark Sites of Agrotechnology Transfer (IBSNAT) Project, they published IBSNAT Technical Report 5, Decision Support System for Agrotechnology Transfer (DSSAT), Documentation for IBSNAT Crop Model Input and Output Files, Version 1.0. Such standardization should facilitate cooperation and exchange of data among different groups of workers and thus aid plant model development because adhering to the standard will minimize the need to continually reformat. As is the case with all attempts to standardize, however, the individuals doing the standardizing were working from within the confines of their collective experience, and they could not foretell future needs. In this case, the IBSNAT standard 1.0 was authored primarily by individuals involved with the development of the CERESMaize (Jones and Kiniry, 1986), CERESWheat (Ritchie and Otter, 1985), and SOYGRO (Wilkerson et al., 1985), which for  tTB example do not use humidity information nor do they consider CO2 concentration to be a weather or climate variable. They also use different parameters to characterize soil hydraulic properties than do the GOSSYM (Baker et al., 1983) and GLYCIM (Acock et al., 1982) family of models. Another important (to us) aspect not addressed by IBSNAT (1986) is that cotton was not one of the original crops. Hence, it was necessary to establish our own formats for the'N0*((@@ plant growth and yield data, but following IBSNAT as much as possible. To overcome some of these particular limitations and achieve a more general standard, Acock (B. Acock, personal communication, 1988) prepared tentative modifications to the IBSNAT 1.0 standard. The dialogue is continuing, and therefore any format used in this report is destined for early obsolescence. However, in order that these data be useful for validating models for as many users as possible, the parameters required to meet the IBSNAT (1986) standard are presented as closely as possible. Then the IBSNAT parameters are followed by those proposed by Acock or by ourselves and then by any additional required by GOSSYM. The variable definitions and the formats used are listed in the following Format Tables. Those we have proposed and followed for cotton are listed in the FILEA and FILEB Tables. PO0*((@@  uS  CTEXPyy.DIR: Directory of Experiment Files  Variable Name FORTRAN Format0 Description   tT Format for line 1 EXPID` `  A8$hh*0Experiment identifier. EXPDES` `  1X,A40hh*0Experiment description. FILE1` `  1X,A12hh*0Daily weather data file name. FILE2` `  1X,A12hh*0Soil profile file name.  tT` Format for line 2 FILE4` `  A12$hh*0Soil organic residue file name. FILE5` `  1X,A12hh*0Soil profile initial ` `  $hh*06conditions file name. FILE6` `  1X,A12hh*0Irrigation management file name. FILE7` `  1X,A12hh*0Fertilizer management file name. FILE8` `  1X,A12hh*0Treatment management file name. FILE9` `  1X,A12hh*0Genetic coefficients file name.  tT0 Format for line 3 FILEA` `  A12$hh*0Measured crop harvest summary ` `  $hh*06file name. FILEB` `  1X,A12hh*0Measured intermediate growth ` `  $hh*06data file name. OUT1` `  1X,A7$hh*0File name for output file 1. OUT2` `  1X,A7$hh*0File name for output file 2. OUT3` `  1X,A7$hh*0File name for output file 3. OUT4` `  1X,A7$hh*0File name for output file 4. Repeat in blocks of 3 lines for additional experiments. For the cotton experiments presented in this manuscript, each plot was considered a separate experiment each year, and therefore each block of three lines in the CTEXPyy.DIR files contains the file names of the data files appropriate for that plot.  uS@  WTHyy.DIR: Directory of Weather Files  Variable Name FORTRAN Format0 Description   tT(# Format for all lines WTHID` `  A4$hh*0Weather station ID. WTHDES` `  1X,A40hh*0Weather station description. BEGDATE` `  A8$hh*0Beginning date in weather file. ENDDATE` `  1X,A8$hh*0Ending date in weather file. FILE1` `  1X,A12hh*0Weather file name.'P0*((@@  uS  FILE1: Daily Weather File  Variable Name FORTRAN Format0 Description   tT  Format for line 1 INSTW` `  A2$hh*0Code for institute ID. STATW` `  A2$hh*0Code for weather station ID. XLAT` `  1X,F6.2hh*0Latitude of station (degrees). XLONG` `  1X,F6.2hh*0Longitude of station (degrees). PARFAC` `  1X,F5.2hh*0Factor to convert total daily solar  tT ` `  $hh*06radiation in MJ/m2 to PAR in  tT` ` `  $hh*06mol/m2. PARDAT` `  1X,I1$hh*0PAR data available (Yes=1, No=0). WNDDAT` `  1X,I1$hh*0Wind data available (Yes=1, No=0). DEWDAT` `  1X,I1$hh*0Dewpoint temperature data available ` `  $hh*06(Yes=1, No=0). STDAT` `  1X,I1$hh*0Soil Temperature data available ` `  $hh*06(Yes=1, No=0). CO2DAT` `  1X,I1$hh*0Daily average carbon dioxide concentration ` `  $hh*06data available (Yes=1, No=0). CO2YR` `  1X,F4.0hh*0Mean seasonal carbon dioxide concentration  tT0 ` `  $hh*06(mol/mol, use 345 if unknown, must ` `  $hh*06 be present if CO2DAT=0). WINDYR` `  $1X,F5.106Average daily wind run for the site ` `  $hh*06(km, must be present if WNDDAT=0).  tT Format for all other lines INSTW` `  A2$hh*0Code for institute ID. STATW` `  A2$hh*0Code for weather station ID. IYR` `  1X,I2$hh*0Year of weather record. JUL` `  1X,I3$hh*0Julian date (day of year).  tT SOLRAD` `  1X,F5.2hh*0Daily total solar radiation (MJ/m2).  tT XTMAX` `  1X,F5.1hh*0Daily maximum air temperature T($C).  tTX XTMIN` `  1X,F5.1hh*0Daily minimum air temperature T($C). XRAIN` `  1X,F5.1hh*0Daily total precipitation (mm).  tT XPAR` `  1X,F6.2hh*0Daily total PAR (mol/m2). WIND` `  1X,F5.1hh*0Daily total wind run (km).  tTx DEWPT` `  1X,F5.1hh*0Daily average dew point temperature ($C). STMAX` `  1X,F5.1hh*0Daily maximum soil temperature at the  tT ` `  $hh*0610 cm depth ($C). STMIN` `  1X,F5.1hh*0Daily minimum soil temperature at the  tT! ` `  $hh*0610 cm depth ($C). CO2` `  1X,F4.0hh*0Daily carbon dioxide concentration  tT(# ` `  $hh*06(mol/mol). A00` `  1X,A7$hh*0Optional date (day3 character month) $Q0*((@@  uS  FILE2: Soil Profile Properties  Variable Name FORTRAN Format0 Description   tT  Format for line 1 IDUMSL` `  1X,I2$hh*0Number assigned to a soil type. PEDON` `  1X,A12hh*0SCS pedon number. TAXON` `  1X,A60hh*0Soil classification.  tT Format for line 2 SALB` `  F6.2$hh*0Bare soil albedo. U` `  1X,F5.2hh*0Upper limit of stage 1 soil evaporation ` `  $hh*06(mm). SWCON` `  1X,F6.2hh*0Soil water drainage constant, fraction ` `  $hh*06drained per day. CN2` `  1X,F6.2hh*0SCS curve number used to calculate ` `  $hh*06daily runoff.  tT TAV` `  1X,F5.1hh*0Annual average ambient temperature ($C). AMP` `  1X,F5.1hh*0Maximum minus minimum mean monthly  tTh ` `  $hh*06temperature ($C). DMOD` `  1X,F3.1hh*0Zerotounity factor which reduces the ` `  $hh*06rate constant for mineralization ` `  $hh*06of the humus pool for soils which ` `  $hh*06are poor mineralizers due to ` `  $hh*06chemical or physical protection of ` `  $hh*06the organic matter (default=1). SWCON1` `  1X,E9.2hh*0Coefficient in the radial flow root uptake  tT ` `  $hh*06equation (cm3/cm rootday, ` `  $hh*06default = 0.267E02). SWCON2` `  1X,F6.1hh*0Coefficient in the radial flow root uptake  tT ` `  $hh*06equation (cm3/cm rootday, ` `  $hh*06default = 58.0). SWCON3` `  1X,F5.2hh*0Coefficient in the radial flow root uptake  tTX ` `  $hh*06equation (cm3/cm rootday, ` `  $hh*06default = 6.68). RWUMX` `  1X,F5.2hh*0Maximum daily root water uptake per unit  tT ` `  $hh*06root length (cm3/cm rootday, ` `  $hh*06default = 0.03). PHFAC3` `  1X,F4.2hh*0Variable to reduce apparent photosynthesis ` `  $hh*06attributed to soil fertility ` `  $hh*06(default = 1.00).  tT`" Format for line 3 through NLAYR+2 DLAYR(L)` `  F6.0$hh*0Thickness of soil layer L, cm. LL(L)` `  1X,F6.3hh*0Lower limit of plantextractable soil  tT% ` `  $hh*06water for soil layer L (m3/m3). DUL(L)` `  1X,F6.3hh*0Drained upper limit soil water content  tT' ` `  $hh*06for soil layer L (m3/m3). SAT(L)` `  1X,F6.3hh*0Saturated water content for soil layer L'R0*((@@Ԍ tT ` `  $hh*06(m3/m3). DSW(L)` `  $1X,F6.306Default soil water content for soil layer  tTX ` `  $hh*06L (m3/m3). WR(L)` `  1X,F6.3hh*0Weighting factor for soil depth L to ` `  $hh*06determine new root growth ` `  $hh*06distribution. BD(L)` `  1X,F5.2hh*0Moist bulk density of soil layer L  tT@  ` `  $hh*06(g/cm3). OC(L)` `  1X,F5.2hh*0Organic carbon concentration in soil layer  tT ` `  $hh*06L (g/cm3). DNH4(L)` `  1X,F4.1hh*0Default soil ammonium in soil layer L ` `  $hh*06(mg elemental N/kg soil). DNO3(L)` `  1X,F4.1hh*0Default soil nitrate in soil layer L ` `  $hh*06(mg elemental N/kg soil). DPH(L)` `  1X,F4.1hh*0Default pH of soil layer L in a 1:1 ` `  $hh*06soil:water slurry.  tT Format for line NLAYR + 3 Enter a 1 to signal endofdata for this block of soil profile data.  tT0 Format for lines NLAYR+4 through 2*NLAYR+4 ROCK(L)` `  1X,F4.1hh*0Coarse rock fragments in soil layer L (%). SAND(L)` `  1X,F4.1hh*0Sand in soil layer L (%). SILT(L)` `  1X,F4.1hh*0Silt in soil layer L (%). CLAY(L)` `  1X,F4.1hh*0Clay in soil layer L (%). SCOND(L)` `  1X,F4.1hh*0Saturated hydraulic conductivity in  tT ` `  $hh*06soil layer L (m/s). CATEXC(L)` `  1X,F5.1hh*0Cation exchange capacity in soil layer L ` `  $hh*06(meq/kg). ALPHA(L)` `  1X,F6.3hh*0Empirical parameter for the van Genuchten  tT ` `  $hh*06soil water retention curve (mé1). XN(L)` `  1X,F6.3hh*0Empirical parameter for the van Genuchten ` `  $hh*06soil water retention curve. VGTHS` `  1X,F6.3hh*0"Saturated" water content from fitting ` `  $hh*06van Genuchten soil water retention  tT ` `  $hh*06curve (m3/m3). VGTHR` `  1X,F6.3hh*0"Residual" water content from fitting ` `  $hh*06van Genuchten soil water retention  tT ` `  $hh*06curve (m3/m3).  tT(# Format for line 2*NLAYR+5 DATAID` `  20A4$hh*0Soil description  tTH& Format for line 2*NLAYR+6 LYRSOL` `  *$hh*0Number of soil layers (GOSSYM model)'S0*((@@Ԍ tT ԙFormat for lines 2*NLAYR+7 through 2*NLAYR+7+LYRSOL DIFF0(L)` `  *$hh*0Diffusivity of soil layer L at low ` `  $hh*06reference water content for Gardner  tT  ` `  $hh*06and Mayhugh equation (m2/s). THETA0(L)` `  *$hh*0Reference water content for Gardner and ` `  $hh*06Mayhugh equation. Usually at wilting  tTx ` `  $hh*06point or 1.5 MPa (m3/m3). BETA(L)` `  *$hh*0Soil characteristic parameter of layer L ` `  $hh*06for GardnerDand Mayhugh equation. GLAYR(L)` `  *$hh*0Thickness of layer L (cm). THETAS(L)` `  *$hh*0Saturation water content of layer L. ` `  $hh*06Assumed to be same as water content  tT( ` `  $hh*06at the bubbling pressure. (m3/m3). FC(L)` `  *$hh*0Water content at "field capacity" of ` `  $hh*0633 kPa or other potential specified  tT ` `  $hh*06by PSISFC. (m3/m3). THETAR(L)` `  *$hh*0"Residual" water content for Brooks and ` `  $hh*06Corey equation, generally is ` `  $hh*06THETA0(L) slightly reduced,  tT  ` `  $hh*06typically by 0.001 (m3/m3).  tTh AIRDR(L)` `  *$hh*0Water content of air dry soil (m3/m3). ETA(L)` `  *$hh*0Soil characteristic parameter of layer L ` `  $hh*06for Brooks and Corey equation.  tT GBD(L)` `  *$hh*0Bulk density of layer L (g/cm3). PSIBUB(L)` `  *$hh*0Bubbling pressure or "saturation" soil ` `  $hh*06water potential of layer L at which ` `  $hh*06air starts to enter the soil as it ` `  $hh*06dries (kPa, positive number).  tTp Format for line 2*NLAYR+7+LYRSOL+1 Same as immediately above for recently cultivated soil, except the layer thickness, GLAYR, is omitted.  tTX Format for line 2*NLAYR+7+LYRSOL+2 Same as immediately above for soil in a wheel track, and again the layer thickness is omitted.  tT@ Format for line 2*NLAYR+7+LYRSOL+3 TD` `  *$hh*0Number of days after first bloom that ` `  $hh*06water content of bottom soil layer ` `  $hh*06begins to change (days). THETAI` `  *$hh*0Initial water content of bottom layer in  tT# ` `  $hh*06soil profile (m3/m3). BDSLOP` `  *$hh*0Rate of change of water content in bottom ` `  $hh*06soil layer after day of first bloom  tTH& ` `  $hh*06((m3/m3)/day). BDRATO` `  *$hh*0Ratio of the initial to the final water ` `  $hh*06water content of the bottom soil'T0*((@@Ԍ` `  $hh*06layer. PSISFC` `  *$hh*0Soil water potential used to define ` `  $hh*06field capacity (kPa, positive ` `  $hh*06number).  tT Format for line 2*NLAYR+7+LYRSOL+4 SNAME` `  20A4$hh*0Description of soil for root impedance ` `  $hh*06data.  tT Format for line 2*NLAYR+7+LYRSOL+5 NCURVE` `  *$hh*0Number of isowater content curves of ` `  $hh*06root impedance versus bulk density.  tT Format for next NCURVE sets of data INRIM( )` `  *$hh*0Number of root impedancebulk density ` `  $hh*06pairs of points on curve. GH2OC( )` `  *$hh*0Water content of a particular root ` `  $hh*06impedancebulk density curve ` `  $hh*06(kg/kg).  tT Format for next NRIM( ) pairs of data TSTBD( , )` `  *$hh*0Bulk density at a point on the root ` `  $hh*06impedancebulk density curve  tT ` `  $hh*06(Mg/m3). TSTIMP( , )` `  *$hh*0Root impedance at point on the root ` `  $hh*06impedancebulk density curve  tTp ` `  $hh*06(kg/cm2). U0*((@@  uS  FILE4: Soil Organic Residue Data INSTS` `  A2$hh*0Code for institute ID. SITES` `  A2$hh*0Code for site ID. YEAR` `  A2$hh*0Year number, last 2 digits. EXPTNO` `  I2$hh*0Experiment number. TRTNO` `  1X,I2$hh*0Treatment number. STRAW` `  1X,F5.0hh*0Weight of organic residue of previous ` `  $hh*06crop and/or added green manure ` `  $hh*06(kg/ha). SDEP` `  1X,F5.0hh*0Depth of surface residue incorporation ` `  $hh*06(cm). SCN` `  1X,F5.0hh*0C:N ratio of surface residue of previous ` `  $hh*06crop (kg C/kg N, default = 75.). ROOT` `  1X,F5.0hh*0Dry weight of root residue of previous ` `  $hh*06crop (kg/ha, default = 500.).  V0*((@@  uS  FILE5: Soil Profile Initial Conditions  tT Format for line 1 TRTNO` `  I2$hh*0Treatment number. INSTS` `  1X,A2$hh*0Code for institute ID. SITES` `  A2$hh*0Code for site ID. YEAR` `  A2$hh*0Year number, last 2 digits. EXPTNO` `  I2$hh*0Experiment number.  tT Format for all other lines DLAYR(L)` `  F6.0$hh*0Depth of layer L (cm, enter 1 to signal ` `  $hh*06endofdata). SW(L)` `  1X,F6.3hh*0Initial soil water content of soil layer L  tT  ` `  $hh*06(m3/m3). NH4(L)` `  1X,F4.1hh*0Initial soil ammonium in layer L ` `  $hh*06(mg elemental N/kg soil). NO3(L)` `  1X,F4.1hh*0Initial soil nitrate in layer L ` `  $hh*06(mg elemental N/kg soil). PH(L)` `  1X,F4.1hh*0Initial pH of soil layer L in a 1:1 ` `  $hh*06soil:water slurry. 0W0*((@@  uS  FILE6: Irrigation Management File Description. This file contains the information about irrigations applied to the crop. The first line contains the usual identification information. Then, each subsequent line contains date, irrigation amount, method, and placement information. The water is assumed to be applied as rectangular line source defined by the upperleft and lowerright coordinates. Using the plant row for  tTx the origin, a flood irrigation would be described by say 0,0 (X1,Z1) and  tT@ 100,0 (X2,Z2) for a 100 cm row spacing. The coordinates for a drip irrigation tube installed at a depth of 30 cm midway between the rows would be 50,30 and 50,30, thus effectively defining a line source. The coordinates for a furrow irrigation that floods approximately half the land surface (with the crop row at the top of the ridge) would be 25,0 and 75,0. A "1" should be placed at the end in the day number column to indicate endofdata. Variable Name FORTRAN Format0 Description   tT Format for line 1 TRTNO` `  1X,I2$hh*0Treatment number. INSTE` `  A2$hh*0Code for institute ID. SITEE` `  A2$hh*0Code for site ID. YEAR` `  I2$hh*0Year number, last 2 digits. EXPTNO` `  I2$hh*0Experiment number.  tT Format for line 2 and subsequent lines JDIRR(J)` `  I3$hh*0Day of Year (Julian) of the Jth ` `  $hh*06irrigation. AMTIRR(J)` `  F5.0$hh*0Effective amount of irrigation added, ` `  $hh*06(mm).pp<This is expressed in terms of ` `  $hh*06unit volume of water applied per  tT ` `  $hh*06unit of land area, ie (m3/m2)*1000. IRRCOD(J)` `  I2$hh*0Code for irrigation method as listed by ` `  $hh*06IBSNAT (1988). X1IRR(J)` `  F6.1$hh*0Distance from plant row to closest edge ` `  $hh*06of water application area (cm). Z1IRR(J)` `  F6.1$hh*0Topmost depth in soil where water is ` `  $hh*06applied (cm). For flood, sprinkler, ` `  $hh*06or furrow irrigation, the ` `  $hh*06appropriateDentry would be zero. X2IRR(J)` `  F6.1$hh*0Distance from plant row to furthest edge ` `  $hh*06of water application area (cm). Z2IRR(J)` `  F6.1$hh*0Lowest depth in soil where water is ` `  $hh*06applied (cm). $X0*((@@  uS  FILE7: Fertilizer Management File Description. This is the file containing the information about fertilizer applications. The first line contains the usual identification information. Then, each subsequent line contains date, fertilizer amount, type, and placement information.  tTx The amount and typesource code for N, P, and K fertilizers and other inoculants and amendments are all supplied. The typesource code is from IBSNAT Technical Report 1, 3rd ed. (1988). The fertilizer is assumed to be applied in a rectangular band defined by the upperleft and lowerright coordinates. Using the plant row for the origin, a  tT( broadcast application would be described by say 0,0 (X1,Z1) and 100,0 (X2,Z2) for a 100 cm row spacing. If this broadcast fertilizer is incorporated into the top 10 cm of soil, then the coordinates would be 0,0 and 100,10. The coordinates for a band 10 cm from the row and 10 cm deep would be 10,10 and 10,10. The coordinates for a sidedress application on the soil surface extending 10 cm from the row to 30 cm from the row would be 10,0 and 30,0. A "1" should be placed at the end in the day number column to indicate endofdata. Variable Name FORTRAN Format0 Description   tT Format for line 1 TRTNO` `  1X,I2$hh*0Treatment number. INSTE` `  A2$hh*0Code for institute ID. SITEE` `  A2$hh*0Code for site ID. YEAR` `  I2$hh*0Year number, last 2 digits. EXPTNO` `  I2$hh*0Experiment number.   tT Format for line 2 and subsequent lines JDFERT(J)` `  I3$hh*0Day of Year (Julian) of the Jth fertilizer ` `  $hh*06application.  tT FERTN(J)` `  F5.1$hh*0Effective amount of nitrogen in ` `  $hh*06fertilizer added, (kg/ha). DFERT(J)` `  F5.1$hh*0Depth of incorporation of fertilizer (cm). INTYPE(J)` `  I3$hh*0Code number for nitrogen fertilizer.  tT FERTP(J)` `  F5.1$hh*0Effective amount of phosphorus in ` `  $hh*06fertilizer added, (kg/ha). IPTYPE(J)` `  I3$hh*0Code number for phosphorus fertilizer.  tT`" FERTK(J)` `  F5.1$hh*0Effective amount of potassium in ` `  $hh*06fertilizer added, (kg/ha). IKTYPE(J)` `  I3$hh*0Code number for potassium fertilizer. FERTIN(J)` `  F5.1$hh*0Effective amount of "other" inoculants ` `  $hh*06and amendments. IITYPE(J)` `  I3$hh*0Code for "other" inoculants and T(#(#Z` `  $hh*06amendments. FERCOD(J)` `  I2$hh*0Code for application placement as 'Y0*((@@Ԍ` `  $hh*06listed by IBSNAT (1988). X1FERT(J)` `  F6.1$hh*0Distance from plant row to closest edge ` `  $hh*06of fertilizer application area (cm). Z1FERT(J)` `  F6.1$hh*0Topmost depth in soil where fertilizer is ` `  $hh*06applied (cm). X2FERT(J)` `  F6.1$hh*0Distance from plant row to furthest edge ` `  $hh*06of fertilizer application area (cm). Z2FERT(J)` `  F6.1$hh*0Lowest depth in soil where fertilizer is ` `  $hh*06applied (cm). Z0*((@@  uS  FILE8: Treatment Management Data  Variable Name FORTRAN Format0 Description   tT  Format for line 1 INSTE` `  A2$hh*0Code for institute ID. SITEE` `  A2$hh*0Code for site ID. YEAR` `  I2$hh*0Year number, last two digits. EXPTNO` `  I2$hh*0Experiment number. TRTNO` `  1X,I2$hh*0Treatment number. TITLET` `  1X,A40hh*0Title of Treatment. ISOILT` `  1X,I4$hh*0Soil number for this treatment. IVARTY` `  1X,I4$hh*0Cultivar number for this treatment.  tT Format for line 2 of each treatment ISIM` `  I4$hh*0Julian date (day of year) simulation ` `  $hh*06begins. ISOW` `  1X,I3$hh*0Sowing date (Julian day of year).  tT PLANTS` `  1X,F6.2hh*0Plant population (plants/m2). ROWSPC` `  1X,F6.3hh*0Row spacing (m). SDEPTH` `  1X,F5.2hh*0Sowing depth (cm). IIRR` `  1X,I2$hh*0Switch describing irrigation (default=1) ` `  $hh*061. no irrigation applied. ` `  $hh*062. irrigation applied according ` `  $hh*06pp<to actual field schedule. ` `  $hh*063. automatically irrigated at ` `  $hh*06pp<threshold soil water. ` `  $hh*064. assume no water stress, water ` `  $hh*06pp<balance not used. ISWNIT` `  1X,I2$hh*0Switch to indicate if nitrogen routines ` `  $hh*06are used (default=0). ` `  $hh*060. nitrogen routines are not used ` `  $hh*06pp<and assume adequate N. ` `  $hh*061. nitrogen routines are used. EFFIRR` `  1X,F6.2hh*0Irrigation efficiency (fraction). DSOIL` `  1X,F5.2hh*0Irrigation management depth (m). THETAC` `  1X,F6.1hh*0Available water triggering irrigation (%). IMERGE` `  1X,I3$hh*0Emergence date (Julian day of year). ISWEED` `  1X,I2$hh*0Switch to indicate weed data are available ` `  $hh*06(0=No). ISWINS` `  1X,I2$hh*0Switch to indicate insect data are ` `  $hh*06available (0=No). ISWNEM` `  1X,I2$hh*0Switch to indicate nematode data are ` `  $hh*06available (0=No). ISWDIS` `  1X,I2$hh*0Switch to indicate disease data are ` `  $hh*06available (0=No).  tTH& Format for line 3 of each treatment HISTRY` `  1X,A76hh*0Brief description of previous cropping ` `  $hh*06history for site.'[0*((@@Ԍ uS  FILEA: Measured Cotton Crop Harvest Summary Data Description. FILEA is a small file which contains a summary of measured field harvest data for the crop for a particular treatment averaged over all replications. The following proposed format for cotton follows the IBSNAT formats for maize, wheat, and soybean as closely as possible. The "X" for the initial letter of variable names indicates experimentally observed data.  Variable Name FORTRAN Format0 Description   tT Format for line 1 INSTE` `  A2$hh*0Code for institute ID. SITEE` `  A2$hh*0Code for site ID. YEAR` `  I2$hh*0Year number, last two digits. EXPTNO` `  I2$hh*0Experiment number. TRTNO` `  1X,I2$hh*0Treatment number. XLTYLD` `  1X,F7.0hh*0Actual fieldmeasured lint yield ` `  $hh*06(dry weight basis, kg/ha). This ` `  $hh*06is the yield attainable by hand ` `  $hh*06harvesting. If the cotton was ` `  $hh*06harvested by machine, then the ` `  $hh*06data must first be corrected for ` `  $hh*06"gin turnout" or "harvest ` `  $hh*06efficiency". XSDYLD` `  1X,F7.0hh*0Actual fieldmeasured seed yield ` `  $hh*06(dry weight basis, kg/ha). XSDWT` `  1X,F7.4hh*0Measured seed dry weight (g/seed).  tT XBLSM` `  1X,F6.0hh*0Fieldmeasured boll number (bolls/m2). XSPB` `  1X,F4.0hh*0Fieldmeasured number of seeds per boll ` `  $hh*06(seeds/boll). XLAIMX` `  1X,F5.2hh*0Maximum leaf area index during season  tT ` `  $hh*06(m2/m2). XBIOM` `  1X,F6.0hh*0Fieldmeasured aboveground dry biomass at ` `  $hh*06maturity (kg/ha). XSTMBR` `  1X,F6.0hh*0Fieldmeasured stem plus burr weight at ` `  $hh*06maturity (kg/ha).  tT Format for line 2 XSDTN` `  F6.2$hh*0Measured nitrogen concentration in seed ` `  $hh*06at maturity (%). XTOTNP` `  1X,F5.1hh*0Measured total crop nitrogen content at ` `  $hh*06maturity (kg/ha). XAPTNP` `  1X,F5.1hh*0Measured nitrogen content of stems plus (#(#Z` `  $hh*06burrs plus lint at maturity (kg/ha). XSDN` `  1X,F5.1hh*0Measured nitrogen content of seed at ` `  $hh*06maturity (kg/ha). %\0*((@@  uS  FILEB: Measured Cotton Intermediate Growth Data (for graphics) Description. FILEB contains observed crop growth data obtained at intermediate days during the growing season. Data for all treatments of an experiment are stored in one file. The first line is header information which identifies the source of the data followed by days of year at which emergence, first square, and first flower occurred. On the second line, the first variable identifies the number of state variables for which there are matching field data, and the rest of the variables on this line are pointers which indicate the state variable for each column of data. Starting with the third line, there is one line of data for each observation date recorded in the minimum data set (MDS). Replication data for each treatment can be included by using a different line of data for each rep. A "1" on the line immediately following a data line indicates the end of data for a specific treatment. The "X" or "J" for the first letter of the variable names indicates experimentally observed data. Variable Name FORTRAN Format0 Description   tT Format for line 1 INSTE` `  * or A2hh*0Code for institute ID. SITEE` `  * or A2hh*0Code for site ID. YEAR` `  * or I2hh*0Year number, last two digits. EXPTNO` `  * or I2hh*0Experiment number. TRTNO` `  * or 1X,I2hh*0Treatment number. JEMRGD` `  * or 1X,I3hh*0Day of emergence which is the day ` `  $hh*06that 50% of the plants ` `  $hh*06emerge from the soil ` `  $hh*06(day of year). JSQRJD` `  * or 1X,I3hh*0Day of first square which is the day ` `  $hh*06of year that 50% of the plants ` `  $hh*06in the field displayed their ` `  $hh*06first square (day of year). JFLRJD` `  * or 1X,I3hh*0Day of first flower which is the ` `  $hh*06day of yearDthat 50% of the ` `  $hh*06plants in the field displayed ` `  $hh*06their first flower (day of ` `  $hh*06year). ]0*((@@  tT Format for line 2 NOVAR` `   * or I3hh*0No. of state variables for which ` `  $hh*06there are matching field data. NV(I),I=1,NOVAR  * or (NOVAR)I30Pointer which indicates state ` `  $hh*06variable number for each column ` `  $hh*06pp<of data  tT@ Format for line 3 and beyond JDOY()` `   * or I3hh*0Day of Year (Julian). NV(1)or XPLTHT()  * or F5.0hh*0Plant height (cm).  tT( NV(2)or XLAI()  * or F6.2hh*0Leaf Area Index (m2/m2).  tT NV(3)or JNNODM()  * or I5hh*0No. nodes/m2.  tT NV(4)or JNSQRM()  * or I4hh*0No. squares/m2.  tT NV(5)or JNFLWM()  * or I4hh*0No. flowers/m2.  tTH NV(6)or JNGBLM()  * or I4hh*0No. green bolls/m2.  tT NV(7)or JNMBLM()  * or I4hh*0No. open mature bolls/m2.  tT NV(8)or JNABSM()  * or I5hh*0No. abscised sites/m2. NV(9)or XWLEFH()  * or F8.0hh*0Dry leaf weight (kg/ha). NV(10) or XWSTMH() * or F8.0hh*0Dry stem weight (kg/ha). NV(11) or XWROTH() * or F8.0hh*0Dry root weight (kg/ha). NV(12) or XWGBLH() * or F8.0hh*0Dry green boll weight (kg/ha). NV(13) or XWMBLH() * or F8.0hh*0Dry mature boll weight (kg/ha). NV(14) or XWLINT() * or F8.0hh*0Dry lint weight (kg/ha). NV(15) or XWSEED() * or F8.0hh*0Dry seed weight (kg/ha). NV(16) or XWBURR() * or F8.0hh*0Dry burr weight (kg/ha). _______________ * Format uses one or more spaces to separate one variable from the next. p^0*((@@  uT n( Treatment Sorting Order ă In the data tables that follow, the presentation sequence follows treatment order, rather than plot number order. As can be seen from a scan of the Table  tT of Contents, the treatments are listed in the order CO2 irrigation  tTB nitrogen replication, with replication cycling fastest and CO2 slowest. We  tT used opentop field chambers whose environment might be affected by the CO2, irrigation, or nitrogen treatment. Therefore, each plot was regarded as a separate site and was outfitted with a psychrometer to measure temperatures, as discussed previously. The sites were identified by plot codes in 1983, as shown in Fig. 1 (2B, 3A, etc.), or sequence numbers in 19841987, as shown in Figs. 23 (01, 02, ..., 16), and these codes or sequence numbers were used for the site ID, SITEE (and SITES for FILE4 and FILE5). Each site had its own weather file each year, as well as treatment management (FILE8), final harvest summary (FILEA), and intermediate growth (FILEB) files. However, the data for some of the files were common to all plots or sites that year. To save space and needless duplication, therefore, in such cases the data for all plots is in the file named for the site with the lowest level of all variables. For example, in 1983 all plots received the same fertilizer  tT management, and these data are in the file for the lowest CO2 level (open field) and rep I, which was plot 4A (Fig. 1). Thus, the fertilizer management data for all plots in 1983 are in file PX4A8301.CT7. Another example is the soil organic residue data for 1986, which also was the same for all plots. In  tT# this case the lowest treatment level would be ambient CO2, dry, N, rep I, which was in plot 07 (or ID3 in Fig. 3), and these soil organic residue data are in file PX078601.CT4. In any event, the experiment file directories'_0*((@@ (CTEXP__.DIR) list the proper group of files needed for each yearplot combination that defines the individual experiments.  uT@  " Spreadsheet Template for Data Entry ă As can be seen from the Format Tables, the IBSNAT (1986) standard prescribes numerous files with detailed naming and format conventions. To aid data entry in the correct sequence and format, and to assist in properly naming all the files, a spreadsheet template was written. It was written using macros for Quatro but should work with Lotus 123 also. It is available upon request from the authors. `0*((@@ # DOCUMENTATION OF PROGRAM RETRVE 3 &FOR READING THE DATA FILES 3  An algorithm named 'RETRVE.FOR' was developed to read in all the data files, i.e., daily weather, soil profile properties, soil organic residue, soil profile initial conditions, irrigation management, fertilizer management, treatment management, measured cotton crop harvest summary, and measured cotton intermediate growth data, which correspond to a particular treatment combination for a given year. A PCbased OTG ANSI FTN77 software package (University of Salford FTN77/386 revision C, P.O. Box 5250, 308 Mulberry Street, Scranton, PA 18505) was used as the developmental tool; therefore, all code complies with the ANSI standard and should be transportable to other FORTRAN compilers. Wherever possible, RETRVE.FOR was structured to adhere to the IBSNAT standard as closely as possible, as discussed in the 'File Structure and Data Formats' section of this document. Structure of RETRVE.FOR: RETRVE.FOR, contains a main routine, MAIN, and three subroutines, LISTID, READDA and VALFAB. MAIN will prompt the user to verify that the include file, 'PATHINP.INC', has been edited correctly. It then reads this file into the appropriate variables. MAIN calls subroutine LISTID which will use the information in the PATHINP.INC file to identify which files are to be retrieved. It also sets the path to those files. Subroutine LISTID then calls subroutine READDA. This routine will open both the files to be read and the output file, RETRVE.OUT. All files which contain variables to be interfaced with the simulation model are read in consecutive order by READDA. Files A and B, which are validation data files, are processed in a separate'a0*((@@ subroutine, VALFAB. All data read in are echoed into an output file, 'RETRVE.OUT'. File B data are loaded into a multidimensional array. An array indexing technique is employed to load the individual state variable arrays according to identification array, NV(I). These single dimensional arrays along with the JDOY array provide real world observed data which can be plotted along with simulated values from a plant growth model to validate the model's predictions. PATHINP.INC file: The 'CTEXPyy.DIR' files contain a list of all the file names associated with the results from a particular year's experiment. These files, therefore are the key to the database. An input file named PATHINP.INC is the control input file used by RETRVE.FOR to retrieve the desired files within the database. As illustrated in Table 2, it contains four lines. Each line corresponds to the following variables in the RETRVE.FOR algorithm. (1) EXPID an 8 digit character string which identifies the experimental data set desired for a particular treatment combination within a year, ex. 'PX158701'. (2) DIRFN sets the path and gives the file name of the desired year 'CTEXPyy.DIR'. For example, if you desire to extract data files from a 1987 experiment using your floppy drive A:, then set DIRFN to 'A:\CTEXTP87.DIR'; otherwise, if you have loaded the data files onto your hard drive, then simply set the path to the proper directory. For example, a subdirectory used to hold your FORTRAN data files could be used 'C:\FTN77\DATA87\CTEXTP87.DIR'. (3) PTODAT sets the path to the directory which contains the actual data files for a given EXPID and year. This path must be set similarly to that in the previous example, i.e. 'A:\' for files located in floppy drive A:, and 'C:\FTN77\DATA87\' for the data files located'b0*((@@ in a \FTN77\DATA87\ subdirectory of your hard drive, C:, (4) NPATH gives the number of characters used in the character string which sets the path in variable PTODAT; for example, if the path is set to the floppy drive A:\, 'A:\', then NPATH is 3. If it is set to 'C:\FTN77\DATA87\', then NPATH is set to 16. It is important that NPATH be set correctly so that the name of each file to be opened can be concatenated to the existing path, thereby enabling those files to be opened and eventually read by RETRVE.FOR. Note!, NPATH must not exceed 45 which limits the character string used to set the paths in both the DIRFN and PTODAT variables. Running RETRVE.FOR: After editing the PATHINP.INC file, and compiling and linking the source code in RETRVE.FOR, RETRVE is ready for execution. During a run it will prompt the user to verify that the PATHINP.INC file is set up correctly. The user should enter '1' if PATHINP.INC is correct; otherwise enter '2' to discontinue and edit PATHINP.INC. During execution, the complete path and file names being processed are echoed to the screen. This will allow verification that the correct files are being loaded.  Output from RETRVE.FOR: Execution RETRVE will generate an output file called RETRVE.OUT. This file is output in the same directory that RETRVE is executed. It consists of nothing more than an echo of the input files read. This is a further check to verify that the actual data is the same as that in the text portion of the document. Interfacing algorithm RETRVE.FOR with a cotton simulation model: TThe'c0*((@@ RETRVE.FOR source code has been developed so that it can be readily interfaced with any cotton simulation model. After each of the input files has been opened in subroutine READDA, the variables within the file are read and can be interfaced directly with the corresponding variables in the cotton model being used, as illustrated by tables after each of the READ statements in READDA. The original source code makes these interface tables nonexecutable by having C's in the first column to make then comments. To use these interface tables, simply insert the variable names used in your model which correspond to the input variable names used in the IBSNAT standard. Check for consistency in units. If the units differ, a conversion factor can be used to scale the variables properly. Lastly, remove the C's so that the comments become executable statements, recompile, link, and run the executable program. Pd0*((@@  ` `  Table 2. example of a PATHINP.INC file ` `  $used to extract data files located ` `  $on a floppy inserted into drive A: (a), ` `  $or on a subdirectory in the ` `  $harddrive C: (b), respectively. ` `  ___________________________________________ ` `  a) ` `  $PX158701 ` `  $A:\CTEXTP87.DIR ` `  $A:\ ` `  $3 ` `  ___________________________________________ ` `  b) ` `  $PX158701 ` `  $C:\FTN77\DATA87\CTEXP87.DIR ` `  $C:\FTN77\DATA87\ ` `  $16 ` `  ___________________________________________ e0*((@@ CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C DOCUMENTATION OF PROGRAM 'RETRVE.FOR' C C C C----------------------------------------------------------------------C C C C DATE: 2/27/92 C C VERSION: 1.0 C C AUTHORS: Mr. Matthew L. Reaves and Dr. Gerard W. Wall C C C C ADDRESS: USDA/ARS C C U.S. Water Conservation Laboratory C C Phoenix, AZ 85040 C C C C----------------------------------------------------------------------C C C C DESCRIPTION: RETRVE.FOR IS AN ALGORITHM DEVELOPED TO C C READ EXPERIMENTAL DATASETS FOUND IN THE C C ENCLOSED DOCUMENT C C C C C C STRUCTURE: RETRVE - MAIN ROUTINE C C LISTID - SUBROUTINE CALLED BY MAIN C C READDA - SUBROUTINE CALLED BY LISTID C C VALFAB - SUBROUTINE CALLED BY READDA C C C C----------------------------------------------------------------------C C C C INPUT: A FOUR LINE FILE CALLED 'PATHINP.INC' WHICH C C CONTAINS THE FOLLOWING VARIABLES: C C C C (1) EXPID - AN 8 DIGIT CHARACTER STRING WHICH C C IDENTIFIES THE DESIRED EXPERIMENTAL C C DATA SET (ex. 'PX158701'). C C C C (2) DIRFN - SETS PATH TO ONE OF THE FOLLOWING C C EXPERIMENTAL DIRECTORY FILES. C C C C 2.1) 'CTEXP83.DIR' C C 2.2) 'CTEXP84.DIR' C C 2.3) 'CTEXP85.DIR' C C 2.4) 'CTEXP86.DIR' C C 2.5) 'CTEXP87.DIR' C C C C (3) PTODAT - SETS PATH TO THE C C DESIRED EXPERIMENTAL C C DATA FILES. C C C C (4) NPATH - SETS THE NUMBER OF CHARACTERS USED C C IN SETTING PATH GIVEN IN VARIABLE C C PTODAT C C C'f0*((@@ԌC----------------------------------------------------------------------C C C C OUTPUT: A DATA FILE 'RETRVE.OUT' WHICH ECHOS THE C C INPUT TO VERIFY CORRECT INPUT C C C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC PROGRAM RETRVE CHARACTER EXPID*8, DIRFN*45, PTODAT*45 INTEGER NPATH WRITE(6,'(//,A60)')'PROGRAM * RETRIEVE * # ','******************************************************' WRITE(6,'(A60)')'PLEASE READ THE FOLLOWING # ','****************************************************** # ',' # ','THIS ALGORITHM WILL RETRIEVE THE REQUESTED DATA FILES # ','ASSOCIATED WITH A SPECIFIC TREATMENT COMBINATION WITHIN # ','A GIVEN YEAR USING INFORMATION SUPPLIED IN THE # ','"PATHINP.INC" FILE. # ',' # ',' ***************************************************** # ',' # ',' NOTE ! # ',' **** # ',' # ',' - THE USER MUST HAVE ALREADY EDITED # ',' THE "PATHINP.INC" FILE . # ',' # ',' - THIS IS A 4 RECORD INPUT FILE TO TELL RETRVE.FOR # ',' WHERE TO LOCATE THE REQUESTED DATA FILES. # ',' # ',' - ENTER A "1" TO CONTINUE OR A "2" TO STOP. # ','****************************************************** # ' READ(6,*)ICHECK IF(ICHECK.NE.1) THEN STOP ENDIF OPEN(38, FILE='PATHINP.INC') CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C 'PATHINP.INC' = PATH LIST INPUT FILE NAME C C CONTAINS 4 RECORDS C C C C EXPID - EXPERIMENTAL IDENTIFIER C C DIRFN - FULL PATH NAME TO CTEXP.DIR FILE C C PTODAT - FULL PATH TO DATA FILES C'g0*((@@ԌC NPATH - NUMBER OF CHARACTERS IN PATH TO DATA FILES C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC READ(38,'(A8)')EXPID READ(38,'(A45)')DIRFN READ(38,'(A45)')PTODAT READ(38,*)NPATH C======================================================================= C CALL SUBROUTINE LISTID C C======================================================================= CALL LISTID(EXPID,DIRFN,PTODAT,NPATH) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C PROGRAM COMPLETE, TELL USER WHERE TO VIEW DATA C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC WRITE(6,'(/)') PRINT*,'**********************************************' PRINT*,'* *' PRINT*,'* PROGRAM HAS EXECUTED CORRECTLY *' PRINT*,'* *' PRINT*,'**********************************************' WRITE(6,'(/)') PRINT*,'**********************************************' PRINT*,'* *' PRINT*,'* EDIT "RETRVE.OUT" TO REVIEW THE DATA *' PRINT*,'* *' PRINT*,'**********************************************' END CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C SUBROUTINE LISTID - THIS ROUTINE WILL USE INFORMATION IN THE C C 'PATHINP.INC' INPUT FILE TO IDENTIFY C C WHICH FILES ARE TO BE RETRIEVED AND TO C C SET THE PATH TO THESE FILES C C C C - CALLED BY MAIN C C C C - CALLS READDA C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC SUBROUTINE LISTID(EXPID,DIRFN,PATH1,NPATH) CHARACTER EXPID*8,VARSEE*8,EXPDES*42,FIL1*12'h0*((@@Ԍ # , FIL2*12, FIL4*12,FIL5*12,FIL6*12,FIL7*12 # , FIL8*12,FIL9*12,FILA*12, FILB*12 # , TEMP*45,FPATH1*45,PATH1*45 # , FPATH2*45,FPATH4*45,FPATH5*45,FPATH6*45 # , FPATH7*45,FPATH8*45,FPATHA*45,FPATHB*45 # , DIRFN*45 INTEGER NPATH C======================================================================= C A CHARACTER VARIABLE FILE NAME IS OPENED C THAT WAS READ FROM THE FILE 'PATHINP.INC' C======================================================================= OPEN(10,FILE= DIRFN) C======================================================================= C THIS LOOP CHECKS EACH 'EXPID' FOR A MATCH C WITH THE CHOSEN TEST C======================================================================= DO 1500 J=1,500 READ(10,'(A8,A42,A12,1X,A12)',END=230)VARSEE # ,EXPDES,FIL1,FIL2 IF (VARSEE.EQ.EXPID)THEN C======================================================================= C ONCE A MATCH IS FOUND ALL FILE NAMES FOR THAT TEST SET C ARE READ INTO CHARACTER VARIABLES C======================================================================= READ(10,'(5(A12,1X),A12)')FIL4,FIL5,FIL6,FIL7,FIL8,FIL9 READ(10,'(A12,1X,A12)')FILA,FILB C====================================================================== C C C FILE NAMES TO BE READ ARE ADDED C C TO PATH GIVEN IN "PATHINP.INC" FILE C C C C====================================================================== WRITE(6,'(/)') TEMP=FIL1 FPATH1=PATH1 FPATH1(NPATH+1:)=TEMP PRINT*,FPATH1 TEMP=FIL2 FPATH2=PATH1 FPATH2(NPATH+1:)=TEMP PRINT*,FPATH2 TEMP=FIL4 FPATH4=PATH1 FPATH4(NPATH+1:)=TEMP PRINT*,FPATH4 TEMP=FIL5'i0*((@@Ԍ FPATH5=PATH1 FPATH5(NPATH+1:)=TEMP PRINT*,FPATH5 TEMP=FIL6 FPATH6=PATH1 FPATH6(NPATH+1:)=TEMP PRINT*,FPATH6 TEMP=FIL7 FPATH7=PATH1 FPATH7(NPATH+1:)=TEMP PRINT*,FPATH7 TEMP=FIL8 FPATH8=PATH1 FPATH8(NPATH+1:)=TEMP PRINT*,FPATH8 TEMP=FILA FPATHA=PATH1 FPATHA(NPATH+1:)=TEMP PRINT*,FPATHA TEMP=FILB FPATHB=PATH1 FPATHB(NPATH+1:)=TEMP PRINT*,FPATHB C======================================================================= C CALL SUBROUTINE READDA C C======================================================================= CALL READDA(EXPID,EXPDES,FPATH1,FPATH2,FPATH4,FPATH5,FPATH6 # ,FPATH7,FPATH8, FPATHA,FPATHB) GOTO 23 ENDIF 1500 CONTINUE 230 PRINT*,'THIS EXPERIMENT ',EXPID,' WAS NOT FOUND' STOP 23 PRINT*,' SUBROUTINE LISTID COMPLETE' RETURN END CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C SUBROUTINE READDA - OPENS THE INPUT FILES TO BE READ C C WHICH WERE DETERMINED IN SUBROUTINE C C LISTID. OPENS THE OUTPUT FILE 'RETRVE.OUT'. C C READS THE INPUT FILES AND ECHOS THE OUTPUT. C C IT IS ALSO THE SUBROUTINE WHERE INTERFACING C C THE IBSNAT STANDARD VARIABLES WITH THOSE C C IN THE COTTON MODEL USED WILL OCCUR. C C C C - CALLED BY LISTID C C C'j0*((@@ԌCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC SUBROUTINE READDA(EXPID,EXPDES,FILE1,FILE2,FILE4,FILE5,FILE6 # ,FILE7,FILE8,FILEA, FILEB) REAL LL,DFERT,NH4,NO3,PH INTEGER CO2,IDUMSL,TRTNO,YEAR # ,PARDAT,WNDDAT, DEWDAT, STDAT, CO2DAT,EXPTNO,FERCOD CHARACTER INSTW*2,STATW*2,A00*7 ,EXPID*8, EXPDES*42 # , INSTS*2,SITES*2, INSTE*2,SITEE*2 # , PEDON*12,TAXON*60 # , DATAID*60 # , SNAME*60 # , HISTRY*76 # , FILE1*45,FILE2*45,FILE4*45,FILE5*45,FILE6*45 # , FILE7*45,FILE8*45,FILEA*45,FILEB*45 # , TITLET*40 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C OPEN INPUT AND OUTPUT FILES C C C C EXAMPLE: OPEN (12,FILE ='D:\FTN77\AVONDALE.CT2') C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC OPEN (3,FILE = 'RETRVE.OUT') OPEN (11,FILE = FILE1) OPEN (12,FILE = FILE2) OPEN (14,FILE = FILE4) OPEN (15,FILE = FILE5) OPEN (16,FILE = FILE6) OPEN (17,FILE = FILE7) OPEN (18,FILE = FILE8) WRITE(3,*)'THE EXPERIMENTAL IDENT. EXPID= ',EXPID WRITE(3,*)'THE EXPERIMENTAL DESCRIPTION= ',EXPDES WRITE(3,*) WRITE(3,'(A72)')'************ FILE 1 ***********************' # , 'THE FOLLOWING WILL BE FROM THE WEATHER FILE = ' # ,FILE1 # ,' *********************************************' WRITE(3,*) PRINT*,'READING FILE 1 ' READ(11,1111)INSTW,STATW,XLAT,XLONG,PARFAC,PARDAT,WNDDAT # ,DEWDAT,STDAT,CO2DAT,CO2YR,WINDYR CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE OF WEATHER DATA - FILE1 C'k0*((@@ԌC C C *********************************************************** C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C *********************************************************** C C C C = INSTW C C = STATW C C = XLAT C C = XLONG C C = PARFAC C C = PARDAT C C = WNDAT C C = DEWDAT C C = STDAT C C = CO2DAT C C = CO2YR C C = WINDYR C C C C ************************************************************ C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC WRITE(3,*)'INSTW/ OPTIONS (PARDAT,WNDDAT,DEWD #AT,STDAT,CO2DAT)' WRITE(3,*)'STATW XLAT XLONG PARFAC CO2YR WINDYR' WRITE(3,1112)INSTW,STATW,XLAT,XLONG,PARFAC,PARDAT,WNDDAT # ,DEWDAT,STDAT,CO2DAT,CO2YR,WINDYR 1111 FORMAT(2(A2),2(1X,F6.2),1X,F5.2,5(1X,I1),1X,F4.0,1X,F5.1) 1112 FORMAT(2(A2),2(1X,F6.2),1X,F5.2,5(1X,I1),1X,F4.0,1X,F5.1,/) WRITE(3,*)'STATW IYR SOLRAD XTMIN XPAR DEWP #T STMIN A00' WRITE(3,*)'INSTW/ JUL XTMAX XRAIN WIND # STMAX CO2 ' CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C FILL APPROPRIATE COTTON MODEL ARRAYS TO END OF FILE C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC DO 678 I=1,360 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C READ EVERYTHING AS IS. AND ENTER INTO APPROPRIATE C C COTTON MODEL ARRAY AFTER EACH LINE IS READ C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC READ(11,1120,END=1999,ERR=193)INSTW,STATW,IYR,JUL,SOLRAD,XTMAX # ,XTMIN,XRAIN,XPAR,WIND,DEWPT,STMAX,STMIN,CO2,A00'l0*((@@Ԍ 1120 FORMAT(2(A2),1X,I2,1X,I3,1X,F5.2,3(1X,F5.1),1X,F6.2 # ,4(1X,F5.1),1X,I4,1X,A7) WRITE(3,1120)INSTW,STATW,IYR,JUL,SOLRAD,XTMAX # ,XTMIN,XRAIN,XPAR,WIND,DEWPT,STMAX,STMIN,CO2,A00 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE OF DAILY WEATHER DATA - FILE1 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C C C = INSTW C C = STATW C C = IYR C C = JUL C C = SOLRAD C C = XTMAX C C = XTMIN C C = XRAIN C C = XPAR C C = WIND C C = DEWPT C C = STMAX C C = STMIN C C = CO2 C C = A00 C C *********************************************************** C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 678 CONTINUE 193 WRITE(3,*)'ERROR PROCESSING FILE 1' WRITE(6,*)'ERROR PROCESSING FILE 1' STOP 1999 PRINT *,'LAST DAY OF DATA WAS ',A00 WRITE(3,'(/,A8,/,A45,/,8X,A25)')'FILE 1= # ',FILE1,'WAS PROCESSED CORRECTLY' WRITE(3,'(//,A50)')'********FILE 2********************** # ' WRITE(3,'(A50)')' THE FOLLOWING WILL BE FROM # ',' THE SOIL PROPERTIES = FILE 2 # ',FILE2 # ,'******************************************* # ' WRITE(3,*)'IDUMSL PEDON TAXON''m0*((@@Ԍ PRINT*,'READING FILE 2 ' READ(12,121)IDUMSL,PEDON,TAXON CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE OF SOIL PROFILE PROPERTIES C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************* C C C C = IDUMSL C C = PEDON C C = TAXON C C ********************************************************** C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC READ(12,122)SALB,U,SWCON,CN2,TAV,AMP,DMOD # ,SWCON1,SWCON2,SWCON3,RWUMX,PHFAC3 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE OF SOIL PROFILE PROPERTIES C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************* C C C C = SALB C C = U C C = SWCON C C = CN2 C C = TAV C C = AMP C C = DMOD C C = SWCON1 C C = SWCON2 C C = SWCON3 C C = RWUMX C C = PHFAC3 C C C C ************************************************************ C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 121 FORMAT(1X,A2,1X,A12,1X,A60) 122 FORMAT(F6.2,1X,F5.2,2(1X,F6.2),2(1X,F5.1),1X,F3.1 # ,1X,E9.2,1X,F6.1,2(1X,F5.2),1X,F4.2)'n0*((@@Ԍ WRITE(3,121)IDUMSL,PEDON,TAXON WRITE(3,*) WRITE(3,*)'SALB U SWCON CN2 TAV AMP DMOD SWCON1 SWCON2 # SWCON3 RWUMX PHFAC3' WRITE(3,126)SALB,U,SWCON,CN2,TAV,AMP,DMOD # ,SWCON1,SWCON2,SWCON3,RWUMX, PHFAC3 126 FORMAT(F5.2,1X,F5.1,F6.2,F5.1,2(1X,F5.1),1X,F3.1 # ,1X,E10.3,1X,F6.1,2(1X,F5.2),1X,F4.2) WRITE(3,*) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INPUT ARRAY DATA FOR LAYERS C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC INLAYS=0 WRITE(3,2129)'DLAYR','LL','DUL','SAT','DSW' # ,'WR','BD','OC','DNH4','DNO3','DPH' 2129 FORMAT(A6,5(1X,A6),2(1X,A5),3(1X,A4)) DO 300 IL = 1,20 INLAYS=INLAYS+1 READ(12,123,ERR=310)DLAYR, LL, DUL, SAT, DSW # ,WR,BD,OC,DNH4 ,DNO3,DPH IF(INT(DLAYR).NE.-1)THEN WRITE(3,123)DLAYR, LL, DUL, SAT, DSW # ,WR,BD,OC,DNH4,DNO3,DPH 123 FORMAT(F6.0,5(1X,F6.3),2(1X,F5.2),3(1X,F4.1)) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE SOIL PROFILE PROPERTIES - FILE 2 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C C C = DLAYR C C = LL C C = DUL C C = SAT C C = DSW C C = WR C C = BD C C = OC C C = DNH4 C C = DNO3 C C = DPH C C ************************************************************* C'o0*((@@ԌC C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC ELSE WRITE(3,'(F6.0)') DLAYR GOTO 333 ENDIF 300 CONTINUE 310 WRITE(6,*)' NLAYERS =',INLAYS,' ERROR' WRITE(3,*)' NLAYERS =',INLAYS,' ERROR' STOP 333 PRINT*,'CORRECT PROCESSING' C======================================================================= C FILE 2 SECOND LOOP C======================================================================= WRITE(3,*) WRITE(3,*)'ROCK SAND SILT CLAY SCOND CATEXC ALPHA XN VGTHS #VGTHR' DO 234 I=1,INLAYS-1 READ(12,124)ROCK,SAND,SILT, CLAY,SCOND,CATEXC # ,ALPHA,XN,VGTHS,VGTHR WRITE(3,124)ROCK,SAND,SILT, CLAY,SCOND,CATEXC # ,ALPHA,XN,VGTHS,VGTHR 124 FORMAT(5(1X,F4.1),1X,F5.1,4(1X,F6.3)) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE SOIL PROFILE PROPERTIES - FILE 2 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C C C = ROCK C C = SAND C C = SILT C C = CLAY C C = SCOND C C = CATEXC C C = ALPHA C C = XN C C = VGTHS C C = VGTHR C C C C *********************************************************** C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC'p0*((@@Ԍ 234 CONTINUE WRITE(3,*) READ(12,'(A60)')DATAID WRITE(3,*)' DATAID' WRITE(3,'(A60,/)')DATAID READ(12,*)LYRSOL WRITE(3,*)'LYRSOL' WRITE(3,'(I4,/)')LYRSOL WRITE(3,*)'DIFF0 THETA0 BETA GLAYR THETAS FC THETAR AIRD #R ETA GBD PSIBUB' CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE SOIL PROFILE PROPERTIES - FILE 2 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************* C C C C = DATAID C C = LYRSOL C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC DO 245 I=1,LYRSOL+2 IF(I.LE.LYRSOL) THEN READ(12,*)DIFF0,THETA0,BETA,GLAYR,THETAS,FC # ,THETAR,AIRDR,ETA,GBD,PSIBUB WRITE(3,342)DIFF0,THETA0,BETA,GLAYR,THETAS,FC # ,THETAR,AIRDR,ETA,GBD,PSIBUB ELSE READ(12,*)DIFF0,THETA0,BETA,THETAS,FC # ,THETAR,AIRDR,ETA,BD,PSIBUB WRITE(3,343)DIFF0,THETA0,BETA,THETAS,FC # ,THETAR,AIRDR,ETA,GBD,PSIBUB ENDIF 342 FORMAT(E8.2,2X,F5.3,2X,F4.1,2X,E8.2,7(2X,F4.2)) 343 FORMAT(E8.2,2X,F5.3,2X,F4.1,2X,8X,7(2X,F4.2)) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE SOIL PROFILE PROPERTIES - FILE 2 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C'q0*((@@ԌC ************************************************************** C C C C = DIFF0 C C = THETA0 C C = BETA C C = GLAYR C C = THETAS C C = FC C C = THETAR C C = AIRDR C C = ETA C C = GBD C C = PSIBUB C C ************************************************************ C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 245 CONTINUE READ(12,*)TD,THETAI,BDSLOP, BDRATO, PSISFC WRITE(3,*) WRITE(3,*)'TD THETAI BDSLOP BDRATO PSISFC' CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE SOIL PROFILE PROPERTIES - FILE 2 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C = TD C C = THETAI C C = BDSLOP C C = BDRATO C C = PSISFC C C ************************************************************* C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC WRITE(3,113)TD,THETAI,BDSLOP, BDRATO, PSISFC 113 FORMAT (F4.1,1X,F5.3,3(2X,F5.2),/) READ(12,'(A60)') SNAME WRITE(3,*)' SNAME' WRITE(3,'(A60)') SNAME READ(12,'(I2)')NCURVE WRITE(3,'(/,A10)')'NCURVE' WRITE(3,'(4X,I2,/)')NCURVE C======================================================================= 'r0*((@@ C SETUP LOOP TO READ IN CURVES C======================================================================= DO 129 I=1,NCURVE READ(12,*)INRIM,GH2OC WRITE(3,*)'INRIM(I) GH2OC(I)' WRITE(3,'(2X,I2,9X,F5.2)')INRIM,GH2OC WRITE(3,*)' TSTBD(I,K)' ,' TSTIMP(I,K)' CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE SOIL PROFILE PROPERTIES - FILE 2 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C C C = INRIM(I) C C = GH2OC(I) C C C C ************************************************************* C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC DO 131 K=1,INRIM READ(12,*)TSTBD ,TSTIMP CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE SOIL PROFILE PROPERTIES - FILE 2 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C C C = TSTBD(I,K) C C = TSTIMP(I,K) C C C C ************************************************************* C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC WRITE(3,'(F7.2,4X,F7.2)')TSTBD ,TSTIMP 131 CONTINUE WRITE(3,*) 129 CONTINUE WRITE(3,'(A8,/,A45,/,8X,A25)')'FILE 2= # ',FILE2,'WAS PROCESSED CORRECTLY' 's0*((@@Ԍ C======================================================================= C (ONE LINE FILE) C C======================================================================= PRINT*,'READING FILE 4 ' WRITE(3,'(//,A50)')'********FILE 4********************** # ' WRITE(3,'(A50)')' THE FOLLOWING WILL BE FROM # ',' THE SOIL RESIDUE = FILE 4 # ',FILE4 # ,'******************************************* # ' IER=0 119 IER=IER+1 IF (IER.EQ.2) THEN WRITE(3,*)'ERROR READING FILE 4',FILE4 WRITE(6,*)'ERROR READING FILE 4',FILE4 STOP ENDIF READ(14,141,ERR=119) INSTS,SITES,YEAR,EXPTNO # ,TRTNO,STRAW,SDEP,SCN,ROOT WRITE(3,*)'INSTS/' WRITE(3,*)'SITES/' WRITE(3,*)' YEAR EXPTNO TRTNO STRAW SDEP SCN ROOT' WRITE(3,142) INSTS,SITES,YEAR,EXPTNO # ,TRTNO,STRAW,SDEP,SCN,ROOT 141 FORMAT(2(A2),2(I2),1X,I2,4(1X,F5.0)) 142 FORMAT(2(A2),2(I2),6X,I2,3X,F6.1,1X,F4.1,1X,F5.1,1X,F6.1) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE ORGANIC RESIDUE - FILE 4 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************* C C C C = INSTS C C = SITES C C = YEAR C C = EXPTNO C C = TRTNO C C = STRAW C C = SDEP C C = SCN C C = ROOT C C C C ************************************************************ C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC't0*((@@Ԍ WRITE(3,'(/,A8,/,A45,/,8X,A25)')'FILE 4= # ',FILE4,'WAS PROCESSED CORRECTLY' C======================================================================= C FILE 5 FOLLOWS,'SOIL PROFILE INITIAL CONDITIONS' C======================================================================= WRITE(3,'(//,A50)')'********FILE 5********************** # ' WRITE(3,'(A50)')' THE FOLLOWING WILL BE FROM # ','SOIL PROFILE INITIAL CONDITIONS = FILE 5 # ',FILE5 # , '******************************************* # ' PRINT*,'READING FILE 5' READ(15,151)TRTNO,INSTS,SITES,YEAR,EXPTNO WRITE(3,*)'TRTNO INSTS/SITES/YEAR EXPTNO' WRITE(3,151)TRTNO,INSTS,SITES,YEAR,EXPTNO WRITE(3,*) 151 FORMAT(I2,1X,2A2,2I2) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE SOIL PROFILE INITIAL CONDITIONS - FILE 5 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C = TRTNO C C = INSTS C C = SITES C C = YEAR C C = EXPTNO C C ************************************************************ C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC WRITE(3,'(5A10)')'DLAYR','SW','NH4','NO3','PH' DO 52 I=1,30 READ(15,152,ERR=159,END=160)DLAYR,SW,NH4,NO3,PH IF (INT(DLAYR).NE.-1) THEN WRITE(3,'(5F10.3)')DLAYR,SW,NH4,NO3,PH 152 FORMAT(F6.0,1X,F6.3,3(1X,F4.1)) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE SOIL PROFILE INITIAL CONDITIONS - FILE 5 C 'u0*((@@ԌC C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************* C C = DLAYR C C = SW C C = NH4 C C = NO3 C C = PH C C C C ************************************************************** C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC ELSE WRITE(3,'(F6.0)')DLAYR ENDIF 52 CONTINUE 159 WRITE(3,'(A30)') 'ERROR READING FILE 5 ARRAY' WRITE(6,'(A30)') 'ERROR READING FILE 5 ARRAY' STOP 160 WRITE(3,'(/,A8,/,A45,/,8X,A25)')'FILE 5= # ',FILE5,'WAS PROCESSED CORRECTLY' WRITE(3,'(//,A50)')'********FILE 6********************** # ' WRITE(3,'(A50)')' THE FOLLOWING WILL BE FROM # ',' THE IRRIGATION MANAGEMENT = FILE 6 # ',FILE6 # ,'******************************************* # ' WRITE(3,*)'TRTNO INSTE/SITEE/YEAR EXPTNO' PRINT*,'READING FILE 6' READ(16,161,ERR=169) TRTNO,INSTE,SITEE,YEAR,EXPTNO 161 FORMAT(I2,1X,2(A2),2(I2)) WRITE(3,161) TRTNO,INSTE,SITEE,YEAR,EXPTNO WRITE(3,*) WRITE(3,*)' AMTIRR X1IRR X2IRR' WRITE(3,*)'JDIRR IRRCOD Z1IRR Z2IRR' CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE OF IRRIGATION MANAGEMENT - FILE 6 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************* C C C C = TRTNO C'v0*((@@ԌC = INSTE C C = SITEE C C = YEAR C C = EXPTNO C C C C ************************************************************* C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC DO 168 I=1,200 READ(16,162,END=170,ERR=169)JDIRR,AMTIRR,IRRCOD,X1IRR # ,Z1IRR,X2IRR,Z2IRR IF(JDIRR.GT.0 )THEN WRITE(3,162)JDIRR,AMTIRR,IRRCOD,X1IRR # ,Z1IRR,X2IRR,Z2IRR 162 FORMAT(I4,F5.0,I2,4(F6.1)) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACING OF IRRIGATION MANAGEMENT C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C C C = JDIRR C C = AMTIRR C C = IRRCOD C C = X1IRR C C = Z1IRR C C = X2IRR C C = Z2IRR C C C C ************************************************************* C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC ELSE IF(JDIRR.EQ.-1) THEN WRITE(3,'(1X,I3)')JDIRR ENDIF ENDIF 168 CONTINUE 169 WRITE(3,'(A30)') 'ERROR READING FILE6 ARRAY' WRITE(6,'(A30)') 'ERROR READING FILE6 ARRAY' STOP 170 WRITE(3,'(/,A8,/,A45,/,8X,A25)')'FILE 6= # ',FILE6,'WAS PROCESSED CORRECTLY' C======================================================================= 'w0*((@@ C FILE 7 , "FERTILIZER MANAGEMENT" C======================================================================= WRITE(3,'(//,A50)')'********FILE 7********************** # ' WRITE(3,'(A50)')' THE FOLLOWING WILL BE FROM # ',' FERTILIZATION MANAGEMENT = FILE 7 # ',FILE7 # ,'******************************************* # ' WRITE(3,*)'TRTNO INSTE/SITEE/YEAR EXPTNO' PRINT*,'READING FILE 7 ' READ(17,171,ERR=179) TRTNO,INSTE,SITEE,YEAR,EXPTNO WRITE(3,171) TRTNO,INSTE,SITEE,YEAR,EXPTNO 171 FORMAT(I2,1X,2(A2),2(I2)) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE OF FERTILIZATION MANAGEMENT - FILE 7 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C = TRTNO C C = INSTE C C = SITEE C C = YEAR C C = EXPTNO C C C C ************************************************************** C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC WRITE(3,*) WRITE(3,*)' FERTN INTYPE IPTYPE IKTYPE IITYPE X1F #ERT X2FERT' WRITE(3,*)'JDFERT DFERT FERTP FERTK FERTIN FERCOD # Z1FERT Z2FERT' DO 178 I=1,200 READ(17,172,END=180,ERR=179)JDFERT,FERTN,DFERT,INTYPE,FERTP # ,IPTYPE,FERTK,IKTYPE,FERTIN,IITYPE,FERCOD,X1FERT,Z1FERT # ,X2FERT,Z2FERT 172 FORMAT(1X,I3,F6.1,4(F6.1,I3),I3,4(F6.1)) IF(JDFERT.NE.-1) THEN WRITE(3,172)JDFERT,FERTN,DFERT,INTYPE,FERTP # ,IPTYPE,FERTK,IKTYPE,FERTIN,IITYPE,FERCOD # ,X1FERT,Z1FERT,X2FERT,Z2FERT CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C'x0*((@@ԌC INTERFACE OF FERTILIZATION MANAGEMENT - FILE 7 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C C C = JDFERT C C = FERTN C C = DFERT C C = INTYPE C C = FERTP C C = IPTYPE C C = FERTK C C = IKTYPE C C = FERTIN C C = IITYPE C C = FERCOD C C = X1FERT C C = Z1FERT C C = X2FERT C C = Z2FERT C C C C *************************************************************** C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC ELSE WRITE(3,'(1X,I3)')JDFERT ENDIF 178 CONTINUE 179 WRITE(3,'(A30)') 'ERROR READING FILE 7 ARRAY' WRITE(6,'(A30)') 'ERROR READING FILE 7 ARRAY' STOP 180 WRITE(3,'(/,A8,/,A45,/,8X,A25)')'FILE 7= # ',FILE7,'WAS PROCESSED CORRECTLY' C======================================================================= C FILE 8 TREATMENT MANAGEMENT C======================================================================= WRITE(3,'(//,A50)')'********FILE 8********************** # ' WRITE(3,'(A50)')' THE FOLLOWING WILL BE FROM #' ,' THE TREATMENT MANAGEMENT = FILE 8 #' ,FILE8 #, '******************************************* #' WRITE(3,*)'INSTE/' WRITE(3,*)'SITEE/ TRTNO' WRITE(3,*)'YEAR EXPTNO TITLET ISO'y0*((@@Ԍ #ILT IVARTY' PRINT*,'READING FILE 8 ' READ(18,181)INSTE, SITEE, YEAR, EXPTNO, TRTNO # ,TITLET, ISOILT, IVARTY WRITE(3,181)INSTE, SITEE, YEAR, EXPTNO, TRTNO # ,TITLET, ISOILT, IVARTY WRITE(3,*) 181 FORMAT(2(A2),2(I2),1X,I2,1X,A40,2(1X,I4)) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE OF TREATMENT MANAGEMENT - FILE 8 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C C C = INSTE C C = SITEE C C = YEAR C C = EXPTNO C C = TRTNO C C = TITLET C C = ISOILT C C = IVARTY C C C C ************************************************************* C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC READ(18,182)ISIM,ISOW,PLANTS,ROWSPC,SDEPTH #, IIRR,ISWNIT,EFFIRR,DSOIL,THETAC,IMERGE,ISWEED,ISWINS #, ISWNEM,ISWDIS WRITE(3,*)' ISOW ROWSPC IIRR EFFIRR THETAC # ISWEED ISWNEM ' WRITE(3,*)'ISIM PLANTS SDEPTH ISWNIT DSOIL IM #ERGE ISWINS ISWDIS' WRITE(3,183)ISIM,ISOW,PLANTS,ROWSPC,SDEPTH #, IIRR,ISWNIT,EFFIRR,DSOIL,THETAC,IMERGE,ISWEED,ISWINS #, ISWNEM,ISWDIS CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE OF TREATMENT MANAGEMENT - FILE 8 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C C C = ISIM C'z0*((@@ԌC = ISOW C C = PLANTS C C = ROWSPC C C = SDEPTH C C = IIRR C C = ISWNIT C C = EFFIRR C C = DSOIL C C = THETAC C C = IMERGE C C = ISWEED C C = ISWINS C C = ISWNEM C C = ISWDIS C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 182 FORMAT(1X,I3,1X,I3,1X,F6.2,1X,F6.3,1X,F5.2,2(1X,I2) #, 1X,F6.2,1X,F5.2,1X,F6.2,1X,I3,4(1X,I2)) 183 FORMAT(1X,I3,1X,I3,1X,F6.2,1X,F6.3,1X,F5.2,2(1X,I2) #, F6.2,1X,F5.2,F6.2,1X,I3,4(2X,I2)) READ(18,'(A76)')HISTRY WRITE(3,*) WRITE(3,*)' HISTORY ' WRITE(3,'(A76)')HISTRY CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C INTERFACE OF TREATMENT MANAGEMENT - FILE 8 C C C C ************************************************************* C C *MODEL* = *EXPERIMENT* X CONVERSION FACTOR C C ____(VAR.NAME)________(var.NAME)_________(must be checked)_ C C ************************************************************** C C = HISTRY C C C C ************************************************************ C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC WRITE(3,'(/,A8,/,A45,/,8X,A25)')'FILE 8= # ',FILE8,'WAS PROCESSED CORRECTLY' CALL VALFAB(FILEA,FILEB) RETURN END #{0*((@@ CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C SUBROUTINE VALFAB - THIS ROUTINE IS USED TO READ THE C C VALIDATION DATA, FILES A AND B C C C C - CALLED BY 'READDA' C C C C - CALLS NO OTHER SUBROUTINEDS C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC SUBROUTINE VALFAB(FILEA,FILEB) REAL JNNODM,JNSQRM,JNFLWM,JNGBLM,JNMBLM,JNABSM INTEGER YEAR,EXPTNO,TRTNO CHARACTER INSTE*2,SITEE*2,FILEB*45,FILEA*45,LABELB*7 DIMENSION JDOYA(100), NV(16),RVDATA(16,100),RV(16) #, XPLTHT(100) , XLAI(100) , JNNODM(100) , JNSQRM (100) #, JNFLWM (100), JNGBLM (100) , JNMBLM (100) , JNABSM (100) #, XWLEFH(100), XWSTMH (100) , XWROTH (100) #, XWGBLH (100) , XWMBLH (100) , XWLINT (100) #, XWSEED (100) , XWBURR (100) #, LABELB(16) C======================================================================= C PROCESS VALIDATION DATA - FILE A C======================================================================= OPEN (19,FILE = FILEA) WRITE(3,'(//,A50)')'********FILE A********************** # ' WRITE(3,'(A50)')' THE FOLLOWING WILL BE FROM # ',' MEASURED HARVEST DATA = FILE A # ',FILEA # ,'******************************************* # ' WRITE(3,*)'INSTE/' WRITE(3,*)'SITEE/' WRITE(3,*)'YEAR EXPTNO XLTYLD XSDWT XSPB # XBIOM' WRITE(3,*)' TRTNO XSDYLD XBLSM XLAIMX # XSTMBR ' PRINT*,'READING FILE A ' READ(19,191)INSTE, SITEE, YEAR, EXPTNO, TRTNO # ,XLTYLD ,XSDYLD,XSDWT,XBLSM,XSPB, XLAIMX,XBIOM ,XSTMBR WRITE(3,191)INSTE, SITEE, YEAR, EXPTNO, TRTNO # ,XLTYLD ,XSDYLD,XSDWT,XBLSM,XSPB, XLAIMX,XBIOM ,XSTMBR 191 FORMAT(2(A2),2(I2),1X,I2,2(1X,F7.0),1X,F7.4,1X,F6.0 # ,1X,F4.0,1X,F5.2,2(1X,F6.0)) READ(19,'(F6.2,3(1X,F5.1))')XSDTN,XTOTNP,XAPTNP,XSDN WRITE(3,*) WRITE(3,*)'XSDTN XTOTNP XAPTNP XSDN' WRITE(3,'(F6.2,3F6.1)')XSDTN,XTOTNP,XAPTNP,XSDN'|0*((@@Ԍ WRITE(3,'(/,A8,/,A45,/,8X,A25)')'FILE A= # ',FILEA,'WAS PROCESSED CORRECTLY' C======================================================================= C FILE B FOLLOWS C======================================================================= OPEN(20,FILE=FILEB) WRITE(3,'(///,A50)')'********FILE B********************** # ' WRITE(3,'(A50)')' THE FOLLOWING WILL BE FROM # ',' MEASURED INTERMEDIATE GROWTH = FILE B # ',FILEB # ,'******************************************* # ' PRINT*,'READING FILE B ' C====================================================================== C INITIALIZE ARRAYS FOR UP TO 100 ROWS OF DATA C====================================================================== DO 129 I=1,16 DO 131 J=1,100 NV(I)=0 RVDATA(I,J)=0.0 RV(I)=0.0 JDOYA(J)=0 131 CONTINUE 129 CONTINUE READ(20,323,ERR=209)INSTE,SITEE,YEAR,EXPTNO # ,TRTNO,JEMRGD,JSQRJD,JFLRJD WRITE(3,*)'INSTE/' WRITE(3,*)'SITEE/' WRITE(3,*)'YEAR EXPTNO TRTNO JEMRGD JSQRJD JFLRJD' WRITE(3,324)INSTE,SITEE,YEAR,EXPTNO # ,TRTNO,JEMRGD,JSQRJD,JFLRJD 323 FORMAT(2(A2),2(I2),1X,I2,3(2X,I3)) 324 FORMAT(2(A2),2(I2),6X,I2,5X,I3,4X,I3,3X,I3) READ(20,*,ERR=209)NOVAR,(NV(I),I=1,NOVAR) WRITE(3,'(/,A10)')' NOVAR' WRITE(3,'(17I3)')NOVAR,(NV(I),I=1,NOVAR) C====================================================================== C LOADING LARGE ARRAY C====================================================================== DO 908 ICOUNT=1,100 READ(20,*,ERR=209,END=488)JDOYA(ICOUNT),(RV(I),I=1,NOVAR)'}0*((@@Ԍ IF(JDOYA(ICOUNT).NE.-1)THEN DO 888 J=1,NOVAR RVDATA(ICOUNT,J)=RV(J) 888 CONTINUE ELSE IF(ICOUNT.EQ.1)THEN WRITE(3,*)'-1 ' WRITE(3,*)'NO VALIDATION DATA ' GOTO 434 ENDIF IROW=ICOUNT-1 ENDIF 908 CONTINUE 209 WRITE(3,*)'ERROR IN READING FILE B' WRITE(6,*)'ERROR IN READING FILE B' STOP C=====================================================================C C INITIALIZE INDIVIDUAL ARRAYS FOR EACH STATE VARIABLE C C=====================================================================C 488 IROW=ICOUNT-1 DO 139 I=1,NOVAR DO 141 J=1,IROW LABELB (I)='EMPTY' XPLTHT(J)=0.0 XLAI(J)=0.0 JNNODM(J)=0.0 JNSQRM (J)=0.0 JNFLWM (J)=0.0 JNGBLM (J)=0.0 JNMBLM (J)=0.0 JNABSM (J)=0.0 XWLEFH(J)=0.0 XWSTMH (J)=0.0 XWROTH(J)=0.0 XWGBLH (J)=0.0 XWMBLH (J)=0.0 XWLINT (J) =0.0 XWSEED (J)=0.0 XWBURR (J)=0.0 141 CONTINUE 139 CONTINUE C====================================================================== C LOAD STATE VARIABLE INTO CORRECT 1-DIMENSIONAL ARRAY NAME C ALSO ASSIGNING CORRECT HEADER TO ARRAY 'LABELB' C====================================================================== DO 263 I = 1, NOVAR DO 265 J = 1, IROW '~0*((@@Ԍ IF(NV(I).EQ.1) THEN LABELB(I)='XPLTHT' XPLTHT(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.2) THEN LABELB(I)='XLAI' XLAI(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.3) THEN LABELB(I)='JNNODM' JNNODM(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.4) THEN LABELB(I)='JNSQRM' JNSQRM(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.5) THEN LABELB(I)='JNFLWM' JNFLWM(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.6) THEN LABELB(I)='JNGBLM' JNGBLM(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.7) THEN LABELB(I)='JNMBLM' JNMBLM(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.8) THEN LABELB(I)='JNABSM' JNABSM(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.9) THEN LABELB(I)='XWLEFH' XWLEFH(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.10)THEN LABELB(I)='XWSTMH' XWSTMH(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.11) THEN LABELB(I)='XWROTH' XWROTH(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.12) THEN LABELB(I)='XWGBLH' XWGBLH(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.13) THEN LABELB(I)='XWMBLH' XWMBLH(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.14) THEN LABELB(I)='XWLINT' XWLINT(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.15) THEN LABELB(I)='XWSEED' XWSEED(J) = RVDATA(J,I) ELSEIF (NV(I).EQ.16) THEN LABELB(I)='XWBURR' XWBURR(J) = RVDATA(J,I) ELSE WRITE(3,*)'ERROR IN INDEXING OF STATE VARIABLES' WRITE(6,*)'ERROR IN INDEXING OF STATE VARIABLES' STOP'0*((@@Ԍ ENDIF 265 CONTINUE 263 CONTINUE C======================================================================= C THIS WRITES-OUT HEADERS FOR INDEXED STATE VARIABLE. C HEADERS NOT USED WILL BE 'EMPTY' IN PRINT-OUT. C======================================================================= IF (NOVAR.LE.8) THEN WRITE(3,*) WRITE(3,'(A4,8(A8))')'JDOY',(LABELB(I),I=1,NOVAR) DO 99 J=1,IROW WRITE(3,'(I3,8(F8.2))')JDOYA(J),(RVDATA(J,L),L=1,NOVAR) 99 CONTINUE WRITE(3,'(I3)')JDOYA(ICOUNT) ELSE WRITE(3,*) WRITE(3,'(A4,8(A8))')'JDOY',(LABELB(I),I=1,8) DO 599 J=1,IROW WRITE(3,'(I3,8(F8.2))')JDOYA(J),(RVDATA(J,L),L=1,8) 599 CONTINUE WRITE(3,'(I3)')JDOYA(ICOUNT) WRITE(3,*) WRITE(3,'(A4,8(A8))')'JDOY',(LABELB(I),I=9,NOVAR) DO 699 J=1,IROW WRITE(3,'(I3,8(F8.2))')JDOYA(J),(RVDATA(J,L),L=9,NOVAR) 699 CONTINUE ENDIF 434 WRITE(3,'(/,A8,/,A45,/,8X,A25)')'FILE B= # ',FILEB,'WAS PROCESSED CORRECTLY' RETURN END 0*((@@