|
Statistics
|   |
|
Statistics
| |
$STATISTICS Call the STATISTICS Block with a '$' in first column.
*============================================================================
* The 'A1' line defines the starting and ending time and date.
* The meaning of the values listed on each line are as follows:
*============================================================================
* A1 Line :
* ISTART : Starting date, 8 digits, year/mo/day (19970925)
* If year is 2 digits, program assumes is 1900.
* TSTART : Starting time, decimal hours
* IEND : Ending date, 6 digits, year/mo/day (19970925)
* If year is 2 digits, program assumes is 1900.
* TEND : Ending time, decimal hours
* INLOG : = 0, calculate only moments of untransformed data.
* = 1, also calculate moments of natural logs of storm
* event data.
* JCUBE : = 0, storm event "volumes" are given in inches [mm].
* = 1, storm event volumes are listed in cubic feet
* [cubic meters].
*
* The following parameter is optional and need not be entered.
* JNEG : = 0, Program processes positive flow events as in
* the original program.
* = 1, Program processes negative flow events. Flow
* less than BASE indicate that you are in an event.
* Events with absolute volumes less than EBASE are
* not considered in the statistical analysis.
*============================================================================
* ISTART TSTART IEND TEND INLOG JCUBE
A1 600101 0.0 601231 0.0 0 0
*============================================================================
* The 'B1' line defines control information for the Statistics Block.
*============================================================================
* B1 Line :
* MIT : Minimum interevent time, decimal hours.
* BASE : Cutoff flow (baseflow), cfs [m3/s]. Flows less
* than BASE are not used. This value is not
* subtracted from flow rate but only serves as a
* threshold for analysis.
* EBASE : Cutoff storm event volume, inches [millimeters].
* Storm events below EBASE are not included in
* the statistical analysis.
* Note: If JCUBE = 1 and flows from interface file
* are being analyzed (LOCRQ > 0), then EBASE has
* units of volume, ft^3 [m^3].
* LOCRQ : ID of flow location requested. Enclose name
* in quotes if using $ANUM option. Enter a -1 or
* '-1' to process flows at all locations contained
* in the interface transfer file.
* LOCRN : Alphanumeric ID of rainfall gage for analysis.
* Enclose in single quotes, or use ' ' to indicate
* no rainfall analysis. [Note peculiarity of Lahey
* Fortran: alphanumeric input apparently still read
* correctly even if values are not in quotes! Thus,
* should be able to use old Statistics Block input
* files with integer LOCRN. But better idea to
* change to value enclosed in quotes.]
* Note: LOCRQ or LOCRN must be zero. Both cannot
* be non-zero.
* NPR : Number of pollutants requested (maximum is limited
* by MQUAL parameter).
* NPOINT : Number of events printed in tables.
* Print only the top NPOINT flows, events etc.
* IF NPOINT is 0 then all events are printed.
* METRIC : Requests type of units for output
* = 0, U.S. customary,
* = 1, Metric.
* LRET : Units of return period,
* = 0, Return period in years,
* = 1, Return period in months.
* A : Plotting postion parameter (see text for
* explanation). Briefly, the return period, T, is
* calculated as:
* T = (n+1-2A)/(m-A) , where
* n = total number of years or months, and
* m = rank. E.g., use A = 0 for Weibull plotting
* position and A = 0.4 for Cunnane's recommended
* plotting position.
*============================================================================
* MIT BASE EBASE LOCRQ LOCRN NPR NPOINT METRIC LRET PLOTTING POS(A)
*============================================================================
B1 6 0.0 0.0 1 ' ' 5 50 0 1 0.4
*B1 6 0.0 0.0 1 0 5 50 0 1 0.4
*============================================================================
* The 'B2' line defines print control information.
*============================================================================
* B2 Line :
* KSEQ : Request to print sequential series
* of flow events? No = 0, Yes = 1.
* KTERM : Code for terminating program if number of
* events exceeds allowable memory space (5000 events).
* = 0, Do not terminate (perform analyses
* on those events already identified).
* = 1, Terminate program (no event analysis performed).
* KTSEQS : Code for printing sequential series if the
* number of events exceed limit and KTERM = 1.
* = 0, Do not print sequential series.
* = 1, Print sequential series of those
* events already identified.
*============================================================================
* KSEQ KTERM KTSEQS
B2 1 0 1
*============================================================================
* The 'B3' line defines the interface file location
* of water quality information.
* Required only if NPR > 0 on Data Group B1.
* If NPR = 0 then skip to Data Group C1.
*
*============================================================================
* B3 Line :
* IPOLRQ(1) : First pollutant requested, identified
* by position on interface file.
* . .
* IPOLRQ(NPR): Last pollutant requested, identified by
* position on interface file.
*============================================================================
* IPOLRQ
B3 1 2 3 4 5
*============================================================================
* Enter the C1 line only if LOCRQ on Data Group B1 is <> 0.
*============================================================================
* The C1 line describes the Statistics Block Options for Flow.
*
* This data group controls the printing or plotting of information on
* magnitude, return period and frequency for each of the five flow parameters.
* In all cases, No = 0, Yes = 1. The control information is entered as a four
* digit integer number with each column controlling a different table or graph.
*
* First column - table of magnitude, return period
* and frequency.
* Second column - print graph of magnitude versus
* return period.
* Third column - print graph of magnitude versus
* frequency.
* Fourth column - print moments.
*
* For example, enter 1111 to print/plot all tables/graphs, 1000 to
* print the table only, 1100 to print the table and the graph of
* magnitude versus return period only, and 0000 to bypass printing/plotting
* of the flow parameter.
*============================================================================
* C1 Line :
* ISFLOW(1,1) : Request for total flow?
* ISFLOW(1,2) : Request for average flow?
* ISFLOW(1,3) : Request for peak flow?
* ISFLOW(1,4) : Request for event duration?
* ISFLOW(1,5) : Request for interevent duration?
*============================================================================
* ISFLOW(1,1) ISFLOW(1,2) ISFLOW(1,3) ISFLOW(1,4) ISFLOW(1,5)
C1 1111 1001 1001 1001 1001
*============================================================================
* Data Group D1 is only required if NPR > 0 on Data Group B1.
*============================================================================
* The D1 line describes the Statistics Block Options for Water Quality.
*
* If NPR > 0 use one D1 Data Group for each pollutant requested, up to
* ten sets of D1 lines, in the order defined by Group B2. The first
* index of ISPOLL(K,I,J) identifies the pollutant. Follow the
* instructions for Data Group C1 in entering data for Data Group D1.
*============================================================================
* D1 Line :
* ISPOLL(1,1,1) : Request for total load?
* ISPOLL(1,1,2) : Request for average load?
* ISPOLL(1,1,3) : Request for peak load?
* ISPOLL(1,1,4) : Request for flow weighted average concentration?
* ISPOLL(1,1,5) : Request for peak concentration?
*============================================================================
* ISPOLL(1,1,1) ISPOLL(1,1,2) ISPOLL(1,1,3) ISPOLL(1,1,4) ISPOLL(1,1,5)
D1 1111 1001 1001 1001 1001
D1 1001 1001 1001 1001 1001
D1 1001 1001 1001 1001 1001
D1 1001 1001 1001 1001 1001
D1 1001 1001 1001 1001 1001
*============================================================================
* Enter the E1 line only if LOCRN on Data Group B1 is > 0.
*============================================================================
* The E1 line describes the Statistics Block Options for Rainfall Analysis.
*
* Follow the instructions for Data Group C1 in entering data
* for Data Group E1.
*============================================================================
* E1 Line :
* ISFLOW(2,1) : Request for total volume?
* ISFLOW(2,2) : Request for average intensity?
* ISFLOW(2,3) : Request for peak intensity?
* ISFLOW(2,4) : Request for event duration?
* ISFLOW(2,5) : Request for interevent duration?
*============================================================================
* ISFLOW(2,1) ISFLOW(2,2) ISFLOW(2,3) ISFLOW(2,4) ISFLOW(2,5)
* E1 1001 1001 1001 1001 1001
* E1 1001 1001 1001 1001 1001
* E1 1001 1001 1001 1001 1001
* E1 1001 1001 1001 1001 1001
* E1 1001 1001 1001 1001 1001
*============================================================================
* End your input data set with a $ENDPROGRAM.
$ENDPROGRAM