#****************************************************************
#      file:  erod.in
#****************************************************************
#
#     +++ PURPOSE +++
#
#     Input file which is read by subroutine erodin.for
#
#     All lines beginning with a "#" character are assumed to
#     be comment lines and are skipped.
#
#     +++ DEFINITIONS +++
#
#     * = inputs NOT presently used by erosion
#         (most of these are expected to be used in the future)
#     All other input values must be correctly specified.
#
#     All comments prior to each line of data input
#     in this template input file have the following format:
#
#     Variable_Name, Var_type, (include file)  Text Definition
#
#     where Var_type is: I = integer L = logical R = real
#
#
# +++ DEBUG STUFF +++
#
#     debugflg - debug flag for providing different levels of debug inof
#                value of 1 will print out and numboer all input lines
#                value of 2 will print out and number all data input lines
#                value of 3 will do both 1 and 2 input line debug output
 0
#
#
# +++ INIT STUFF +++
#
#     am0eif, L, (m1flag.inc) EROSION initialization flag
 .TRUE.
#     am0efl, L, (m1flag.inc) EROSION "print" flag
 1
#
# +++ SIMULATION REGION +++
#
#     amxsim(x,y), R, (m1geo.inc) Simulation Region diag. coordinates (m)
#                     Input x,y coordinates in this form: x1,y1  x2,y2
 0.0  0.0  180.0  180.0
#
#     amasim, R, (m1geo.inc) Simulation Region orientation angle (deg)
 0.0
#
# +++ ACCOUNTING REGIONS +++
#
#     nacctr, I, (m1geo.inc) Number of accounting regions (must be 1 for now)
 1
#
#     amxar(x,y,a), R, (m1geo.inc) Accounting Region diag. coordinates (m)
#                     Input x,y coordinates in this form: x1,y1  x2,y2
#                     for each accounting region specified (nacctr)
 0.0, 0.0  180.0, 180.0
#
# +++ BARRIERS +++
#
#     nbr, I, (m1geo.inc) Number of barriers (0-5)
 0
#     NOTE: Remaining BARRIER inputs are repeated for each barrier specified
#     If no barriers specified (nbr=0), then no BARRIER inputs will be here
#
#     amxbr(x,y,b), R, (m1geo.inc) Accounting Region diag. coordinates (m)
#                     Input x,y coordinates in this form: x1,y1  x2,y2
#                     for each barrier specified (nbr)
#    0.0, 0.0  0.0  100.0
#
#       amzbr(b), R, (m1geo.inc) Barrier height (m)
#       ampbr(b), R, (m1geo.inc) Barrier porosity (m^2/m^2)
#       amxbrw(b), R, (m1geo.inc) Barrier width (m)
#  1.2   0.50   2.0
#
# +++ SUBREGION REGIONS +++
#
#     nsubr, I, (m1subr.inc) Number of subregions (1-5)
 1
#
#     NOTE: Remaining SUBREGION inputs (BIOMASS, SOIL, and HYDROLOGY,
#     ie. variables defined by subregion) are repeated for nsubr"
#     subregions specified
#
#     amxsr(x,y,s), R, (m1subr.inc) Subregion diag. coordinates (m)
#                     Input x,y coordinates in this form: x1,y1  x2,y2
#                     for each subregion specified (subr)
 0.0, 0.0  180.0, 180.0
#
#     +++ BIOMASS +++
#
#       abzht(s), R, (b1glob.inc) Overall biomass height (m)
 0.0
#
#       aczht(s), R, (c1glob.inc) Crop height (m)
 0.0
#
#       acrsai(s), R, (c1glob.inc) Crop stem area index (m^2/m^2)
#       acrlai(s), R, (c1glob.inc) Crop leaf area index (m^2/m^2)
 0.0    0.0
#
#       adrsaitot(s), R, (d1glob.inc) Residue stem area index (m^2/m^2)
#       adrlaitot(s), R, (d1glob.inc) Residue leaf area index (m^2/m^2)
 0.0   0.0
#
#       acxrow(s) Crop row spacing (m) and pla
#       ac0rg(s)  Crop seed placement (0 - furrow, 1 - ridge)
 0.33 0
