#**************************************************************** # erod.in.template Template INPUT DATA FILE # Updated January 2006 - LEW # Updated February, 2006 - JT #**************************************************************** # # +++ PURPOSE +++ # # File for input to standalone erosion submodel program (sweep) # # All lines beginning with a "#" character are assumed to # be comment lines and are skipped. # # +++ DEFINITIONS +++ # # All comments prior to each line of data input # in this template input file have the following format: # # Variable_Name, Var_type, Text Definition # # where Var_type is: I = integer L = logical R = real # # +++ VERSION HEADER +++ # included to allow running old erosion file formats with the multiple subregion enabled version # this must be the first non-comment line in the file # Version: 100 multiple subregions # # # +++ DEBUG FLAG +++ # # debugflg - debug flag for providing different levels of debug info # currently useful to debug/check input file data format # # value of 0 will print no debug information # value of 1 will print out and number all input file lines # value of 2 will print out and number all data input lines # value of 3 will do both 1 and 2 0 # # # +++ INITIALIZATION +++ # # am0eif, L, EROSION "initialization" flag # Must be set to .TRUE. for standalone erosion runs .TRUE. # # # +++ SIMULATION REGION +++ # # amxsim(x,y), R, Simulation Region diagonal coordinates (meters) # Input (x,y) coordinates in this form: x1,y1 x2,y2 # If orientation in 0 degrees, then x axis is East-West # and y axis is North-South. # Typical Range: 10.0 to 1600.0 # # NOTE: Accounting region and Subregion coordinates # must also be set to the same values # 0.0, 0.0 1000.0, 200.0 # # amasim, R, Simulation Region orientation angle # clockwise is positive with 0=north 0.0 # # +++ NUMBER OF GRID POINTS +++ # # xgdpt, ygdpt, I, Number of grid points in simulation region x and y directions # A value of 0 for both values uses the default internal grid. # A value of 0 for only one generates and error. 0 0 # # +++ ACCOUNTING REGIONS +++ # # nacctr, I, Number of accounting regions 1 # # accounting region polygon, count, xy pairs (subregions_mod) # polygons can be: # - open (the last point connects to the first point) # - closed (last point is same as first point) # - complex (multiple closed polygons entered as one) 4 0.0 0.0 1000.0 0.0 1000.0 200.0 0.0 200.0 # # +++ BARRIERS +++ # # nbr, I, Number of barriers (0-5) 2 # # NOTE: Remaining BARRIER inputs are repeated for each barrier specified # If no barriers specified (nbr=0), then no BARRIER inputs will be here # # barrier(b)%np, I, number of points in barrier polyline # Inputs are repeated for each point specified # barrier(b)%points(n)%x, barrier(b)%points(n)%y, R, x,y coordinate pair # Example: 0.0 500.0 # barrier(b)%points(n)%amzbr, R, Barrier height (m) # barrier(b)%points(n)%amxbrw, R, Barrier width (m) # barrier(b)%points(n)%ampbr, R, Barrier porosity (m^2/m^2) # Example: 1.2 2.0 0.50 # # # barrier(b)%np, I, number of points in barrier polyline 2 # barrier(b)%points(n)%x, barrier(b)%points(n)%y, R, x,y coordinate pair 0.0 0.0 # amzbr(b), amxbrw(b), ampbr(b) 0.2 15.0 0.5 # barrier(b)%points(n)%x, barrier(b)%points(n)%y, R, x,y coordinate pair 0.0 200.0 # amzbr(b), amxbrw(b), ampbr(b) 0.2 15.0 0.5 # # Barrier #2 number of points 3 # Barrier #2 coordinates (x1,y1) 0.0 0.0 # Barrier #2 height, width and porosity 1 0.2 15.0 0.5 # # Barrier #2 coordinates (x2,y2) 500.0 10.0 # Barrier #2 height, width and porosity 2 0.3 10.0 0.4 # # Barrier #2 coordinates (x3,y3) 1000.0 0.0 # Barrier #2 height, width and porosity 3 0.2 15.0 0.5 # # +++ SUBREGIONS +++ # # nsubr, I, Number of subregions (at least 1) # NOTE: together all must cover simulation region 1 # # NOTE: Remaining SUBREGION inputs (BIOMASS, SOIL, and HYDROLOGY, # ie. variables defined by subregion) are repeated for "nsubr" # subregions specified # # subregion polygon, count, xy pairs (subregions_mod) # polygons can be: # - open (the last point connects to the first point) # - closed (last point is same as first point) # - complex (multiple closed polygons entered as one) 4 0.0 0.0 1000.0 0.0 1000.0 200.0 0.0 200.0 # # +++ BIOMASS +++ # # # adzht_ave(s), R, Height of standing residue (meteres) # WEPS generated input files will provide # "SAI weighted" average residue height # across all residue pools. # Typical Range: 0.0 to 3.0 0.21 # # aczht(s), R, Average height of growing crop (meters) 0.0 # # acrsai(s), R, Growing crop stem area index (m^2/m^2) # Typical Range: 0.0 to 3.0 # acrlai(s), R, Growing crop leaf area index (m^2/m^2) # Typical Range: 0.0 to 8.0 0.0 0.0 # # adrsaitot(s), R, Residue stem area index (m^2/m^2) # adrlaitot(s), R, Residue leaf area index (m^2/m^2) # WEPS generated input files will provide # total "SAI" and "LAI" values # across all residue pools. 0.02 0.00 # # acxrow(s), R, Growing crop row spacing (meters) # Use value of 0.0 if not planted in rows, # e.g. broadcast seeded # ac0rg(s) , I, Specify seed location (0=furrow,1=ridge) # Value doesn't matter if no ridges exist 0.3, 0 # # abffcv(s), R, Flat biomass cover (m^2/m^2) 0.0 # # # +++ SOIL +++ # # nslay(s), I, Number of soil layers (1-100) # NOTE: Only surface soil layer necessary 1 # # NOTE: Remaining SOIL inputs are repeated on each input line # for each layer specified # # aszlyt(l,s), R, Thickness (mm) 1000.0 # # asdblk(l,s), R, Bulk density of soil layer (Mg/m^3) # Typical Range: >0.0 to 10.0 1.8 # asfsan(l,s), R, Fraction of sand content in soil layer (Mg/Mg) # Range: 0.0 to 1.0 (sand + silt + clay = 1.0) 0.90 # asfvfs(l,s), R Fraction of very fine sand in soil layer (Mg/Mg) # Range: 0.0 to 1.0 (fraction of total soil < 2.0 mm) 0.21 # asfsil(l,s), R, Fraction of silt content in soil layer (Mg/Mg) # Range: 0.0 to 1.0 (sand + silt + clay = 1.0) 0.08 # asfcla(l,s), R, Fraction of clay content in soil layer (Mg/Mg) # Range: 0.0 to 1.0 (sand + silt + clay = 1.0) 0.02 # # asvroc(l,s), R, Rock volume in soil layer (m^3/m^3) # Range: 0.0 to 1.0 0.30 # # asdagd(l,s), R, Average aggregate density for the soil layer (Mg/m^3) # Typical Range: 0.5 to 2.5 1.8 # aseags(l,s), R, Average dry aggregate stability for the soil layer [ln(J/kg)] # Typical Range: 0.1 to 7.0 2.50 # # ---- Size distribution of soil aggregates ---- # GMD - Geometric Mean Diameter of aggregates # GSD - Geometric Mean Standard Deviation of aggregates # # aslagm(l,s), R, GMD of aggregate sizes in soil layer (mm) # Typical Range: 0.03 to 30.0 0.47 # aslagn(l,s), R, Minimum aggregate size in soil layer (mm) # Typical Range: 0.001 to 5.0 0.043 # aslagx(l,s), R, Maximum aggregate size in soil layer (mm) # Typical Range: 1.0 to 1000.0 89.8 # as0ags(l,s), R, GSD of aggregate sizes in soil layer (mm/mm) # Typical Range: 1.0 to 40.0 12.0 # # +++ SOIL SURFACE +++ # # asfcr(s), R, Surface crust fraction (m^2/m^2) # Range: 0.0 to 1.0 # aszcr(s), R, Surface crust thickness (mm) # Typical Range: 0.0 to 23.0 # asflos(s), R, Fraction of crusted surface with loose material on top of crust (m^2/m^2) # Range: 0.0 to 1.0 # asmlos(s), R, Mass of loose material on top of crust (kg/m^2) # Typical Range: 0.0 to 3.0 # asdcr(s), R, Density of soil crust (Mg/m^3) # Typical Range: 0.6 to 2.0 # asecr(s), R, Dry crust stability [ln(J/kg)] # Typical Range: 0.0 to 7.0 0.6 7.0 0.2 0.4 0.1 1.0 # # aslrr(s), R, Allmaras random roughness (mm) # Typical Range: 1.0 to 60.0 5.0 # aszrgh(s), R, Ridge height (mm) # Typical Range: 0.0 to 500.0 # asxrgs(s), R, Ridge spacing (mm) # Typical Range: 0.0 to 2000.0 # asxrgw(s), R, Ridge width (mm) # Typical Range: 0.0 to 4000.0 # asxrgo(s), R, Ridge orientation (degrees) # Range: 0.0 to 179.99 # NOTE: If no ridges, then specify 0.0 for height, width and spacing 0.0 0.0 0.0 0 # asxdks(s), R, Dike spacing (mm) # Typical Range: 0.0 to 1000.0 # NOTE: If no dikes, then specify 0.0 0.0 # # +++ HYDROLOGY +++ # # ahzsnd(s), R, Snow depth (mm) # Typical Range: 0.0 to 1000.0 0.0 # # ahrwcw(l,s), R, Wilting point water content of soil layer (Mg/Mg) # Typical Range: 0.0 to 0.25 0.077 # # ahrwca(l,s), R, Current water content of soil layer (Mg/Mg) # Typical Range: 0.0 to 0.50 0.0 # # # ahrwc0(h,s), R, Surface layer water content (Mg/Mg) # Typical Range: 0.0 to 0.50 # 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 +++ # # awzypt, R, Average annual precipitation (mm) # 300.0 # # awdair, R, Air density (kg/m^3) # Typical Range: 0.7 to 1.5 1.2 # # awadir, R, Wind direction (degrees) measured clockwise from North # Typical Range: 0.0 to 359.9 270.0 # # ntstep, I, Number of intervals/day to run EROSION # Range: 24 to 96 # # NOTE: ntstep = 24 means hourly updates # ntstep = 48 means 30 minute updating # ntstep = 96 means 15 minute updating 24 # # anemht, R anemometer height (m) # Typical Range: 0.5 to 30.0 # awzzo, R aerodynamic roughness at anemometer site (mm) # Typical Range: 0.5 to 2000.0 # wzoflg, I (global variable) zo location flag # (flag =0 - at weather station location - zo is a constant) # (flag = 1 - on field location - zo varies based on field surface) 10.0, 10.00 0 # # wflg, I, Wind/Weibull flag # (0 - read in Weibull parameters, 1 - read in wind speeds) 1 # # NOTE: This is only present when (wflg=0) # wfcalm, R, Fraction of time winds are calm (hr/hr) # Range: 0.0 to 1.0 # wuc, R, Weibull "c" factor (m/s) # Typical Range: >0.0 to 30.0 # w0k, R, Weibull "k" factor # Typical Range: 1.0 to 3.0 # 0.217 7.125 2.971 <--- Example data line for wind expressed as Weibull parameters # # NOTE: The remaining data is only present when (wflg=1) # awu(i), R, Wind speed for (ntstep) intervals (m/s) # Typical Range: 0.0 to 30.0 # # NOTE: We can read multiple lines with 6 values per line # Wind data should be AVERAGES for the period. # Hourly averages will often under estimate wind erosion. # 30 minute averages or shorter time interval is more suitable. 8.181 4.068 4.068 4.426 5.052 5.052 4.739 4.292 4.515 3.353 3.621 2.280 5.275 6.750 7.242 7.868 9.835 13.814 17.211 12.651 11.712 12.964 10.014 8.583 #