SWMM4

SWMM 4.45A(beta) RUNOFF DATA FILE


This is the data input file for a CDM/CDMSmith Pittsburgh versionSWMM4 model from 2000 to 2004. It has additional capabilities beyond SWMM 4.4 and was made as a bridge between SWMM 4 and SWMM 5, thus the 4.5 number.

 

* <<<<<<<<< SWMM 4.45A(beta) RUNOFF DATA FILE >>>>>>>>>
*
*
* This is an input data file to the SWMM 4.4 Runoff Block for
* modeling watershed quantity and quality. All lines with an
* asterisk in column 1 are comment lines and are ignored
* by the program.
*
* Input data are free format and may be up to 230 columns wide.
* You must have a value for every data column even if the program
* will not actually use a given value. A slash (/) may be used
* to indicate that remaining fields should be filled with “no data
* entry” or null-entry. This almost always means those input
* parameters will be zero. (See example for data group L1.) A very
* common data input error is to accidentally omit required parameters
* at the end of a data group. There must be at least one space
* (or comma) between every input value.
*
* Caution! Data lines that are “wrapped around” (continued on
* two or more lines) should have a blank in column 1, unless a
* card identifier is needed.
*
* Alphanumeric data ($ANUM option) should be enclosed in single
* quotes. These include all references to subcatchment and
* channel/pipe/inlet names.
*
* In general, SWMM parameters with names that begin with the letters
* I,J,K,L,M,N are integers (e.g., NSCRAT() ), following the usual
* Fortran convention, and entered values must not include a
* decimal point.
*
* To avoid literary quotes being printed in output, use $NOQUOTE
* after MM line.
*
* SWMM uses both U.S. customary units and metric units. The
* examples use feet, cfs, acres, inches, inches/hour, and miles/hr.
* If metric is specified substitute meters, cms, hectares,
* millimeters, millimeters/hour, and km/hr.
*============================================================================
* The SW card sets up the interface files to be used or created.
* There is one output file (#9) that will contain the time series
* of flows and pollutant loads for subsequent blocks.
*============================================================================
* NBLOCK JIN(1) JOUT(1)
SW 1 0 9
*============================================================================
* The MM card opens the scratch files to be used by different subroutines.
* Up to 8 scratch files are required by the Runoff Block.
*============================================================================
* NITCH NSCRAT(1) NSCRAT(2) NSCRAT(3) NSCRAT(4) NSCRAT(5) NSCRAT(6) NSCRAT(7)
MM 8 1 2 3 10 11 12 13
* NSCRAT(8)
14
*============================================================================
* The @ command is used to permanently save an interface or
* scratch file. This line should be placed before the first SWMM
* block call. The format of the @ command is as follows:
*============================================================================
*Column 1 Unit number of the Name of the interface
* interface file saved file (any valid DOS filename),
* or utilized including optional path.
*
@ 9 ‘RUNOFF.DNT’
*============================================================================
*Column 1
* $ANUM ==> Use alphanumeric labels for subcatchment and channel/pipe
* labels — WHEREVER ENCOUNTERED AND IN ALL SUBSEQUENT BLOCKS.
* Names (IDs) must be enclosed in single quotes. A maximum length
* of 6 characters for a label is recommended. Longer names (max
* of 8 characters) may not print out correctly but will
* be input OK.
*============================================================================
*Column 1
* $NOQUOTE ==> Omit on-screen and printed literary quotations in SWMM output.
*============================================================================
$RUNOFF Call the RUNOFF block with a ‘$’ in first column.
*============================================================================
* Create title lines for the simulation. There are two title lines
* for the Runoff Block. Titles are enclosed in single quotes.
*============================================================================
* A1 Line :
* Title : Two lines (both with A1 identifier) with heading
* to be printed on output.
* Each line has format A76 (76 characters, maximum).
*============================================================================
A1 ‘RUNOFF example: parabolic channels, groundwater, water quality’
A1 ‘Lake modeled as a wide parabolic channel with outlet weir’
*============================================================================
* The ‘B’ lines are for program control purposes.
*============================================================================
* B1 Line :
* METRIC : Metric input-output.
* = 0, Use U.S. customary units
* = 1, Use metric units. Metric input indicated
* in brackets [] in remainder of this table.
* ISNOW : Snowmelt parameter.
* = 0, Snowmelt not simulated.
* = 1, Single event snowmelt simulation.
* = 2, Continuous snowmelt simulation.
* NRGAG : Number of hyetographs (rain gages),
* Maximum is limited by MAXRG parameter in Tapes.inc
* INFILM : Choice of infiltration equation
* = 0, Horton equation used.
* = 1, Green-Ampt equation used.
* = 2, Horton equation with maximum infiltration
* volume limiting infiltration.
* In this version, the available infiltration volume for the
* Horton option will recover during dry periods.
* = 3, Green-Ampt equation with maximum infiltration
* volume limiting infiltration.
* DON’T USE INFILM = 3 TEMPORARILY. WCH, 6/10/97.
* In this version, 3 for G-A can only be used for single
* event simulation as the infiltration volume is not
* regenerated during dry periods.
* KWALTY : Quality (or erosion) simulated?
* = 0, No.
* = 1, Yes.
* IVAP : Evaporation parameter
* = 0, Evaporation data not read in,
* default rate used of 0.1 in/day [3.0 mm/day].
* = 1, Read monthly evaporation data in Group F1
* in units of inch/day [mm/day].
* = 2, Read monthly evaporation data in Group F1
* in units of inch/month [mm/month].
* = 3, Read monthly evaporation data on lines
* F1 and F2,
* in units of inch/month [mm/month].
* = 4, Read evaporation time series on NSCRAT(3)
* file as created by the TEMP Block of SWMM.
* NOTE! If it is desired to have no (zero) evaporation
* during time steps when it is raining or snowing, input
* IVAP as a negative number, i.e., IVAP = -1, -2, -3
* or -4 instead of a positive number. This option can
* only be used if IVAP not equal to 0. If IVAP < 0,
* there will be zero surface or subsurface evaporation
* during any time step with rain or snow on that
* subcatchment and zero evaporation from all channel/
* pipes if there is rain or snow on any subcatchment.
* Normal evaporation continues when precipitation = 0.
*
* NHR : Hour of day of start of storm (24 hour clock,
* midnight = 00).
* NMN : Minute of hour of start of storm (0 – 59).
* NDAY : Day of month of start of simulation ( 1 – 31).
* MONTH : Month of start of simulation (1 – 12).
* IYRSTR : Year of start of simulation (4 digits).
* If less than 4 digits are entered, then
* program assumes 1900.
* Optional input to control evaporation on
* channels but does not need to be entered.
* If not entered or 0, then the default is to allow
* evaporation as controlled by IVAP.
* IF 1 then evaporation is never allowed from channels.
*
* IVCHAN : 0 – Allow evaporation from channels.
* : 1 – Don’t allow evaporation from channels.
*============================================================================
* METRIC ISNOW NRGAG INFILM KWALTY IVAP NHR NMN NDAY MONTH IYRSTR [IVCHAN]
B1 0 0 1 1 1 1 00 0 1 10 1989
*============================================================================
* B2 Line :
* IPRN(1) : Print control for SWMM input.
* = 0, Print all input data.
* = 1, Do not print channel/pipe, snowmelt,
* subcatchment, or quality data, only control
* information is printed.
* = K, where K equals possible combinations of
* channel/pipe(2), snowmelt(3), subcatchment(4),
* or water quality(5). For example:
* Channel/pipe + subcatchment would be 24,
* Channel/pipe + subcatchment + quality would be 245.
* IPRN(2) : Print control for Runoff Block graphs.
* = 0, Plot all graphs.
* = 1, Do not plot hyetograph(s) (for each gage),
* or inlet hydrograph (sum of all inlets).
* IPRN(3) : Print control for output of SWMM. ‘Totals’
* below refer to precipitation, runoff and all
* quality parameters. Done for each inlet. Daily,
* monthly, and yearly printouts only function if
* simulation is long enough.
* = 0, Do not print daily, monthly, or yearly totals.
* = 1, Monthly and annual totals only, one year
* per page.
* = 2, Daily, monthly and annual totals, two months
* per page. Daily totals are printed whenever
* there is non-zero precipitation and/or runoff.
*
* The following parameter is truly optional and may
* be omitted from line B2 without an error.
*
* IRPNGW = 0, Print up to 10,000 ground water routine error
* messages.
* > 0, Print limit of IPRNGW ground water routine
* error messages.
*============================================================================
* IPRN(1) IPRN(2) IPRN(3) IRPNGW
B2 0 0 1
*============================================================================
* The B3 line contains time step and duration-of-run parameters.
* The program starts at date/time indicated on line B1. It then
* uses time steps WET, WETDRY and DRY to simulate to an ending date/time
* specified by parameter LONG.
*
* B3 Line :
* WET : Wet time step (seconds). WET must be => 1 second.
* Typical: 60-300-900 sec for event simulation; 900
* or 3600 sec for continuous simulation. WET time
* step is used only during time steps with precip.
* WETDRY : Transition (no rain but water on surface or in
* channels) between wet and dry time step (seconds).
* WETDRY is used during 1) residual overland flow
* (no precipitation), 2) residual channel/pipe flow,
* 3) snowmelt, 4) groundwater outflow to channel/
* pipes. WETDRY should be greater than or equal
* to WET and less than or equal to DRY.
* Typical: = WET for event simulation; 3600 – 7200
* for continuous simulation.
* Note, decrease WETDRY toward WET for better
* resolution and lower continuity errors, but at
* the expense of greater computer time during
* continuous simulation.
* DRY : Dry time step (seconds). DRY must be greater
* than or equal to WET. Typical: = WET for event
* simulation; 7200 – 86400 sec for continuous
* simulation. DRY time step principally affects
* groundwater ET and deep percolation and residual
* surface evaporation and infiltration.
*
* Note: DRY and WETDRY time steps are only approximated during
* time intervals with no precipitation. Thus, print-outs may
* occur at intervals that do not correspond exactly to DRY
* or WETDRY.
*
* LUNIT : Units of LONG (simulation length).
* = 0, seconds. = 1, minutes.
* = 2, hours. = 3, days.
* = 4, ending date, a eight figure number
* (year/mo/dy), e.g. 19870730.
* If year is two digits, program assumes 1900.
* LONG : Simulation length (units from LUNIT). A real
* number, not an integer.
*============================================================================
* SIMULATION LENGTH OF 6 DAYS
* WET WET/DRY DRY LUNIT LONG
B3 600. 1200.0 7200. 3 6.0
*============================================================================
* B4 is an optional data group. The B4 data group is used only when the
* user desires to modify one of SWMM’s subcatchment default parameters.
*============================================================================
* B4 Line :
* PCTZER : Percent of impervious area with zero detention
* (immediate runoff). Default = 25%.
* REGEN : For continuous SWMM, infiltration capacity is
* regenerated using a Horton type exponential rate
* constant equal to REGEN*DECAY, where DECAY is the
* Horton rate constant read in for each subcatchment
* in group H1. Default = 0.01. Not required for
* Green-Ampt infiltration.
*============================================================================
* Use line C1 to input general snow input data.
* If ISNOW = 0 in group B1, skip to group D1.
*============================================================================
* C1 Line :
* ELEV : Average watershed elevation, ft, msl [m, msl].
* FWFRAC(1) : Ratio of free water holding capacity to snow depth
* (in. or mm w.e.= water equivalent) on snow
* covered impervious area.
* FWFRAC(2) : Ratio of free water holding capacity to snow depth
* (in. or mm w.e.) on snow covered pervious area.
*============================================================================
* Note: The following parameters are required only for ISNOW=2.
*============================================================================
* FWFRAC(3) : Ratio of free water holding capacity to snow depth
* (in. or mm w.e.) for snow on normally bare
* impervious area.
* SNOTMP : Dividing temperature between snow and rain,
* F [C]. Precipitation occurring at air
* temperatures above this value will be rain,
* at or below will be snow.
* SCF : Snow gage catch correction factor.
* Snow depths computed from NWS precipitation tape
* will be multiplied by this value.
* TIPM : Weight used to compute antecedent temperature index,
* 0 <= TIPM <= 1.0. Low values (e.g., 0.1) give
* more weight to past temperatures. Values > 0.5
* essentially give weight to temperatures only
* during the past day.
* RNM : Ratio of negative melt coefficient to melt
* coefficient. “Negative melt coefficient” is used
* when snow is warming or cooling below the base melt
* temperature without producing liquid melt. RNM is
* usually <= 1.0 with a typical value of 0.6.
* ANGLAT : Average latitude of watershed, degrees north.
* DTLONG : Longitude correction, standard time minus
* mean solar time, minutes (of time).
*============================================================================
* Use line C2 to input average Monthly Wind Speeds. Enter pairs of values
* (month number, wind speed) only for months with potential snow
* melt. Enter values for months in any order. Months not entered
* are assumed to have zero wind.
*============================================================================
* C2 Line :
* NUMB : Enter number of months with wind speed data.
* (Maximum = 12)
* [NOTE. Option on page 69 of User’s Manual to
* set NUMB = 999 to indicate NOAA wind data is
* not valid. Use ISNOW=2 to indicate use of
* NOAA data on NSCRAT(3) from Temp Block.]
* MONTH : Integer number of first month.
* WIND(MONTH): Average wind speed for first month, mi/hr [km/hr].
* . .
* MONTH : Integer number of last month.
* WIND(MONTH): Average wind speed for last month, mi/hr [km/hr].
*============================================================================
* Use line C3 to input Areal Depletion Curve for Impervious Area.
* IF ISNOW=1 IN GROUP B1, SKIP TO DATA GROUP C5.
*============================================================================
* C3 Line :
* ADCI(1) : Fraction of area covered by snow (ASC) at “zero+”
* ratio of snow depth to depth at 100 percent
* cover (AWESI).
* ADCI(2) : Value of ASC for AWESI = 0.1.
* ADCI(3) : Value of ASC for AWESI = 0.2.
* . .
* ADCI(9) : Value of ASC for AWESI = 0.8.
* ADCI(10) : Value of ASC for AWESI = 0.9.
* Note: Program automatically assigns value of ADCI=1.0 when AWESI = 1.0.
*============================================================================
* Use the C4 line to define an Areal Depletion Curve for Pervious Area.
*============================================================================
* C4 Line :
* ADCP(1) : Fraction of area covered by snow (ASC) at “zero+”
* ratio of snow depth to depth at 100 percent cover
* (AWESI).
* ADCP(2) : Value of ASC for AWESI = 0.1.
* ADCP(3) : Value of ASC for AWESI = 0.2.
* . .
* ADCP(9) : Value of ASC for AWESI = 0.8.
* ADCP(10) : Value of ASC for AWESI = 0.9.
* Note: Program automatically assigns value of ADCP = 1.0 when AWESI = 1.0.
*============================================================================
* READ GROUP C5 ONLY IF ISNOW = 1. SKIP TO GROUP D1 IF ISNOW = 2.
*
* For ISNOW = 2 (continuous SWMM), air temperatures are entered
* in the Temp Block. For ISNOW = 1, read an air temperature for each
* time interval DTAIR, for a total of NAIRT values. (Maximum number
* of values = 200. If more are needed, use ISNOW = 2 option.) DTAIR,
* the time step of air temperatures, is not necessarily equal to the
* time steps entered on data group B1. Air temperatures are considered
* constant over the air time step.
*============================================================================
* C5 Line :
* DTAIR : Time interval for input of air temperatures,
* hours. First line only.
* NAIRT : Number of air temperatures read. First line only.
* TAIR(1) : Air temperature during time interval 1, F [C].
* . .
* TAIR(NAIRT): Air temperature during time interval NAIRT, F [C].
*============================================================================
* Line D1 is the first rainfall control line.
*============================================================================
* D1 Line :
* ROPT : Precipitation input option.
* = 0, Read NRGAG hyetographs on E1, E2 and E3
* data groups. (Rain data can be saved permanently
* on NSCRAT(1) using the @ function.)
* = 1, Read processed precipitation file on NSCRAT(1)
* file [not JIN!]. This file is either from the Rain
* Block (earlier saved JOUT file) or from a previous
* run of the Runoff Block (earlier saved NSCRAT(1)
* file). Unless blocks are run as part of a single
* overall SWMM run, access to earlier saved files is
* through the @ function described at the beginning
* of this file.
*============================================================================
* ROPT
D1 0
*============================================================================
* Line E1 is the second rainfall control line.
*============================================================================
* E1 Line :
* KTYPE : Type of precipitation input. Precipitation
* is in units of in./hr [mm/hr] for THISTO minutes or
* hours. Use variable KTIME to select units of time.
* = 0, Read KINC precipitation values per line.
* = 1, Read KINC time and precipitation pairs per line.
* = 2, Read time and NRGAG precipitation values per line.
* KINC : Number of precipitation values or time/precipitation
* pairs per line. Enter any number if KTYPE = 2.
* KPRINT : Print control for precipitation input.
* = 0, Print all precipitation input.
* = 1, Suppress all but summary of precipitation input.
* KTHIS : Variable THISTO option. Data input on E2 lines.
* = 0, precipitation interval (THISTO) is constant.
* = K, where K is the number of variable precipitation
* intervals entered on the E2 data group lines.
* Precipitation values outside the time frame
* of any variable rainfall interval uses THISTO
* as the precipitation interval.
* KTIME : Precipitation time units.
* = 0, time in minutes.
* = 1, time in hours.
* KPREP : Precipitation unit type.
* = 0, intensity, in./hr [mm/hr].
* = 1, total precipitation volume over
* the interval, in. [mm]
* NHISTO : Number of data points for each hyetograph.
* THISTO : Time interval between values (and duration of
* precipitation value), units of KTIME.
* TZRAIN : Initial time of day of precipitation input, units
* of KTIME, or off-set time added to times entered
* in groups E2 and E3. (If first time entered in
* groups E2 and/or E3 is 0.0, TZRAIN will ordinarily
* correspond to time of start of storm entered on
* group B1.)
* Caution. When precipitation times are not included
* with rainfall values, TZERO will usually correspond
* to time of day of start of storm entered on line B1
* or else there is a danger that rainfall times may
* not overlap with simulation times and zero runoff
* will result.
*============================================================================
* KTYPE KINC KPRINT KTHIS KTIME KPREP NHISTO THISTO TZRAIN
E1 1 1 0 0 1 1 24 1.0 0.0
*============================================================================
* Line E2 lists the variable rainfall interval information.
* Required only if KTHIS > 0. Enter variable precipitation intervals,
* for a total of KTHIS intervals. Do not repeat the E2 line identifier
* after the first line. (Wrap around, leaving at least the first column
* blank in each succeeding row.) This data group is used
* to interleave rainfall records of differing intervals, for example, a
* period of 5 minute rainfall between periods of 15 minute rainfall.
*============================================================================
* E2 Line :
* WTHIS(1,1) : Start time for first variable precipitation
* interval. Units of KTIME.
* WTHIS(1,2) : End time for first variable precipitation
* interval. Units of KTIME.
* WTHIS(1,3) : Length of THISTO for the first precipitation
* interval. Units of KTIME.
* . .
* WTHIS(KTHIS,1): Start time for last variable precipitation
* interval. Units of KTIME.
* WTHIS(KTHIS,2): End time for last variable precipitation
* interval. Units of KTIME.
* WTHIS(KTHIS,3): Length of THISTO for the last precipitation
* interval. Units of KTIME.
*============================================================================
* Use line E3 to input precipitation input. Input is a function
* of the parameter KTYPE on data group E1.
*
* Note: If ISNOW = 1, snowfall during a time step may be entered as
* a negative value. Units are in. [mm] water equivalent/hr.
*============================================================================
* Precipitation input if KTYPE = 0.
*
* KINC precipitation values per line, up to NHISTO values.
* Repeat group E3 for each hyetograph, up to NRGAG times.
*
* Note, you must include the E3 identifier at the beginning of each
* group of KINC rainfall entries. An individual line of KINC entries
* may be “wrapped around,” but each new line of KINC entries must include
* the E3 identifier.
*============================================================================
* E3 Line :
* RAIN(1) : Rainfall intensity, first interval, in./hr [mm/hr].
* . .
* RAIN(KINC) : Rainfall intensity, last interval per line,
* in./hr [mm/hr].
*============================================================================
* Precipitation input if KTYPE = 1.
*
* Read KINC pairs per line, up to NHISTO values.
* Repeat group E3 for each hyetograph, up to NRGAG times.
*============================================================================
* E3 Line :
* REIN(1) : Time of first precipitation. Units of KTIME.
* REIN(2) : Precipitation in./hr [mm/hr], for first interval.
* . .
* REIN(2*KINC-1): Time of last precipitation. Units of KTIME.
* REIN(2*KINC) : Precipitation for last interval, in./hr [mm/hr].
*============================================================================
* Precipitation input if KTYPE = 2.
*
* Read NRGAG precipitation values per line. Repeat NHISTO times.
*============================================================================
* E3 Line :
* REIN(1) : Time of precipitation. Units of KTIME.
* REIN(2) : Precipitation, first raingage, in./hr [mm/hr].
* . .
* REIN(NRGAG+1) : Precipitation, last raingage, in./hr [mm/hr].
*============================================================================
* STEP-FUNCTION HYETOGRAPH
* TIME=REIN(1) RAIN=REIN(2)
E3 0.0 0.5
E3 1.0 1.0
E3 2.0 0.2
E3 10.0 0.4
E3 11.0 0.2
E3 20.0 0.1
E3 21.0 0.2
E3 22.0 0.3
E3 23.0 0.1
E3 30.0 1.5
E3 40.0 0.2
E3 41.0 0.2
E3 42.0 0.2
E3 50.0 2.0
E3 60.0 0.5
E3 61.0 0.4
E3 63.0 0.2
E3 64.0 0.1
E3 120.0 0.2
E3 121.0 0.3
E3 122.0 0.4
E3 123.0 0.2
E3 124.0 0.1
E3 125.0 0.1
*============================================================================
* Use the F1 line to input evaporation data if IVAP >= 1 on group B1.
*============================================================================
* F1 Line :
* Note, units depend on value of IVAP. This example is for IVAP = 1.
* VAP(1) : Evaporation rate for month 1 (January)
* in./day [mm/day].
* . .
* VAP(12) : Evaporation rate for month 12 (December)
* in./day [mm/day].
*============================================================================
* EVAPORATION DATA
F1 0.1 0.2 0.2 0.2 0.3 0.4 0.5 0.5 0.4 0.3 0.2 0.1
*============================================================================
* Use the F1 and F2 lines to input evaporation data if
* IVAP = 3 on data group B1.
*============================================================================
* F1 line :
* Only for IVAP = 3:
* NVAP(1) : Start year of evaporation data (4-digit integer)
* If less than 4 digits are entered program assumes 1900.
* NVAP(2) : Number of months of evaporation data to be
* entered (maximum = 600). First month must be for January.
* NVAP(1) NVAP(2)
*F1 80 120
*============================================================================
*
* Line F2: Read 12 monthly values per line to a maximum of 600 values.
* JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC
*F2 1.1 2.6 2.5 4.0 5.1 6.7 6.0 4.9 4.0 3.1 2.3 0.9
*F2 0.9 2.1 2.8 4.5 4.8 6.3 6.0 5.4 3.8 3.5 1.5 0.5
*F2 1.2 1.3 2.1 3.4 3.5 7.0 5.5 5.6 3.2 2.8 2.2 1.1
*F2 0.5 0.9 3.2 4.1 4.6 6.1 6.2 5.8 2.9 3.2 1.8 0.7
*F2 0.8 0.9 2.4 4.2 4.7 5.7 5.8 5.9 4.4 3.1 1.7 0.8
*F2 1.3 2.4 2.5 4.0 5.1 6.7 6.0 4.5 4.0 3.1 2.3 0.9
*F2 0.7 2.1 2.8 4.2 4.8 6.3 6.0 5.4 3.8 3.5 1.5 0.6
*F2 1.2 1.3 2.1 3.4 3.5 7.0 5.5 5.6 3.2 2.6 2.2 1.1
*F2 0.6 0.9 3.2 4.1 4.4 6.1 6.2 5.8 2.9 3.2 1.8 0.6
*F2 0.9 1.6 2.4 4.2 4.7 5.7 5.2 5.9 4.4 3.1 1.7 0.5
*===========================================================================
* Rainfall-Dependent Infiltration/Inflow (RDII or I/I) Data, Lines F3 and F4.
* New, 9/4/93. Chuck Moore
* Camp, Dresser and McKee, Inc., Annandale, VA
*
* These lines, plus H5 lines, define triangular unit hydrographs (UH)
* to compute subcatchment I/I response from rainfall record on NSCRAT(1).
* The response is computed before the time step simulation and stored on
* NSCRAT(8) (required if this I/I procedure is used).
* An initial abstraction of up to DSTORE in [mm] is subtracted from
* rainfall before computing rainfall excess at each time step. The
* initial abstraction is regenerated during dry weather at a rate of
* DREC in/day [mm/day].
* Up to five sets of three triangular UHs may be input in the F4 lines.
* Any subcatchment may use a fraction (defined on H5 lines) of flow
* produced from each of three UHs selected from any of the up to
* five sets input on the F4 lines.
*
* The triangles are defined by TP = time to peak and K = ratio of recession
* limb to TP, so that the time base = TP*(1+K). Peak flow, Qp, is then
* calculated in the program so that the volume (area) of the triangle =
* 1 cfs/ac-in or 1 cms/ha-mm.
*
* Three triangles may be used (all starting at the same time) so that
* one may define a fast response, one a delayed response, and one a
* lengthy response, if desired.
*
* The time step TSTEP (UH duration) must equal the rainfall time
* step (THISTO). Rainfall can be input on the E-lines as usual, or
* rainfall already stored on NSCRAT(1) may be used, i.e., from prior
* Runoff or Rain Block runs.
* Time step TSTEP2 is the time step used for computation of the UH
* response and should be an integer fraction of (or equal to) TSTEP.
* TSTEP2 = WET time step in the Runoff Block is a good idea.
* There is a limit of 300 UH points for a given rainfall increment. An
* error message is printed if this is exceeded, which can be caused by
* too long a time base and too short a value of TSTEP2.
* The RDII routine is designed to be compatible with the variable time
* step used in the Runoff Block. Values of I/I are linearly interpolated
* from the file on NSCRAT(8) at each time during the simulation.
*
* H5 lines may be entered only for desired subcatchments.
* If quality is simulated, constant concentrations for I/I are entered
* on data line J6.
*
*
*==========================================================================
* Lines F3 and F4 are optional and may be omitted.
*
* F3 Line : Line identifier
* IIRDII : = 0, compute new I/I response from rainfall and
* store on NSCRAT(8).
* = 1, use I/I response already calculated during previous
* run. (NSCRAT(8) must be defined on @-line)
* TSTEP : Time interval for rainfall, hr. (Must equal value from
* E-lines or from Rain Block.) TSTEP is the duration of
* each of the three UHs.
* TSTEP2 : Time step for computation of I/I response, hr. Must be
* equal to or integer fraction of TSTEP.
*==========================================================================
* IIRDII RAINFALL TIME STEP (hr) RDII TIME STEP (hr)
F3 0 1.0 0.25
*==========================================================================
* Enter up to five F4 lines, with identifier.
*
* F4 Line : Line identifier
* NRDHYET : Number of hyetograph (rain gage) to use for these
* triangular UHs.
* RDIIT : Time to peak of triangular UH, hr.
* RDIIK : Ratio of recession limb to time to peak.
* DSTORE : Maximum initial abstraction, to compute rainfall excess for
* UH, in. [mm].
* STORAGE : Initial storage (max = DSTORE), in. [mm].
* DREC : Recovery rate for storage (initial abstraction) during dry
* time steps, in./day [mm/day].
*
* IF NRDHYET is entered as a negative number, the values on the
* first F4 card are for January. The program then reads
* 11 additional F4 cards with only the parameters RDIIT, RDIIK,
* DSTORE, STORAGE, and DREC for each month. Storage is
* used only for the first timestep.
*=========================================================================
*REPEAT F4 DATA LINES FOR UP TO 5 SETS OF I/I BASE RESPONSE CURVES
* NNRDHYET RDIIK STORAGE etc.
* RDIIT DSTORE DREC etc.
* show input for monthly factors
F4 -1 1.0 2.0 0.1 0.1 0.05 3.0 2.0 0.1 0.1 0.05 10.0 2.0 1.1 0.1 0.05
F4 1.0 2.0 0.1 0.1 0.05 3.0 2.0 0.1 0.1 0.05 10.0 2.0 1.1 0.1 0.05
F4 1.0 2.0 0.1 0.1 0.05 3.0 2.0 0.1 0.1 0.05 10.0 2.0 1.1 0.1 0.05
F4 1.0 2.0 0.1 0.1 0.05 3.0 2.0 0.1 0.1 0.05 10.0 2.0 1.1 0.1 0.05
F4 1.0 2.0 0.1 0.1 0.05 3.0 2.0 0.1 0.1 0.05 10.0 2.0 1.1 0.1 0.05
F4 1.0 2.0 0.2 0.1 0.05 3.0 2.0 0.2 0.1 0.05 10.0 2.0 1.5 0.1 0.05
F4 1.0 2.0 0.3 0.1 0.05 3.0 2.0 0.3 0.1 0.05 10.0 2.0 1.9 0.1 0.05
F4 1.0 2.0 0.3 0.1 0.05 3.0 2.0 0.3 0.1 0.05 10.0 2.0 1.9 0.1 0.05
F4 1.0 2.0 0.3 0.1 0.05 3.0 2.0 0.3 0.1 0.05 10.0 2.0 1.5 0.1 0.05
F4 1.0 2.0 0.2 0.1 0.05 3.0 2.0 0.2 0.1 0.05 10.0 2.0 1.1 0.1 0.05
F4 1.0 2.0 0.1 0.1 0.05 3.0 2.0 0.1 0.1 0.05 10.0 2.0 1.1 0.1 0.05
F4 1.0 2.0 0.1 0.1 0.05 3.0 2.0 0.1 0.1 0.05 10.0 2.0 1.1 0.1 0.05
* slower response
F4 1 2.0 2.0 0.1 0.1 0.05 4.0 2.0 0.1 0.1 0.05 15.0 2.0 1.1 0.1 0.05
* even slower response
F4 1 4.0 2.0 0.1 0.1 0.05 5.0 2.0 0.1 0.1 0.05 20.0 2.0 1.1 0.1 0.05
*============================================================================
* Enter Channel/Pipe data on line G1.
*
* Channel/pipe data: one line per channel/pipe (if none, leave out).
* Maximum number of channels or pipes plus inlets is defined by parameter
* NG on the ‘TAPES.INC’ common. An inlet is any location identified by NGTO
* on a G1 or H1 line that is not listed in group G1 as a channel or pipe.
* All inlets (and only these inlets) are saved on the output interface
* file, if JOUT > 0.
*
* Note: Variables with asterisks can be modified using the Default/Ratio
* option. A -1 entered for NAMEG means non-zero entries for data with
* asterisks are ratios, by which subsequent entries for those parameters
* will be multiplied. E.g., enter a ratio of 1.2 to increase all
* parameters by 20%. A -2 entered for NAMEG means non-zero entries for
* data with asterisks are default values. If a data line has a zero for
* a parameter for which a default value has been defined, the parameter
* is assigned this value. As many ratio and default lines may be inserted
* as desired within a data group. Of course, if the alphanumeric option
* is being used, the -1 or -2 should be enclosed in single quotes, e.g.,
* ‘-1’ or ‘-2’.
*============================================================================
* G1 Line :
* NAMEG : Channel/pipe number or name.
* NGTO : Channel/pipe or inlet number or name for drainage.
* NPG=NP : Type of channel or pipe.
* = 1 for trapezoidal channel,
* = 2 for circular pipe,
* = 3 for dummy channel/pipe, inflow = outflow,
* = 4 for parabolic channel,
* = 5 for trapezoidal channel with weir or orifice
* (follow with G2 data group),
* = 6 for circular pipe with weir or orifice
* (follow with G2 data group), and
* = 7 for parabolic channel with weir or orifice
* (follow with G2 data group).
*============================================================================
* The following parameters are not used if NP = 3.
*============================================================================
* GWIDTH* : Bottom width of trapezoidal channel, diameter.
* of pipe, or top width of parabolic channel, ft [m].
* GLEN* : Length of channel/pipe, ft [m].
* G3* : Invert slope, ft/ft (dimensionless).
* GS1 : Left-hand side slope, ft/ft. (Slope = horiz./vert.).
* GS2 : Right-hand side slope, ft/ft.
* G6* : Manning’s roughness coefficient.
* DFULL* : Depth of channel when full, ft [m].
* (N.R. if NP equals 2, 3, or 6)
* GDEPTH* : Starting depth of pipe/channel, ft [m].
*============================================================================
* CHANNEL 20 – CIRCULAR PIPE
* CHANNEL 30 – TRAPEZOIDAL CHANNEL
* CHANNEL 1 – PARABOLIC CHANNEL
* CHANNEL 2 – LAKE (PARABOLIC) CHANNEL WITH WEIR
* NAMEG NGTO NPG GWIDTH GLEN G3 GS1 GS2 G6 DFULL GDEPTH
*G1 -1 0 0 0 0.8 0 0 0 0 0 0
*G1 -2 0 0 0 0 0 0 0 0.020 0 0
G1 20 2 2 3.0 1000. 0.01 0.00 0.00 0.020 0.0 0.0
G1 30 2 1 10.0 500. 0.01 5.00 4.00 0.030 8.0 0.0
G1 1 2 4 20.0 5000. 0.01 0.00 0.00 0.020 10.0 0.0
*Note, initial water level below weir crest.
G1 2 3 7 500.0 2000. 0.01 0.00 0.00 0.020 10.0 2.5
*============================================================================
* Enter control structure data on line G2.
*
* Note: A G2 data group must follow a G1 line if NPG is greater than 4.
* Then continue with further G1 lines and/or G1-G2 pairs.
*============================================================================
* G2 Line :
* WTYPE : Type of weir/orifice,
* = 0, Broad or narrow crested weir,
* = 1, V-notched weir, or
* = 2, Orifice.
* WELEV : Elevation of weir (bottom of notch for V-notch) or
* of orifice centerline, referenced to bottom of
* channel/pipe, ft [m].
* WDIS : Discharge coefficient of the weir or orifice
* (parameter C in equations 4-5, 4-6, 4-7). Units
* for equations 4-5 or 4-6: ft1/2/sec [m1/2/sec].
* Parameter Cd in equation 4-7 is dimensionless.
* SPILL : Weir length (e.g., width of spillway) for a broad
* or narrow crested weir, ft [m]. The angle (degrees)
* of the notch for a V-notch weir. The cross
* sectional area of the outflow orifice, ft2 [m2].
*============================================================================
* WTYPE WELEV WDIS SPILL
G2 0 3.0 3.3 10.0
*============================================================================
* Enter Subcatchment Data on line H1. Repeat for each subcatchment
* (Maximum of NW different subcatchments).
*
* Note: Variables with asterisks can be modified using
* the Default/Ratio option. If any of H2-H5 lines follow for
* this subcatchment, must have a non-ratio/default line follow
* this ratio/default line. That is, cannot have H1 ratio/default
* line followed immediately by H2 or H5 line.
*============================================================================
* H1 Line :
* JK : Hyetograph number (based on the order
* in which they are input, in Group E3).
* NAMEW : Subcatchment number or name.
* NGTO : Channel/pipe or inlet (manhole) number for drainage.
* WW(1)* : Width of subcatchment, ft [m].
* This term actually refers to the physical width of
* overland flow in the subcatchment and may be estimated
* as illustrated in the text or by ratio of subcatchment
* area to average length of overland flow.
* WAREA* : Area of subcatchment, acres [ha].
* WW(3)* : Percent imperviousness of subcatchment,
* (percent hydraulically effective or directly
* connected impervious area).
* WSLOPE* : Ground slope, ft/ft (dimensionless).
* WW(5)* : Impervious area Manning’s roughness.
* WW(6)* : Pervious area Manning’s roughness.
* WSTORE1* : Impervious area depression storage, in. [mm].
* WSTORE2* : Pervious area depression storage, in. [mm].
*============================================================================
* Last three parameters on line H1 if Horton equation
* is used, INFILM = 0 on data group B1.
*============================================================================
* WLMAX* : Maximum initial infiltration rate, in./hr [mm/hr].
* WLMIN* : Minimum (asymptotic) infiltration rate, in./hr [mm/hr].
* DECAY* : Decay rate of infiltration in Horton’s equation, 1/sec.
*============================================================================
* Last three parameters on line H1 if Green-Ampt equation
* is used, INFILM = 1 on group B1.
*============================================================================
* SUCT* : Average capillary suction, in. [mm] of water.
* HYDCON* : Saturated hydraulic conductivity of soil,
* in./hr [mm/hr].
* SMDMAX* : Initial moisture deficit for soil, volume
* air/volume voids (fraction).
* IF INFILM = 2 or 3, the maximum infiltration volume is entered
* as the last variable on the H1 line.
* Only needed for INFILM = 2 or 3.
* RMAXINF : Maximum infiltration volume, in. [mm] of water.
*============================================================================
* =====> SURFACE WATER DATA
* JK NAMEW NGTO WIDTH AREA %IMP SLP IMPN PERVN IDS PDS SUCT HYDCON SMDMAX RMAXINF
H1 1 100 1 100.0 300. 20.0 .001 .04 .30 .05 .10 4.00 1.00 .34 5.0
*============================================================================
* Input Groundwater Subcatchment Data on lines H2, H3 and H4.
*
* Data groups H2, H3, and H4 describe the groundwater portion of the
* subcatchment. They should follow the correct H1 data group line. There
* are a maximum of NGW subcatchments with groundwater simulation allowed.
*
* Note: Variables with asterisks can be modified using the
* the Default/Ratio option. Indicator variable is NMSUB. Be sure
* to include all three H2-H4 lines if entering Default/Ratio data.
*============================================================================
* H2 Line :
* NMSUB : Subsurface subcatchment indicator variable,
* must be same as preceding NAMEW on H1 line.
* NGWGW : Number or name of inlet, channel or pipe for
* subsurface drainage. Does not have to be the
* same as preceding NGTO for surface runoff.
* ISFPF : Indicator variable for saving soil moisture,
* water table elevation and outflow for printing.
* = 0, do not save subsurface information, or
* = 1, save subsurface information for printing.
* ISFGF : Indicator variable for saving soil moisture,
* water table elvation and outflow for graphing.
* = 0, do not save subsurface information, or
* = 1, save subsurface information for graphing.
*============================================================================
* Note: See Figure X-1 for definition of elevation variables.
*============================================================================
* BELEV* : Elevation of bottom of water table aquifer, ft [m].
* GRELEV* : Elevation of ground surface, ft [m].
* STG* : Elevation of initial water table stage, ft [m].
* BC* : Elevation of channel bottom or threshold stage
* for groundwater flow, ft [m].
* TW* : Channel water influence parameter
* >= BC, average elevation of water in channel
* or pipe over run, ft [m] or,
* < 0, (e.g., -1) channel water influence will be
* determined by depth in channel or pipe at the end
* of the previous time step.
*============================================================================
* GROUNDWATER DATA
* NMSUB NGWGW ISFPF ISFGF BELEV GRELEV STG BC TW
H2 100 1 1 0 0.0 20.0 5.00 5.00 5.00
*============================================================================
* Input Groundwater Flow Coefficients And Exponents from
* (Equations X-24 and X-25) on line H3.
*============================================================================
* H3 Line :
* A1* : Groundwater flow coefficient, in/hr-ft^B1 [mm/hr-m^B1].
* B1* : Groundwater flow exponent, dimensionless.
* A2* : Coefficient for channel water influence,
* in/hr-ft^B2 [mm/hr-m^B2].
* B2* : Exponent for channel water influence, dimensionless.
* A3* : Coefficient for the cross product between groundwater
* flow and channel water, in/hr-ft^2 [mm/hr-m^2].
* POR* : Porosity expressed as a fraction.
* WP* : Wilting point expressed as a fraction.
* FC* : Field capacity expressed as a fraction.
* HKSAT* : Saturated hydraulic conductivity, in./hr [mm/hr].
* TH1* : Initial upper zone moisture expressed as a fraction.
*============================================================================
* A1 B1 A2 B2 A3 POR WP FC HKSAT TH1
H3 4.5E-5 2.6 0.0 1.0 0.0 .46 .15 .30 5.0 .301
*============================================================================
* Input more groundwater parameters on line H4.
*============================================================================
* H4 Line :
* HCO* : Hydraulic conductivity vs. moisture content
* curve-fitting parameter (Eqn. X-21), dimensionless.
* PCO* : Average slope of tension versus soil
* soil moisture curve (see Figures X-2, X-3 and
* X-4), ft/fraction [m/fraction].
* CET* : Fraction of maximum ET rate assigned to the upper zone.
* DP* : Coefficient for unquantified losses,
* (Eqn. X-23), in./hr [mm/hr].
* DET* : Maximum depth over which significant lower zone
* transpiration occurs, ft [m].
*============================================================================
* HCO PCO CET DP DET
H4 10. 15. 0.35 2.E-03 14.0
*============================================================================
* Define subcatchment response to infiltration/inflow on H5 line.
* See explanation of F3 and F4 lines for additional detail.
*
* Any desired subcatchment may define response parameters. Generated I/I
* will enter in channel/pipe or inlet NGTO. The sewered area defined on
* the H5 line is used, not the subcatchment area. The response is defined
* as a fraction RDIIR from each of the three triangular UHs defined on
* one of the F4 lines. The fractions do not have to sum to 1.0.
*============================================================================
* Input H5 line only for a subcatchment for which I/I response is desired.
* Lines H1-H5 follow in groups. Do not “cluster” all H5 lines together.
*
* H5 Line : Line identifer
* SEWAREA : Sewered area, ac [ha]
* RDIIR(1) : Fraction of first UH response toward total I/I response.
* RDIIR(2) : Fraction of second UH response toward total I/I response.
* RDIIR(3) : Fraction of third UH response toward total I/I response.
* ICURVE : Indicator for which set of three UHs to use, from sequence
* of F4 lines. (E.g., if ICURVE = 2, use second set for
* this subcatchment.)
*
* If SEWAREA is entered as a negative number, then values on first H5
* line is for the month of January. Eleven additional H5 lines will then
* be read for each month of the year with the following parameters only:
* RDIIR(1),RDIIR(2),RDIIR(3)
*============================================================================
* SEWAREA RDIIR1 RDIIR2 RDIIR3 ICURVE
H5 -100 0.03333 0.03333 0.0333 1 * January
H5 0.03333 0.03333 0.0333 * February
H5 0.03333 0.03333 0.0333 * March
H5 0.03333 0.03333 0.0333 * April
H5 0.03333 0.02222 0.0222 * May
H5 0.03333 0.02222 0.0111 * June
H5 0.03333 0.02222 0.0 * July
H5 0.03333 0.01111 0.0 * August
H5 0.03333 0.02222 0.0 * September
H5 0.03333 0.02222 0.0 * October
H5 0.03333 0.03333 0.0111 * November
H5 0.03333 0.03333 0.0222 * December
*============================================================================
* Additional subcatchments with no subsurface data or I/I data.
* JK NAMEW NGTO WIDTH AREA %IMP SLP IMPN PERVN IDS PDS SUCT HYDCON SMDMAX RMAXINF
H1 1 200 20 50.0 100. 30.0 .001 .04 .30 .05 .10 3.00 1.00 .34 5.0
H1 1 300 30 150.0 400. 10.0 .002 .04 .30 .05 .10 4.00 0.50 .34 5.0
*============================================================================
* Enter Subcatchment Snow Input Data on data groups I1 and I2 (if modeled).
*============================================================================
*
* Note: If ISNOW = 0, skip to group J1.
* If ISNOW = 1, read only group I1.
* If ISNOW = 2, read both groups I1 and I2, in pairs.
*
* Order of subcatchments must be same as in group H1, and there
* must be snow data group(s) for each H1 line. All snow-depth related
* parameters refer to depth of snow water equivalent (w.e.).
*
* Note: Variables with asterisks can be modified
* using the Default/Ratio option.
*============================================================================
* I1 Line :
* JK1 : Subcatchment number or name. Must correspond to NAMEW
* entered in Group H1.
* SNN1 : Fraction of impervious area with 100 percent
* snow cover (ISNOW = 1) or subject to areal
* depletion curve (ISNOW = 2).
* SNCP(N) : Fraction of pervious area subject to 100 percent
* snow cover (ISNOW = 1). N.R. if ISNOW = 2.
* WSNOW(N,1) : Initial snow depth of impervious area that is
* normally snow covered, in. water equivalent
* [mm w.e.]
* WSNOW(N,2) : Initial snow depth on pervious area,
* in. w.e. [mm w.e.].
* FW(N,1) : Initial free water on snow covered impervious
* area, in. [mm].
* FW(N,2) : Initial free water on snow covered pervious
* area, in. [mm].
* DHMAX(N,1)* : Melt coefficient (ISNOW = 1) or maximum melt
* coefficient, occurring on June 21 (ISNOW = 2)
* for snow covered impervious area,
* in. w.e./hr-F [mm w.e./hr-C].
* DHMAX(N,2)* : Melt coefficient (ISNOW = 1) or maximum melt
* coefficient, occurring on June 21 (ISNOW = 2)
* for snow covered pervious area,
* in. w.e./hr-F [mm w.e./hr-C].
* TBASE(N,1)* : Snow melt base temperature for snow covered
* impervious area, F [C].
* TBASE N,2)* : Snow melt base temperature for snow covered
* pervious area, F [C].
*============================================================================
* Enter Subcatchment Snow Input Data on data group I2 if ISNOW = 2.
*============================================================================
* I2 Line :
* JK2 : Subcatchment number or name. Must correspond to JK1
* on Line I1 and NAMEW in Group H1.
* WSNOW(N,3) : Initial snow depth on impervious area that is
* normally bare, in. [mm].
* FW(N,3) : Initial free water on impervious area that is
* normally bare, in. [mm].
* DHMAX(N,3)* : Maximum melt coefficient occurring on June 21,
* for snow on normally bare impervious area,
* in. w.e./hr-F [mm w.e./hr-C].
* TBASE(N,3)* : Snow melt base temperature for normally bare
* impervious area, F [C].
* DHMIN(N,1)* : Minimum melt coefficient occurring on December 21
* for snow covered impervious area, in. w.e./hr-F
* [mm w.e./hr-C].
* DHMIN(N,2)* : Minimum melt coefficient occurring on December 21
* for snow covered pervious area, in. w.e./hr-F
* [mm w.e./hr-C].
* DHMIN(N,3)* : Minimum melt coefficient occurring on December 21
* for snow on normally bare impervious area,
* in. w.e./hr-F [mm w.e./hr-C].
* SI(N,1)* : Snow depth above which there is 100 percent cover
* on snow covered impervious areas, in. [mm] w.e.
* SI(N,2)* : Snow depth above which there is 100 percent cover
* on snow covered pervious areas, in. [mm] w.e.
* WEPLOW(N) : Redistribution (plowing) depth on normally bare
* impervious area, in. [mm] w.e. Snow above this
* depth redistributed according to fractions below.
*
* Note: Redistribution (plowing) fractions (see Figure 4-25). Snow above
* WEPLOW in. [mm] w.e. on normally bare impervious area will be
* transferred to area(s) indicated below. The five fractions should
* sum to 1.0.
*
* SFRAC(N,1) : Fraction transferred to snow covered impervious area.
* SFRAC(N,2) : Fraction transferred to snow covered pervious area.
* SFRAC(N,3) : Fraction transferred to snow covered pervious area
* in last catchment.
* SFRAC(N,4) : Fraction transferred out of watershed.
* SFRAC(N,5) : Fraction converted to immediate melt on
* normally bare impervious area.
*============================================================================
* IF KWALTY = 0 on data group B1 skip to data group M1.
*###########################################################################
* Optional input of multiple land uses per subcatchment.
*
* IMUL is a variable to trigger multiple land uses
* per subcatchment. Any value > 0 will cause the
* model to use JLAND land uses per subcatchment.
* CAUTION: IMUL > 0 requires input on L2 lines
* for each subcatchment, even if JLAND = 1.
* ALSO, if IMUL > 0 and JLAND > 1, read JLAND J3 lines
* for each constituent.
* IMUL
JJ 1
*############################################################################
* Enter General Quality Control on data group J1.
*============================================================================
* J1 Line :
* NQS : Number of quality constituents. Maximum is controlled
* by parameter statement but should generally be limited to
* 20. NQS must be one less than maximum if erosion is simulated
* (IROS = 1).
* JLAND : Number of land uses (Maximum controlled by NLU paramter
* in TAPES.INC).
* IROS : Erosion simulation parameter
* = 0, Erosion not simulated.
* = 1, Erosion of suspended solids simulated using
* the Universal Soil Loss Equation. Parameters input
* in Group K1. Output will be last quality constituent
* (i.e., constituent NQS+1).
* IROSAD : Option to add erosion constituent to constituent
* number IROSAD. E.g., if IROSAD = 3, erosion will
* be added to constituent 3 (perhaps suspended solids).
* No addition if IROSAD = 0. N.R. if IROS = 0.
* DRYDAY : Number of dry days prior to start of storm.
* CBVOL : Average individual catchbasin storage volume, ft3 [m3].
* DRYBSN : Dry days required to recharge catchbasin concentrations
* to initial values (CBFACT on group J3). Must be > 0.
* RAINIT : For erosion, highest average 30-minute rainfall
* intensity during the year (continuous SWMM) or during
* the storm (single event), in./hr [mm/hr].
* N.R. if IROS = 0.
*
* The next three parameters are for modeling street sweeping.
* KLNBGN and KLNEND are only used if the simulation is greater
* than one month.
*
* REFFDD : Street sweeping efficiency (removal)
* fraction) for “dust and dirt.”
* KLNBGN : Day of year on which street sweeping
* begins (e.g. March 1 = 60).
* KLNEND : Day of year on which street sweeping
* stops (e.g. Nov. 30 = 334)
*============================================================================
* NQS JLAND IROS IROSAD DRYDRY CBVOL DRYBSN RAINIT REFFDD KLNBGN KLNEND
J1 4 2 1 0 5.00 2.0 1.00 0.30 0.50 0 0
*============================================================================
* Enter JLAND (from data group J1) Land Use data lines. One line for
* each land use. Land use 1 will be that of first group, land use 2 will be
* that of the second group etc.
*
* Note: Variables with asterisks can be modified
* using the Default/Ratio option.
*============================================================================
* J2 Line :
* LNAME(J) : Name of Land use.
* METHOD(J) : Buildup equation type for ‘dust and dirt'(see text).
* = -2, New default values,
* = -1, New ratios,
* = 0, Power-linear,
* = 1, Exponential,
* = 2, Michaelis – Menten.
* JACGUT(J) : Functional dependence of buildup parameters.
* = 0, Function of subcatchment gutter length,
* = 1, Function of subcatchment area,
* = 2, Constant.
*
* Following are up to three buildup parameters. (See Table 4-16).
*
* DDLIM(J)* : Limiting buildup quantity.
* DDPOW(J)* : Power or exponent.
* DDFACT(J)* : Coefficient.
*
* Following are three street sweeping parameters.
*
* CLFREQ(J)* : Cleaning interval, days.
* AVSWP(J)* : Availability factor, fraction
* DSLCL(J)* : Days since last cleaning, DSLCL <= CLFREQ
*============================================================================
* LNAME METHOD JACGUT DDLIM DDPOW DDFACT CLFREQ AVSWP DSLCL
J2 ‘SINGLE’ 0 0 1.E04 1.0 10.0 30.0 0.80 15.0
J2 ‘MULTPL’ 0 1 5.E04 1.5 50.0 7.0 0.80 5.0
*============================================================================
* Enter data for quality constituent(s) on data group J3. Repeat for
* each constituent, total of NQS groups. Constituent 1 will be that of the
* first line constituent 2 that of the second line, etc.
*
* If IMUL > 0 (line JJ) and JLAND > 1, then read JLAND J3 lines for
* each constituent. I.e., read JLAND lines for constituent 1, followed by
* JLAND lines for constituent 2, etc. In this case, each constituent
* can have different J3 parameters for each land use. These different
* parameters will be used on each land use fraction for each subcatchment,
* as defined in data group L2.
*
* Note: Variables with asterisks can be modified
* using the Default/Ratio option.
*============================================================================
* J3 Line :
* PNAME(K) : Constituent name.
* PUNIT(K) : Constituent units.
* NDIM(K) : Type of units.
* = 0, mg/l
* = 1, “Other” per liter, e.g., MPN/l or ug/l
* = 2, Other concentration units, e.g., pH, JTU
* KALC(K) : Type of buildup calculation.
* = 0, Buildup is fraction of “dust and dirt”
* for each land use.
* = 1, Power-linear constituent buildup
* = 2, Exponential constituent buildup
* = 3, Michaelis-Menten constituent buildup
* = 4, No buildup required (with KWASH = 1)
* KWASH(K) : Type of washoff calculation
* = 0, Power-exponential
* = 1, Rating curve, no upper limit (see note, below)
* = 2, Rating curve, upper limit by buildup equation
* KACGUT(K) : Functional dependence of buildup
* parameters. N.R. for KALC = 0 or 4.
* = 0, Function of subcatchment gutter length
* = 1, Function of subcatchment area
* = 2, Constant
* LINKUP(K) : Linkage to snowmelt. N.R. if ISNOW = 0 or KALC = 4.
* = 0, No linkage to snow parameters
* = 1, Constituent buildup during dry weather only when
* snow is present on impervious surface of subcatchment.
*
* Following are up to five buildup parameters
* (see text and Tables 4-17, 4-18).
*
* QFACT(1,K)* : First buildup parameter, e.g., limit.
* QFACT(2,K)* : Second buildup parameter, e.g., power or exponent.
* QFACT(3,K)* : Third buildup parameter, e.g. coefficient.
* QFACT(4,K)* : Fourth buildup parameter, N.R. if KALC > 0
* or JLAND < 4.
* QFACT(5,K)* : Fifth buildup parameter, N.R., if KALC > 0
* or JLAND < 5.
*
* Following are two washoff or rating curve parameters.
*
* WASHPO(K)* : Power (exponent) for runoff rate.
* RCOEF(K)* : Coefficient.
*
* CBFACT(K)* : Initial catchbasin concentration.
* (units according to NDIM).
* CONCRN(K)* : Concentration in precipitation.
* (units according to NDIM).
* REFF(K)* : Street sweeping efficiency (removal fraction)
* for this constituent.
*
*============================================================================
* For rating curve, equation is Load (mg/sec) = RCOEF*FLOW^WASHPO
* where FLOW is in cfs (METRIC=0) or cms (METRIC=1). If WASHPO = 1,
* this equation can be used to get constant concentration = EMC,
* and RCOEF must include conversion coefficient of 28.316 L/ft3 or 1000 L/m3.
* Then RCOEF = EMC*conversion. See example for TN, below.
* Another way to get a constant concentration is to set rainfall
* concentration to desired EMC and zero-out buildup-washoff parameters.
* Caution for constant concentration: dilution can result from inflows of
* I/I and groundwater (but both may be set to non-zero concentrations) and
* from initial water stored in channel/pipes. Cannot set non-zero
* concentrations for latter.
*============================================================================
* PNAME PUNIT NDIM KALC KWASH KACGUT LINKUP QFACT1 QFACT2 QFACT3 QFACT4 QFACT5 WASHPO RCOEF CBFACT CONCRN REFF
* Land use 1:
J3 ‘TOT.SOL’ ‘MG/L’ 0 2 0 0 0 900.0 2.0 0.0 0.0 0.0 2.0 1.5 100.0 2.0 0.7
* Land use 2:
J3 ‘TOT.SOL’ ‘MG/L’ 0 0 0 1 0 1000. 200. 0.0 0.0 0.0 2.0 1.5 150.0 2.0 0.7
* Land use 1:
J3 ‘ BOD5 ‘ ‘MG/L’ 0 1 0 0 0 60.0 1.5 0.3 0.0 0.0 2.0 1.2 20.0 0.1 0.5
* Land use 2:
J3 ‘ BOD5 ‘ ‘MG/L’ 0 0 0 1 0 200. 70. 0.0 0.0 0.0 2.0 1.2 30.0 0.1 0.5
*Simulate Total-N by rating curve. Want constant concentration = 25 mg/L.
*Use RCOEF = 25 mg/L * 28.316 L/ft3 = 707.9, and WASHPO = 1.0
* Land use 1:
J3 ‘ TOT-N ‘ ‘MG/L’ 0 4 1 0 0 0.0 0.0 0.0 0.0 0.0 1.0 707.9 0.0 0.0 0.0
* Land use 2:
J3 ‘ TOT-N ‘ ‘MG/L’ 0 4 1 0 0 0.0 0.0 0.0 0.0 0.0 1.0 707.9 0.0 0.0 0.0
* Use rating curve data from User’s Manual Fig. 4-37(e).
* Slope ~log(11/2.7)/log(1/0.1) = 0.61 = WASHPO
* At flow = 0.1 cfs, load ~ 2.7 mg/L ==>, RCOEF = 2.7/[0.1^(1/0.61)] = 118
* Land use 1:
J3 ‘NO2+NO3’ ‘MG/L’ 0 4 1 0 0 0.0 0.0 0.0 0.0 0.0 0.61 118. 0.0 0.0 0.0
* Land use 2:
J3 ‘NO2+NO3’ ‘MG/L’ 0 4 1 0 0 0.0 0.0 0.0 0.0 0.0 0.61 118. 0.0 0.0 0.0
*============================================================================
* Enter data for fractional contributions from other constituents
* on data group J4. Repeat until all desired fractions are entered.
*============================================================================
* J4 Line :
* KTO : Number (from order in Group J3) of constituent to
* which fraction will be added.
* KFROM : Number of constituent from which fraction is computed.
* F1(KTO,KFROM) : Fraction of constituent KFROM to be added
* to constituent KTO.
* [Note, these fractions will be applied for all
* land use segments if multiple land use option
* is used.]
*============================================================================
* KTO KFROM F1
J4 2 1 0.02
*============================================================================
* In data group J5, enter a constant groundwater concentration for every
* water quality constituent. Same units as NDIM in data group J3.
*============================================================================
* TS BOD5 TN NO2+NO3
J5 5.0 0.1 25.0 0.0
*Note, constant concentration of 25 mg/L will be maintained here, but
*concentration of NO2+NO3 will be diluted.
*============================================================================
* In data group J6, enter a constant infiltration/inflow concentration for
* every water quality constituent. Same units as NDIM in data group J3.
* Required only if I/I option used (lines F3, F4, H5). If omitted, I/I
* is assumed to have zero concentrations.
*============================================================================
* TS BOD5 TN NO2+NO3
J6 50.0 0.05 25.0 0.0
*Note, constant concentration of 25 mg/L will be maintained here, but
*concentration of NO2+NO3 will be diluted.
*============================================================================
* Enter Erosion Data on data group K1.
*
* If IROS = 0 on data group J1, skip to group L1.
*
* Note: Repeat group K1 ONLY for each subcatchment that is subject
* to erosion computations. The order of lines is arbitrary,
* but a match must be found of subcatchment number/name with a
* value of NAMEW used in group H1.
*
* Note: Variables with asterisks can be modified
* using the Default/Ratio option.
*============================================================================
* K1 Line :
* N=NAMEW : Subcatchment number or name matched with H1 line.
* ERODAR* : Area of subcatchment subject to erosion, acres [ha].
* ERLEN* : Flow distance in feet [meters] from point of
* origin of overland flow over erodible area to
* point at which runoff enters channel/pipe or inlet.
* SOILF* : Soil factor ‘K’.
* CROPMF* : Cropping management factor ‘C’.
* CONTPF* : Control practice factor ‘P’.
*============================================================================
* EROSION DATA
* NAMEW ERODAR ERLEN SOILF CROPMF CONTPF
K1 100 30.0 300.0 0.43 1.0 1.0
K1 200 4.0 200.0 0.43 1.0 1.0
K1 300 20.0 300.0 0.33 1.0 1.0
*============================================================================
* Enter Subcatchment Surface Quality data on data group L1.
*
* One line for each subcatchment is required. The order
* is arbitrary, but a match must be found for each subcatchment
* number (NAMEW) used earlier in group H1.
*
* Note: Variables with asterisks can be modified
* using the Default/Ratio option.
*============================================================================
* L1 Line :
* N=NAMEW : Subcatchment number or name.
* KL : Land use classification. 0 < KL < 5. Numbers
* correspond to input sequence of Group J2.
* Note: default value for KL = 1. If L2 lines are
* being used, still must enter a “representative”
* value for KL. Fraction data on L2 line will
* over-ride value of KL.
* BASINS(N)* : Number of catchbasins in subcatchment.
* GQLEN(N)* : Total curb length within subcatchment hundreds
* of feet [km]. May not be required depending on
* method used to calculate constituent loadings
* (Groups J2 and J3).
*
* The following initial constituent loading values may be input as an
* alternative to computation of loadings via methods specified in groups
* J2 and J3 (for initial conditions only). For any non-zero values
* read in, initial constituent loadings will be calculated simply by
* multiplication of the value by the subcatchment area (or fractional area
* if IMUL > 0 on line JJ). (I.e., if a loading value is entered on
* line L1, it will be apportioned over land uses with non-zero fractions.)
* “Load” has units depending on value of NDIM (Group J3),
* according to the following table:
*
* NDIM LOAD
* 0 pounds [kg]
* 1 10^6 x quantity, e.g. 10^6 MPN
* 2 10^6 x quantity x ft3,
* e.g. 10^6 pH-ft3.
*
* PSHED(1,N) : Initial loading, first constituent,
* load/acre [load/ha].
* . .
* PSHED(10,N) : Initial loading, tenth constituent,
* load/acre [load/ha].
*############################################################################
* Note, line L1 below illustrates use of slash to fill in any needed
* remaining zeros automatically.
*############################################################################
* NAMEW KL BA GQ PSHED(1) PSHED(2)
L1 100 1 12.0 20.0 /
*============================================================================
* L2 line :
*
* If required, this line immediately follows each L1 line.
* Enter land use fractions for subcatchment on previous L1 line.
* Not required if IMUL = 0 on line JJ (or line JJ omitted).
*
* If IMUL > 0, an L2 line must follow each L1 line, even if JLAND = 1.
* The fractions on line L2 must total to exactly 1.0 (tolerance = 0.001).
*
* These fractions (PLAND) are used to multiply the subcatchment area
* or curb length for all buildup calculations, depending on
* parameters JACGUT (line J2) and/or KACGUT (line J3).
* Similarly, if initial loads are input on line L1, the
* quantity/area values will be multiplied by the fraction area
* for each land use. Thus, the quantity/area values in line L1 cannot
* differentiate between different land uses.
* Throughout the simulation, separate buildup and washoff parameters
* will be used for each land use fraction, as input in the multiple
* J3 lines.
* If catchbasin quality information is entered, the total catchbasin
* load for a subcatchment is the sum over the number of land uses of:
* CBFACT*BASINS*CBVOL*PLAND.
*
* Enter JLAND fractions. Value 1 corresponds to land use 1, etc.
*
* PLAND(1,N) : Fraction of subcatchment N consisting of land use 1.
* PLAND(2,N) : Fraction of subcatchment N consisting of land use 2.
* Etc. Read JLAND values of PLAND.
*
*============================================================================
* Single Multiple
* PLAND1 PLAND2
L2 0.7 0.3
*============================================================================
* Additional L1/L2 lines. Must enter L1/L2 pairs for each subcatchment.
L1 200 2 20.0 30.0 /
L2 0.2 0.8
L1 300 1 6.0 10.0 /
L2 1.0 0.0
*============================================================================
* Enter data for Channel/Inlet Print Control on data group M1.
*============================================================================
* M1 Line :
* NPRNT : Total number of channels/pipes/inlets for which
* non-zero flows (and concentrations) are
* to be printed (maximum = NG).
* NPRNT < 0 means skip the M2 and M3 lines
* NPRNT = 12345 means print every inlet
* INTERV : Print Control.
* = 0, Print statistical summary only.
* = 1, Print every time step.
* = K, Print every K time steps.
*============================================================================
* NPRNT INTERV
M1 5 6
*============================================================================
* IF NPRNT = 0 on line M1 SKIP groups M2 and M3.
*
* Enter Print Period information on data group M2.
*============================================================================
* M2 Line :
* NDET : Number of detailed printout periods.
* (Maximum of 10 periods.)
*
* Note: If NDET = 1 and STARTP(1) = 0 and STOPPR(1) = 0 then the
* total simulation period will be printed as a default.
*
* STARTP(1) : First starting printout date, year, month,
* day, e.g., October 2, 1949 = 19491002.
* If year is entered as 2 digits, program assumes 1900.
* STOPPR(1) : First stopping printout date.
* . .
* STARTP(NDET): Last starting date.
* STOPPR(NDET): Last stopping date.
*============================================================================
* NDET STARTP(1) STOPPR(1)
M2 1 0 0
*============================================================================
* Enter channel/inlet printout locations on data group M3.
*============================================================================
* M3 Line :
* IPRNT(1) : First channel/inlet numbers or name for which flows
* and concentrations are to be printed.
* . .
* IPRNT(NPRNT): First channel/inlet numbers or name for which flows
* and concentrations are to be printed.
*
* Note: INflows to channel/pipes are printed as the default option. To
* print the OUTflow from a channel/pipe, give the ID as a negative
* number. A channel/pipe may be listed with both a positive and negative
* number. The negative option is not available for alphanumeric labeling.
*============================================================================
* IPRNT(1) … IPRNT(NRPNT)
M3 1 2 3 20 30
*============================================================================
* Enter channel/pipe depth output locations on data group M4.
* ****NOTE: M4 line is optional and may be omitted.****
*============================================================================
* M4 Line :
* MDEEP : Number of depth locations for printout (max = NG).
* KDEEP(1) : First conduit selected.
* . .
* KDEEP(MDEEP): Last conduit selected.
*============================================================================
* Select two conduits for depth printout
* # of conduits …conduits
* MDEEP KDEEP(1)..KDEEP(2)
M4 2 1 2
*============================================================================
* End your input data set with a $ENDPROGRAM.
$ENDPROGRAM

 

Categories: SWMM4

Tagged as:

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s