#**************************************************************** # april1l.txt: INPUT DATA FOR STORM, BURLINGTON, April 11, 2001 #**************************************************************** # # +++ 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, 600.0, 415.0 # 0.0, 0.0, 100.0, 100.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, 600.0, 415.0 # 0.0, 0.0, 100.0, 100.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) # # 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) # # +++ 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, 600.0, 415.0 # 0.0, 0.0, 100.0, 100.0 # # +++ BIOMASS +++ # # abzht(s), R, (b1glob.inc) Biomass height (m) # 0.149 with wheat crop # 0.010 0.0 no wheat crop # # abrsai(s), R, (b1geom.inc) Biomass stem area index (m^2/m^2) # abrlai(s), R, (b1geom.inc) Biomass leaf area index (m^2/m^2) # 0.11 1.03 with wheat crop # 0.01 0.1 0.0 0.0 no wheat crop # # 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^3/m^3) # abfscv(s), R, (b1geom.inc) Standing biomass cover (m^3/m^3) # abftcv(s), R, (b1geom.inc) Total biomass cover (m^3/m^3) # (should be 3 values here when abffcv(s) and abfscv(s) get used in the future) # 0.098 with residue 0.0 no residue # # +++ SOIL +++ # # nslay(s), I, (s1layr.inc) Number of soil layers (3-10) 1 # # NOTE: Remaining SOIL inputs are repeated for each layer specified # # aszlyt(l,s), R, (s1layr.inc) Soil layer thickness (mm) 1000.0 # # asdblk(l,s), R, (s1phys.inc) Soil layer bulk density (Mg/m^3) 1.8 # # asfsan(l,s), R, (s1dbh.inc) Soil layer sand content (Mg/Mg) 0.31 # asfvfs(l,s), R (s1dbh.inc) Soil layer very fine sand(Mg/Mg) 0.15 # asfsil(l,s), R, (s1dbh.inc) Soil layer silt content (Mg/Mg) 0.53 # asfcla(l,s), R, (s1dbh.inc) Soil layer clay content (Mg/Mg) 0.16 # # asvroc(l,s), R, (s1dbh.inc) Soil layer rock volume (m^3/m^3) 0.0 # # asdagd(l,s), R, (s1agg.inc) Soil layer agg density (Mg/m^3) 1.8 # aseags(l,s), R, (s1agg.inc) Soil layer agg stability ln(J/kg) 2.51 # aslagm(l,s), R, (s1agg.inc) Soil layer GMD (mm) 0.75 # aslagn(l,s), R, (s1agg.inc) Soil layer minimum agg size (mm) 0.02 # aslagx(l,s), R, (s1agg.inc) Soil layer maximum agg size (mm) 44.4 # as0ags(l,s), R, (s1agg.inc) Soil layer GSD (mm/mm) 14.7 # # # 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) 1.0 21.0 0.0 0.0 1.5 2.51 # # # # aslrr(s), R, (s1sgeo.inc) Allmaras random roughness (mm) 5.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) 34 305 0.0 90 # 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.077 # # ahrwca(l,s), R, (h1db1.inc) Soil layer water content (Mg/Mg) 0.0 # # 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 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 # 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 # 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 # 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 # 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 # 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 # 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 0.07 # 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 # 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 0.08 # 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 # 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.10 # 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 top 8 mm measurement # 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 0.12 # # # 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) 324 # # ntstep, I, (local variable) Number of intervals/day to run EROSION 96 # # 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, 3.36 0 # wflg, I, (local variable) Wind/Weibull flag # (0 - read in Weibull parameters, 1 - read in wind speeds) 1 # # 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) # # NOTE: The remaining data is only present when (wflg=1) # # 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 ( Note: no blank lines allowed - LH) # 7.6 6.0 6.1 6.2 6.8 6.6 6.4 6.7 7.2 7.6 8.2 8.5 8.4 8.3 8.9 8.0 8.5 7.9 7.0 6.6 5.7 5.4 5.8 6.0 5.5 5.0 3.5 3.9 2.4 1.9 3.0 3.9 3.7 3.6 3.5 3.4 4.2 4.5 5.7 9.6 13.8 14.4 15.3 16.2 16.8 17.0 17.2 17.0 17.7 18.1 18.3 17.1 17.2 17.8 17.2 17.6 17.8 16.4 17.6 20.9 20.9 19.2 17.0 15.7 16.3 14.7 14.5 14.4 15.6 15.2 14.8 14.3 12.9 13.4 12.3 11.0 11.1 10.0 9.2 8.3 7.5 7.4 6.2 5.6 5.1 5.8 4.9 3.8 4.0 4.0 4.1 4.5 4.8 4.6 4.6 5.2