# day mon yr 15 3 20 #***************************************************************** # file: erod.in #***************************************************************** # # +++ PURPOSE +++ # # Input file which is read by stand alone erosion model # # +++ DEFINITIONS +++ # # Lines beginning with a "#" character are comments # # * = inputs NOT presently used by erosion # # All other input values must be correctly specified. # # Comments prior to each line of data have the following format # # Variable_Name, Var_type, (inc file) Text Definition # # where Var_type is: I = integer, L = logical, R = real # # +++ DEBUG FLAG +++ # # debugflg - debug flag for providing different levels of debug output # value of 0 will print no debug output # value of 1 will print out and number 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 # # +++ INITIALIZATIONS +++ # # 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 coordinates (m) # Input x,y coordinates in this form: x1,y1 x2,y2 0.00000000E+00 0.00000000E+00 705.609985 705.609985 # # amxsim(x,y), R, (m1geo.inc) Simulation Region orientation angle # clockwise with 0=north 0.00000000E+00 # # +++ 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 coordinates (m) # Input x,y coordinates in this form: x1,y1 x2,y2 # for each accounting region specified (nacctr) 0.00000000E+00 0.00000000E+00 705.609985 705.609985 # # +++ 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 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 # # # +++ SUBREGIONS +++ # # 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 coordinates (m) # Input x,y coordinates in this form: x1,y1 x2,y2 # for each subregion specified (subr) 0.00000000E+00 0.00000000E+00 705.609985 705.609985 # # +++ BIOMASS +++ # # adzht_ave(s), R, (d1glob.inc) Average residue height (m) 0.764339924 # # aczht(s), R, (c1glob.inc) Crop height (m) 0.00000000E+00 # # 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.00000000E+00 0.00000000E+00 # # 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) 1.78759263E-04 3.43357387E-04 # # acxrow(s) Crop row spacing (m) # ac0rg(s) Crop seed placement (0 - furrow, 1 - ridge) 0.203200415 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) are used) 7.79799893E-02 # # +++ SOIL +++ # # nslay(s), I, (s1layr.inc) Number of soil layers (3-10) 27 # # NOTE: Remaining SOIL inputs are repeated for each layer specified # # aszlyt(l,s), R, (s1layr.inc) Soil layer thickness (mm) 2.25385237 2.81731796 3.52164578 4.40162992 5.50202751 6.87754059 8.59592438 10.7436314 13.4249296 16.7773876 20.9715958 26.2144432 30.4377422 44.5573730 55.4666786 74.1902618 64.0634689 80.0799561 88.1016693 110.127556 137.659378 172.073746 215.093079 268.865234 336.081146 353.295288 353.295288 # # asdblk(l,s), R, (s1phys.inc) Soil layer bulk density (Mg/m^3 1.09146523 1.09146523 1.09146523 1.09156930 1.09156930 1.09156930 1.09169888 1.09182715 1.09220445 1.09244990 1.09244990 1.09244990 1.17609429 1.00425744 1.00842249 0.942409754 1.36421740 1.36421740 1.54999399 1.54999399 1.54999399 1.54999399 1.54999399 1.54999399 1.54999399 1.54999399 1.54999399 # # asfsan(l,s),R,(s1dbh.inc) Soil layer sand content (Mg/Mg) 0.630000114 0.630000174 0.630000114 0.630000114 0.630000114 0.630000114 0.630000114 0.630000114 0.630000114 0.630000114 0.630000114 0.630000114 0.630000114 0.630000114 0.630000174 0.630000114 0.630000114 0.630000114 0.629999995 0.629999995 0.629999995 0.629999995 0.629999995 0.629999995 0.629999995 0.629999995 0.629999995 # # asfvfs(l,s), R, (s1dbh.inc) Soil layer very fine sand (Mg/Mg) 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999976 0.149999961 0.149999976 0.149999991 0.149999991 0.150000006 0.150000006 0.150000006 0.150000006 0.150000006 0.150000006 0.150000006 0.150000006 0.150000006 # # asfsil(l,s),R,(s1dbh.inc) Soil layer silt content (Mg/Mg) 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000070 0.270000011 0.270000011 0.270000011 0.270000011 0.270000011 0.270000011 0.270000011 0.270000011 0.270000011 # # asfcla(l,s),R,(s1dbh.inc) Soil layer clay content (Mg/Mg) 9.99999791E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999791E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 9.99999866E-02 # # asvroc(l,s), R, (s1dbh.inc) Soil layer rock volume (m^3/m^3) 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 # # asdagd(l,s),R,(s1agg.inc) Soil layer agg density (Mg/m^3) 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 1.79999995 # # aseags(l,s), R, (s1agg.inc) Soil layer agg stability ln(J/kg) 2.06865215 2.07169437 2.07450843 2.07828903 2.08216619 2.08827281 2.09681535 2.12306881 2.15287089 2.12216759 2.03156281 2.01851749 2.00698280 2.05899954 2.05005240 2.05664492 2.29534745 2.28258538 2.50804782 2.75096917 2.55720448 2.56031418 2.56119537 2.56141210 2.51703143 2.50755548 2.48583508 # # aslagm(l,s), R, (s1agg.inc) Soil layer GMD (mm) 0.980590522 0.986580193 0.992808521 1.00206757 1.00991833 1.02313817 1.04281640 1.10328031 1.17268944 1.16871572 1.15512419 1.15682638 1.78477991 1.45615876 1.39278185 1.06871045 2.55711150 2.53949594 2.81950855 2.97761631 2.85584664 2.85400295 2.85686302 2.86027527 2.82935333 2.82227087 2.80364299 # # aslagn(l,s), R, (s1agg.inc) Soil layer minimum agg size (mm) 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 9.99999978E-03 # # aslagx(l,s), R, (s1agg.inc) Soil layer maximum agg size (mm) 34.2809372 34.2749176 34.2688904 34.2603111 34.2534180 34.2425423 34.2280159 34.1948242 34.1753273 34.1759529 34.1785583 34.1781921 34.5014038 34.2433243 34.2119255 34.5131836 35.4157333 35.3926315 35.7656097 35.9800835 35.8147202 35.8122177 35.8160858 35.8207054 35.7789001 35.7693329 35.7442093 # # as0ags(l,s), R, (s1agg.inc) Soil layer GSD (mm/mm) 10.1532612 10.1761236 10.1997852 10.2347527 10.2642069 10.3134098 10.3857460 10.6015291 10.8379250 10.8247004 10.7791891 10.7849131 12.5198727 11.6962948 11.5178041 9.19020271 13.9540358 13.9274111 14.3233700 14.5231056 14.3707151 14.3683338 14.3720264 14.3764248 14.3362837 14.3269997 14.3024206 # # 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.190630719 7.59999943 0.164768413 0.222776875 1.42840958 2.06865215 # # aslrr(s), R, (s1sgeo.inc) Allmaras random roughness (mm) 15.2698593 # # aszrgh(s), R, (s1sgeo.inc) Ridge height (mm) # asxrgs(s), R, (s1sgeo.inc) Ridge spacing (mm) # asxrgw(s), R, (s1sgeo.inc) Ridge width (mm) # asargo(s), R, (s1sgeo.inc) Ridge orientation (deg) 46.6173782 228.600464 0.00000000E+00 0.00000000E+00 # # asxdks(s), R, (s1sgeo.inc) Dike spacing (mm) 0.00000000E+00 # # +++ HYDROLOGY +++ # # ahzsnd(s), R, (s1sgeo.inc) Snow depth (mm) 0.00000000E+00 # # ahrwcw(l,s), R, (h1db1.inc) Soil layer wilting point water content (Mg/Mg) 9.21605974E-02 9.21605974E-02 9.21605974E-02 9.21405926E-02 9.21405926E-02 9.21405926E-02 9.21157449E-02 9.20911655E-02 9.20189247E-02 9.19720158E-02 9.19720158E-02 9.19720158E-02 7.43870884E-02 8.95260051E-02 8.90982151E-02 9.57716927E-02 6.22166544E-02 6.22166544E-02 5.32381572E-02 5.32381572E-02 5.32381572E-02 5.32381572E-02 5.32381572E-02 5.32381572E-02 5.32381572E-02 5.32381572E-02 5.32381572E-02 # # ahrwca(l,s), R, (h1db1.inc) Soil layer water content (Mg/Mg) 0.111949265 0.111949250 0.111949265 0.111913644 0.111913644 0.111913644 0.111869365 0.111825563 0.111696906 0.111613400 0.111613400 0.111613400 8.14211816E-02 0.107279144 0.106525391 0.118426971 6.19299561E-02 6.19299561E-02 4.82975505E-02 4.82975505E-02 4.82975505E-02 4.82975505E-02 4.82975505E-02 4.82975505E-02 4.82975505E-02 4.82975505E-02 4.82975505E-02 # # 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.0472277 0.0479427 0.0486578 0.0493729 0.0500879 0.0508030 0.0515181 0.0522331 0.0529482 0.0536633 0.0543783 0.0550934 0.0558085 0.0565235 0.0572386 0.0579537 0.0586687 0.0593838 0.0600989 0.0608139 0.0615290 0.0622441 0.0629591 0.0636742 # # NOTE: This is the end of the SUBREGION variables # # +++ WEATHER +++ # # awdair, R, (w1pavg.inc) Air density (kg/m^3) 1.19958723 # # awadir, R, (w1wind.inc) Wind direction (deg) 225.000000 # # ntstep, I, (local variable) Number of intervals/day to run EROSION 24 # # 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) 10.0000000 25.0000000 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 the above (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.263 5.856 1.720# # NOTE: The remaining data is only present when (wflg=1) # 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. # 1.60000002 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.00000000 2.09999990 2.79999995 3.29999995 4.00000000 5.59999990 6.59999990 12.6999998 14.5000000 19.3999996 27.7999992 14.8000002 14.3000002 7.90000010 5.90000010 4.50000000 3.59999990 3.00000000 2.40000010 # # + + + 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 = plot indep.variables with 1 flag) -1