#
# These are not implemented within EROSION
#       abrsaz(h,s), R, (b1geom.inc) Biomass stem area index by ht (1/m)
# (should be 5 values here when used)
#       abrlaz(h,s), R, (b1geom.inc) Biomass leaf area index by ht (1/m)
# (should be 5 values here when used)
#
# Only abffcv(s) is currently implemented within EROSION
#       abffcv(s), R, (b1geom.inc) Flat biomass cover (m^2/m^2)
#       abfscv(s), R, (b1geom.inc) Standing biomass cover (m^2/m^2)
#       abftcv(s), R, (b1geom.inc) Total biomass cover (m^2/m^2)
# (should be 3 values here when abffcv(s) and abfscv(s) get used in the future)
 0.03
#
#     +++ SOIL +++
#
#     nslay(s), I, (s1layr.inc) Number of soil layers (3-10)
 3
#
#     NOTE: Remaining SOIL inputs are repeated for each layer specified
#
#     aszlyt(l,s), R, (s1layr.inc) Soil layer thickness (mm)
 10.0 40.0 200.0
#
#     asdblk(l,s), R, (s1phys.inc) Soil layer bulk density (Mg/m^3)
 1.4 1.4 1.4
#
#     asfsan(l,s), R, (s1dbh.inc) Soil layer sand content (Mg/Mg)
 0.83 0.83 0.83
#     asfvfs(l,s), R  (s1dbh.inc) Soil layer very fine sand(Mg/Mg)
 0.20   0.20    0.20
#     asfsil(l,s), R, (s1dbh.inc) Soil layer silt content (Mg/Mg)
 0.09 0.09 0.09
#     asfcla(l,s), R, (s1dbh.inc) Soil layer clay content (Mg/Mg)
 0.08 0.08 0.08
#
#     asvroc(l,s), R, (s1dbh.inc) Soil layer rock volume (m^3/m^3)
 0.0 0.0 0.0
#
#     asdagd(l,s), R, (s1agg.inc) Soil layer agg density (Mg/m^3)
 1.8 1.8 1.8
#     aseags(l,s), R, (s1agg.inc) Soil layer agg stability ln(J/kg)
 1.93   2.65   2.65
#     aslagm(l,s), R, (s1agg.inc) Soil layer GMD (mm)
 0.38   0.98  0.98
#     aslagn(l,s), R, (s1agg.inc) Soil layer minimum agg size (mm)
 0.001 0.01 0.01
#     aslagx(l,s), R, (s1agg.inc) Soil layer maximum agg size (mm)
 18.16  14.8 48.8
#     as0ags(l,s), R, (s1agg.inc) Soil layer GSD (mm/mm)
 26.03    20.4   20.4
#
#
#     asfcr(s), R, (s1surf.inc) Surface crust fraction (m^2/m^2)
#     aszcr(s), R, (s1surf.inc) Surface crust thickness (mm)
#     asflos(s), R, (s1surf.inc) Fraction of loose material on surface (m^2/m^2)
#     asmlos(s), R, (s1surf.inc) Mass of loose material on crust (kg/m^2)
#     asdcr(s), R, (s1surf.inc) Soil crust density (Mg/m^3)
#     asecr(s), R, (s1surf.inc) Soil crust stability ln(J/kg)
 0.4    5.30    0.28    0.255   1.2   1.93
#
#     aslrr(s), R, (s1sgeo.inc) Allmaras random roughness (mm)
  2.0
#     aszrgh(s), R, (s1sgeo.inc) Ridge height (mm)
#     asxrgs(s), R, (s1sgeo.inc) Ridge spacing (mm)
#     asxrgw(s), R, (s1sgeo.inc) Ridge width (mm)
#     asxrgo(s), R, (s1sgeo.inc) Ridge orientation (deg)
 00.0  00.0    0.0   0.0
#     asxdks(s), R, (s1sgeo.inc) Dike spacing (mm)
 0.0
#
#     +++ HYDROLOGY +++
#
#     ahzsnd(s), R, (s1sgeo.inc) Snow depth (mm)
 0.0
#
#     ahrwcw(l,s), R, (h1db1.inc) Soil layer wilting point water content (Mg/Mg)
 0.05 0.05 0.05
#
#     ahrwca(l,s), R, (h1db1.inc) Soil layer water content (Mg/Mg)
 0.01  0.05 0.05
#
#     ahrwc0(h,s), R, (h1db1.inc) Surface layer water content (Mg/Mg)
#                  NOTE: the near surface water content is specified on an
#                        hourly basis.  We read in the hrly water content
#                        on two lines, with 12 values in each line.
 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
#
#
# NOTE: This is the end of the SUBREGION variables
#
#     +++ WEATHER +++
#
#     awdair, R, (w1pavg.inc) Air density (kg/m^3)
 1.2
#
#     awadir, R, (w1wind.inc) Wind direction (deg)
 270.0
#
#     ntstep, I, (local variable) Number of intervals/day to run EROSION
 48
#
#     anemht, R  anemometer height (m)
#     awzzo,  R  aerodynamic roughness at anemometer site (mm)
#     wzoflg, I (global variable) zo location flag
#               (flag =0 - zo fixed at wx sta. location)
#               (flag = 1 - zo variable at field location)
  2.0   25.0    1
#     wflg, I, (local variable) Wind/Weibull flag
#              (0 - read in Weibull parameters, 1 - read in wind speeds)
 0
#
# NOTE: This is only present when (wflg=0)
#     wfcalm, R, (local variable) Fraction of time winds are calm (hr/hr)
#     wuc, R, (local variable) Weibull "c" factor (m/s)
#     w0k, R, (local variable) Weibull "k" factor (fraction)
     0.142  6.500  4.767
#
# NOTE: The remaining data is only present when (wflg=1)
#       condition wflg=1 uses standard input from windgen in WEPS.
#    awu(i), R, (w1wind) Wind speed for (ntstep) intervals (m/s)
# I think I can read multiple lines with variable number of values
# We will try and see - LEW  Must use 6 values per line LH.
#
#    12.0      12.0       12.0        12.0    12.0       12.0
#    12.0      12.0       12.0        12.0    12.0       12.0
#    0.0       12.0       12.0        12.0    12.0       12.0
#  14.0       15.0       12.00       12.0    12.0       10.0
#
#    + + + DATA TO PLOT + + +
#
#    names and values to input for plot
# place 1 flag in 1st line after #-name line for variables to include in plot
#
#  initial xplot value,I, (-1=no plot, 0 = output indep.variables with 1 flag)
  -1
####
#    xin(i), R, (field length)
   1
  Length(m)
# abzht, R, (biomass ht.(m))
   0
 bio_ht(m)
# abrsai, R (stem area index)
  1
 stem_area
#   abrlai(s), R, (b1geom.inc) Biomass leaf area index (m^2/m^2)
  1
 lai_area
# abffcv, R,( biomass flat fraction cover)
  0
 flat_cov
#     asfvfs(1,s), R, (soil fraction very fine sand in layer 1)
  0
  vfsand
#     asfsan(1,s), R, (soil fraction sand in layer 1)
  1
   sand
#     asfsil(1,s), R (soil fraction silt in layer 1)
  0
  silt
#     asfcla(1,s), R (soil fraction clay in layer 1)
  0
  clay
#     asvoc(1,s), R (soil volume roc in layer 1)(m^3/m^3)
  0
  rock_vol
#     aseags(1,s), R (soil aggregate stability) (ln J/m^3)
  0
  ag_stab
#     aslagm(1,s), R (soil aggregate geom. mean dia.) (mm)
  0
  ag_gmd
#     aslagn(1,s), R (soil aggregate min. dia.) (mm)
  0
  ag_min
#     aslagx(1,s), R (soil aggregate max. dia.) (mm)
  0
  ag_max
#     as0ags(1,s), R (soil aggregate geo. std. dev.)
  0
  ag_std
#     asfcr(s), R, (s1surf.inc) Surface crust fraction (m^2/m^2)
  0
  crust_cv
#     aszcr(s), R, (s1surf.inc) Surface crust thickness (mm)
  0
 crust_z(mm)
#     asflos(s), R, (s1surf.inc) Fraction of loose material on surface (m^2/m^2)
  0
  los_cv
#     asmlos(s), R, (s1surf.inc) Mass of loose material on crust (kg/m^2)
  0
  los(kg/m^2)
#     asdcr(s), R, (s1surf.inc) Soil crust density (Mg/m^3)
  0
  cr_den(Mg/m^3)
#     asecr(s), R, (s1surf.inc) Soil crust stability ln(J/kg)
  0
  cr_se
#     aslrr(s), R, (s1sgeo.inc) Allmaras random roughness (mm)
  0
  rr(mm)
#     aszrgh(s), R, (s1sgeo.inc) Ridge height (mm)
  0
  z_rgh(mm)
#     asxrgs(s), R, (s1sgeo.inc) Ridge spacing (mm)
  0
  x_rgs(mm)
#     asxrgw(s), R, (s1sgeo.inc) Ridge width (mm)
  0
  x_rgw(mm)
#     asxrgo(s), R, (s1sgeo.inc) Ridge orientation (deg)
  0
  a_rgo(deg)
#


