ADA204562
ADA204562
0 ti
Ltl
US Army Corps
0 of Engineers
The Hydrologic
Engineering Center
DTIC
IIIELECTE
FEB 2
S .c--
HEC-I
Flood Hydrograph Package
Users Manual
rImamTIKUATEN=T
*n
S "i n Pub*i ...
.A . Vk~mb~udm UnlIWitd
September 1981
t.89 2 24
30272-101
REPORT DOCUMENTATION 1. REPORT NO. 2. 3. Recipient's Acceson No.
PAGE DOD/SW/DK-89/004a
4. Title and Stbtltle 5. Report Date 37
HEC-1, Flood Hydrograph Package Users Manual Sep 81 (Rev.
6.
Department of Army NA
The Hydrologic Engineering Center 11. Contract(C) or Grant(G) No.
Corps of Engineers ' r (c) N/A
609 Second Street (G)
Davis. Californiai 95616
12. Sponsoring Orgization Name and Address f13. Type of Report & Period Covered
NotesRelated Technical Papers: #82 ew HEC-1 klood Hydrograph Package (5-81); #95
15. SuPplementaryt
Infiltration and Soil Moisture Redistribut n in HEC-1 \(1-84); #106 Flood-Runoff
Forecasting with HEC1F (5-85); 1116 The H 's Activities in Watershed Modeling; HEC-l
Flood Hydrzaph Package pom~uter Implemntion Guide (5-85). Copies available from HEC.
16. Abstract (L it: 200 words) Four- D skette e
,'! HEC- is ae watershed -designed to simulate the surface runoff response
of a river basin to precipitation by representing the basias an interconnected system of
hydrologic and hydraulic components. Each component models an aspect of the
precipitation-runoff pfoCess within a portion of the basin. A component may represent a
surface runoff entity, a stream channel or a reservoir, defined by parameters specifying
its characteristics and the mathematical relations describing the physical processes. The
result of the modeling process is computation of streamflow hydrographs at desired
locations in the river basin. All ordinary flood hydrograph computations associated with
a single recorded or hypothetical storm can be accomplished with this package.
Capabilities include rainfall-snowfall-snowmelt determinations, computations of basin
precipitation, unit hydrographs. kinematic wave transforms, and hydrograph; routing by
reservoir, storage-lag, multiple storage, straddle-stagger, Tatum, Irlurskingum, and
kinematic wave methods; and complete stream system hydrograph combining and routing.
Best-fit unit hydrograph, loss-rate, snowmelt, base freezing temperatures and routing
coefficients can be derived automatically. Automatic printer plot routines are also
provided. HEC-l may also be used to simulate flow over and through breached dams.
Expected annual flood damage can be also be computed for any location in the river basin.
b. Identlflers/Open-Ended Terms
c. COSATI Field/Group
18. Availability Statement 19. Security Class (This Report) 21. No. of Pages
Unclassified 314
Release Unlimited 20. Security Class (This Page) 22. Price
HEC-1
Users Manual
Acceslon For
NTIS CRA&I
DTtC TAB E3
Unannoariced 0
Justification
(916) 551-1748
(FTS) 460-1748
TABLE OF CONTENTS
Section Page
List of Figures vi
Foreword ix
Introduction
2 Model Components
3 Rainfall-Runoff Simulation
3.1 Precipitation 8
3.1.1 Precipitation Hyetograph 8
3.1.2 Historical Storms 9
3.1.3 Synthetic Storms 10
3.1.4 Snowfall and Snowmelt 13
3.2 Interception/Infiltration 15
3.2.1 Initial and Uniform Loss Rate 16
3.2.2 Exponential Loss Rate 16
3.2.3 SCS Curve Number 18
3.2.4 Holtan Loss Rate 19
3.3 Unit Hydrograph 20
3.3.1 Basic Methodology 20
3.3.2 Synthetic Unit Hydrographs 20
3.4 Kinematic Wave 23
3.4.1 Basic Concepts 23
3.4.2 Solution Procedure 24
3.4.3 Element Application 27
3.5 Base Flow 29
3.6 Flood Routing 31
3.6.1 Channel Infiltration 31
3.6.2 Muskingum 32
3.6.3 Modified Puls 32
3.6.4 Working R & D 35
iii
TABLE OF CONTENTS (Contd)
Section Pale
3.7 Diversions 41
3.8 Pumping Plants 42
A Parameter Calibration
5 Multiplan-Multiflood Analysis 51
iv 4
TABLE OF CONTENTS (Cont'd)
Section page
11 Program Output
12 Example Problems
13 Computer Requirements
14 References 189
Index
v
LIST OF FIGURES
vi
LIST OF TABLES
No. Title
3.1 Distribution of Maximum 6-Hour SPS or PMP 11
3.2 Distribution of Maximum 1-Hour SPS or PMP 11
3.3 Partial-Duration Series to Equivalent-Annual Series
Conversion Factors 13
3.4 Adjustment of Point Rainfall Conversion Factors 14
3.5 Effective Roughness Parameters for Overland Flow 28
3.6 Typical Kinematic Wave Data 30
3.7 Spillway Rating Coefficients 38
3.8 Submergence Coefficients 40
4.1 Constraints on Unit Graph and Loss Rate Parameters 45
4.2 HEC-1 Unit Hydrograph and Loss Rate Optimization Output 46
4.3 HEC-1 Default Initial Estimates for Unit Hydrograph
and Loss Rate Parameters 47
10.1 HEC-1 Input Data Identification Scheme 73
10.2 Subdivisions of Simulation Data 76
10.3 Data Repetition Options 77
10.4 Precipitation Data Input Options 81
10.5 Hydrograph Derivation Input Options 81
10.6 Runoff and Routing Optimization Input Data Options 82
10.7 Channel and Reservoir Routing Methods Input Data Options 82
10.8 Spillway Routing, Dam Overtopping and Dam Failure
Input Data Options 83
10.9 Net Benefit Analysis Input Data Options 83
10.10 Flood Control Project Optimization Input Data Options 84
10.11 Hydrograph Transformation, Comparison and I/O 84
11.1 HEC-l Error Messages 87
12.la Red River Watershed: Rainfall and Observed Hydrograph Data 92
12.lb Subbasin Physical Parameters (Test 1) 92
12.lc Channel Storage Routing and Diversion Data 93
12.ld Example Problem 1: Input and Output 94
12.2a Subbasin Characteristics 107
12.2b Channel Data (Test 2) 107
12.2c Precipitation Data 108
12.2d Example Problems 2: Input and Output 109
12.3a Snowmelt Data 115
12.3b Example Problem 3: Input and Output 115
12.4 Example Problem 4: Input and Output 120
12.5 Example Problem 5: Input and Output 126
12.6a Depth Area Simulation Data 129
12.6b Example Problem 6: Input and Output 130
12.7a Reservoir Data 135
12.7b Example Problem 7: Input and Output 136
12.8a Dam Failure Analysis Results 140
12.8b Example Problem 8: Input and Output 140
12.9 Example Problem 9: Input and Output 155
12.10a Multiplan Analysis - Rockbed Watershed Flood Control Data 160
12.10b Example Problem 10: Input and Output 161
12.11a Flood Damage Reduction Analyses Economic Data 170
12.11b Example Problem 11: Input and Output 171
12.12a Flood Control System Optimization Data 177
12.12b Example Problem 12: Input and Output 178
13.1 I/O and Scratch Files 187
13.2 Computer Memory and Time Requirements 188
vii
FOREWORD
In 1973, the 1970 version of the program underwent a major revision. The
computational methods used by the program remained basically unchanged;
however, the input and output formats were almost completely restructured.
These changes were made in order to simplify input requirements and to make
the program output more meaningful and readable.
The present program again represents a major revision of the 1973 version
of the program. The program input and output formats have been completely
revised and the computational capabilities of the dam-break HEC-1DB), project
optimization (HEC-IGS) and kinematic wave (HEC-lKW) special versions of HEC-1
have been combined in the one program. The new program gives the powerful
analysis features available in all the previous programs, together with some
additional capabilities, in a single easy to use package.
Up-to-date information and copies of source code for the programs are
available from the Center. While the Government is not responsible for the
results obtained when using the programs, assistance in resolving malfunctions
in the programs will be furnished to the extent that time and funds are
available. It is desired that users notify the Center of inadequacies in, or
desirable modifications to, the program.
ix
Section 1
INTRODUCTION
1
Section 9: The optimal size of a flood control system can be estimated
using an optimization procedure provided by HEC-1. The option utilizes data
provided for the economic assessment option together with data on flood
control project costs to determine a system which maximizes net benefits with
or without a specified degree of protection level for the components.
Program Usage
Section 10: The data input conventions are discussed, emphasizing the
data card groups used for the various program options.
Section 11: Program output capabilities and error messages are explained.
Section 12: Test examples are displayed, including example input data and
computed output generated by the program.
Section 13: The computer hardware requirements are discussed, and computer
run times for the example problems are given. A programers supplement
provides detailed information about the operational characteristics of the
computer program.
Appendix A: The input description details the use of each data card and
input variable in the program. The input description is contained in under
seperate cover.
2
1.4 Computer Requirements
The HEC-1 program requires 371,000 octal words (130,000 decimal) of core
storage. Disk storage is needed for the 16 output and scratch files used by
the program. For further information on the program's computer requirements,
see Section 13 and the Programmers Supplement.
1.5 Acknowledgments
This manual was written by Messrs. David Goldman and Paul Ely. Mr. Ely
was also responsible for the design and implementation of the new computer
code. Mr. John Tracy modified the code for use on microcomputers. Messrs.
John Peters, Darryl Davis and Arthur Pabst made many excellent contributions
to the development of the modeling concepts and the documentation. The
development of this new version of HHC-1 was managed by Mr. Arlen D. Feldman,
Chief of the HEC Research Branch. Mr. Bill S. Eichert was the Director of the
HEC during this time. The word processing for this document was performed by
Ms. Cathy Lewis.
3
Section 2
MODEL COMPONENTS
The stream network model is the foundation capability of the HEC-i program.
All other program computation options build on this option's capability to
calculate flood hydrographs at desired locations in a river basin. Section
2.1 discusses the conceptual aspects of using the HEC-1 program to formulate a
stream network model from river basin data. Section 2.2 discusses the model
formulation as a step-by-step process, where the physical characteristics of
the river basin are systematically represented by an interconnected group of
HEC-1 model components. Sections 2.3 - 2.8 discuss the functions of each
component in representing individual characteristics of the river basin.
!~2 1- Ito
.3 2050
A COMPONENT
- - - CHANNEL ROUTINW
Figure 2.1 ETsmple River Basin Figure 2.2 Example River Basin Schematic
4
(2) Segmentation of the basin into a number of subbasins determines the
number and types of stream network components to be used in the model. Two
factors impact on the basin segmentation: the study purpose and the
hydrometeorological variability throughout thp basin. First, the study
purpose defines the areas of interest in the basin, and hence, the points
where subbasin boundaries should occur.
(4) The subbasins and their components are linked together to represent
the connectivity of the river basin. HEC-1 has available a number of methods
for combining or linking together outflow from different components. This
step finalizes the basin schematic.
The subbasin land surface runoff component, such as subbasins 10. 20, 30,
etc. in Fig. 2.1 or equivalently as element 10 in Fig. 2.2, is used to
represent the movement of water over the land surface and in stream channels.
The input to this component is a precipitation hyotograph. Precipitation
excess is computed by subtracting infiltration and detention losses based on a
soil water infiltration rate function. Note that the rainfall and infiltration
are assumed to be uniform over the subbasin. The resulting rainfall excesses
are then routed by the unit hydrograph or kinematic wave techniques to the
outlet of the subbasin producing a runoff hydrograph. The unit hydrograph
technique produces a runoff hydrograph at the most downstream point in the
subbasin. If that location for the runoff computation is not appropriate, it
may be necessary to further subdivide the subbasin or use the kinematic wave
method to distribute the local inflow.
5
2.3 River Routing Component
Consider the use of subbasin runoff components 10 and 20 and river routing
reach 1020 in Fig. 2.2 and the corresponding subbasins 10 and 20 in Fig. 2.1
The runoff from component 10 is calculated and routed to control point 20 via
routing reach 1020. The runoff hydrograph at analysis point 20 can be
calculated by methods employing either the unit hydrograph or kinematic wave
techniques. In the case that the unit hydrograph technique is employed,
runoff from component 10 is calculated and routed to contrl point 20 via
routing reach 1020. Runoff from subbasin 20 is calculated and combined with
the outflow hydrograph from reach 1020 at analysis point 20. Alternatively,
runoff from subbasins 10 and 20 can be combined before routing in the case
that the lateral inflows from subarea 20 are concentrated near the upstream
end of reach 1020. In the case, that the kinematic wave technique is
employed, the runoff from subbasin 20 is modeled as a uniformly distributed
lateral inflow to reach 1020. The runoff from subbasin 10 is routed in
combination with this lateral inflow via reach 1020 to analysis point 20.
6
2.6 Diversion Component
The pump component can be used to simulate action of pumping plants used
to lift runoff out of low lying ponding areas such as behind levees. Pump
operation data describes the number of pumps, their capacities, and "on" and
"off" elevations. Inflow to the pump station comes from the river channel.
7
Section 3
RAINFALL-RUNOFF SIMULATION
3.1 Precipitation
INTENSITY
(in /hr)
to tt ta to 14 to
TIME (w)
8
3.1.2 Historical Storms
n
E PRCPN(J) * WTN(J)
Jul
PRCPA =. ........ . . . ................ (3.1)
n
E WTN(J)
J-l
n
I PRCPN(J) WTN(J)
PRCPA = SNAP *... . .... ............... (3.2)
n
E ANAPN(J) * WTN(J)
J=l
where ANAPN is the station normal annual precipitation, and SNAP is the
subbasin-average normal annual precipitation. Use of this option may be
desirable in cases where precipitation measurements are known to be biased.
For example, data obtained from a gage located on the floor of a valley may
consistently underestimate subbasin average precipitaticn for higher
elevations. ANAPN may be used to adjust for this bias.
n
E PRCPR(I,J) * WTR(J)
PRCP(I) = Jul..... ................ (3.3)
n
I WTR(J)
J=1
where PRCP(I) is the basin-average precipitation for the Ith time interval,
PRCPR(I, J) is the recording station precipitation for the Ith time interval,
and WTR(J) is the relative weight for gage J.
9
The subbasin-average hyetograph is computed using the temporal pattern,
Synthetic storms are frequently used for planning and design studies.
Criteria for synthetic storms are generally based on a detailed analysis of
long term precipitation data for a region. There are three methods in HKC-l
for generating synthetic storm distributions:
(I) Standard Project Storm. The procedure for computing Standard Project
Storms, SPS, programmed in HKC-l is applicable to basins of area 10 to 1,000
square miles located east of 105" longitude. The SP5 is determined by
specifying an index precipitation, SPFE, a storm reduction coefficient, TRSPC,
and the area over which the storm occurs, TRSDA. SPFE and TRSPC are
determined by referring to manual SH-1110-2-1411 (Corps of Engineers, 1952).
A total storm depth is determined and distributed over a 96-hour duration
based on the following formulas which were derived from design charts in the
referenced manual.
R24HR(2) = 15.5
R24HR(4) = 6.0
where R24HRM() is the percent of the index precipitation occurring during the
Ith 24-hour period.
Each 24-hour period is divided into four 6-hour periods. The ratio of the
24-hour precipitation occurring during each 6-hour period is calculated as
361133)
- 13.42
11.0). 9 3
R6HR(3) - (SPrZ +13.2 ....... ................... (3.5)
'
R6HR(2) = 0.055 * (SPF1 - 6.0)0 ...... ................. .(3.6)
where R6HR(I) is the ratio of 24-hour precipitation occurring during the Ith
6-hour period and SPF3 ts the index precipitation in inches.
The precipitation for each time interval, except during the peak 6-hour
period, is computed as
10
The peak 6-hour precipitation of each day is distributed according to the
percentages in Table 3.1 If time intervals less than one hour ase used, the
peak 1-hour precipitation is distributed according to the percentages in Table
3.2. The time interval must divide evenly into one hour. When the time
interval is larger than shown in Tables 3.1 and 3.2, the percentage for the
peak time interval is the sum of the highest percentages; e.g. for a 2-hour
time interval, the values are (14+12)%, (38+15)%, and (11+10)%. The interval
with the largest percentage is preceded by the second largest and followed by
the third largest. The second largest percentage is preceded by the fourth
largest, the third largest percentage is followed by the fifth largest, etc.
TABLE 3.1
1 10 4
2 12 8
3 15 19
4 38 50
5 14 11
6 11 8
TABLE 3.2
5 3 3
10 4 7
15 5 12
20 6 18
25 9 27
30 17 44
35 25 69
40 11 80
45 8 88
50 5 93
55 4 97
60 3 100
11
(ii) Probable Maximum Precipitation. Current Probable Maximum Precipita-
tion. PMP. computation methods are not available in HEC-1. The PEP must be
determined according to the National Weather Service's Hydrometeorological
Reports Nos. 36, 43, 49, 51, 52, or 55, depending upon geographic location.
Computer program RMM52 (REC, 1984) is available to assist with PuP and Probable
Maximum Storm determination for the eastern United States. The PHP computed
from 1MR52 or any other method may be input to fEC-1 to calculate runoff.
112 - 16
R6HR(2) - 124...... ... ....................... .(3.8b)
R24
0.3008
TRSPC = 1 - TRSDA 0.17718......................(3.9)
12
Weather Service, 1961) data are used, the program will automatically make the
partial-to-annual series conversion using the factors in Table 3.3 (which is
table 2 of TP-40) if desired.
TABLE 3.3
Conversion
Return Period Frequency Factor C
Depths for 10-minute and 30-minute durations are interpolated from 5-, 15-,
and 60-minute depths using the following equations from NYDRO-35 (National
Weather Service, 1977):
13
TABLE 3.4
Point-to-Areal Rainfall Conversion Factors
0.5 .48
1 .35
3 .22
6 .17
24 .09
48 .068
96 .055
168 .049
240 .044
variations from 320F (OC) might be warranted considering both spatial and
temporal fluctuations of temperature within the zone.
where SNWKT is the melt in inches (m) per day in the elevation zone, TNPR is
the air temperature in OF or *C lapsed to the midpoint of the elevation zone,
FRZTP is the temperature in OF or *C at which snow melts, and COEF is the melt
coefficient in inches (mm) per degree-day (OF or *C).
SNW T = COEF (.09 + (.029 + .00504 WIND + .007 RAIN) (TMPR - FRZTP)) . (3.14)
14
SNWKT = COEF (.002 SOL (1-ALBDO) + (.0011 WIND + .0145) +
(TMPR - FRZTP) + .0039 WIND (DKWPT - FRZPT)) .. ....... .(3.15)
where SNWtT is the melt in inches per day in the elevation zone, TMPR is the
air temperature in OF lapsed at the rate TLAPS to midpoint of the elevation
zone, DRWPT is the dewpoint temperature in OF lapsed at a rate 0.2 TLAPS to
the midpoint of the elevation zone. A discussion of the decrease in dewpoint
temperature with higher elevations is found in (Miller, 1970). FRZTP is the
freezing temperature in 0F, COEF is the dimensionless coefficient to account
for variation from the general snowmelt equation referenced in EM 1110-2-1406.
RAIN is the rainfall in inches per day, SOL is the solar radiation in langleys
per day, ALBDO is the albedo of snow, .75/(D2), constrained above .4, D is
the days since last snowfall, and WIND is the wind speed in miles per hour, 50
feet above the snow.
3.2 Interception/Infiltration
The precipitation loss computations can be used with either the unit
hydrograph or kinematic wave model components. In the case of the unit
hydropraph component, the precipitation loss is considered to be a subbasin
average (uniformly distributed over an entire subbasin). On the other hand,
separate precipitation losses can be specified for each overland flow plane
(if two are used) in the kinematic wave component. The losses are assumed to
be uniformly distributed over each overland flow plane.
There are four methods that can be used to calculate the precipitation
loss. Using any one of the methods, an average precipitation loss is
determined for a computation interval and subtracted from the
rainfall/snowmelt hyetograph as shown in Fig. 3.2. The resulting
precipitation excess is used to compute an outflow hydrograph for a subbasin.
A percent imperviousness factor can be used with any of the loss rate methods
to guarantee 100 runoff from that portion of the basin.
15
A percent Impervious factor can be used with any of the loss rate methods;
it guarantees 100% runoff from that percent of the subbasin.
PRECIPIrATION EXCES
INTENSITY
(in/kr)
or AVERASE INFILTRATION
(bOVER PERIOD
tt to to t4 to to
TIME (1I)
An Initial loss, STRTL (units of depth), and a constant loss rate, CKSTL
(units of depth/hour), are specified for this method. All rainfall is lost
until the volume of initial loss is satisfied. After the initial loss is
satisfied, rainfall is lost at the constant rate, CWSTL.
16
where ALOSS is the potential loss rate in inches (mu) per hour during the time
interval, AK is the loss rate coefficient at the beginning of the time
interval, and DLTK is the incremental increase in the loss rate coefficient
during the first DLTKR inches (mr) of accumulated loss, CUNL. The accumulated
loss, CUML, is determined by summing the actual losses computed for each time
interval. Note that there is not a direct conversion between metric and
English units for coefficients of this method, consequently separate
calibrations to rainfall data are necessary to derive the coefficients for
both units of measure.
STRKR ---- 1
e tO /
It-
Arithmetic Scale
Figure 3.3 General HEC Loss Rate Function for Snow-free Ground
- •=• I i l i l | ||I
DLTKR is the amount of initial accumulated rain lose during which the loss
rate coefficient is increased. This parameter is considered to be a function
primarily of antecedent soil moisutre deficiency and is usually storm
dependent. STRKR is the starting value of loss coefficient on exponential
recession curve for rain losses (snow-free ground). The starting value is
considered a function of infiltration capacity and thus depends on such basin
characteristics as soil type, land use and vegetal cover.
17
function that veflects the influence of precipitation rate on basin-average
loss characteristics. It reflects the manner in which storms occur within an
area and may be considered a characteristic of a particular region. ERAIN
varies from 0.0 to 1.0.
AIS=ACRAN
AACRA -- IA + S................................. (3.17)
IA+S
In the case that the user does not wish to specify IA, a default value is
computed as
Since the SCS method gives total excess for a storm, the incremental
excess (the difference between rainfall and precipitation loss) for a time
period is computed as the difference between the accumulated excess at the end
of the current period and the accumulated excess at the end of the previous
period.
18
3.2.4 Holtan Loss Rate
Holtan et al. (1975) compute loss rate based on the infiltration capacity
given by the formula:
where f is the infiltration capacity in inches per hour, GIA is the product of
GI a "growth index" representing the relative maturity of the ground cover and
A the infiltration capacity in inches per hour (inch 1" 4 of available
storage), SA is the equivalent depth in inches of pore space in the surface
layer of the soil which is available for storage of infiltrated water, PC is
the constant rate of percolation of water through the soil profile below the
surface layer, and BRIP is an empirical exponent, typically taken equal to 1.4.
Holtan et al. (1975) have made estimates of the value of "A" for several
vegetation types. Their estimates were evaluated at plant maturity as the
percent of the ground surface occupied by plant stems or root crowns.
Estimates of PC can be based on the hydrologic soil group given in the SCS
Handbook (1972 and 1975). Musgrave (1955) has given the following values of
PC in inches per hour for the four hydrologic soil groups: A, 0.45 to 0.30;
B, 0.30 to 0.15; C, 0.15 to 0.05; D, 0.05 or less.
19
3.3 Unit Hydrozraph
I
Q(i) = Z U(j)*X(i-j+l) ............................ (3.25)
J=l
where Q(i) is the subbasin outflow at the end of computation interval I, U(j)
is the jth ordinate of the unit hydrograph, X(i) is the average rainfall
excess for computation interval I.
The parameters for the synthetic unit hydrograph can be determined from
gage data by employing the parameter optimization option described in Section
4. Otherwise, these parameters can be determined from regional studies or
from guidelines given in references for each synthetic technique. There are
three synthetic unit hydrograph methods available in the model.
(I) Clark Unit Hydrograph. The Clark method (1945) requires three
parameters to calculate a unit hydrograph: TC, the time of concentration for
the basin, R, a storage coefficient, and a time-area curve. A time-area curve
defines the cumulative area of the watershed contributing runoff to the
subbasin outlet as a function of time (expressed as a proportion of TC).
In the case that a time area curve is not supplied, the program utilizes a
dimensionless time area curve:
20
where Al is the cumulative area as a fraction of total subbasin area and T is
the fraction of time of concentration. The ordinates of the time-area curve
are converted to volume of runoff per second for unit excess and interpolated
to the given time interval. The resulting translation hydrograph is then
routed through a linear reservoir to simulate the storage effects of the
basin; and the resulting unit hydrograph for instantaneous excess is averaged
to produce the hydrograph for unit excess occurring in the given time interval.
(ii) Snyder Unit Hydrograph The Snyder method (1938) determines the unit
graph peak discharge, time to peak, and widths of the unit graph at 50 and 75%
of the peak discharge. The method does not produce the complete unit graph
required by HEC-1. Thus, HEC-1 uses the Clark method to affect a Snyder unit
graph. The initial Clark parameters are estimated from the given Snyder's
parameters, Tp and Cp. A unit hydrograph is computed using Clark's method and
Snyder parameters are computed from the resulting unit hydrograph by the
following equations:
CPTKP = QMAX * Tpeak - 0.5* At
C*A..................(3.32)
*
ALAG = 1.048 (Tpeak - 0.75* At) ..... ............... . (3.33)
where CPTKP is Snyder's Cp for computed unit hydrograph, QKAX is the maximum
ordinate of unit hydrograph, Tpeak is the time when QNAX occurs, in hours,
At is the duration of excess, in hours, A is the subbasin area in square
miles (sq km), C is a conversion factor, and ALAG is Snyder's standard Lag, Tp
for the computed unit hydrograph. Snyder's standard Lag i for a unit
hydrograph which has a duration of excess equal to Tp/5.5. The coefficient,
1.048, in equation results from converting the duration of excess to the given
time interval.
21
Clark's TC and R are adjusted to compensate for differences between values
of Tp and Cp calculated by equations 3.32 and 3.33 and the given values. A
new unit hydrograph is computed using these adjusted values. This procedure
continues through 20 iterations or until the differences between computed and
given values of Tp and Cp are less than one percent of the given values.
(iii) SCS Dimensionless Unit Hydrograph. Input data for the Soil
Conservation Service, SCS, dimensionless unit hydrograph method (1972) consists
of a single parameter, TLAG, which is equal to the lag (hrs) between the
center of mass of rainfall excess and the peak of the unit hydrograph. Peak
flow and time to peak are computed as:
t.0
i 4
.2
0i
0 1-
0 1 2 3 4
DIMENSIONLESS TIM, tpT'pM.E
22
The selection of the program computation interval, which is also the
duration of the unit hydrograph, is based on the relationship At = .2*TPKAK
(SCS, 1972, Chapters 15, 16). There is some latitude allowed in this
relationship; however, the duration of the unit graph should not exceed
At < .25*Tpeak. These relations are based on an empirical relationship,
TLAG = .6*Tc, and 1.7*TPEAK - At+Tc where Te is the time of concentration
of the watershed. Using these relationships, along with equation (3.34) it is
found that the duration should not be greater than At < .29*TL&G.
Sf = SO . . . . . . . . . . . . . . . . . . . . . . . .. .. ... . . . . . . . . . . .. (3.36)
where Sf is the friction slope and So is the channel bed slope. Thus flow
at any point in the channel can be computed from Manning's formula.
1/2 2/3
Q = 1.486 S R A ......... ....................... .(3.37)
n
where a and m are related to flow geometry and surface roughness. Figure 3.6
gives relations for a and m for channel shapes used in NEC-1. Note that flow
depths greater than the diameter of the circular channel shape are possible, which
only approximates the storage characteristics of a pipe or culvert.
Since the momentum equation has been reduced to a simple functional relation
between area and discharge, the movement of a flood wave is described solely by the
continuity equation
23
3.4.2 Solution Procedure
The governing equations for either overland flow
solved In the same manner. The method assumes that or channel routing are
inflows, whether It be
rainfall excess or lateral inflows, are constant
within a time step and
uniformly
0 EC~o
Flows from
overland flow
elements
L
INFLOW
oe Collector channel$
SUBSASIN OUTFLOW
24
distributed along the element. By combining equations 3.38 and 3.39, the
governing equation is obtained as:
CIRCULAR
1
.804 S1/2 D /6
(APPROXIMATE)
m = 5/4
TRIANGULAR
0,94 s1/2 Z >I/3
- -
m =4f3
SQUARE
T CL = _.7 S11 2
n
Y:w
- 4/ 3
-m
RECTANGULAR
11 2 2/3
1.49 s W-
S-m = 5/3
TRAPEZOIDAL
25
A is the only dependent variable In the equation; a and m are considered
constant. The standard form of the finite difference approximation to this
equation is developed as:
At 2
Al -l) - iq.... q ......... (3.41)
t
t
j 0
A - A
At
tj- 1 0
j indicates the current time step of the solution scheme. i-1, J-1 indicate,
respectively, positions and times removed a value Ax and At from the
current position of the solution scheme. The only unknown value in the
equation is the current value A(jj). All other values are known from
either a solution of the equation at a previous position 1-1 and time J-i, or
from a boundary condition. Solving for the unknown:
A(ij ) = q At + A(i,j-l)
26
Once A(j,j) is known, the flow can be computed as:
+ A(ilj)]. (3.46)
Q(i,j) - Q(i-1,j) q Ax - [a(iNlj) - .....
The space increment Ax and time step At are chosen by the program to
insure the scheme stability and convergence, and are based on experience in
the use of the scheme. Ax is computed by the program to fall between the
limits:
LREACH LREACH
so0 50 < Ax -< 2..................................
2 (3.48)
where LREACH is the length of a channel reach. At has a minimum value based
on variable array requirements which are defined as:
27
TABLE 3.5
Surface N
Two overland flow elements may be used for each subbasin. The total
discharge, Q, from each element is computed as
Q = q * AREA
Q~~ L......................................... (3.50)
where q is the discharge per unit width from each overland flow element
computed from equations (3.44) or (3.46), AREA is the area represented by each
element, and L is the overland flow length.
(ii) Channel Elements. Flow from the overland flow elements travels
to the subbasin outlet through one or two successive channel elements, Fig.
3.5. A channel is defined by length, slope, roughness, shape, width or
diameter, and side slope, Fig. 3.6. The last channel in a subbasin is called
the main channel, and any intermediate channels between the overland flow
elements and the main channel are called collector channels. Use of a
collector channel is optional.
Lateral inflow into a channel element from overland flow is the sum of the
total discharge computed by equation (3.50) for both elements divided by the
channel length. If the channel is a collector, the area used in equation
(3.50) is the area serviced by the collector. Lateral inflow, q, from a
collector channel is computed as
Q * AREA2 3.
AR.Al L........................(3.51)
where Q is the discharge from the collector, AREAl is a typical area served
by this collector, AREA2 is the area served by the channel receiving flow
28
from the collector, and L is the length of the receiving channel. If the
receiving channel is the main channel, ARZA2 is the subbasin area.
The model requires that at least one overland flow plane and one main
channel be used in kinematic wave applications. In the above example, fewer
elements might have been used depending on the level of detail required for
the hydrologic analysis.
29
TABLE 3.6
Typical Kinematic Wave Data
Average
Identifica- Overland Flow Slope Roughness Percentage of
tion Length (ft) (ft/ft) Coefficient Subbasin Area
Channel Data
DISCHARGE TRcuslN
DIRECT THaIIN0L
SURFACE RUNOFF
rRAWFALL/SNOWM.LT1
I TRANSFORMED TO I
STARTING.
S SUSBASIN RUNOFF J
I "
0
/O
F L OW
TIME
30
Q = Qo (RTIOR)-n a......... ........................ . (3.52)
where Qo is STRTQ or QRCSN, and nat is the time in hours since recession
was inftiated. QRCSN and RTIOR can be obtained by plotting the log of
observed flows versus time. The point at which the recession limb fits a
straight line defines QRCSN and the slope of the straight line is used to
define RTIOR. Alternatively, QRCSN can be specified as a ratio of the peak
flow. For example, the user can specify that the exponential recession is to
begin when the "falling limb" discharge drops to 0.1 of the calculated peak
discharge.
The rising limb of the streamflow hydrograph is adjusted for base flow by
adding the recessed starting flow to the computed direct runoff flows. The
falling limb is determined in the same manner until the computed flow is
determined to be less than QRCSN. At this point, the time at which the value
of QRCSN is reached is estimated from the computed hydrograph. From this time
on, the streamflow hydrograph is computed using the recession equation unless
the computed flow rises above the base flow recession. This is the case of a
double peaked streamflow hydrograph where a rising limb of the second peak is
computed by combining the starting flow recessed from the beginning of the
simulation and the direct runoff.
Storage routing methods in HKC-1 are those methods which require data that
define the storage characteristics of a routing reach or reservoir. These
methods are: modified Puls, working R and D, and level-pool reservoir routing.
There are also two routing methods in HEC-1 which are based on lagging
averaged hydrograph ordinates. These methods are not based on reservoir
storage characteristics, but have been used on several rivers with good
results.
31
A second methods computes channel loss during storage routing based on a
constant channel loss (cfs/acre) per unit area and the surface area of channel
flow. The surface area of channel flow is computed as:
where Q(I) is the routed outflow and QO(I) is the flow adjusted for the
constant channel loss rate PERCRT (cfs/acre).
3.6.2 Muskingum
At - 2
AMSKK +X
*
CB = 2 * AMSKK * (l-X) + At...................(3.59)
where QIN is the inflow to the routing reach in cfs (m3 /sec), QOUT is the
outflow from the routing reach in cfs (m3 /sec), AMSKK is the travel time
through the reach in hours, and X is the Muskingum weighting factor
(0 < X < .5). The routing procedure may be repeated for several subreaches
(designated as NSTPS) so the total travel time through the reach is AMSKK. To
insure the method's computational stability and the accuracy of computed
hydrograph, the routing reach should be chosen so that:
1 AMSKK 1
2(l-X) NSTPS*dt 1 2x.......................(3.60)
The modified Puls routing method (Chow, 1964) is a variation of the storage
routing method described by Henderson (1966). It is applicable to both channel
and reservoir routing. Caution must be used when applying this method to
channel routing. The degree of attenuation introduced in the routed flood wave
32
varies depending on the river reach lengths chosen, or alternatively, on the
number of routing steps specified for a single reach. The number of routing
steps (variable NSTPS) is a calibration parameter for the storage routing
methods; it can be varied to produce desired routed hydrographs. A storage
indication function is computed from given storage and outflow data.
STRII) - C * STOR(I)
at
____ ++ OTLI.........................
2 (3.61)
where STRI is the storage indication in cfs (m/sec), STOR is the storage
in the routing reach for a given outflow in acre-ft (1000 m), OUTFL is
the outflow from routing reach in cfs (m/sec), C is the conversion factor
from acre-ft/hr to cfs (1000 ms/hr to m3/sec), At is the time
interval in hours, and I is a subscript indicating corresponding values of
storage and outflow. Storage indication at the end of each time interval is
given by
where QIN is the average inflow in cfs (m/sec), and Q is the outflow in
cfs (ms/sec), and subscripts 1 and 2 indicate beginning and end of the
current time interval.
The outflow at the end of the time interval is interpolated from a table
of storage indication (STRI) versus outflow (OUTFL). Storage (STR) is then
computed from
SR= (STRI - 2 ) At
STR = C. . .......................... (3.63)
When stage data are given, stages are interpolated for computed storages.
(i) Given Storage versus Outflow Relationship. The modified Puls routing
may be accomplished by providing a storage versus outflow relationship as
direct imput to HEC-1. Such a relationship can be derived from water surface
profile studies or other hydraulic analyses of rivers or reservoirs.
(ii) Normal-Depth Storage and Outflow. Storage and outflow data for use
in modified Puls or working R&D (see next subsection) routing may be computed
from channel characteristics. The program uses an 8-point cross section which
is representative of the routing reach (Fig. 3.9). Outflows are computed for
normal depth using Manning's equation. Storage is cross-sectional area times
reach length. Storage and outflow values are computed for 20 evenly-spaced
stages beginning at the lowest point on the cross section to a specified
maximum stage. The cross section is extended vertically at each end to the
maximum stage.
As shown in Fig. 3.9, the input variables to the program are the hydraulic
and geometric data: ANL, ANCH, ANR, RLNTH, SEL, ELMAX, and (X,Y) coordinates.
ANL, ANCH, ANR are Manning's n values for left overbank, main channel, and
33
right overbank, respectively. RLNTH is routing reach length in feet (meters).
SEL is the energy gradient used for computing outflows. (X,Y) are coordinates
of an 8-point cross section.
Storage and outflow should not be calculated from normal depth when the
storage limits and conveyance limits are significantly different. Also, if
the cross section is "representative" for a reach that is not uniform, the
stages will not be applicable to any specific location. Generally, the stages
produced by the method are of limited value because downstream effects are not
taken into account.
RSTIVE
ROUTING REAC
Reach Leneth feet) e iRLeTH
Finery ra Line Slope ( Strgtt) aetl
pt. I Pt. 0I
X(I),Y(I} X(a), Y(s) )
ANCH .
X (6), Y(G)
Pt.4
X (4), Y(4)
x(s),V(s)
34
3.6.4 Working R and D
St 2
S2 - S2 _ Il + 12 - 01 + 02 ......................... (3.68)
At 2 2
35
where Y is the continuity error for the estimated elevation. The estimated
elevation is adjusted until Y is within ± 1 cfs (m/sec).
AV12 (A1 + +
WMTS
-• oseto bosel,
elevation given, even if the surface area is greater than zero at that point.
36
(ii) Orifice and Weir Flow. This option is often used in spillway
adequacy investigations of dam safety, see Example Problems, Sections 12.7 and
12.8.
If pumps or dam breaks are not being simulated, an outflow rating curve is
computed for 20 elevations which span the range of elevations given for
storage data. Storages are computed for those elevations. The routing is
then accomplished by the modified Puls method using the derived storage-
outflow relation. For level-pool reservoir routing with pumping or dam-break
simulation, outflows are computed for the orifice and weir equations for each
time interval.
Approach
to
4-Reservoir Pool
OG Spillway
37
The available energy head HE for flow over the spillway is computed as
HEAD.
HE = HEAD - [APLOSS *HAD.I ...... ................... (3.72)
where APLOSS is the approach loss at design head, HEAD is the water surface
elevation minus spillway crest elevation, and DESHD is the design head.
Design head is the difference between the normal maximum pool elevation and
the spillway crest elevation.
Pier and abutment energy losses are computed by interpolation of the data
shown in Table 3.7 based on HE/DESHD.
TABLE 3.7
Specific:
Energy/: : Approach
Design : : Depth : Pier : Abutment Contraction
Head, : Discharge : Adjustment : Contraction : Coefficients, KA
HE : Coefficient,: Exponent, : Coefficients,:
DKSHD : CC : EC : KP (3) Concrete (1) Earth (2)
(3) Pier contraction coefficients for type 3 piers are from Plate 7 of EA
1110-2-1603 (Corps of Engineers, 1965).
38
Effective length of the spillway crest ZEFFL is computed as
where SPUID is the spillway crest length, N is the number of piers, KP is the
pier contraction coefficient, and KA is the abutment contraction coefficient.
Hv = V .......................... (3.74)
COF=
COFQ = CC *
~c PDPTH K
--
T E . . . . . . . . . . . . . . . . .. . (3.75)
DESHD
where PDPTH is the approach depth to spillway, and CC and EC are interpolated
from Table 3.5 based on HE/DESHD. The spillway discharge QFREE assuming no
tailwater submergence is
's . . . . . . . . . . . . . . ... .. . . . . . . ...
QFREE = COFQ * ZEFFL * HE (3.76)
where het is the specific energy at toe of spillway, HE is the specific energy
at crest of spillway, ELSPI is the spillway crest elevation, and APEL is the
spillway apron (toe) elevation. Tailwater depth is then computed by trial and
error until
and aHE
D =E..................................(3.80)
HE + ELSPI - APEL - D
HE HE
39
The corrected flow is then
where QCORR is the spillway discharge corrected for tailwater submergence, and
SUBQ is the submergence coefficient in percent. A new corrected discharge is
assumed, and tailwater and submergence correction is computed until the change
in QCORR is less than one percent.
TABLE 3.8
Submergence Coefficients
(HE + D)/U MD/ME
1.07 1.10 1.15 1.20 1.30 1.40 1.50 1.60 1.70 1.80 1.90 2.00 2.25 2.50 3.00 3.50 4.00 4.50
PECENT SUINKIOGICK
100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 .00
55.0 54.0 52.0 49.0 45.0 42.0 40.0 39.0 38.0 38.0 37.5 39.0 40.5 43.0 53.0 58.0 60.0 60.0 .05
36.5 35.0 33.0 31.0 27.0 23.5 21.0 19.0 18.5 18.0 18.785 18.88 19.52 21.15 26.25 29.0 31.0 32.0 .10
27.5 25.0 22.0 19.5 17.5 15.5 14.0 13.5 13.0 12.5 12.45 12.21 12.63 13.44 15.0 17.0 18.3 21.0 .15
21.0 18.0 17.0 15.0 13.0 11.3 9.8 9.0 8.5 8.2 8.0 8.0 8.19 8.56 9.41 11.2 12.0 13.0 .20
18.0 15.5 13.5 12.0 10.0 8.4 7.2 6.0 5.4 5.0 4.9 4.914 5.375 5.688 7.0 7.85 8.5 9.0 .25
16.0 13.5 12.0 10.5 8.0 6.1 4.3 3.7 3.3 3.1 3.00 3.02 3.333 3.82 5.123 6.08 6.66 7.0 .30
15.0 13.0 10.0 8.0 5.5 3.6 2.5 1.8 1.7 1.5 1.450 1.438 1.625 1.88 2.717 3.73 4.19 4.5 .40
15.0 13.0 10.0 8.0 5.5 3.3 2.0 1.2 .96 .87 .857 .842 .853 .933 1.62 2.24 2.70 2.9 .50
15.0 13.0 10.0 8.0 5.5 3.3 2.0 1.1 .90 .75 .525 .515 .562 .600 .860 1.27 1.65 1.8 .60
15.0 13.0 10.0 8.0 5.5 3.3 2.0 1.1 .80 .50 .475 .450 .390 .385 .470 .69 0.93 t.0 .70
15.0 13.0 10.0 8.0 5.5 3.3 2.0 1.1 .70 .49 .450 .413 .323 .250 .110 .20 0.34 0.3 .80
15.0 13.0 10.0 8.0 5.5 3.3 2.0 1.1 .70 .49 .445 .410 .310 .220 .030 0.0 0.0 0.0 .85
15.0 13.0 10.0 8.0 5.5 3.3 2.0 1.1 .70 .49 .445 .400 .300 .200 0.0 0.0 0.0 0.0 .90
where CQFREE is the conduit discharge for unsubmerged outlet, COQL is the
discharge coefficient, CARKA is the conduit cross-sectional area, EL is the
reservoir water surface elevation, and ELEVL is the center elevation of the
conduit outlet. Tailwater elevation is interpolated from the tailwater rating
table and the corrected conduit flow is computed from
40
3.6.6 Average-Lam
NSTDL
1+ 2
QOUT(I) = 2
NSTDL NSTDL....................(3.86)
2
where LAG is the number of time intervals to lag inflow hydrograph, NSTDL is
the number of ordinates to average to compute the outflow, QIN is the inflow
hydrograph ordinate, Q is the lagged hydrograph ordinate, and QOUT is the
outflow hydrograph ordinate.
where QIN is the inflow hydrograph ordinate, and Q is the routed hydrograph
ordinate. This averaging is repeated NSTPS times to produce the outflow
hydrograph.
3.7 Diversions
41
3.8 Pumling Plants
Pumped flow may be retrieved at any point downstream of the pump location
in the same manner as a diverted hydrograph.
42
Section 4
PARAMETER CALIBRATION
The runoff parameters that can be determined in the calibration are the
unit hydrograph parameters of the Snyder, Clark and SCS methods and the loss
rate parameters of the exponential, Holtan, SCS, and initial/constant
methods. The melt rate and threshold melt temperature can also be optimized
for snow hydrology studies. If the Snyder method is employed, the Clark
coefficients will be determined and converted to the Snyder parameters.
n 1/2
STDER = [ E (QOBSi - QCOIPi) 2 * WTi/n]. .... ................ (4.1)
i=l
where QCOMPi is the runoff hydrograph ordinate for time period i computed by
HEC-1, QOBSi is the observed runoff hydrograph ordinate i, n is the total
number of hydrograph ordinates, and WTi is the weight for the hydrograph
ordinate i computed from the following equation.
43
Vt = (QOBS i + QAVE)/(2 * QAVE) ................. (4.2)
Limits of Optimitatiew.---mi
The optimization procedure does not guarantee that a "global" optimum (or
a global minimum of the objective function) will be found for the runoff
parameter; a local minimum of the objective function might be found by the
procedure. To help assess the results of the optimization, HEC-1 provides
graphical and statistical comparisons of the observed and computed hydrographs.
From this, the user can then judge the accuracy of the optimization result.
It is possible that the computed hydrograph will not meet with the criteria
44
TABLE 4.1
TC I 1.034t
R 1 .52
At = Computation Interval
Exponential SCS
Snowmelt
RTIOK 4 1.0
- l.116C 4 FRZTP 4 3.336C
Uniform Holtan
STRTL k 0 FC k 0
CNSTL 4 0 GIA k 1.0
BEXP a 0
1. Initial values are assigned for all parameters. These values may be
assigned by the user or program-assigned default values, Table 4.3, may
be used. In the example output, four parameters are optimized: unit
hydrograph parameters TC and R, and exponential loss infiltration
parameters STRKR and DLTKR (ERAIN and RTIOL are constant). In this case,
initial values were chosen by the user, STRKR = .20, etc. Note that the
unit hydrograph parameters TC, R are displayed as the sum (TC+R) and
ratio R/(TC+R) which are adjusted by the program during the optimization
process.
45
TABLE 4.2
OPTIMIZATION RESULTS
• WC01ARISON
OF COMPUTED AND OBSERVED RIDRDGRAPBS
TIME TO LAG
• SOF Boiv MEAN CENTER C.M. TO PEAR TIME OF
•FLOWS FLI OFOMASS
or C.K. FLOW PEAK
* PRECIPITATION EXCESS 0.937 4.13 *
• :
C OMUTED HYTROGRAPB 84787. 0.867 1390. 8.51 4.38 3621. 7.25
OBSERVED HYDROGRAPH 04787. 0.867 1390. 8.16 4.03 3540. 7.25
• DIFFERc 0. 0.000 0. 0.35 0.35 01. 0.00
P]CEN DIFFRECE 0.00 8.68 2.28
* STANDARD ERROR 270. AVERAGE ABSOLUTE ERROR 208.
SOBJECTIVE FUNCTIi 284. AVERAGE PERCUIT ABSOLUTE ERROR 27.27
46
TABLE 4.3
Initial
Parameter Value
Unit Graph
Clark TC+R (TAREA)I1/ 2
R/(TC+R) 0.5
Loss Rates
Holtan FC 0.01
GIA 0.5
SA 1.0
BEXP 1.4
47
equally-spaced values of the para-,ter with all other parameters held
constant. The "best" value of the parameter is then estimated using
Newton's method. This is demonstrated in the example by the asterisked
values of each of the optimization variables (e.g., TC+R = 6.895*, R/(TC+R)
= .522*, etc.). A parameter which does not improve the objective function
under this procedure is maintained at its original value. This is
indicated by a plus (+) in place of an asterisk (M) in the computed output;
this circumstance does not occur in the example.
5. Step 3 is then repeated for the parameter that most improved the value
of the objective function in its last change. This is continued until no
single change in any parameter yields a reduction of the objective function
of more than one percent. In the example this leads to changes to STRKR
and DLTKR.
The final results of the optimization are also summarized in Table 4.2,
TC = 3.16, R = 3.88, etc. Additional information is displayed comparing
computed and observed hydrograph statistics, which are defined as follows:
Average Percent
Absolute Error the average of absolute value of percent
difference between computed and observed
hydrograph ordinates.
48
4.1.3 Application of the Calibration Capability (from Ford et al., 1980)
Due to the varying quantity and form of data available for precipitation-
runoff analysis, the exact sequence of steps in application of the automatic
calibration capability of HEC-1 varies from study to study. An often-used
strategy employs the following steps when using the complete exponential loss
rate equation:
1. For each storm selected, determine the base flow and recession
parameters that are event dependent. These are not included in the
set of parameters that can be estimated automatically. These
parameters are the recession flow for antecedent runoff (STRTQ), the
discharge at which recession flow begins (QRCSN), and the recession
coefficient that is the ratio of flow at some time to the flow one
hour later (RTIOR).
2. For each storm at each gage, determine the optimal estimates of all
unknown unit hydrograph and loss rate parameters using automatic
calibration.
5. With ERAIN and RTIOL fixed, estimate the remaining unknown parameters
using the optimization scheme. Select a value of STRKR for each
storm being used for calibration. If parameter values for adjacent
basins have been determined, check the selected value for regional
consistency.
6. With ERAIN, RTIOL, and STRKR fixed, use the parameter estimation
algorithm to compute all remaining unknown parameters. DLTKR can be
generalized and fixed if desired at this point, although this
parameter is considered to be relatively event-dependent.
7. Using the calibration capability of HEC-1, determine values of TC+R
and R/(TC+R). Select appropriate values of TC+R for each gage. In
order to determine TC and R, an average value of R/(TC+R) is
typically selected for the region.
8. Once all parameters have been selected, the values should be verified
by simulating the response of the gaged basins to other events not
included in the calibration process.
49
certain hydrologic routing techniques. Criteria can be derived for the Tatum,
straddle-stagger and Muskingum routing methods only.
Inputs to this method are observed inflow and outflow hydrographs and a
pattern local inflow hydrograph for the river reach. The pattern hydrograph
is used to compensate for the difference between observed inflow and outflow
hydrographs. The assumed pattern hydrograph can have a significant effect on
the optimized routing criteria.
50
Section 5
KULTIPLAN-MULTIFLOOD ANALYSIS
The input conventions for the use of this option are described in the
input description. Section 10 gives specific examples on the use of data set
update techniques for the multiplan option. Example problems 9 and 10,
Section 12, illustrate the use of this HKC-1 option.
51
ra 10 raio
tt
Figre.1raItifood autiplnHdrga
520
Section 6
The dam safety analysis capability was added to the HEC-1 model to assist
in studies required for the National Non-Federal Dam Safety Inspection
Program. This option uses simplified hydraulic techniques to estimate the
potential for and consequences of dam overtopping or structural failures on
downstream areas in a river basin. Subsequent paragraphs describe dam
overtopping analysis, dam-break model formulation, the methodology used to
simulate dam failures, and the limitations of the method. An example of dam
overtopping analysis with HEC-1 is given in example problem 7, Section 12.
Example problem 8 simulates dam failures.
The dam safety simulation differs from the previously described reservoir
routing in that the elevation-outflow relation is computed by determining the
flow over the top of the dam (dam overtopping) and/or through the dam breach
(dam break) as well as through other reservoir outlet works. The elevation-
outflow characteristics are then combined with the level-pool storage routing
(see Section 3) to simulate a dam failure.
The discharge over the top of the dam is computed by the weir flow equation
Where h, is the depth of water over the top of dam, COQW is the weir
discharge coefficient, DANWID is the effective width of top-of-dam weir
overflow, and EXPD is the exponent of head. These variables are illustrated
in Fig. 6.1. The top-of-dam weir crest length, DANWID, must not include the
spillway. Spillway discharges continue to be computed by the spillway
equation (see Section 3) even as the water surface elevation exceeds the top
of the dam. The weir flow for dam overtopping is added to the spillway and
low level-outlet discharges.
Critical flow over a non-level dam crest is computed from crest length
53
E
011=00 * OAMWO * h XPO
whore
-
hI current Wm TePEL
d=c 2 1/ DA WI aX . + L2 . . . . . . . . . 6 3
wher
where ue yi the choaeo in ueain
cXrt o s aiacrossteto
i y. Z e t t scton
l f o (E(in
v r t l) t - pEsaThVWis)
e6.- am i
FloreaA is cmputed ascfo
rectangul
taeziar sections and ase 1/2wT(dcs
N2 CEL
awoatWOE- LL CO * ARE *'lg4
Thepflothrough h sectioni optdfo
5it
d................................................... (6.2)
an eeaidt. A da crs suc assoni i.6.ai rnfre
where gm is acceeaionade spogait y Thctotaflo oe the top of mi
Figure 6.2a Non-Level Damn Crest
55
then the sum of flows through each section (Fig. 6.2e). When a dam is being
breached the width of the breach is subtracted from the crest length beginning
at the lowest portion of the dam (Fig. 6.2f).
Dam breaks are simulated using the methodology proposed by Fread (National
Weather Service, 1979) with the exception that no reduction in the breach
discharge is made for submergence caused by downstream flow controls.
Structural failures are modeled by assuming certain geometrical shapes for the
dam breach. The variables used in the analysis, as well as the dam breach
shapes available in the program, are shown in Fig. 6.3.
Ieletloe of Soft e
~\
ELSM 1a
ofTOefE Mtsz
,- TOPEL TOPEL
I a I
56
Flow Q through a dam breach is computed as
s
Q = Cl * BWID * (WSEL - BREL)' + C2 * (WSEL - BREL) 2 .S .......(6.5)
where WSEL is the reservoir water surface elevation, BREL is the elevation at
base of breach, BRWID is the breach width, Cl is the broad-crested rectangular
weir coefficient (3.08), and C2 is the V-notch weir coefficient (2.44Z) with
side slope Z, horizontal to vertical.
The breach is initiated when the water surface in the reservoir reaches a
given elevation (FAILEL). The breach begins at the top of the dam and expands
linearly to the bottom elevation of the breach (ELB) and to its full width in
a given time (TFAIL). Note that the top-of-dam elevation must be specified to
fully determine the breach geometry.
6.3 Limitations
The dam-break simulation assumes that the dam-break hydrograph will not
be affected by tailwater constraints i.e., no correction for submergence of
the weir outflows is made. Also, the reservoir pool remains level. Also,
HEC-1 hydrologic routing methods are assumed appropriate for the dynamic flood
wave. Under the appropriate conditionp, these assumptions will be
approximately true and the analysis will give answers which are sufficiently
accurate for the purpose of the study. However, care should be taken in
interpreting the results of the dam-break analysis. If a higher order of
accuracy is needed, then an unsteady flow model, such as the National Weather
Service's DAMBRK (1979), should be used.
57
Section 7
PRECIPITATION
DEPTH-AREA RELATIONSHIP SIMULATION
58
A(A13 +B1)
A1 BAO+B
(A8+,8)
Precipitation Depth-Drainage
Area Function Legend
59
from subbasins A and B and are labeled (A15 + B15), (A13 + B13), etc., to so
indicate. The consistent hydrograph for the confluence of A and B is then
determined by interpolating between the two combined index hydrographs that
bracket the sum of drainage areas A and B, as shown on the Fig.7.1.
A2 A2 Ax A2
Q = Ql* (log Ax /log A2-)+ (Q2* log Al /log A2. ........ (7.1)
The interpolation formula would be exact if the loss function applied was
uniform and if the precipitation depth-drainage area relationship was in fact
a straight line on semilogarithmic paper. Because the interpolation formula
is not exact, the computer program insures that the peak of the interpolated
hydrographs below all confluences are not smaller than any of the interpolated
hydrographs above the confluence.
Operation of HEC-1 for the depth-area computation requires that the basin
be modeled (Section 2) and that the desired precipitation depth- drainage area
relationship be defined by up to five pairs of values that include the range
of tributary areas to be encountered. A different temporal pattern may be
specified for each depth-area point. Successive runs of the depth-area
feature with and without a proposed project will provide a balanced evaluation
of that project on downstream flood hydrographs. A single run will provide a
set of hydrographs at all locations within the basin that conform consistently
with the precipitation depth-drainage area function.
60
500 50m 103
1000 8I
LegeBd
* I~ouA8
611
Section 8
The basic technique used in the HAD analysis is to form the damage
frequency curve by combining damage versus flow (stage) and flow (stage)
versus frequency relations which are characteristic of the area that the
damage reach represents. The damage versus flow (stage) relation ascribes a
dollar damage that occurs in an area to a level of flood flow. The flow
(stage) versus exceedence frequency relation ascribes an exceedence frequency
to the magnitude of flood flow. By combining this Information, the damage
versus frequency curve and, hence, the HAD for a reach can be determined.
In the flood damage analysis, the conceptual model of the river basin
developed for a multiplan-multiflood analysis (example problems 9 and 10,
Section 12) is extended to include damage computations. Damage reaches are
designated by providing economic data, consisting of flow (stage) versus
frequency and flow (stage) versus damage data, for each damage reach in the
multiplan-multiflood model.
62
The development of the conceptual model for the flood damage analysis is
based on the'interrelated requirements for the stream network and damage
calculations. This relationship is shown on Fig. 8.1 where subbasins,
routing reaches, and damage reaches are delineated for an example river
basin. The definition of the subbasins and routing reaches for the stream
network calculations is determined in part by criteria outlined in Section 2,
and in part by the requirements of the damage calculations.
Note that the damage reach may encompass parts of a number of routing
reaches. The flows used in the damage calculation are based on the outflows
from the most downstream of these routing reaches. The flows combined with
damage data for the index location result in the appropriate damage for the
entire damage reach.
11 ACH
Divisioan
INDEX
LOCATION
63
8.3 Damage Reach Data
64
0 FLOW FREQUENCY DATA
Q1 GENERATED FROM
MULTIFLOOO ANALYSI
-Z.-CUIC SPLINE
F?? OF OMT
Q&
DISCHARGE
Q
Q
,- I
too f' o fI I
GENERATED
FLOW-FREQUENCY
QE CURVE
DISCHARGE
Q
Q1
- I I I
100 ft 10 f, I
65
MODIFIED
CURVE
DISCHARGE
(0) Q!
@ 2 - -
- I EXISTING
I
G~ ~
_ _ _ I _
ADOPTEO
fi fa
FREQUENCY (f)
flow-damage data for each PLAN and damage reach (HEC, 1979a) using the
following methodology.
In the case that stages are used, the procedure is the same except
that the stages for generated frequencies are determined using a
linear interpolation procedure. If stages are specified for the
damage data and flows for the frequency data (or vice versa), a rating
curve is used to relate the stages and flows before determining the
appropriate damage.
66
as shown in Fig. 8.4. This extrapolation will not severely affect the
accuracy of the result if.the peak flows generated result in a
relatively small exceedence frequency.
3) The Integral of the damage-frequency curve is the KAD for the reach.
This area is computed using a three point Gaussian Quadrature formula.
4) If more than one damage category is specified for a reach, the above
steps are repeated for each category. The KAD is suined for all the
categories to produce the WD for the reach.
DISCHARGE DISCHARGE
01- Qa-
SENERATED TRPOLAE
DAMAGE
"RUENC'r
flft ft D1 Di Da 4m
AME
FLOW FREQUENCY CURVE FLOW DAMAGE CURVE
~EXTRAPOLATED
hL fn 0
EXCEEDENCE FREQUENCY (FRACTION) (FREQUENCY OF MAXIMUM
SIMULATED FLOW
67
Section 9
68
Uniform Level of Protection. A flood control plan may require that, as
part of the flood control system, levees (local protection projects) provide
the same level or a uniform level of protection at a number of locations
(damage reaches). In this instance, the level of protection refers to the
flood exceedence frequency at which the capacity of the project is surpassed.
The flood control system optimization option can be used to determine the
uniform level of protection that, in conjunction with the structural flood
control components, leads to the maximum net flood loss reduction benefits in
the river basin.
Cost and Capacity Data. Two types of data are required to calculate the
total annual cost of a flood control component. First, capacity versus
capital cost tables are required to determine the capital cost for any
capacity of the flood control component. A capital recovery factor is also
required so that equivalent annual costs for the capital investments can be
computed. Second, operation and maintenance costs are computed as a
proportion of the capital cost. For pumping plants, average annual power
costs for various pump capacities are required. Pump operation costs are
computed in proportion to the volume pumped. Capital and operating costs for
non-optimized components of the system may also be considered.
69
The optimization procedure can be generally described as follows:
(2) The model then uses the univariate search procedure to find a minimum
value for the objective function. (The optimization algorithm is the
same as used for parameter optimization, Section 4.) The procedure
finds a minimum by systematically altering flood control component
capacities in order to calculate various values of the system
objective function. Each time a flood control system capacity is
changed, stream network calculation and EAD calculations are
performed giving a value for the system objective function.
(3) Once the optimization procedure is completed, the costs, damage and
net benefits accrued to the optimized system are computed and
displayed.
where TANCST is the flood control system total annual cost, ANDMG is the river
basin expected annual damage, ODEV is the sum of the weighted deviations from
the target flow or stage, and CONST is a term representing the importance of
the target penalty (default value equal to 1.0). As CONST increases, the
target penalty has less importance in determining STDER.
where ANFCST is the sum of the equivalent annual capital costs for the flood
70
control components, ANOMPR is the sum of the annual operation, maintenance,
power and replacement costs for the flood control components, FDCNT is the
equivalent annual capital cost for non-optimized components, and FAN is the
annual operation, maintenance, power and replacement cost for non-optimized
components.
The annualized capital and operation and maintenance costs are computed
as follows.
where CAPCST is the capital cost of a flood control project, CRF is the
capital recovery factor for a specified project life and interest rate, and
FCAP is the total capital cost of the non-optimized components of the system.
FDCNT may be computed as shown above or the equivalent annual capital cost may
be specified as direct input.
where TRGT is the target flow specified by the user for a given exceedence
frequency, and TMP is the computed flow for the given exceedence frequency with
the flood control projects in operation, i.e., PLAN 2. The exceedence
frequency specified for the target penalty is used to interpolate a value of
TMP from the PLAN 2 flow-frequency curve computed fqr a reach. The
interpolation is accomplished by using the cubic-spline fit procedure.
The penalty, PEN, for deviations from the target conditions are calculated
for stages as:
The sum of the penalties for all reaches is equal to the deviation
penalty ODEV in equation 9.1. The factors CONST (equation 9.1) and ANORM can
be adjusted by the user (ANORM should be greater than or equal to .02) until
satisfactory compliance with the target constraints are met by the
optimization procedure. The default values for these parameters should
suffice for most purposes.
71
Section 10
This section describes: the general organization of the input data, special
features for specifying data, and groupings of data to accomplish specific
simulation options. A detailed descripiton of the individual input data
records and their contents is given in the Appendix A: Input Description.
There are two general types of data records for HEC-l: input control and
river basin simulation data. The input control records tell the program the
format of the river basin data as well as controlling certain diagnostic output.
All input control records begin with an asterisk (0) in column one follwed by a
command. These input controls are discussed in the next subsection and a
detailed explanation is given in Appendix A.
The river basin simulation data are all identified by a unique two-
character alphabetic code in columns one and two of each record. These codes
serve two functions: they identify the data to be read from the record; and
they activate various simulation options. The first character of the code
identifies the general category and the second character identifies a specific
type of data within a category. An overview of these data categories and codes
is shown in Table 10.1. The flood damage data, beginning with the EC record is
placed at the end of the river basin simulation data. These data are not all
labeled as E records because the record code and format were taken from the
Expected Annual Flood Damage (HEC, 1979) program. Thus these same data records
may be used directly in both programs.
The river basin simulation data records are structured by the user to
reflect the topology of the basin. The sequence of the input data prescribes
how the river basin is simulated. There are three general subdivisions of these
data as shown in Table 10.2: job control; hydrology and hydraulics; and
economics. Example input data for a simple river basin are shown in Fig. 10.1.
The data model of a river basin can be thought of as a series of building
blocks, each block beginning with a KK record. The data following each KK
record identifies the type of operation to be performed, e.g., BA signifies
subbasin runoff and R_ signifies a routing. Section 12 gives examples of input
data structures to accomplish various program options.
There are six input control commands: *FREE, *FIX, *LIST, *NOLIST,
*MESSAGE, and *DIAGRAM. Data can be input to the HEC-1 model in a fixed and/or
free format as noted in the Input Data Description. The traditional HEC fixed-
format input structure (ten 8-column fields) is the default option of the
program. The program now provides the capability to enter data in a free
format. All records following a *FREE record in the data will be considered as
being in free format. Free format data fields are separated by commas or one
or more spaces, and successive commas represent blank fields. The fixed format
can be returned to at any point in the data set by providing a *FIX record.
The *FIX will be in control until another *FREE record is encountered, etc.
72
TABLE 10.1
Data Record
Category Identification Description of Data
73
TABLE 10.1: HEC-1 Input Data Identification Scheme (Cont'd)
Data Record
Category Identification Description of Data
74
TABLE 10.1: HEC-1 Input Data Identification Scheme (Cont'd)
Data Record
Zatexory Identification Description of Data
75
TABLE 10.2
Hydrology Economics
Job Control & Hydraulics & End of Job
I_, Job Initialization K_, Job step control E_, etc., Economics,
V, Variable Output Summary H_, Hydrograph trans- data
0_, Optimization formation ZZ, End of job
J_, Job Type Q, Hydrograph data
B_, Basin data
P, Precipitation data
L_, Loss (infiltration) data
U_, Unitgraph data
N_, Melt data
R_, Routing data
S_, Storage data
D_, Diversion data
W_, Pump Withdrawal data
The IN record allows the user to enter time-series data, either hyetographs
or hydrographs, at time steps other than the computation interval specified on
the IT record. This option is convenient when entering data generated by
another program or in a separate HEC-1 simulation. Note that if direct input
unit hydrograph ordinates is used (UI record), they must be at the same time
step as the simulation computation interval and cannot be input with the IN
record.
Data groups for subbasin runoff simulation which need not be repeated (if
they are the same as input for the previous subbasin) are shown in Table 10.3.
HEC-1 automatically uses the previous subbasin's input data for these data
types unless new data are provided for the current subbasin. The source of
76
the data used as identified by tie input record number is printed in the left
hand margin. If a zero is printed as the input record number, this means no
data records have been provided, up to that point, which contain the required
Information. Great care should be taken to verify that the input data used
was so intended. No data are repeatable for routing reaches.
TABLE 10.3
Data Repetition Options
Rainfall P
Infiltration L
Base Flow BF
Snowmelt M
*Unit Hydrograph US, UC, UD
*Kinematic Wave " UK, RK
In the multiplan analysis, data may be supplied for a number of plans for
the same subbasin. Data need not be repeated for each plan by following two
conventions:
(1) Plans not specified in the data set by a KP record are assumed to be
the same as the first plan in the KK record group. (Data for a
particular plan follows a KP record in the data set.)
77
computation as shown in Fig. 10.2 The gage data consists of PG records for
nonrecording gages and PG and PI or PC records for recording gages. These
data are usually grouped toward the beginning of the data set before the first
KK-record runoff computation. Within each KK-record group, the (PR, PW) and
(PT, PW) records are used to specify which gages and corresponding weightings
are to be used for computation of that subbasin's average precipitation. Note
that a recording gage can be used as both a storm total and a recording &age
station. This is indicated by using gage WEST of PT and PW records in Fig.
10.2. If the storm total value is not specified on the PG record for the
recording station (as is the case for the Fig. 10.2 example), the program sums
the incremental values on the PI records to compute that value.
Table 10.8 Spillway Routing, Dam Overtopping and Dam Failure Input
Data Options
10.4 Input Data Retrieval from the HEC Data Storage System (DSS)
The HEC Data Storage System , DSS (HEC, 1984), may be used to supply
certain catchment characteristics and time-series data to the HEC-1 input data
set. Those data are runoff parameters stored by program HYDPAR (Corps of
Engineers, 1978), cumulative and incremental precipitation (PC and PI data),
and streamflows (QI and QO data). The input connections used to retrieve data
are given in the overview of HEC-1 usage with DSS in Appendix B. Access to
DSS is limited to HEC-supported computers, and requires a special version of
HEC-1 and associated DSS software.
78
A Subbasin
Runoff
Combining
Routing C
Card ID Description
ID Title
IT Time interval and beginning time
10 Output control option for whole job
KK Subbasin A
Runoff BA Area
from BF Base flow
Subbasin P Select one precipitation method, use IN if necessary
A L Select one loss rate method
U_ Select one rainfall excess transformation method
KK
Subbasin B BA Similar to above for Subbasin A
runoff BF
P,L_,u
KK
Subbasin C BA Similar to above for Subbasin A
runoff BF
, P_,L_,U_
79
lines delineating
areas pertaining
to rainfall gages
gage WEST
0 gage EAST
A .8 sq i Symbols
A - Recording gage
% o - Nonrecording gage
g0 OTotal Rainfall
EAST - 4.8"
WEST - 3.5"
SOUTH - 4.2"
DATA INPUT
Card
ID Data
ID
IT
(PG EAST 4.8
PG SOUTH 4.2
Rainfall gage data PG WEST
I .02 .05 etc. recording gage
readings for storm
KK 3-gage basin
BA 2.0
Gage weightings for (PT WEST EAST SOUTH
basin-average total jPW .4 .4 .2
Gage weightings for PR WEST
basin-average recorder PW 1
L_
U_
80
TABLE 10.4
TABLE 10.5
Basin Area BA BR BA BA
Base Flow BF BF
81
TABLE 10.6
Optimization Control OU OR
Pattern Hydrograph QP
TABLE 10.7
Modified Puls
Given Storage Normal-Depth
Type of Data Muskingum Outflow Storage Outflow Kinematic Wave
Routing Control RM RS RS RK
82
TABLE 10.8
TABLE 10.9
Damage Reach ID KK
Optional records
83
TABLE 10.10
Optimization OS
Target Penalty 00
Discount Factor + WO SO DO LO
Size Constraint
Degree of Protection DP
Table 10.11
Combination HC
Punch *KO, KF
The use of these options must be in combination with some other hydrograph
computation
84
Section 11
PROGRAM OUTPUT
A large variety and degree of detail in the printer output are available
from HEC-1. This section describes the output in terms of input data feedback,
intermediate simulation results, summary results, and error messages. The
degree of detail of virtually all of the program output can be controlled by
the user.
Several of the sunuary outputs are printed from scratch files generated
during the simulation. If the user desires to save these scratch files for
use in other jobs (say, for a plotting device), their location can be found in
the definition of Input/Output Fortran logical units in Table 13.1 of Section
13.
The input data file for each job are read and copied to a working file.
As the data are copied to the working file they are converted from free format
to fixed format (see Section 10.2.1) and a sequence number is assigned to each
line. The reformatted data are printed so the user can see the data which are
going into the main part of the program.
85
plot. For snowmelt calculations, separate values cf loss and excess are
printed for rainfall and snowmelt. For storage routines, storage and stage
(if stage data are given) are printed/plotted along with discharge.
For optimization jobs (unit graph and loss rate, routing, or flood
control project sizing), the program prints values for the variables and
objective function for each iteration of the process. This output should be
carefully reviewed to understand why changes are being made in the variables
and to verify (using engineering judgment and comparison with similar
results) that the results are reasonable.
The HEC Data Storage System, DSS (HEC, 1984), may be used to save HEC-1
output information for use in another HEC-1 simualtion or by other HEC
computer programs. Time-series data, streamflow or stage, as well as paired-
function data, flow-frequency curves, can be output to DSS. The means by
which this data can be stored is given in the overview of HEC-1 usage with
DSS in Appendix B. Access to DSS is limited to HEC-supported computers, and
requires a special version of HEC-1 and DSS software.
Table 11.1 lists error messages (in capital letters) which HEC-I will
print along with an explanation of the message. Some errors will not cause
the program to stop execution, so the user should always check the output for
possible errors or warnings. The array dimensions listed in Table 11.1 are
those used by HEC-1 on a mainframe computer.
The computer operating system may also print error messages. When an
error occurs, the user should first ascertain if it is generated by HEC-1 or
by the system. If it is generated by HEC-1, i.e., in the format given in
Table 11.1, that table should be referred to and the indicated actions
taken. If the error is system generated, the computer center user service
and/or the in-house computer systems personnel should be contacted to
ascertain the meaning of the error. These errors may be due to incorrectly
input or read data or errors in HEC-1 or the computer system. If these
system errors cannot be resolved in-house or if there is an error in the
HEC-1 program, the HEC should be contacted.
86
TABLE 11. 1
NEC-I Error Messages
87
TABLE 11.1:
HEC-1 Error Messages (Cont'd)
88
TABLE 11.1:
HEC-1 Error Ressages (Cont'd)
99
Section 12
EXAMPLE PROBLEMS
A stream network model was developed for the Red River watershed shown in
Fig. 12.1. The development of this type of model for a watershed is basic to
the use of the HEC-1 program. The example demonstrates the following features
of the program:
Tables 12.1a-12.lc display data for the watershed model; note that the
data record identifiers used to input each type of data are also indicated in
the tables. Important points to note about the stream network model data are
as follows:
(a) Both recording and non-recording gage stations can be used as total-
storm stations for a subbasin as specified on the PT, PW cards. (The total
depth associated with incremental or cumulative rainfall data is automatically
calculated for each recording gage.) In this example, gage 400 is used only
for the temporal pattern. The subbasin storm pattern is calculated as a
weighted average of the recording gage storm patterns indicated on the PR, PW
cards.
(b) The various unit hydrograph options available can be used with any of
the loss rate options. The data in the appropriate HEC-1 format and the
results of the computer simulation are displayed in the Table 12.1d computer
output.
90
LO$T
DR
CHANNEL RoITt 101030
LOSTm OWERT
Figur Strem
12. Nework odelschmai
WEST ED
91
TABLE 12. la
TABLE 12.Ib
SUBSASIN PHYSICAL PARAMETERS (Test 1)
PRECIPITAT ION
SUBBASIN BASIN AREA GAGE WEIGHTS BASE FLOW PARAMETERS
NAME (SQ MI) (PT, PW CARD) LOSS RATE UNIT GRAPH (BF CARD)
(KK CARD) (BA CARD) GAGE 0 NT. (METHOD) (CARD) (METHOD) (CARD) STRTQ QRCSN RTIOR
RED RI .82 400 1 SCS LS SCS U9 10.0 -2.5 1.2
60 .75 N-B=0 E=-.47
61 .25
EASTIO .66 400 1 EXPON. LE SNYDER US 10.0 -.25 1.2
61 .6 9rW-o. 6 1pF:
62 .3 DLTKRrl.0 CP--O.8
63 .1 RTIOR=1.0
ERAIN .0
LOSTBR .36 400 1 UNIFORM LU CLARK UC 10.0 -. 25 1.2
62 .5 STR--0.3 T-8
63 .5 CNSTL-.04 R=1.2
92
TABLE 12.lc
Reach: lOto2O KK
VOLUME-OUTFLOW DATA RS
Reach: 20to30 KK
VOLUME-OUTFLOW DATA RS
Reservoir: E.DAM KK
SPILLWAY SS
Crest elevation - 856 0.8.1.
Width = 60 feet
Weir coefficient 2.7
Head exponent = 1.5
VOLUME-ELEVATION DATA
DIVERSION DATA
Location: EAST1O KK
DIVERSION DESIGNATION DT
Diverted flows labeled: DIVERT
93
TABLE 12.ld
X X XXXXXXX XXXXX X
X X X I X XX
X X X X X
XXXXXXX XXXX X XXXXX X
X X X X X
X XX X X X
X X XXXXXXX XXXXX XXX
THIS PROGRAM REPLACES ALL PREVIOUS VERSIONS OF HEC-1 KXNMN AS HECI (JAN 73), NWIUGS. HIEC1DB, AND HEC1KW.
THE DEFINITIONS OF VARIABLES -RTIMP- AND -RTIOR- HAVE CIANED FROM THOSE USED WITH THE 1973-STYLE INPUT STRUCTURE.
THE DEFINITION OF -A.4SKK- ON RN-CARD) WAS CHANGED WITH REVISIONS DATED 28 SEP 81. THE VERSION RELEASED 313AN85
CONTAINS HEN OPTIONS ON RL AND BA RECORDS, AND ADDS THE RL RECORD. SEE JANUARY 19S5 INPUT
DESCRIPTION FOR NEW DEFINITIONS.
LINE ID ....... 1 ....... 2 ....... 3 ....... 4 ....... 5 ....... 6 ....... 7 ....... 8 ....... 9 ...... 10
13 XX REDRI
14 KO 4
15 KM SCS RNOFr COMPUTATION
16 BA .82
17 aF 10.0 -. 25 1.2
18 PR 400
19 Pw 1
20 PT 60 61
21 P" .75 .25
22 LS 90
23 UD 1.47
24 XI EAST10
25 KO 4
26 KM SNYDER UNIT GRAPE COMPUTATIOW-EXPONDITIAL LOSS RATE
27 BA .66
28 BF 10.0 -. 25 1.2
29 PR 400
30 Pw 1
31 PT 61 62 63
32 Pw .6 .3 .1
33 LE .6 1.0 1.0 0
34 US 1.3 .8
35 KI EASTIO
36 KM DIVERT rLOM TO LSTAR
37 DT DIVERT
3 DI 0 100 300 600 900
39 DQ 0 25 100 180 270
40 I RED10
41 KM4 COMBINE HYDROGRAPHS FROM SUBBNSINS EAST10 AND RED RI
42 sc 2
43 KR 10T020
44 g0 4
45 KM ROOTM FLOWS rR)M STATION RLD01 TO RED 20
46 RS 1 PLO* -1
47 SV 0 18 36 54 84 110 138 174 228 444
48 SQ 0 500 1000 1500 2150 2600 3000 3450 4000 6000
94
IS-1 INMT PAGE 2
LIZI 10 ....... 1 ....... 2 ....... 3 ....... 4 ....... 5 ....... 6 ....... 7 ....... 0 ....... 9 ...... 10
0 IK LOSTBR
so KN RETRIEVE DIVERSION FROM EAST10
51 DR DIVERT
52 KK L06T
53 K4 CIAR4 UNIT GRAPH COMPUTATION-INITIAL AND UNIFORM LOSS RATES
54 SA .36
55 8? 10.0 -. 25 1.2
56 PR 400
57 p1 1
58 PT 62 63
59 P .5 .5
60 LU .3 .04
61 uc .80 1.2
62 KR LOSTBR
63 KM COMINE RUINOPF FROM LOSTOR VITO DIVRTE FLOW
64 BC 2
65 KR S.D4
66 KM ROUTE FLOWS TIHOUGH DAN
67 RS 1 EXV 651.2
68 SV 21 100 205 325 955
69 Se 850 651.5 853.3 856.5 858.0
70 St, 851.2 12 .6 .5
71 SS 856 60 2.7 1.5
72 KR WEST20 P
73 Km SCS RUMOF COMPEUTATION-EOLTAN LOSS RATE
74 gO 1 2
75 SA .80
76 B? 10.0 -. 25 1.2
77 PR 400
78 P11 1
79 PT 63 64
80 Pw .6 .4
81 LX .04 .4 .3 1.4
82 UD .94
83 RK RED20
84 KK COMBINE RUN" FROM WSST20,OU/PLOW FROM K.DAM AND REACH 1020
85 HC 3
86 KR 20TO30
87 K4 ROUTE PLOWS FROM RED20 TO R,130
88 RS 1 FLOW -1
89 5V 0 17 42 67 100 184 274 386 620
90 SQ 0 500 1000 1500 2000 3000 4000 5000 7000
91 KR RED30
92 KM RUNOFF BY THE SCS pITND
93 BA .19
94 BF 10.0 -. 25 1.2
95 PR 400
96 PW 1
97 PT 64 63
98 P" .65 .35
99 LB 79
100 UD 1.03
101 91 RED30
102 KM COMBINE RUNOFF FROM RED30 AND OUTLOW FROM REACH 20TO30
103 HC 2
104 KK GAGE
105 sO 1
106 KM AND OBSERVED YDNP=RAPO
COMPARE COMPUTIED S AT RED30
107 IN 15 12JUN68 715
108 O 10 13 16 20 25 30 51 92 159 241
109 O 332 399 412 393 348 291 255 229 235 321
110 0O 472 705 921 1120 1255 1345 1373 1314 1228 1122
111 QO 996 900 817 742 668 614 549 500 444 409
112 gO 388 372 359 348 330 328 321 310 300 291
113 QO 282 274 267 277 252 240 231 224
114 zz
95
SCHMATIC DIAGRAM Of STIRAN METWOMS
INPUT
LINE (V) ROUTNG (--) DIVERSION
NO. (.) coWNEt (4--) RETURN O DrIV tTD F.
13 R RI
24 . FAST1O
37 ------- DIVERT
35 ZAST0
40 RE10 ............
V
V
43 101020
51 ..-------. DIVERT
49 LOSTBR
52 LOSTBR
62 LOSTa ............
V
V
65 E.DAK
72 WEST20
83 RPD20 ........................
V
V
of 20T030
91 RED30
Yb
U.S. ARMYCORPS OF ENGINEERS
FLOW NYROGRAP PACKAGE (N1C-9)8* 1E NpROLOGIC ENGINEERING CENTER
T
REIED 195* 609 SECAID STREET
J a• DAVIS, CALIFORNIA 95616
9XAMPLZ1PROL, 1O. 1
STRWIA NtIORK MOD
ENGLISH UNITS
DRAINAGE AREA SQUARE MILES
PRECIPITATION DEPTH INCHES
LENGTH , ELEVATION FEET
FLOW CUBIC FEET PER SECOND
STORAGE VULME ACRE-PEET
SURFACE AREA ACRES
TEMPERATURE DEGRES FARENHEIT
*he.* *** act ac e hat ata tee .. aat *a te teee ate a.. ae• tea at. en ta t•a *a eta aeee tat*eaa* .at Ceeeeta .e. te**
atatetee Ca~e
13K * RED RI
* a
at...........
16 BA SUBBASIN CHARACTERISTICS
TARFA 0.82 SUBBASIN AREA
97
17 sp BABS FLi CHARACTERISTICS
STRTQ 10.00 INITIAL FL0
QRMN -0.25 BEGIN BASE FLOW RBESbION
RTIOR 1.20000 RCES ION CONSTANT
PRICIPITATION DATA
TEMPORAL DISTRIBUTIONS
UNIT HYDROGRAPH
31 EWD-OP-PRRIOO ORDINATES
17. 51. 105. 175. 226. 247. 246. 222. 191. 149.
110. 85. 66. 52. 40. 31. 24. 19. 14. 11.
9. 7. 5. 4. 3. 3. 2. 2. 1. 1.
0.
*.e. ... ... .* .. l ... e... ... . . .* *.. ** *.. *.. **. ... .*. ... te. *.. *t* ..
24 KX EAST10
Z
27 BA SUBSASIN CHARACTERISTICS
TAREA 0.66 SUBBASIN AREA
98
28 BP BASE PLOWCHARACTERISTICS
STRTQ 10.00 INITIAL PLOW
QRCSN -0.25 BEGIN BASE FLOW RECESSION
RTIOR 1.20000 RECESSION CONSTANT
PRECIPITATION DATA
34 US SNYDER UNITGRAPB
TP 1.30 LAG
CP 0.80 PEAKING COEFFICIENT
TEMPORAL DISTRIBUTIONS
UNIT HYDROGRAPE
16 ENITOF-PERIOD ORDINATES
23. 79. 146. 210. 252. 261. 238. 181. 116. 73.
46. 29. 18. 12. 7. 5.
43 R * 10T020 *
99
II WIU I iW t ,k ,L
46 RS STORAGE ROUTING
NSTPS I NUMDER OP SUBREACHES
ITYP FLOW TYPE OP INITIAL CONDITION
RSVRIC -1.00 INITIAL CONDITION
X 0.00 NORKING R AND D COEFICIENT
47 SV STORAGE 0.0 18.0 36.0 54.0 84.0 110.0 138.0 174.0 228.0 444.0
48 SQ DISCHARGE 0. 500. 1000. 1500. 2150. 2600. 3000. 3450. 4000. 6000.
*4* *4* 44* *44 *** ** *4*, 4*e *** *4*44 ** 4*4*4e* *4* 4*4 *4* *4*o 4*t ** *44 e* */ 4* 444 ** *44 4*4 4*4 44, *4*e *44
* *
72 KK * WEST20 *
75 BA SUDDASIN CHARACTERISTICS
TARER 0.80 SUBBASIN AREA
PRBCIPITATION DATA
100
APILFLLTATION bTATIUN DATA
TZKPORAL DISTRIBUTIONS
DA MN BRi ORB RAIN LOSS EXCESS OM * DA MON BUS O11 RAIN LOSS EXCSS Com
12 JUN 0715 1 0.00 0.00 0.00 10. * 12 J 1430 30 0.01 0.01 0.00 211.
12 JUN 0730 2 0.03 0.02 0.01 10. • 12 JUN 1445 31 0.01 0.01 0.00 200.
12 JUN 0745 3 0.03 0.02 0.01 11. * 12 JUN 1500 32 0.01 0.01 0.00 192.
12 JUN 0800 4 0.03 0.02 0.01 15. • 12 JUN 1515 33 0.01 0.01 0.00 1e.
12 JUM 0815 5 0.03 0.02 0.01 19. * 12 JUN 1530 34 0.03 0.01 0.01 175.
12 JUN 0830 6 0.30 0.02 0.28 36. * 12 JUN 1545 35 0.03 0.01 0.01 167.
12 JUN 0845 7 0.30 0.02 0.28 81. * 12 JUN 1600 36 0.03 0.01 0.01 160.
12 JUN 0900 8 0.30 0.02 0.28 164. * 12 JUN 1615 37 u.03 0.01 0.01 152.
12 JUN 0915 9 0.30 0.02 0.29 263. * 12 JUN 1630 38 0.01 0.01 0.00 146.
12 JUN 0930 10 0.01 0.01 0.00 344. * 12 JUN 1645 39 0.01 0.01 0.00 139.
12 JUN 0945 11 0.01 0.01 0.00 377. * 12 JUN 1700 40 0.01 0.01 0.00 133.
12 JUN 1000 12 0.01 0.01 0.00 343. * 12 JUN 1715 41 0.01 0.01 0.00 127.
12 JUN 1015 13 0.01 0.01 0.00 275. * 12 JUN 1730 42 0.00 0.00 0.00 121.
12 JUN 1030 14 0.03 0.02 0.01 201. * 12 JUN 1745 43 0.00 0.00 0.00 116.
12 JUN 1045 15 0.03 0.02 0.01 139. * 12 JUN 1800 44 0.00 0.00 0.00 Ill.
12 JUN 1100 16 0.03 0.02 0.01 99. * 12 JUN 1815 45 0.00 0.00 0.00 i06.
12 JUN 1115 17 0.03 0.02 0.01 91. * 12 JUN 1830 46 0.00 0.00 0.00 101.
12 JUN 1130 18 0.63 0.02 0.62 87. * 12 JUN 1845 47 0.00 0.00 0.00 7.
12 JUN 1145 19 0.63 0.01 0.62 160. * 12 JUN 1900 40 0.00 0.00 0.00
12 JUN 1200 20 0.63 0.01 0.62 343. * 12 JUN 1915 49 0.00 0.00 0.00 to
12 JUN 1215 21 0.63 0.01 0.62 556. * 12 JUN 1930 50 0.00 0.00 0.00 8
12 JUN 1230 22 0.18 0.01 0.17 740. * 12 JUN 1945 51 0.00 0.00 0.00 8
12 JUN 1245 23 0.18 0.01 0.17 839. * 12 JUN 2000 52 0.00 0.00 0.00 77
12 JUN 1300 24 0.18 0.01 0.17 817. * 12 JUN 2015 53 0.00 0.00 0.00 74.
12 JUN 1315 25 0.18 0.01 0.17 729. * 12 JUN 2030 54 0.00 0.00 0.00 70.
12 JUN 1330 26 0.02 0.02 0.00 619. * 12 JUN 2045 55 0.00 0.00 0.00 67.
12 JUN 1345 27 0.02 0.01 0.01 503. * 12 JUN 2100 56 0.00 0.00 0.00 64.
12 JUN 1400 28 0.02 0.01 0.01 393. * 12 JUN 2115 57 0.00 0.00 0.00 61.
12 JUN 1415 29 0.02 0.01 0.01 294. * 12 JUN 2130 58 0.00 0.00 0.00 59.
SSUN 4.98 0.55 4.43
101
09 I xx I.i a
S-.
~~. o 0 , .
o .
0 .o~. . .. . . . . . . . . . , , . .. o , o . . . . . . . .o
* 0
o ~, ~
t o . ~ ° o
.0 . .
o .
0 .
, 0 o 0 0 o o o,0,
.. . o . .... ....
. . . . ... .. . . . . .0 . . . . . . . . . . . . . . . . . . . . . .
.00 0
.. . . . . . . . . .. . . . . . . . . . .
0
0 0
o , 0 , 0 0 0 . , , , , 0I, * 0 0 , 0 , o 0 , 0 , 0 0 0 0 °
*; 0.
. . .. ....... .............. .
...........
. . . . . . . . . ..
S0
0
:t.0
0
S0 0 o0 - 4 00
I 0 .. 0
I. .. 0000
000 000
o*00
°
8000
0t
102
... *0* *0* et e• 0. 0*0 *0* 0*0 *04 *** *0* * *0*
•• *0* *** 0** *0* 4*0 *0 • ** *0* *0* *•* *0* *0* *0* 4*4eoe *4* 4 .*4*
* 0
104 xg GAGE •
* •
* TIME TO LAG0
• o
fU OP OOXY MEA CUSTER C.N. TO pIA TINE OF
PLO" OEM rF or MASS C.N. PI flr[
•• STANDARD
O.JUCTIVE AOR
PmOTIOS 21.
22. AV33m AVERAGE ANSOWE
PUnGT &66icx im om '6
18.56 •
*
103
I I
I * ° ° *
.'° . . . . o o o o o , . o * o o , o * * . , * o o I
. . . . . . . . . . . . . . . .
S . . . . . . . . . . . . . . . . . . . .i
. . . . . . . . . .i
. . . . . . . . . .i
. . . . . . . . . .i
......:
. . . . . . . . . . . . . . . . . . .
*o
I
..... . . . .
. . . . . . . . .
o 0"
. .I . .
J . . . . . .. . . . . . .. . . . . .. . .. .. .. . . . . .
0 0
* 0
0
.................... ......... . . 0 ........
... . . . .
0
I
. . . . .. ... . .. .
o , o . o1
I *' I
RKW?
CUICAM
SWOIAR
Fur Iu Bo6RD
TINS IN IN-ltI,
FLOW O NLU
U
O
SIU30G 1S AT uUM 20 639. 5.50 367. 210. 210. 0,0
A?
i 31020 1370. 6.00 733. 444. 444. 2.64
3 CS3O
201030 1224. 6.50 729. 439. 439. 2.64
loIwD to
AT 530 143. 5.75 61. 34. 34. 0.19
SToGwIs1
lf)30 1331. 6.50 709. 473. 473. 2.53
2 Cm I M
105
12.2 Example Problem #2: Kinematic Wave Watershed Model
The input data for the watershed are displayed in Tables 12.2a -
12.2c. The HEC-l data model for the basin is shown in Table 12.2d. There
are a number of important points to note about the data:
(1) Each subbasin has data for two overland flow elements (only
one is required) which is specified on the UK card. The two
elements represent separately the impervious and pervious
areas of a subbasin.
(2) Collector channel and main channel data are specified on the
RK card for each subbasin. As many as two collector channels
can be specified for each subbasin, however, only one
collector channel was used In this example.
The simulation results are displayed in Table 12.2d following the input
listing.
SUB 4
SUB 4
106
TABLE 12.2a
Subbasin Characteristics
OVERLAND FLOW PLANE DATA
(UK RECORD) LOSS RATE
DRAIN (S RECORD)
SUBBASIN O.F. LENGTH O.F. SLOPE MANNING S SUBBASIN AREA SCS CIRVE
DATA (ft.) (ft/ft) N AREA (sc. mi.) NUMBER
SUBI
Iip Catchment 100 .03 .24 15 1.43 go
Perv Catchment 190 .02 .35 85 85
SUB 2
Ip Catchmnt 100 .05 .24 15 .67 98
Perv Catchment 190 .03 .35 85 85
SUB 3
Imp Catchment 100 .05 .24 15 .56 98
Perv Catctwent 190 .03 .35 85 85
SUB 4
Imp Catchment 100 .03 .24 15 1.83 98
Perv Catchment 190 .015 .35 85 85
SUB 5
I Catclment 100 .05 .24 20 .67 98
Perv Catchment 220 .028 .35 80 85
SUB 6
Imp Catchment 100 .03 .24 20 1.43 98
Perv Catchment 200 .02 .35 80 85
SUB7
Ip Catchment 100 .06 .24 15 .96 98
Perv Catchnmnt 190 .03 .35 85 85
TABLE 12.2b
CWNEL DATA (Test 2)
(RK RECORD)
107
TABLE 12.2c
Precipitation Data
1 1.96
2 1.68
3 2.73
4 2.56
5 2.52
SUB 1 10 1 1 .75
3 .25
SUB 2 10 1 1 .75
2 .25
SUB 3 20 1 2 1
SUB 4 40 1 2 .05
3 .40
4 .50
5 .05
SUB 5 40 1 2 .2
4 .8
SUB 6 40 1 4 1
SUB 7 50 1 5 1
108
TABLE 12.2d
LINE to ....... 1 ....... 2 ....... 3 ....... 4 ....... ....... 66 ....... 7 ....... S ....... 9 ...... 10
56 KR SulI
57 KJ4 RUNOFF FROM SU.BASf'4 I
so KO 1
59 PR 10
60 tW 1
61 PT 1 3
62 P" .75 .25
63 BA 1.43
LS 98 85
64
65 U1K 100 .03 .24 15
66 UK 190 .02 .35 85
.008 .02 .45 TwAP 0 1
67 RI 2000
.004 .08 TAP 2 2
68 RI 13500
109
aw-1 INPUT PAGE 2
LIN ID ....... 1 ....... 2 ....... 3 ....... 4 ....... 5....... 6 ....... 7 ....... a ....... 9 ...... 10
69 1! SUB2
70 1K RUNOFF FROM SUBBASNI 2
71 PR 10
72 FN 1
73 PT 1 2
74 PW .75 .25
75 BA .67
76 UK 100 .05 .24 15
77 UK 190 .03 .35 85
78 RK 2400 .01 .02 .39 TRAP 0
79 R 6500 .008 .08 TRAP 2 2
80 K1 S0B3
81 Dt RUNOFF FRO SUBBASI- 3
82 BA .56
83 PR 20
84 9w 1
85 PT 2
86 :" 1
87 UK 100 .05 .24 15
88 UK 190 .03 .35 85
89 Rl 1600 .019 .02 .35 TRAP 0
s0 RK 6500 .012 .08 TRAP 2 2
91 K1 SUB4
92 KM COMBINE RU?4OFF FRO SUB1, SUB2 AND SUB3
93 HC 3
94 UK SUB4
95 K! RUNOFF FROM SUBBASIN 4
96 PR 40
97 PW 1
98 PT 2 3 4 5
99 PW .05 .40 S0 .05
100 BA 1.83
101 UK 100 .03 .24 15
102 UK 190 .015 .35 85
103 RK 2500 .01 .02 .79 TRAP 0
104 AN 12000 .007 .05 TRAP 50 2 YES
105 K SUBS
106 KH RUNOFF FRO SUBBASIrN 5
107 PR 40
100 9W 1
109 PT 2 4
110 9W .2 .8
111 BA .67
112 UK 100 .05 .24 20
113 UK 220 .028 .35 80
114 RK 2000 .013 .02 .42 TRAP 0
115 RK 8000 .01 .0S TRAP 3
116 UK SUBS
117 KM RUNOFF FROM SUBDASIN 6
118 PR 40
119 Pw 1
120 PT 4
121 PW 1
122 BA 1.43
123 UK 100 .03 .24 20
124 UK 200 .02 .35 80
125 RK 2200 .011 .02 .55 TRAP 0
126 RK 14000 .005 .09 TRAP 2 2
127 9K SUB7
128 KM COMBINE RUNOFF FROM SUB4, SUDS, AND SUB6
129 HC 3
110
*taaa~o **ttt*a~tt
tttta tawtlatateit~iI tt tttttttt *a
eaaaaaaaaaaa~aaaaa
eaaaae
EXAMPLE
"MUM NO. 2
KI'NU4AflC WAVI WATIEIm b4V0/
UIGLISR UNITS
DRAINAGE ARM SQUARE MILES
PRECIPITATION ORS XNC
LENGTU, ELEVATION PEET
now CUBIC PR1T PER SEDND
STORMEN VOLUR Am-FE
SURFACE AM ACRES
TEMPERATURE DnuRm FAoREuIT
,et tat tat tat at. eta tea tat tat a.. tea a.. eta taa tea ea e eattt te at. tat at, tat tatta tat at. ate Cra. at ee
56EK SUI •
S*eatat t aaa
63 BA SUEBSMIN CRARACTERISTICS
TARIA 1.43 SUBBSALN AREA
111
." UAla ,',l'l " L;A'A
59 PR RECORDING STATIONS 10
60 PW WEIGHTS 1.00
KINEMATIC NAVE
65 UK OVERLAND-FLOW EWEIT NO. 1
L 100. OVERLAND FLOW LENGTE
S 0.0300 SLOPE
N 0.240 MOUGHWBSS COEFICIENT
PA 15.0 PERCIT OF SUBEASIN
66 UK OVERLAIID-FLOW, ELDI NO. 2
L 190. OVERLAND FLO LZM=
8 0.0200 SLOPE
N 0.350 RUG3IESS CIEICIENT
PA 85.0 PERCENT OP SUeASIN
67 RK COLLECTOR CHANNEL
L 2000. CHANNEL LZSTH
S 0.0080 SLOPE
N 0.020 CHANNEL ROUONNU COEFFICIfNT
CA 0.45 CONTRIBUTING AE
SHAPE TRAP CHANNEL SHAPE
NO 0.00 00TroN WIDTH OR DIAMETER
2 1.00 SIDE SLOPE
68 RI MAIN CHANNL,
L 13500. CHANNEL LDEGTH
S 0,0040 SLOPE
N 0.080 CHANNEL ROUGHNESS COEFFICIENT
CA 1.43 CONTRItNBUTN AREA
SHAPE TRAP CHANNEL SHAPE
NO 2.00 907" Wti)m OR DIAMETER
$ 2.00 SIDE SLOPE
RUPSTQ NO ROUTE UPSTREAM HY'DROGRAM
112
HYDROGRAP AT STATION SUBI
DA NON HPRN ORD RAIN LOSS EXCESS COP Q * DA MN EM ORD RAIN LOSS EXCESS COMP Q
RUNOFF SIDHARy
FLOW IN CUBIC FEET PER SECOND
TIME IN HOURS, AREA IN SQUARE MILES
PEAK TIME OF AVERAGE FLOW FOR MAXIUM PERIOD BASIN MAXIMUM TONE OF
OPERATION STATION PLO; PEAK 6-HOUR 24-HOUR 72-HOUR AREA STAGE MAX STAGE
113
12.3 Example Problem #3: Snowmelt Runoff Simulation
The input data and results of the analysis are displayed in the computer
printout in Table 12.3b.
114
TABLE 12.3a
Snowmelt Data
INFILTRATION CARDS
SJ1R - 0. 24
DLTKR - 0.00
RTIOL - 1.00
ERAIN - 0.70
S7RKS - 0.24
RTIOK - 1.00
UNIT HYDROGRAPH UC
TC - 46
R - 183
ZOE DATA NA
1 1,000 7.5
2 500 6.2
3 370 8.4
MELT COEFFICIENTS MC
TLAPS - 3.3
COEF - .08
FRZTP - 33
TABLE 12.3b
Example Problem #3: Input and Output
H&C-i INPUT PAGE
LINE ID ....... I ....... 2 ....... 3 ....... 4 ....... 5 ....... 6 ....... 7 ....... 6 ....... 9 ...... 10
I ID EXAMPLE PROSLM 1. 3
2 10 SNUOIULT RIOF SIMULATION
12 IX 7
13 K ,MINWNESOTA RIVER BASIN
14 BA 1870
15 or 8 1500 1.0022
16 PT 100
17 PM 1.0
18 PR 100
19 P" 1.0
20 tc 46 183
21 .E .24 0 1.0 .7
22 1M .24 1.0
* * * 11 DATA (L(MEST ZONE)
23 KA 1000 7.5
115
* ..... DATA FOR ZOtEs AT HIGHER ELEVATIONS (1000 rT INCREMENTS)
24 MA% 500 6.2
25 MA 370 8.4
26 BC 3.3 .08 33
27 IN 1440 04APR75 0800
28 HT 18 30 35 31 27 22 32 14 0 2
29 HT 17 37 28 37 38 34 37 48 51 47
30 MT 42 45 55 60 54 53 52 47 45 50
31 MT 55 51 50 49 50 60 55 50 41 46
32 MT 54 57 57 54 64 65 63 58 52 47
33 zz
E4GLISH UNITS
DRAINAGE AREA SQUARE MILES
PRECIPITATION DEPTH INCHES
LENGTH, ELEVATION FEET
FLOW CUBIC FEET PER SECOND
STORAGE VOLUME ACRE-FEET
SURFACE AREA ACRES
TEMPERATURE DEGREES FHREIT
*• 444
* * * 44 *t *4 • • e •eeeJ * t ••••tee•o et e e ,••tte e o * • e
* *
12 KK 7 *
* *
116
IUDA ItMM I II
14 DA BUNAiU CmA~aCTs*ISTICS
TAM 1670.00 SWhWNIAWAMNA
15 w BUZ3 V CNAECTRIXSTICU
fmg 6.00 INITIAL FLOW
OF"N 1500.00 srxB EASE FAI/ EM8IO
ITIOr 1.00220 1A[BSION CONiTAX
U PIT DATA
r'JIO"
IC 81K)MLT DATA
TLApS 3.30 T13MRAT LAPSE PATE
C0lp 0.06 SWONMT CWOTICIUIT
P'11? 33.00 MT T20MRATURZ
IT TUWZRAIURZ DATA
16.0 24.0 30.0 32.5 35.0 33.0 31.0 29.0 27.0 24.5
22.0 27.0 32.0 23.0 14.0 7.0 0.0 1.0 2.0 9.5
17.0 27.0 37.0 32.5 26.0 32.5 37.0 37.5 38.0 36.0
34.0 3S.5 37.0 42.5 66.0 49.5 51.0 49.0 47.0 44.5
42.0 43.5 45.0 50.0 55.0 57.5 60.0 57.0 54.0 53.5
53.0 52.5 52.0 49.5 47.0 46.0 45.0 47.S SO.0 52.5
55.0 53.0 51.0 50.5 50.0 49.5 49.0 49.5 SO.0 55.0
60.0 57.5 55.0 52.5 50.0 45.5 41.0 43.5 46.0 50.0
54.0 55.5 57.0 57.0 57.0 55.5 54.0 59.0 64.0 64.5
20 UC CLRK UNIGPAPU
TC 46.00 TIM OF CONCNfTRATION
R 183.00 STORAGE QOEICIEIT
ISSORAL DISTRIBUTIONS
117
UNIT HYDROGRAPH
83 I"D-OF-PRIOD ORDINATES
602. 2263. 4252. 5476. 5586. 5231. 4899. 4588. 4296. 4024.
3768. 3529. 3305. 3095. 2899. 2715. 2542. 2381. 2230. 2088.
1955. 1831. 1715. 1606. 1504. 1409. 1319. 1235. 1157. 1084.
1015. 950. 890. 833. 781. 731. 685. 641. 600. 562.
527. 493. 462. 433. 405. 379. 355. 333. 312. 292.
273. 256. 240. 224. 210. 197. 184. 173. 162. 151.
142. 133. 124. 116. 109. 102. 96. 90. 84. 79.
74. 69. 65. 60. 57. 53. 50. 46. 44. 41.
38. 36. 33.
HYDR0GRAPH AT STATION 7
DA MON HR84N ORD PRECIP TOM SHcEELT SOLOSS SWOEXCS RAIN RAINLOS RAIMEXS SM0+RAIN LOSS EXCESS COMP Q
4 APR 0800 1 0.00 16.0 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 8.
4 APR 2000 2 0.00 24.0 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 8.
118
STATION 7
(0) OUTFLOW
0. 2000. 4000. 6000. $000. 10000. 12000. 14000. 16000. O. 0. 0. 0.
(T) TMIERATURE
O. 0. 0. 0. 0. 0. 0. 20. 40. 60. 80. 0. 0.
(L) PRECIP, (X) EXCESS
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.2 0.8 0.4 0.0
'.qM PER
42000 20 T
50800 30 T
52000 40 T
60800 so T
62000 60 T
70800 70 T
72000 00 T
119
12.4 Example Problem #4: Unit Graph and Loss Rate Parameter Optimization
TABLE 12.4
LINER 2 . 3 . - .
IN.1 5 . . 7 . S 9 . 1
LIMB 10 ....... I ....... 2 ....... 3 ....... 4 ....... 5 ....... 6 ....... 7 ....... I ....... 0 ...... 10
1 ID ZXMLlE TEST 0. 4
2 D UNIT GRAM AM 1,4058 RATE OPTIXUZATIOU
3 IT 15 67k]G27 1145 61
4 10 1 2
5 O0
6 PC 467042 2.39 1.00
7 PG 100
a ix .00 .00 .03 .06 .45 .42 .29 .14 .08 .04
9 PI .03 . .0 .02 .02 .01 .01 .01 .01 .01 .02
10 PI .01 .01 .02 .01 .01 .01 .01 .01 .01 .00
11 Pr .01 .01 .00 .00 .01 .02 .01 .01 .01 .00
12 PI .01 .01 .01 .01 .00 .00 .00 .01 .00 .00
13 PG 300
14 Pi .00 .00 .00 .00 .00 .00 .00 .00 .00 .00
15 P1 .00 .00 .00 .10 .45 1.45 .73 .02 .80 .50
16 PI .25 .05 .00 .00 .00 .00 .00 .00 .00 .00
17 PZ .00 .00 .00 .00 .00 .00 .00 .00 .00 .00
18 PI .00 .00 .00 .00 .00 .00 .00 .00 .00 .00
19 PG 5000
20 P1 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00
21 P2 .00 .00 .00 .00 .00 .00 .00 .00 .00 .00
22 PI .00 .00 .00 .00 .00 .00 .00 .00 .04 .23
23 PI .39 .18 .56 .00 .00 .00 .19 .06 .20 .20
24 P1 .11 .03 .00 .00 .00 .00 .00 .00 .00 .00
25 KR 467042
26 go 57 57 59 61 63 65 67 69 71 73
27 QO 130 250 370 520 720 920 1170 1470 1720 1900
28 QO 2060 2250 2400 2570 2720 2860 3090 3390 3540 3520
29 QO 3480 3330 3290 3230 3100 2900 2720 2520 2270 2050
30 QO 1800 1570 1430 1300 1200 1100 980 890 800 745
31 QO 690 650 610 570 540 510 490 475 460 445
32 90 430 0 0 0 0 0 0 0 0 0
33 PT 467042
34 PE 1.00 0. 0. 0. 0.
35 PR 100 300 5000 0 0
36 13 .45 .45 .10 0. 0.
37 BA 37.90 0.-
30 9F 57. -.25 1.3195
39 OC -1.00 -1.00
40 La -1. -1. 1. .5
41 Z2
120
ft * * ft
DIGLIS UNITS
DRAINAGE AREA SQUARE NILES
PRECIPITATION DEPTB INCHES
LENGTH ELEVATION TEST
FLOW CUBIC FEET PR
STORAGE VOLUME ACl-FnT
SUPAC AREA ACRES
TZPERR DEGREES FAHRENHIT
25 X 467042
37 BA SUBAMIN CRARACTERISTICS
TARRA 37.90 SUlBSIN AEA
38 BY BASS FLOW CEARACTERISTICS
STRITO 57.00 INITIAL FLOW
121
qiotm -0.1) OmN BASE FW tMMSION
31100 1.31950 RECESSION CONSTANT
PRlIPI7ATIOW DATA
39 OC CLARK UNITG3AI3
IC -1.00 TIM OF CONCRTRATION
3 -1.00 STORAGZ COEFICIENT
TJ3KPAL DISTRIBUTIONS
122
3.10123402 7.101* 0.551 0.452 0.542 1.000 0.500
3.10123402 7.101 0.SS1 0.42 0.542 1.000 0.500
•
3.05772402 7.101 0.551 0.465 0.542 1.000 0.500
2.9360Z+02 7.101 0.551 0.465 0.36 2* 1.000 0.500
* OPTIIZATION 3ESLTS
................................................:
~eeat~
.*S............. t...t. ta t *ttta*aatte~ e a. ~..t.taa. t.t0.att.Ot
* TIM TO LAG
• SON Or BOv MEAN CENTE C.K. TO PEAS TIME Or
• FLOWS DETH FLOW OF "An C.,. FLOW PEAK
123
UIT YDROGRAPE
89 ZWD-OF-PZIMOO ORDINATES
96. 361. 741. 1191. 1690. 2227. 2779. 3285. 3700. 4017.
4232. 4332. 4265. 4050. 3797. 3560. 3338. 3130. 2935. 2752.
2580. 2419. 2268. 2127. 1994. 1870. 1753. 1644. 1541. 1445.
1355. 1270. 1191. 1117. 1047. 982. 920. 863. 809. 759.
711. 667. 625. 586. 550. 515. 483. 453. 425. 398.
374. 350. 328. 308. 289. 271. 254. 238. 223. 209.
196. 184. 172. 162. 152. 142. 133. 125. 117. 110.
103. 97. 91. 85. 80. 75. 70. 66. 62. 58.
54. 51. 48. 45. 42. 39. 37. 34. 32.
1 SEP 1145 1 0.00 0.00 0.00 57. 57. * 1 SEP 1930 32 0.03 0.03 0.00 3311. 333'
1 SEP 1200 2 0.00 0.00 0.00 53. 57. 1 SEP 1945 33 0.02 0.02 0.00 3135. 329-
1 SEP 1215 3 0.00 0.00 0.00 50. 59. * 1 SEP 2000 34 0.04 0.04 0.00 2947. 322'
1 SEP 1230 4 0.01 0.01 0.00 46. 61. * 1 SEP 2015 35 0.00 0.00 0.00 2764. 317'
1 SEP 1245 5 0.02 0.02 0.00 43. 63. * 1 SEP 2030 36 0.00 0.00 0.00 2592. 293
1 SEP 1300 6 0.16 0.10 0.06 46. 65. * 1 SEP 2045 37 0.01 0.01 0.00 2430. 27:
1 SEP 1315 7 0.15 0.09 0.05 64. 67. * 1 SEP 2100 38 0.02 0.02 0.00 2276. 2t2
1 SEP 1330 8 0.10 0.08 0.02 100. 69. * 1 SEP 2119 39 0.01 0.01 0.00 2136. 227,
1 SEP 1345 9 0.0S 0.05 0.00 151. 71. * 1 SEP 2130 40 0.02 0.02 0.00 2003. 205
1 SEP 1400 10 0.03 0.03 0.00 212. 73. 1 SEP 2145 41 0.02 0.02 0.00 1878. 18C:
1 SEP 1415 13 0.01 0.01 0.00 280. 130. * 1 SEP 2200 42 0.01 0.01 0.00 1761. 15:.
1 SEP 1430 12 0.01 0.01 0.00 352. 250. * L SEP 2215 43 0.01 0.01 0.00 1651. 143,
1 SEP 1445 13 0.01 0.01 0.00 423. 370. 1 SEP 2230 44 0.00 0.00 0.00 1546. 130G
1 SEP 1500 14 0.01 0.01 0.00 487. 520. * 1 SEP 2245 45 0.00 0.00 0.00 1451. 121"
1 SEP 1515 15 0.04 0.04 0.00 539. 720. * 1 SEP 2300 46 0.00 0.00 0.00 1361. 11-
1 SEP 1530 16 0.16 0.10 0.06 584. 920. * 1 SEP 2315 47 0.00 0.00 0.00 1276. 98
1 SEP 1545 17 0.52 0.18 0.34 658. 1170. * 1 SEP 2330 48 0.00 0.00 0.00 1196. 89:
1 SEP 1600 18 0.26 0.12 0.14 792. 1470. 1 SEP 2345 49 0.00 0.00 0.00 1122. 81
1 SEP 1615 19 0.01 0.01 0.00 973. 1720. * 2 SEP 0000 50 0.00 0.00 0.00 1052. 74
1 SEP 1630 20 0.29 0.13 0.16 1198. 1900. * 2 SEP 0015 51 0.00 0.00 0.00 986. 6?
1 SEP 1645 21 0.18 0.10 0.08 1481. 2060. * 2 SEP 0030 52 0.00 0.00 0.00 925. 65-
1 SEP 1700 22 0.09 0.07 0.02 1818. 2250. * 2 SEP 0045 53 0.00 0.00 0.00 867. 6,.
i ;zLr ill ;.U2 0.02 O.UU 2MO. 24"~. L
WID A
S4 i I. .;U ULS. Diu.
I SEP 1730 24 0.01 0.01 0.00 2557. 2570. * 2 SEP 0115 55 0.00 0.00 0.00 762. 540.
'SEP 1745 25 0.00 0.00 0.00 2894. 2720. * 2 SEP 0130 56 0.00 0.00 0.00 714. 510.
1 SEP 1800 26 0.00 0.00 0.00 3187. 2660. * 2 SEP 0145 57 0.00 0.00 0.00 670. 490.
1 SEP 1815 27 0.00 0.00 0.00 3420. 3090. * 2 SEP 0200 58 0.00 0.00 0.00 628. 475.
1 SEP 1830 28 0.00 0.00 0.00 3573. 3390. * 2 SEP 0215 59 0.00 0.00 0.00 589. 460.
1 SEP 1845 29 0.00 0.00 0.00 3621. 3540. * 2 SEP 0230 60 0.00 0.00 0.00 552. 445.
1 SEP 1900 30 0.01 0.01 0.00 3568. 3520. * 2 SEP 0245 61 0.00 0.00 0.00 518. 430.
1 SEP 1915 31 0.02 0.02 0.00 3457. 3480. *
124
.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . .
*
.3 3
."
04. I
0 . .
o*I 0.0 0 . . . . .
O• •
.. ; . . . . . ... . . . . . . ..
0
* 0
"0
* * 0
o
. . 0 . .
S. . . . .
- . ..
l • "0
o • I
.0 • • • • • • • • • ° • • • • • • • • • • • • • • . • , . 0 0
0 4.4C Mo M
. . . .. . . . . . . 0
.. .. . . . . . . . . . . . . . . .. .. 00..
0. . . . . . . . . . . . . .. . . . . . .
o
"°
o.
OATE
h .. PZM f ERRAR.... . . .I.I.... 0.. Tl... T. 1.0...1
OU S T 8.....................
DA MEN R V VOL LAG ER IC A T. R/(ITC+R) T1 Cl QFP STAEA DLTKR RTIOL EAJ
67 AUG 27 27.3 0.0 8.7 2.3 3.16 3.88 7.0S 0.55 2.99 0.52 111. 0.49 0.14 1.00 0.50
125
12.5 Example Problem #5: Routing Parameter Optimization
TABLE 12.5
1 D EXAMPLE PROBLDM N. S
2 ID STREANELW ROUTING OPTIMIZATIOM
3 ID MUSKINGUM METHOD
4 IT 720 600000 0 is
5 10 1 2
6 OR 2
7K 1
S QF 2000 2000 7000 11700 16500 24000 29100 28400 23800 19400
9 QP 15300 11200 8200 6400 5200 4600 0 0 0 0
10 QI. 2200 2200 14500 28400 31800 29700 25300 20400 16300 12600
11 QI 9300 6700 5000 4100 3600 2400 0 0 0 0
12 QO 2000 2000 7000 11700 16500 24000 29100 28400 23800 19400
13 QO 15300 11200 8200 6400 5200 4600 0 0 0 0
14 RL 0. 0.
15 RN -1 -1.00 -1.00
16 21
ENGLISH UNITS
DRAINAGE AREA SQUARE MILES
PRECIPITATION DEPTH INCHES
LENGTH. ELEVATION FEET
FLOw CUBIC FEET PER SECOND
STORAGE VOLUME ACR-FrEET
SURFACE ARA ACRES
TEqPERATURE DEGREES FAHRENHEIT
126
KYDROGNAP8 ROUING DATA
14 RL ROUTING as83
is IN N6KIWA ADUTING
I -1.00 HSKINUN I
127
DERIVED COEFFICIENTS
NSTPS NSTDL LAG ANSKK X TSK
2 0 0 21.61 0.04 0.00
STATION 1
(I IW IE. (0) OUTFLOW, (*) OBSRVRD FLOW
0. 4000. 6600. 12000 16000. 28000. 32000. 20000.
0. 24000.
0. 0. 0.
60000 P3*1 .
DAIR
61200 2. . .* 1- . . . . . ........
-.......... . . . ..................................................................................
. . . . . . . . . . . . . . . . . . . . . . . . . . .
70000 3. 0. *
71200 4. o * .1 .
90000 5 .0
81200 6. •0
90000 7 . 0
91200 1. . 0 .* •
100000 9.• 0* •
101200 10. *
110000 11.......... ..... . ... .. e.. ...... ...........................
111200 12 1 . .0. ..
120000 13. 1 . 0
121200 14. 1 0
131200 16.- - -x--.o*,* - --- . . . -. .. . ----.
130000 15. 1. . . --
. ...
--..
--.-.--
.. ..- -. --
........
-- .........
- -- - -- -..
- .. ... .. .
(-) O OPTIMIZATION
128
12.6 Example Problem #6: Precipitation Depth-Area Simulation
TABLE 12.6a
Depth-Area Simulation Data
TRRNSPOSITIaN STORM JD
AREA DEPTH
(sq MI) (in)
1000 9.08
3000 8.93
5000 8.70
7000 8.57
9000 8.43
Please see data input listing for Dattern hyetograph (PI cards)
4606 4306
129
TABLE 12.6b
130
1 3D T STOM NO. 2
STEW 8.93 PRECIPIAION DIM
TWA 3000.00 TRUASPOSITION OPAXSZ AM
a PI PRIPIATIOM PATTm
0.00 0.00 0.00 0.00 0.00 0.00 0.01 0.00 0.00 0.01
0.01 0.01 0.02 0.02 0.02 0.01 0.01 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.01 0.02 0.02 0.02 0.00 0.01
0.00 0.00 0.01 0.00 0.01 0.02 0.00 0.00 0.03 0.00
0.00 0.02 0.02 0.02 0.02 0.01 0.01 0.02 0.01 0.01
0.01 0.01 0.02 0.01 0.01 0.01 0.01 0.01 0.01 0.02
0.02 0.02 0.02 0.03 0.03 0.07 0.02 0.02 0.02 0.02
0.01 0.01 0.01 0.01 0.00 0.01 0.00 0.00 0.00 0.00
0.01 0.00 0.01 0.03 0.00 0.03 0.01 0.00 0.00 0.01
19 JD INDEX STORM NO. 3
STUN 8.70 PRCIPITATION OE
TLUA 5000.00 TUASPOSI'TION DRAINA AMEA
0 Pl P3WIPTATION PATTIM
0.00 0.00 0.00 0.00 0.00 0.00 0.01 0.00 0.00o 0.01
01
0.00
001
0.0
0.02
0.00
0.02
0
S: 8:8
:2 8:81 8:8 8:8 0 :?
0.00 0.00 0.01 0.00 0.01 0.02 0.00 0.00 0.03 0.00
0.00 0.02 0.02 0.02 0.02 0.01 0.01 0.02 0.01 0.01
0.01 0.01 0.02 0.01 0.01 0.01 0.01 0.01 0.01 0.02
0.02 0.02 0.02 0.03 0.03 0.07 0.02 0.02 0.02 0.02
0.01 O.Oi . .6i U. 0.5i
.G 0.00 0.00 .. UU
0.01 0.00 0.01 0.03 0.00 0.03 0.01 0.00 0.00 0.01
22KB * 486S 0
131
DTUG3AM AT STATION 4865
T1MESIosTIOm AM 1000.0 Q INI
CUILATIVE AM - 3503.00 SQ NI
CUMULATIVE AR - 3503.00 SQ KI
RUDOFF SUMMARY
FLOW IN CUBIC PET PER SEoD
TIME IN HOURS, AREA IN SQUARE MILES
PEAR TIlE OF AVERAGE FLOW FOR MAXIMUM PERIOD BASIN NAXIMUN TIME OF
OPERATION STATION FLOM PEAK 6-DOER 24-WiuR 72-OUR AREA STAGE MAR STAGE
132
12.7 Example Problem #7: Dam Safety Analysis
Problem Description
Develop the PMP for area above the reservoir from input index
rainfall parameters.
Determine average basin loss rates and probable maximum
rainfall excess.
Develop a unit hydrograph using the Snyder method.
Generate the runoff hydrograph and add base flow to get
probable maximum inflow hydrograph to the reservoir.
Apply ratios to the PMF to obtain a series of proportional
inflow hydrographs.
Develop reservoir storage-outflow functions from elevation-area
relationship and characteristics of reservoir outlet works and
dam.
Route hydrographs through the reservoir and determine the ratio
of the PMP that overtops the dam.
The input data and output from the HEC-I program are shown in
Table 12.7b.
Discussion of Results
The last page of the HEC-I output (Table 12.7b) provides a "SUMMARY
OF DAM OVERTOPPING/BREACH ANALYSIS FOR STATION DAM" which illustrates the
potential of the dam to overtop as a ratio of the PMF. Also, data on
duration of overtopping and maximum water surface elevations, for use in
determining possible dam failure due to erosion are shown. Interpolation of
the information provided in that summary indicates that a flood of about
thirty percent of the PMF would overtop the dam.
133
BEAR GREEK RESERVOIR
*reale Area at
114148te-72.1 eq. 01.
Rootling Itch I1
Rooting MC1M a
Istaftee 12000 feet
aStop* .0024 Wt/t
134
TABLE 12.7a
Reservoir Data
RECORD IDENTIFIER
Outflow characteristics of the Bear Creek Reservoir:
Diameter = 4 feet
Coefficient of discharge = .7
Downstream centerline elevation of outlet = 380.0 m.s.l.
Exponent of head = .5
Spillway SS
Dam ST
Elevation Area
(m.s.l.) (acres)
340 0
380 100
410 250
420 300
424 320
428 350
432 380
436 410
440 450
444 500
135
TABLE 12.7b
Example Problem #7: Input and Output
1-1 ism PAGE 1
LIIZ ID ....... I....... 2....... 3 ....... 4 ....... 5 ....... 6 ....... 7....... ....... 110
......
1 ID XUSNA PMOBLXN M. 7
2 ID DAN SAFEY ANALYSIS
3 IO ANALYSIS OF DAN OVERTOPPING USING PATIOS OF He
4 IT 15 0 0 121
5 1 5
6 31 FLOW .20 .35 .50 .65 .80 1.0
7 I DIFLO IFLO TO WAR RESVOIR
IA 72.1
9 ip -1.0 -. 05 2.0
10 PM 25 0 0 02 97 110
11 LU 1.0 .04
12 us 4.8 .60
13 xx DAM BEAR CRUX DAN
14 Rs 1 STOP 10000
15 8A 0 100 250 300 320 350 380 410 450 S00
16 BE 340 380 410 420 424 420 432 436 440 444
17 as 420 200 3.1 1.5
18 IL 380 12.6 .7 .5
19 BT 432 900 3.1 1.5
20 U a1
21 K1 RBMZ OUTFLOWM RS IMaR
T00 FIRST CHNNL MRZC DONST EAM
22 RS 1 STOR -1
23 IN .04 .05 .04 15000 .0033 312
24 Ix 0 S00 1400 1425 1450 1475 2500 3000
25 RY 400 350 290 280 280 290 350 400
26 S2
ANWPOLPOLM no. 7
DAMSHAM ANALYBS
ANALYSIS OF DANOVERTPPING USING RATIOS OF N
5 10 OUTPUT CONTROL VARIABLES
IPIMT 5 PRINT CONTROL
IPLOT 0 PLOT ON1 ROL
QSCAL 0. TDGIRAP1 PLT SCALE
ONIG YES PRINT DIAGNOSTIC SAOAUB
IT HYDROGRAPHTIC DATA
115116 15 MNUTES IN COMPUTTION rNTZRVAL
IDATE 1 0 STARTING DATS
INE 0000 STAiIG ThU
no 121 NUMBR OP RTDF40RAPS ORDINATZS
NDOATE 2 0 ESIN0 DATE
NOTIIM 0600 EI31NG TIMR
COMUTATION INTERVAL 0.25 ROURS
TOTKL TIME SARK 30.00 HOURS
ENGLISH UNITS
DRAINAGE ARER SQUARE KILES
PRIPITATION DWPTH
LIN". ZT
MATION FE
FLOW CUBIC FEET PER SlN
STOAGE VOLUME A=R-P?
SURE AM ARE
TEMPEATURE ONORMS FARIUISRIT
3P ISDTI-PIAN OPTION
NPtAN 1 NUMBER O PLANS
3R NULTI-RATIO OPTION
RATIOS 07 M" rF
0.20 0.35 0.50 0.65 0.80 1.00
136
FORM4ULTIPLEPLAN-RATIO ECONOMIC
PEAKFLOW AND STAGE(END-OF-PERIOD)SUMMARY COMPUTATIONS
FLOWSIN CUBIC FEETPER SECOND, AREA IN SUJAREMILES
TIMETO PEAK IN HOURS
HYDROGRAPH AT
+ INFLOU 72.10 1 FLOW 19032. 33305. 47579. 61853, 76126. 95158.
lNK 2c.?- 20.25 .0.25 20.25 20.25 20.25
ROUTEO TO
* DAN 72.10 1 FLOu 17040. 31622. 46856. 61272. 75545. 94537.
TIME 21.50 21.00 20.75 20.50 20.50 20.50
ROUTED TO
* RCHI 72.10 1 FLOW 16886. 31224. 46515. 60842. 75115. 94074.
TIE 22.00 21.50 21.00 20.7 20.75 20.75
137
12.8 Example Problem #8: Dam Failure Analysis
Discussitn of Results
The plots illustrate how well the hydrograph depicted by normal time
steps represents the breach hydrograph generated using smaller time steps.
PLAN I has a volume gain of 2330 acre-feet from the peak portion of the
hydrograph indicating that a smaller time step should be used. The plot for
PLAN 3 indicates that the peak flow from the dam occurs after thu breach is
fully formed. Characterization of the outflow hydrograph and peak discharge
will depend on the specified time step as in a standard storage routing.
138
TYI CAL CRVOSS SECTION 1
mooch me" I
1000' 2t00',
1251 050ms,
Figure~~8
126BerCeek Don1rem0rssSetin
TYICL nO139CIO
TABLE 12.8a
Dam Failure Analysis Results
Plan No. Reservoir RCHI RCH2
and Peak Peak Peak Peak Peak Peak
Breach Criteria %W.S. El. Q W.S. El. Q W.S. El.
Plan 1
fail time - 15 sin. 1,244,000 433.5 610,000 334.1 422,000 280.7
total dam
Plan 2
fail time - 3 hrs. 209.000 434.1 197,000 317.7 184,000 276.0
total dam
Plan 3
fail time - 3 hra 135,000 435.4 127,000 312.9 122,000 274.2
breach depth - 50 ft.
b.w. - 50 ft.
@.a. - 2:1
Plan 4
fail time - 3 hr. 180,000 434.6 175,000 316.3 171,000 275.6
breach depth - 70 ft.
b.v. - 200 ft.
s.a. - 2:1
Plan 5
fail time - 10 hra 109,000 436.6 109,000 311.4 108,000 273.7
breach depth = 70 ft.
b.v. - 200 ft.
a.s. - 2:1
TABLE 12.8b
Example Problem #8: Input and Output
L=C- I23 .79 PAGE 1
L~INE ID ....... I ....... 2 ....... 3 ....... 4 ....... 5 ....... 6 ....... 7 ....... a ....... 9 ...... 10
1 ID MI P LEM NO. 8
2 ID DAM PFILRN ANALYSIS
e FIX ***
3 IT is 140
4 10 4
5 JP 5
6 VS Will 1 SMI H2 ITCR2 K2 BCH2
7 W 2.11 2.51 7.11 7.51 2.11 2.51 7.11 7.51
140
35 a aE1
36 KO 5
37 IN COMM R0TING 31UM 2-3
38 Its 1 SJOR 0
39 MC .04 .05 .04 ISOO .0033 33S
40 BE 0 500 1400 1425 1450 1475 2S00 3000
41 BY 400 350 290 280 200 290 350 400
42 [ 3l2
43 10 5
44 UO CAME ROUTING ZMIC 3-4
45 RS I STOR 0
46 PC .03 .04 .03 12000 .0025 200
47 RX 0 500 1000 10S0 1150 1200 3S00 5000
48 RY 400 275 26S 255 255 26S 275 400
49 Zz
ENGLISH UNITS
DRAINAGE AM SQUARE MILES
PRECIPITATION DiE nius
tasT, .EVATION FET
FLOW CUBIC FEET PER SECOND
STORM!GRVOLUME ACRE-FEET
SURFACE ARM ACRES
TEMPERATURE DEGEE FAHRENEIT
U ER-DiINED OUTPUT SPECIPICATIONS
TABLE 1
VS STATION RCH1 Nal ami CRCi RCH2 K22 MC12 MC82
W VARIABLE ODE 2.11 2.51 7.11 7.51 2.31 2.51 7.11 7.51 0.00 0.00
JP IIILTI-PIAN OPTION
UPIAN 5 NUMBER OF PLANS
JR NULTI-RATIO OPTION
RATIOs OF RUOFF
1.*00
e* *9* 9*9 9*9 99* * *9* 9* *e999 *9 e e*9 9*9 *99te*9* *999* *9*ee*9*9* 99* e* ,99* 99* *9* *** *9* 9ee *•* e* 99* 9*
141
CALCUATION OF INFLOW TO BEAR CRACK RESEVSIR
10 BA SUBBASrN CHARACTERISTICS
TAM 72.10 SUASrSN AREA
PRECIPITATION DATA
11 P PODABLE MAXIMN ST RM
PSI 25.00 INDEX PRECIPITATION
TRSPC 0.86 TRANSPOSITION CONICIEIT
TRSDA 72.10 TRANSPOSITION ARIA
SWD SO OSX SWD DISTRIBUTION
13 US SNYDER UNITORAPE
TP 4.80 LAG
CP 0.60 PEAlING COEFICIENT
UNIT SYDROGRAPH
115 ZRD-O-PERIO0 ORDINATES
70. 265. 547. 883. 1259. 1667. 2099. 2552. 3020. 3501.
3984. 4437. 4834. 5174. 5457. 5683. 5851. 5958. 60C, 5963.
5810. 5554. 5277. 5013. 4762. 4524. 4298. 4083. 3879. 3685.
3501. 3326. 3160. 3002. 2851. 2709. 2574. 2445. 2323. 2207.
2096. 1991. 1892. 1797. 1707. 1622. 1541. 1464. 1391. 1321.
1255. 1192. 1133. 1076. 1022. 971. 923. 877. 833. 791.
752. 714. 678. 644. 612. 582. 552. 525. 499. 474.
450. 428. 406. 386. 367. 348. 331. 314. 299. 284.
269. 256. 243. 231. 219. 209. 190. 188. 179. 170.
161. 153. 146. 138. 131. 125. 119. 113. 107. 102.
97. 92. 87. 83. 79. 75. 71. 67. 64. 61.
58. 55. 52. 50. 47.
*** CCC 0Cctt *C0 CCC **C CCC dCe Ccc *** *** CCC CCC CCC CCC
CCC CCC tee 0CC Ccc ccc ccc ene ti ee rCC CCC en *ee cc eee eec
eo ee **c ***
ee *** *** e**
*** *** *** *e* ** *** *0
ccc en* eec Ccc *c* CCC ehe e eec eec tee *n eec Ccc ccc
•** * Ccc C eec ene eec ccc eec tee eec eec *c cc cc ere ***
e * ***
ec eec eec to*
ee ***
*t* CCC *to ec ccc c.e et
142
* t
Is8 * OT 3m C
a Kmb
18 I P 1 P*
FAN TATION OUT am CHEC DAN
NI0Ymom WIUING DATA
i9 11 8aOlUd10
PAGE no
momP 1 NUaE oF SOSREACBM
ITYlP RIM TrM ( INITIAL -bDITIOi
MoVRC 420.00 INITIAL ' DXTION
1 0.00 WORKING R AND D COEFFICIENT
20 fA ARIA 0.0 100.0 250.0 300.0 320.0 350.0 300.0 410.0 450.0 S00.0
213 H wIY/A IO 340.00 380.00 410.00 420.00 424.00 428.00 432.00 436.00 440.00 444.00
24 8L LAO-Lirm OULET
E.I, 380.00 ELVATION AT CEINTE OF OUTLET
CAM 12.60 CS-SBCWTONAL AMRA
0.70 COEFICIMT
O.S0 Iapomm 0F READ
22 98 SPILMMAY
cm 420.00 8PILLMAY CR=T ALVATION
813TD 200.00 SPILIMAY WIDTH
ODR 3.10 33I3 COEFICIENT
1.50 R
a" Or EAD
23 38 TOP or DN
Tom 432.00 X, IVATION AT TOP UP DAN
DAWInD 900.00 DAN W1DTH
Q3.10 MIR COEFICINT
1.50 110001101T OF EAD
25 HE BUACB DATA
w.i 340.00 33VATION AT 907109 OP BREAC
BIM[D 900.00 WIDTH OF IREACI 901701
z 0.00 RSrACBBIDS SLOPE
TAIL 0.25 TITll FOR BRlAeI TQ oVELOP
FAILS, 433.00 W.. E,3 I TO PIGGER FAILURE
STORAGE-ELEVATION
COMPUTED DATA
STORAGE 0.00 1333.33 6414.47 9160.68 10400.46 11740.01 13199.60 14779.2 16498.60 18397.72
ELEVATION 340.00 380.00 10.00 420.00 424.00 428.00 432.00 436.00 440.00 444.00
COMPUTED DATA
OUTFLORI-ELEVATION
(EXCLUD)ING OVERDAN)
FLOWJ
OUTFLOW 0.00 0.00 102.19 114.85 131.09 t52.68 182.78 227-66 301.76 417.38
ELEVATION 340.00 380.00 382.09 382.64 383.43 384.66 386.68 390.36 398.20 420.00
OUTFLOW 524.09 1044.57 2444.54 5159.65 9625.53 16277.80 25552.06 37W3.90 53708.9 T3462.72
ELEVATION 420.25 420.97 422.17 423.85 426.01 428.65 431.77 435.37 439.45 444,.00
(INCLUING FLOWOVERDAN)
STORAGE 0.00 1333.33 1550.59 1610.64 1700.12 1842.50 2091.06 2590.89 3870.13 6414.47
OUTFLOU 0.00 0.00 102.19 114.85 131.09 152.68 182.78 227.66 301.76 387.44
ELEVATION 340.00 380.00 382.09 382.64 383.43 384.66 386.68 390.36 398.20 410.00
STORAGE 9160.68 9234.42 9453.83 9824.07 10353.84 10400.46 11060.17 11740.01 11970.54 13113.30
OUTFLOW 447.38 524.09 1044.57 2444.54 5159.65 5429.21 9625.53 14519.08 16277.80 25552.06
ELEVATION 420.00 420.25 420.97 422.17 423.85 424.00 426.01 428.00 428.65 431.77
143
fYD0OGRAPB AT STATION Our
FW 1. PATIO - 1.00
144
THE DAMBREACH HYOORAPS WS OEVEWPED USING A TIME INTERVAL OF 0.*005 HOURSWRING BREACH FORMATION.
OWNSTREAM CALCULATIOSNS WILL USE A TIME INTERVAL OF 0.250 SOURS.
THIS TABLE GOPARMS ?HE H1DNGAPH FOR DrAESTHEAK CALCULATIONS WITH THE COMPUTED BRERCH HTI3NOGAl.
INfTREIATE FLOWS ARE INITERPOLATED FROM DUD-0P-FERRIQ VALOM.
Tim nHom IVFRPOLATIND COMPUTED
TimE EGINING BREACH - BREACH Z RA ACCUMULATED
ACCUULATED
OF BRER4H HIDROGRAPH HIDROGRAPH ERROR ERROR
(HOURS) (HOURS) (CPS) (CPS) (CPS) (Cps) (AC-FrT)
16.750 0.000 36297. 36297. 0. 0. -0.
16.755 0.005 60454. 36651. 23803. 23803.
16.760 0.010 04612. 37634. 46978. 70781.
16.765 0.01 108769. 39427. 69342. 140123. 58.
16.770 0.020 132927. 42174. 90753. 230876. 95.
16.77S 0.025 157004, 45994. 111090. 341967. 141.
16.784 0.030 181241. 50991. 130251. 472217. its.
16.785 0.035 205399. 57253. 148146. 620363. 256.
16.790 0.040 229556. 6485S. 164698. 785062. 324.
16.795 0.045 253713. 73874. 179840. 964901. 399.
16.800 0.050 277871. 84359. 193512. 1158414. 479.
16.605 0.055 302020. 96362. 205666. 1364080. 564.
16.810 0.060 326186. 109922. 216263. 1580343. 653.
16.815 0.065 350343. 12507S. 225266. 1805611. 746.
16.020 0.070 374500. 141868. 232653. 2038264. 842.
16.82S 0.075 398658. 160258. 236400. 2276663. 941.
16.830 0.080 422815. 180320. 242495. 2519158. 1041.
16.835 0.065 446973. 202040. 244932. 2764091. 1142.
16.840 0.090 471130. 225421. 245709. 3009800. 1244.
16.845 0.095 49S287. 250461. 244826. 3254626. 134S.
MGM85 0.100 519445. 277147. 242277. 3496903. 1445.
16.855 0.105 S43602. 305594. 238008. 3734911. 1543.
16.860 0.110 567759. 335715. 232044. 3966955. 1639.
16.865 0.115 591917. 367271. 224646. 4191601. 1732.
16.070 0.120 616074. 400189. 215885. 4407466. 1821.
16.87S 0.125 640232. 434384. 205847. 4613333. 1906.
16.880 0.130 664389. 469763. 194626. 4807959. 1987.
16.66S 0.135 608546. 506221. 182325. 4990204. 2062.
16.890 0.140 712704. 543644. 169060. 5159344. 2132.
16.895 015 73601 581904. 154957. 5314301. 2196.
16.900 0.150
IS 7'6019. 620868. 140151. 5454452. 2254.
16.905 0.155 75176. 660181. 124795. 5579246. 2305.
16.910 0.160 009333. 700283. 109051. 5688297. 2351.
16.915 0.165 833491. 740403. 93088. 5701385. 2309.
16.920 0.170 857648. 780562. 77066. S56470. 2421.
16.925 0.17S 881805. 820563. 61222. 5919693. 2446.
16.930 0.160 905963. 860405. 45558. 5965250. 246S.
16.935 0.185 930120. 099949. 30171. 5995422. 2477.
16.940 0.190 954278. 939176. 15102. 6010524. 2484.
16.945 0.195 978435. 978312. 123. 6010647. 2484.
16.950 0.200 1002592. 1016305. -13713. 5996934. 2470.
16.955 0.205 1026750. 1052711. -25961. 5970973. 2467.
16.960 0.210 1050907. 1087237. -36329. 5934643. 2452.
16.965 0.215 1075065. 1119568. -44502. 5890142. 2434.
16.970 0.220 1099222. 1149390. -50168. 5839974. 2413.
16.975 0.225 1123379. 1174368. -5299. 5786985. 2391.
16.980 0.230 1147537. 1199788. -52251. 5734734. 2370.
16.985 0.235 1171694. 1218913. -47219. 5607514. 2350.
36.990 0.240 1195051. 1233145. -37293. 5650221. 2335.
16.995 0.245 1220009. 1241813. -21804. 5628417. 2326.
17.000 0.250 1244166. 1244166. -0. 5628417. 2326.
145
,, ,o ,, .. ° , .o . ,o ~ , o ° , o , , ° o o .. . . . . . oo~
. . ..
. • a 0
t I, ..... . .. .. .. .. .. .. .. . .. . . .. . . ..
.. .. . . . .° ., . .. . .. . ., . ,. . .. . . ,. . . e o , o
.. ....
t0
. 0 a
L .
f1.. . . . . . . . . . ...... ..
........... 0....... ......
00. ..... . . .
... .....
.. ..
101 .
'A
a
0ll .. . .
R ..oo..
° 00000a a
C;u
(1 0 0 0 3 4 al go,24.
16 0 26 a~ 622 2.3
I I( D I l G-HR I a
24-M 72-HR 34.75-HR
147MI6. 11t 13104. so... 35124. 4|73.
M3 PA427.OS
1 401.28 363.63 i$3.63
11l715A R ARMi - 72.10 SQK
146
26 KS PLAN 2 POR STATION OUT SUAR CREEK DAN
27 KO OUTT COWZRL VARIARLES
P136T 1 PRINT COWNL
IPW1T 0 PL0 CO WTOL
QSCAL 0. NYORIGRAPH PLC" SCALK
21 BE ELEVATION 340.00 380.00 410.00 420.00 424.00 428.00 432.00 436.00 440.00 444.00
24 EL Li4W-LEVEL OUTLET
9EE9 380.00 ELEVATION AT CENTER 0r OU"UE
CARERA 12.60 CitoSS-seCTromAL AREA
COQ, 0.70 COEVICIVET
Rx?. 0.S0 EXPONENT OF ERAD
22 SE SPILLWAY
CREL 420.00 SPILIMAY CRUET ELEVATION
5566D 2pO.00 BPZ1LU.AY WID"R
COwM 3.10 IR COEFICIET
WXPW 1.50 ERtPOEI? o RnAD
23 ST TOP O DAN
TPELt 432.00 ELEVATION AT ToP OF DAn
DMI96D 900.00 DAN IDTH
COQn 3.10 '63R COEPVICTI T
3151 1.50 EXPORENT Or READ
28 SI BREACH DATA
WLEN 340.00 ELEVATION AT SOTTON OF BRACR
810410 900.00 WIDTH 0F BREACH T
2 0.00 3R11CR SIO SLOPE
TFAIL 3.00 TINE FOR BREACH TO DEVELOP
FAILEL 433.00 M.S. ELEVATION TOD TRIGGER PAIWR
STORAGE 0.00 1333.33 1550.59 1610.64 1700.12 1342.50 2091.06 2590.89 3870.13 9160.68
oUTLOW 0.00 0.00 102.19 114.85 131.09 152.68 182.78 227.66 301.76 447.38
ELEVATION 340.00 380.00 382.09 382.64 383.43 384.66 386.68 390.36 398.20 420.00
STORAGE 9234.42 9453.83 9824.07 10353.84 11060.17 11970.54 13113.30 14522.20 16250.56 18397.72
OuTPLON 524.09 1044.57 2444.54 5159.65 9625.53 16277.80 25552.06 37883.09 53708.91 73462.71
ELE1VFATION 420.25 420.97 422.17 423.85 426.01 428.65 431.77 435.37 439.45 444.00
DA RON HRH ORD OUTFLOW STORAGE STAGS * DA MON ROIN ORD OUTVLON STORAGE STAGE DA NON RIN OR) OUTPW STORAGE STAGS
1 0000 I 447. 9160.7 420.0 * 1 1145 48 4989. 10324.0 423.6 * 1 2330 95 62759. 10.7 348.0
1 0015 2 447. 9152.9 420.0 * 1 1200 49 5714. 10449.0 424.2 a 1 2345 96 60177. 9.8 347.8
1 0030 3 447. 9145.0 419.9 a 1 1215 50 6450. 10571.1 424.5 * 2 0000 97 57723. - 9.0 347.6
1 0045 4 447. 9137.0 419.9 * 1 1230 51 7193. 10690.5 424.9 * 2 001S 98 SS388. 8.3 347.4
1 0100 5 447. 9128.9 419.9 1 1245 52 7945. 10008.0 425.3 * 2 0030 99 53158. 7.7 347.2
1 oils 6 447. 9120.8 419.9 * 1 1300 53 8712. 10924.6 425.6 * 2 004S 100 51023. 7.1 347.0
1 0130 7 446. 9112.6 419.8 * 1 1315 54 9501. 11041.9 426.0 * 2 0100 101 48972. 6.5 346.8
1 014S 8 446. 9104.3 419.8 * 1 1330 S5 10323. 11161.4 426.3 * 2 011s 102 46999. 6.0 346.6
1 0200 9 446. 9096.0 419.8 * 1 1345 S6 11190. 11285.0 426.7 * 2 0130 103 45097. 5.5 346.4
1 0215 10 446. 9087.6 419.9 * 1 1400 57 12117. 11414.6 427.1 * 2 0145 104 43259. S.1 346.2
147
1 0230 11 446. 9079.2 419.7 * 1415 58 13119. 11SS2.0
11699.1
427.S
427.9
*
*
2
2
0200
021S
105
106
41479.
39753.
4.7
4.3
346.1
345.9
024S 12 446. 9070.7 419.7 * 1430 59 14211.
0300 13 446. 9062.1 419.7 1 1445 60 15413. 11656.0 428.3 * 2 0230 107 38076. 3.9 345.7
0315 14 445. 053.6 419.6 e1 1S00 f1 16742. 12030.S 428.8 * 2 024S 108 36443. 3.6 345.6
0330 1 44S. 9044.9 419.6 1 1515 62 18213. 12219.0 429.3 * 2 0300 109 34853. 3.3 345.4
034S 16 445. 9036.3 419.6 * 1530 63 19842. 12422.0 429.9 * 2 0315 110 33305. 3.0 345.2
0400 17 445. 9027.6 419.6 1 154S 64 21666. 12646.6 430.5 2 0330 ti1 31800. 2.7 345.1
0415 16 445. 9016.9 419.3 * 1600 65 23747. 12898.3 431.2 * 2 0345 112 30338. 2.5 344.9
0430 19 445. 9010.1 419.5 I 1615 66 26139. 13162.6 432.0 * 2 0400 113 28919. 2.3 344.6
0445 20 444. 001.4 419.5 1 1630 67 30405. 13482.9 432.7 * 2 0415 114 27544. 2.1 344.6
1 500 21 444. 6992.6 419.4 1 1645 68 36297. 13773.4 433.5 2 0430 115 26213. 1.9 344.S
051S 22 444. 8963.6 419.4 * 1700 69 47771. 13997.9 434.1 * 2 0445 116 24929. 1.7 344.3
0530 23 444. 0974.9 419.4 * 1715 70 72262. 13951.6 433.9 * 2 0500 117 23694. 1.5 344.2
0545 24 444. 8966.2 419.3 1 1730 71 107456. 1338.3 432.S * 2 0515 118 22511. 1.4 344.0
0600 25 444. 8957.6 419.3 1 1745 72 151366. 12082.6 429.0 * 2 0530 119 21386. 1.2 343.9
0615 26 443. 8949.S 419.3 1 1800 73 165889. 10046.0 422.9 * 2 0545 120 20316. 1.1 343.8
0630 27 443. 8942.4 419.3 1 1615 74 208217. 7513.2 414.2 * 2 0600 121 19300. 1.0 343.6
0645 20 443. 8937.0 419.2 1 1830 7S 202406. 4679.2 403.3 * 2 0615 122 16335. 0.9 343.S
0700 29 443. 6934.1 419.2 1 1845 76 172793. 2721.1 391.2 * 2 0630 123 17419. 0.8 343.4
071S 30 443. 8934.7 419.2 * 1900 77 141211. 1298.4 379.6 * 2 0645 124 16548. 0.7 343.3
0730 31 443. 8939.9 419.3 1 1915 78 117067. 509.6 369.0 * 2 0700 125 15720. 0.7 343.2
0745 32 443. 8950.7 419.3 1 1930 79 103749. 1SS.0 319.5 * 2 0715 126 14934. 0.6 343.1
1 0600 33 444. 8968.3 419.4 1 1945 80 97536. 25.8 350.7 * 2 0730 127 14168. 0.5 343.0
1 0615 34 444. 6994.0 419.4 * 2000 81 95630. 24.6 350.6 * 2 0745 128 13476. 0.S 342.9
0630 31 445. 9026.6 419.6 1 2015 62 95238. 24.6 350.6 * 2 0800 129 12605. ,.4 342.8
0845 36 446. 9074.0 419.7 1 2030 83 94506. 24.2 350.5 2 0815 130 12164. 0.4 342.7
0900 37 447. 9130.9 419.9 * 204S 84 93048. 23.5 350.4 * 2 0830 131 11556. 0.4 342.6
0915 38 477. 9199.9 420.1 1 2100 85 90939. 22.4 350.2 * 2 0845 132 10979. 0.3 342.5
0930 39 605. 9280.2 420.4 * 2115 86 88478. 21.2 350.1 * 2 0900 133 10430. 0.3 342.4
1 0945 40 613. 9371.2 420.7 1 2130 87 65830. 20.0 349.9 * 2 0915 134 9908. 0.3 342.3
1 1000 41 1099. 9471.7 421.0 * 1 2145 88 83044. 18.7 349.6 - 2 0930 135 9413. 0.2 342.3
1 OIS 42 1464. 9560.5 421.4 - 1 2200 89 80154. 17.4 349.4 - 2 0945 136 8942. 0.2 342.2
1 1030 43 1904. 9696.1 421.8 1 2215 40 77201. 16.2 349.2 * 2 1000 137 8495. 0.2 342.1
1 1045 44 2414. 9817.2 422.2 * 2230 91 74224. 14.9 349.0 * 2 1015 138 8070. 0.2 342.0
1 1100 45 2987. 9942.1 422.6 1 2245 92 71254. 13.8 348.7 * 2 1030 139 7665. 0.2 342.0
1 1115 46 3615. 10069.2 423.0 * 2300 93 68323. 12.7 346.5 * 2 1045 140 7280. 0.1 341.9
1 1130 47 4286. 10197.0 423.4 * 2315 94 65478. 11.6 346.2
THE DAM BREACH HYDROGRAPH WAS DEVR.IPED USINC. A TIME INTERVAL. OF 0.062 HOURS IRNG BREACH
rORMATION.
nmwSTREAM CAtCULATIONS WILL 11S9 A TIME INTERVAL OF 0.2"0 HOURS.
THIS TARLE COMPARES TE HYI R(3GRAPH FOR DOWNSTREAM CALrI-LATTONS WITH THE COMPUTED BREACH HYOIROGRAPI.
INN4cDtATz rLmS ARE INTERPOLATZD FROM END-OF-PERI0O VALuES.
TIME FROM IHTZRPOTAT?.fn COMTTE!)
7IME REGIMNING BRECH
8REACH = ERROR ACCUMULATE ACCUMILATED
Or BREACH HYnROGRAPH HYDOWGRAPH ERROR ERROR
(HOURS) (HOURS) (CPS) (Cr5) (CRS) (CPS) (AC-FT)
148
C ..
I . . . . . . .. . . . . C
. .~NI. . . . . . . . . . . . . . .
01m
ofo
=Cot
M el
InI
cooI
000
I . . ...
. . . . . . . . .. . . . . . . . . . .
.00 ... ..
(Noe:C1an RC2ouptmied
~0149
PRAX PLOW AND STAGE (Da-OF-PBRIOD) SUfMARY FOR MULTIPLE PLAN-RATIO ECONOMIC CONPUTATIONS
FLOW IN CUIC Y PER SECOND, ARM IN SQOAR MILES
TIM TO PEAX IN HOURS
m
RATIOS APPLIED TO FLOWS
STATIOm AREA PLAN RATIO 1
OPnATIOM
1.00
iYtmo~MAPi AT IN 72.10 1 rLOW 955158.
Tim 20.25
2 FLOW 95158.
TIME 20.25
3 FLOW 95158.
TIMe 20.25
4 FLOW 95158.
TINN 20.25
5 rLm 95158.
TINS 20.25
BaD TO OUT 72.10 1 Fum 1244166.
TIme 17.00
2 FU 208217.
TIME 18.25
3 FLOW 135679.
TImm 19.75
4 FI 179533.
TIME 19.50
5 nLOw 109207.
TIM 21.75
PEAK gTAGES IN rZXT
1 STAGE 433.54
TIM 16.78
2 STAGE 434.13
TIME 17.12
3 STAGE 435.35
TIN 17.94
4 STAGK 434.61
TIME 17.37
5 STAGS 436.63
TIME 19.00
ROUTED TO "OC1 72.10 1 FLOW 610857.
TriM 17.25
2 FLOW 197055.
TIME 18.50
3 FLO 127387.
TIME 20.00
4 FLOW 175177.
TIME 19.75
s Iow 108707.
TIN 22.00
150
SULARY Or DAN oVNlTOPPING/10RBAU ANALYSIS FOR STATION OUT
............... IITIAL
TI VALUE SPILLWAY CRT 'oP OF DAN
PLEVATIO 420.00 420.00 432.00
STORAGE 9161. 9161. 13200.
OUTFLO 447. 447. 26263.
151
ITASLE 1 STATION acul Dtll
*CHI ECNI aw ACN2 CN2 TCa?
FLOW FLOW STAGE STAGE FLOW FLOW STAGE STAGE
PLAN 1 5 1 5 1 5 1 5
RATIO 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
125 2 0700 17234.57 18132.50 298.13 298.39 19302.40 20286.18 267.35 267.53
126 2 0715 16374.46 17229.65 297.89 298.13 18359.13 19302.42 267.18 267.35
127 2 0730 15556.60 16370.82 297.65 297.89 17455.23 18357.52 267.02 267.18
128 2 0745 14779.22 15554.40 297.43 297.65 16591.47 17453.13 266.87 267.02
129 2 0800 14093.57 14778.59 297.15 297.43 15757.18 16589.63 266.64 266.87
130 2 0815 13439.12 14094.50 296.86 297.15 14982.55 15756.06 266.41 266.4
131 2 0830 12797.01 13"1.65 296.5? 296.86 14257.08 14982.53 266.20 266.41
132 2 0845 12174.88 12801.22 296.30 296.58 13567.93 14258.38 266. 00 266.20
133 2 0900 11576.66 12180.84 296.03 296.30 12909.51 13570.71 265.81 266.00
134 2 0915 11004.08 11584.40 295.78 296.04 12279.51 12913.87 265.62 265.81
135 2 0930 10457.59 11013.59 295.53 295.78 11583.02 12285.53 265.36 265.63
136 2 0945 9936.92 10468.86 295.30 295.54 10856.00 11593.43 265.02 265.36
137 2 1000 9441.39 9949.93 295.08 295.31 10243.P4 10867.32 264.74 265.03
138 2 1015 8969.93 9456.04 294.87 295.09 9699.18 10256.47 264.48 264.74
139 2 1030 8521.30 8986.14 294.67 294.88 9198.91 9713.28 264.25 264.49
140 2 1045 8094.50 8539.10 294.48 294.68 8737.39 9214.53 263.99 264.26
152
12.9 Example Problem #9: Multif lood. Analysis
The input data and program output are shown in Table 12.9a. In this
case, the runoff from the design flood is input directly; these data would
have been obtained from previous rainfall-runoff simulations. The RCH1
channel routing data are for the modified Puls method in which previous water
surface profile studies have determined the storage-outflow characteristics
of the reach. The RCH2 routing is from the ponding area, through the levee
culvert, and into the main river. Two important points should be made about
the input and output for this example:
(1) The mu~tiflood analysis data deck differs from a stream network
data deck by the addition oF a JR card (see problem 1 for an
example of a stream network analysis).
(2) The resulting peak flows and stages for each ratio of the design
flood are displayed in the summary output at the end of the
exhibited printout.
153
OS SUUUASI
flOUTIAGG MEAC I
r
urene 12.7N
Fimg Rocbe RieIai
154p
TABLE 12.9
Example Problem #9: Input and Output
aw-1 T PAGE 1
LINE ID ....... 1 ....... 2 ....... 3 ....... 4 ....... S ....... 6 ....... 7 ....... a ....... 9 ......10
7 K 100
8 KH DESIGN FLOOD FOR SIUASIN 100
9 A 35.1
10 QI 24 24 24 26 33 50 86 189 376 516
11 Q1 594 657 710 760 801 839 910 1044 1287 1921
12 QI 2995 3953 4599 5077 5363 5374 5099 4603 3980 3325
13 QI 2719 2200 1844 1540 1251 994 777 605 471 365
14 QI 281 0 0 0 0 0 0 0 0 0
15 EK RCHI
16 KM LOCATION OF EXISTING FLOOD HAZARD
17 RS 1 STOR -1. 0.
18 SV 0. 50. 475. 940. 2135. 3080. 0. 0. 0. 0.
19 SQ 0. 200. 1020. 2050. 6100. 10250. 0. 0. 0. 0.
20 !K 300
21 KM RUIOFF FR4 StU9AIN 300
22 A 49.1
23 QI 32 32 32 35 44 67 114 252 501 688
24 QI 789 877 940 1013 1068 1119 1214 1392 1717 2561
25 QI 3993 4273 6139 6727 7163 7179 6789 6137 5308 4433
26 QI 3622 2930 2458 2053 1665 1325 1032 806 628 487
27 QI 374
28 RK 300
29 KM COMINED UPSTREAM INFLOWS
30 aC 2
31 KR F "2
32 KM DAMAGE =U CATE
OR IN THIS REAC, LONLAND FLOODING
33 RS I STOR -1. 0.
34 SV 0. 400. 30000. 35000. 40000.
35 SE 840 845 855 857 859
36 SQ 0 1250 1500 1800 2000
37 Zz
7 100
V
V
15 DCN1
20 * 300
28 300 ...........
V
V
31 DCH2
I I I II I I
(***) 1DINOF7 ALS O 00IPUTD AT THIS LOCATION
155
EXAMPLE PRDBL04 NO. 9
MULTIFLOOD AIALYSIS
CqKBED WATERSHED
ENGLISH UNITS
DRAINAGE AREA SQUARE MILES
PRECIPITATION DEPTH INCHES
LENGTH, ELEVATION FEET
FLOW CUBIC FEET PER SECOND
STORAGE VOLUME ACRE-FEET
SURFACE ARM ACRES
TEMPERATURE DEGREES FAHRENHEIT
iP MLTI-PLAN OPTION
NPLAN I NUMBER OF PLANS
JR MULTI-RATIO OPTION
RATIOS OF RUNOFF
0.11 0.26 0.45 0.65 0.86 1.00 1.20 1.40 1.50
• *4 ,44 44 * *44 *4* 444 *4*te e **4* t* *4 * 44 *44 44* 444 44* 44* 444 4*4 444 *44444 444 444 *4*4*4*e **
7 K * 100 *
BA SUBBAS IN CHARACTERISTICS
TAREA 35.10 SUBBASIN AREA
KR RCHI
RS STORAGE ROUTING
NSTPS 1 NUMBEROF SUBREACHES
ITYP STOR TYPE OF INITIAL CONDITION
RSVRIC -1.00 INITIAL CONDITION
X 0.00 WORKING R AND D COEFFICIENT
156
20 9 • 300
*•• •• *••
22 BA 58BASIN CHARACTERISTICS
TAM 49.10 SIBJEBI" AA
28 9 * 300 •
30 HC BYDROGRAPH COMBINATION
ICUOI 2 NIUER OF NYDROGRAPU TO COMINE
KK RCH2 *
RS S70RAGE ROUTING
NSTPS I NUMBER OF SUDREACEES
ITYP STOR TYPE OF INITIAL CONDITION
RSVRIC -1.00 INITIAL CONDITION
X 0.00 YORxING R AND 0 0o3WICIDT
PEAK FLOW AND STAGE (3D-OF-PERIOD) SUMMARYFOR MULTIPLE PLAN-RATIO ECONOMIC COIPUTATIONS
FLOWS IN CUBIC FEET PER SRCOD, AREA IN SQUARE MILES
TINE TO PEAK IN HORS
ENTDOGRAPE AT 100 35.10 1 FLOW 591. 1397. 2418. 3493. 4622. 5374. 6449. 7524. 8061.
TinE 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
POUTED TOA cE 35.10 1 FLOW 429. 978. 1742. 2680. 3668. 4313. 5232. 6156. 6701.
TIME 28.00 28.00 28.00 28.00 28.00 27.00 27.00 27.00 27.00
E!RGRAPH AT 300 49.10 1 FLOW 790. 1867. 3231. 4666. 6174. 7179. 8615. 10051. 10768.
TIME 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
2 COUNINED AT 300 84.20 1 FLow 1162. 2688. 4687. 6892. 9339. 10959. 132S0. 15529. 16663.
TINS 25.00 25.00 25.00 26.00 25.00 25.00 25.00 25.00 25.00
POUTED TO IC82 84.20 1 FLOW 964. 1257. 1273. 1291. 1312. 1326. 1347. 1369. 1379.
TIME 28.00 33.00 37.00 39.00 40.00 41.00 43.00 45.00 46.00
00 PEAK STAGES IN FET **
1 STAGE 843.86 845.27 845.90 846.65 847.48 848.05 848.89 849.74 850.17
TI E 28.00 33.00 37.00 39.00 40.00 42.00 43.00 45.00 46.00
* NORMAL IM OF UC-1*
157
12.10 Example Problem #10: Multiplan, Multiflood Analysis
Preparation of the multiplan data for input into the required HEC-l
format can be simplified by following input conventions described in Section
10. Examples which demonstrate these conventions in the problem 10 data deck
are as follows:
(1) Inflows from subareas 100 and 300 are only specified once for all
three plans; same as for problem 9. Because the rainfall-runoff
response is assumed constant in all three PLANS, ratios are taken
of the runoff.
(2) Routing reach RCH2 specifies data for a storage routing in PLAN 1;
a KP card specifying PLAN 2 updates the storage routing with pump
information; and lastly, a KP card specifying PLAN 3, not followed
by any data, indicates PLAN 2 and PLAN 3 data for reach RCH2 are
equivalent.
(3) Note the use of the RN card for routing reach 2OO. In the existing
plan, PLAN 1, a reservoir is not included, and this is indicated
with an RN card. The PLAN 2 flood control scenario includes a
reservoir at station 200, which is indicated by the appropriate KP
card and routing data. There is no data specified for PLAN 3 in
this case (the KP card is absent) and hence the program defaults to
the PLAN 1 data and prints a message to that effect. This is
appropriate since there is no reservoir at station 200 for PLAN 3.
(4) Only PLAN 3 calls for a diversion as part of the flood control
system. However, diversion data are included in all three plans.
By program input convention, the data for PLANS 1 and 2 specify a
diversion of zero capacity which has the intended effect of
omitting a diversion for these plans.
158
too
RCH
1
300 300
RCH
a
100
RCH
1
RC"
2
159
Multiplan Analysis Results:
TABLE 12.lOa
STORAGE ROUTING RS
N STPS - 1
ITYP - SlUR
RSVRIC . -1
LOW-LEVEL OUTLET SL
Invert Elevation - 975 (m.s.l.)
Cross section Area - 35 (sq-ft.)
Discharge Coefficient = .7
Exponent of Head - .5
SPILLWAY SS
Crest Elevation = 1105 (m.s.1.)
Width = 35 (ft.)
Weir Coefficient - 2.8
Exponent of Head = 1.5
VOLUME-ELEVATION DATA
VOLUME: 0, 2500, 4000, 5200, 6800, 9000, 11500, 15500, 21000, 30000 SV
ELEVATION: 965, 1000, 1015, 1030, 1045, 1060, 1075, 1090, 1105, 1120 SE
DIVERSION FLOW: 0, 1400, 2000, 3400, 4800, 8000, 12200, 16200, 18550, DT
20000
PUMP DATA WP
Threshold Reservoir Elevation = 843.5 (ft.)
Pump Capacity - 3000 (cfs)
160
TABLE 12.10b
Example Problem #10: Input and Output
uM-1 M PAM 1
LINE ID ....... I ....... 2 ....... 3 ....... 4 ....... S ....... 6 ....... 7 ....... a ....... 9 ...... 10
8 II 100
9 04 POIENTIAL IMVOIR I'FLOW
10 SA 35.1
11 QI 24 24 24 26 33 50 86 189 376 516
12 QI 594 657 710 766 801 839 910 1044 1287 1921
13 QI 2995 3953 4599 5077 5363 5374 5099 4603 3980 3325
14 QI 2719 2200 1844 1540 1251 994 777 605 471 365
is QI 281 0 0 0 0 0 0 0 0 0
* ****e.e ******* * ****** * POPOSET) REUSEOIR DATA
16 KI 200
17 1I PWPOSEDMR RVOIR
18 MR
19 KP 2
20 RS 1 STOR -1. 0.
21 SL 975 35 .7 .5
22 SS 1105 35 2.8 1.5
23 SV 0 2500 4000 5200 6800 9000 11500 15500 21000 30000
24 99 965 1000 1015 1030 1045 1060 1075 1090 1105 1120
•* 1NO RESEVOIR PLAN 3
25 9P 3
26 RN
27 KR 325
28 14 DIVERT FLOW PLAN 3
* *************** b**********e**a*e* DIO" DTVMION
29 DT FLOW 20000
30 DI 0 2300 4100 6300 8800 1430, 20200 30400 33250
31 DQ
32 RP 2
* ~ DDIIY DIVERSION
33 D1 FLO 20000
34 DI 0 2300 4100 6300 6800 14300 20200 30400 33250
35 Do
36 U 3
37 DT FLOW 20000
38 DI 0 2300 4100 6300 8800 14300 20200 30400 33250 38000
39 DQ 0 1400 2000 3400 4000 8000 12200 16200 1a550 20000
40 K HCH1
41 I" POTENTIAL CEANNEL IVDIFICATION REAC
42 RS 1 STOR -1. 0.
43 SV 0. 50. 475. 940. 2135. 3080. 0. 0. 0. 0.
44 SQ 0. 200. 1020. 2050. 6100. 10250. 0. 0. 0. 0.
161
IC-1 =PUT PA 2
45 99 300
46 US RUN=l r1 sl0N Ur 300
47 IA 49.1
32 35 44 67 114 252 501 688
48 QI 32 32
1013 1068 1119 1214 1392 1717 2561
49 01 789 877 940
6727 7163 7179 6789 6137 5308 4433
so QI 3993 4273 6139
1665 1325 1032 806 628 487
51 Q' 3622 2930 2456 2053
52 01 374
53 99 300
54 10n cowmnw o101 IUm1os
55 BC 2
56 I9 350
57 13 3it VX DrVI Tdw rLON
58 OR ?LO1
59 RU 400
60 K" C0 a' I DIrVEM?
UPSTRMAND nIILNS
61 BC 2
62 RI ]CB2
63 I pwVosD PolwrNG PIAjT aITE
64 its 1 1CR -1. 0.
.65 SV 0. 400. 30000. 35000. 40000.
SE 840 845 855 857 859
66
67 SQ 0 1250 1500 1600 2000
* *.*.*.......**.**.*...*****ee*@e****e* PLAX 2 PWU DATA
6 RI 2
69 vI $43.5 3000
70 HP 3
71 z
8 100
16 200
29 .--- LUI
27 325
V
V
40 am
45 300
53 300 ............
58 .. ,- 1
56 350
59 400 ............
V
V
62 MIC12
162
* U.S. ARMY CORPS OF ENGINEERS *FLOW YDROGRAPH PACKAGE (NEC-I)
STNE NYDIROLOGIC ENGINEERING CENTER * * L FEBROAN 1981 N
609 SECOND STREET *FEBRUARY 1981
* REVISED 14 JUN 85
' DAVIS, CALIFORNIA 95616 •
' (91., 0-32.85 OR (FTS) 8"3285 • RUN DATE 2 JUL 8S TINE 13:45:17 *
EGLISH UNITS
DRAINAGE AREA SARE MILES
PRECIPITATION DIPTM INCHES
Ia.G'ru* ELEVATION TEXT
FLO CUBIC FEET PER SCOKND
STORAGE VOLWE ACRE-FUT
SUI ACE ARi ACRES
TOAETUR DOGRERS FAHREEIT
JP MTLTI-PLhN OPTION
HPLAN 3 NUMBER OP PLANS
8 K * 10
a ioo
* a
PLAN 2 INPUT DATA FOR STATION 100 ARE SAME AS FOR PLAN I
•*, a.. a.. a.. aa. a.. a.. a.. a.. a.. *aa a.. a.. a.. a.. a.. a..
PLAN 3 INPUT DATA FOR STATION 100 ARE SAME AS FOR PLAN I
.. a..
e, a.. ata a.. e*a
.. a a.. aa etaett a..
a. an at. a.. . tea a.. aaa a.. aa
a..tet aa a.. *a a.. tat a.. *.a a..
tee a a..
163
16K1 * 200 *
* *
PI4POSXD RUURYOIR
IS RN No ROUTING
20 RS STORAGE ROOTING
KSTPS 1 NU1MM OF SODRnj
ITYP STOR T OF INITIAL CONDITION
RSVRIC -1.00 INITIAL CONDITION
x 0.00 RiKING R AND 0 OEFICIENT
23 SV STORAGE 0.0 2500.0 4000.0 5200.0 6800.0 9000.0 11500.0 15500.0 21000.0 30000.0
24 SE ELEVATION 965.00 1000.00 1015.00 1030.00 1045.00 1060.00 1075.00 1090.00 1105.00 1120.00
21 SL O-IxBVJM OUTLET
KzV 975.00 S.EVATIOU AT CNE OF OUTLET
CAMS 35.00 CN6-S3RMIOL AIA
COQL 0.70 )IWtICIRMT
1PL 0.50 or mAD
22 eS SPILLWAY
CARL 1105.00 8PILIMYT CRT ESZVTION
SPwID 35.00 IPILLIMY WID
COON 2.80 WEIR COTVFICIEr
1.50 MO T OFRAD
26 RN NO ROOTING
164
27 I 32S *
DZVERT FLOW 1M 3
DT DIVRSiONi
ZSTAD FL DIVRION grmcw IDENTIFICAlON
DSUN 20000.00 IMUM POLUM TO BE DIVlETE
DT DIVERSION
ST&D FLOW DIVESION Rf10000311 IDENTIFICATION
DSTRIX 20000.00 NAZIWN[ VOLUME T0 BE DIVERTED
DI INFLOW 0.00 2300.00 4100.00 6300.00 8300.00 14300.00 20200.00 30400.00 33250.00
DO DIVEI i0W 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
DT DIVERSION
19TAD PLOW DIVERSION 3100003113 IDENTIFICATION
DBTRUX 20000.00 NAXINUI4 VOLUM 10 BE DIVEIDD
DI INLww 0.00 2300.00 4100.00 6300.00 8800.00 14300.00 30200.00 30400.00 33250.00 38000.00
DO DIVERTED FLOW 0.00 1400.00 2000.00 3400.00 4800.00 8000.00 12200.00 16200.00 18S50.00 20000.00
40 0*
42 RS STORUE ROUTING
USTPS 1 NUMER OF SOBiRBCE
u
ITYP 8T0R TYPE O INITIAL CONDITION
RSVRIC -1.00 INITIAL CONDITION
x 0.00 WOW1ING R ANDD OEFICIENT
*0* *00L L* *0* 0*0 00* 0*0 000 00* 0*0 *00 0*0 000 000 0*0 000 00*
165
* 0
45 * 300 0
* *
uBAm In o DATA
47 SA SSBASN CNACTS3ISTICS
TARSA 49.10 SUBSASI AM
0*0 *0* *0 0*0 0*0 000 00* *00 0*0 e**0 0e* 000 00 * 000 0*0 t* 00* 0*0 0*0 000 00* 0*0 000 0*0 000 **0 000 *0 00*
0.000**.oooo*0
* *
53 K * 300 *
0 0
4*000*000*4*00
55 sC UTOROmtAP OmINATION
IoIp 2 NOUIR OP UTDOGRAPBS TO OMBINE
SEn
K 350
59 0 400 0
61 aC ITDROGRAPH comINATION
XCOP 2 lMONSI OF E GRAPHS TO COMBINE
166
62 I *• C2 •
64 RS STOG NOUTING
"PS 1 UHM OF SURZ3S
ITP S1OR TPZ OF INITIAL CONDITION
3SVR C -1.00 INITIAL CONDITION
X 0.00 WMiING R ANDD OIPCrI3T
64 RS 97RAGM WOOTING
momP I 3533* Or SUB*MACU
ITTP 8YO* TP Of INITIAL CONDITION
38VU1C -1.00 INITIAL ONDITION
x 0.00 W*ZINGR AND P CIDIFICIENT
69 VP IMING DATA
64 RS STORAGE WraTNG
xSTPS 1 3MUMB OF SUBRErHmE
ITYP M1CR TTPX OP INITIAL CMDITION
RSVMIC -1.00 INITIAL ONDITION
X 0.00 NODIfIG R AND D 0WFICIENTO
69 P PUMPIN DATA
167
PEAK FLOW AND STAGE (100-Of-PERI5) SSYVr FOR NULTIPLEPAN.RATIO ECMIC CMWWUTATIOS
0aic
FLUM Ill FEET PEE SCOND, AREA 1II 0888E MILES
TIE TO PEAK 108 85R
EIOROAfl AT
100 35.10 1 FLOW 591. 1397. 2418. 3493. 62. 5374. 6449. 7524. Sl.
TIE 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
2 FLU 591. 1397. 2418. 3493. 4622. S374. 6449. 75A4. 8061.
TIE 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 2S.00
3 FLOW 591. 1397. 241R. 3493. 4622. 5374. 6449. 7524. 8061.
TINE 25.00 25.00 25.00 25.00 25.00 25.00 25.00 2S.00 25.00
ROUTEDTO
200 35.10 FLO 91. 139?. 2418. 3493. 4622. 5374. 649. 75214. 8061.
TIE 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
2 FLOW 367. 617. 864. 1052. 1206. 1317. 1467. 1573. 1627.
TIE 29.00 31.00 32.00 33.00 33.00 34.00 34.00 34.00 3S.00
3 FLOU 591. 139?. 2418. 3493. 4622. 5374. 6449. 7524. 8061.
TIME 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
DIVERSION TO
FLOW 35.10 1 FLOt 0. 0. 0. 0. 0. 0. 0 0. 0.
TIME 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
2 FLOW 0. 0. 0. 0. 0. 0. 0. 0. 0.
TIE 1.00 1.00 1.00 1.00 1.00 1,00 1.00 1.00 1.00
3 FLtO 360. 800. 1439. T98. 2332. 2811. 3483. 4085. 4386.
TINE 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
MTOROGRAPHAT
325 35.10 1 FLOW 91. 1397. 2418. 3493. 4622. 5374. 6449. 7524. 8061.
TINE 25.00 25.0 25.00 25.00 25.00 25.00 25.00 25.00 25.00
2 FLOW 367. 617. 86. 1052. 1206. 1317. 1467. 1573. 1627.
TIME 29.00 31.00 32.00 33.00 33.00 34.00 34.00 34.00 35.00
3 FLOt 231. 547. 979. 1695. 2290. 2563. 2965. 3438. 3675.
TIME 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
RUTED TO
RaT 35.10 1 FLW 429. 978. 1742. 2680. 3668. 4313. 5232. 6156. 6701.
TINE 28.0 U.N 28.00 28.00 28.00 27.00 27.00 27.00 27.00
2 FLUS 305. 551. 784. 910. 1135. 1241. 1389. 150s. 1557.
TIE 34.00 8.00 39.00 41.00 41.00 41.00 42.00 43.00 43.00
3 FLUS 199. 399. 675. 1129. 1626. 1868. 222. 2646. 2853.
TIME 27.00 28.00 28.00 28.00 28.00 28.00 28.00 28.00 28.00
T'W*OGAPN
AT
300 49.10 FLOW 790. 1867. 3231. 466. 6174. 717. 8615. 10051. 107..
TINE 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
2 FLUS 790. 1867. 3231. 4666. 6174. 7179. 8615. 10051. 10768.
TINE 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
3 FLUS 790. 1867. 3M1. 466. 6174. 7179. 8615. 10051. 10768.
file 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
2 CSMINED AT
300 84.20 1 Ftt 1162. 2688. 4687. 6892. 9339. 10059. 13250. 15529. 16663.
TINE 25.00 25.00 25.00 26.00 25.00 25.00 25.00 25.00 25.00
2 PLUS 979. 2176. 3649. 5181. 6177. 7333. 9332. 1I0S. 11571.
TINE 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
3 FLt 974. 2215. 3805. 5597. 750. 8712. 10420. 12175. 13108.
TIME 25.00 25.00 25.00 25.80 25.00 25.00 25.00 25.00 25.00
- PEAKSTAGESIN FIT
1 STAGE 843.86 845.27 845.90 866.65 847.48 868.05 848.89 849.74 850.17
TIME 28.00 33.00 37.00 39.00 40.00 42.00 43.00 45.00 46.00
2 STA0 843.21 864.51 845.02 345.09 845.41 845.62 845.97 846.36 846.60
TIME 28.00 25.00 24.00 28.00 30.00 30.00 32.00 33.00 33.00
3 ST0G 043.74 845.01 845.89 845.51 846.1 4 W6.57 847.23 847.93 848.32
168
12.11 Example Problem #11: Flood Damage Analysis
The data for the flood damage analysis is shown in Table 12.11a. The
listing of the input data deck and a summary of the analysis results is given
in Table 12.11b. Note that the economic data (beginning with the EC card) is
added at the end of the multiplan-multiflood data deck (no changes are made
to the multiplan-multiflood data).
Discussion of Results
169
TABLE 12.11a
Flood Damage Reduction Analysis Economic Data
1. RECORD IDENTIFIERS
LAND USE CATEGORIES: CN
1. 95 843.6 1. 845.0 0 0
2. 81 844.8 2. 845.5 720 10.5
3. 60 846.6 3. 847.0 1380 15.0
4. 45 846.0 4. 847.6 2710 52.5
5. 25 846.6 5. 848.3 5200 105.0
6. 11 847.3 6. 849.0 8000 202.5
7. 5 857.9 7. 849.8 10050 540.0
8. 2.5 848.4 8. 851.0 11250 585.0
9. 1 849.1
10. .5 849.5
11. .2 850.0
12. .1 850.3
170
TABLE 12.11b
Example Problem #11: Input and Output
LINE ID ........1 ....... 2 ....... 3 ....... 4 ....... S ....... 6 ....... 7 ....... 8 ....... 9 ...... 10
r
1 ID EXAMPL PROBLE NO. 11
2 ID F100D DAMAGEAZALTSIS
3 ID ROCKBED
WATRESWED
*DIAGRAM
4 IT 60 0 0 130
5 10 5
6 **M*********************•****** WLTI PLAN AND RATIO DATA
f 3P 3
7 JR FIAN .11 .26 .45 .65 .86 1.00 1.20 1.40 1.50
8 Xg 100
9 RN DEIGN FLOOD SUSBASIN 100
10 BA 35.1
11 QI 24 24 24 26 33 50 66 189 376 516
12 QI 594 657 710 760 801 839 910 1044 1287 1921
13 QI 2995 3953 4599 5077 5363 5374 5099 4603 3980 3325
14 0! 2719 220C 1844 1540 1251 994 777 605 471 365
15 QI 281 0 0 0 0 0 0 0 0 0
*t PRPOSED RESERVOIR DATA
16 PK 200
17 K PROPOSED RZSERVOIR
18 RN
19 K? 2
20 RS 1 STOR -1 0.
21 SL 975 35 .7 .5
22 SS 1105 35 2.8 1.5
23 SV 0 2500 4000 5200 6800 9000 11500 15500 21000 30000
24 SE 965 1000 1015 1030 1045 1060 1075 1090 1105 1120
f ******f*O***fft***** ISORSRVIR PLAN 3
25 KK 325
26 DI DIVERT FLOW PLAN 3
* ********************* KY DIVERSION
27 OT FLOW 20000
28 DI 0 2300 4100 6300 8600 14300 20200 30400 33250
29 DQ
30 K? 2
* DUMIYf DIVERSION
31 DT FLOW 20000
32 DI 0 2300 4100 6300 8800 14300 20200 30400 33250
33 DO
34 Kr 3
35 OT FLOW 20000
36 DI 0 2300 4100 6300 8800 14300 20200 30400 33250 38000
37 DQ 0 1400 2000 3400 4800 8000 12200 16200 18550 20000
38 KR RCH1
39 KM LOCAL PROTECTION P2OJUCT PMDJWI FOERAC RC81
40 RS 1 O -1. 0.
41 SV 0. 50. 475. 940. 2135. 3080. 0. 0. 0. 0.
42 SQ 0. 200. 1020. 2050. 6100. 10250. 0. 0. 0. 0.
43 KR 300
44 KI DaIN FLOOD SUBSASIN 300
45 KM R8R? ?M SUDESIN 300
46 BA 49.1
47 QI 32 32 32 35 44 67 114 252 501 88
48 QI 789 877 940 1013 1068 1119 1214 1392 1717 2561
49 0 3993 4273 6139 6727 7163 7179 6789 6137 5308 4433
50 01 3622 2930 2458 2053 1665 1325 1032 806 628 487
51 QI 374
52 KR 300
53 KR CowUID UPSTREAM INFO
54 9C 2
55 IN 350
56 IN RETRIEVI DIVr T FLOW
57 OR FLO
58 K 400
59 KRM COlM E UPSTREAM AND OIrTB rNI.
60 vC 2
171
61 IU AC]2
62 Al DaINGE RMUM LlAND FLOWING PEDBL
63 10 PK)POSD PUMPING PIANT SITE
64 RS 1 IFOW -1. 0.
65 SV 0. 400. 30000. 35000. 40000.
66 OR 840 845 855 857 059
67 SQ 0 1250 1500 1800 2000
* *..*.*.***~t*.i******************o** PLAN 2 PUMP DATA
68 EP 2
69 VP 843.5 3000
70 RP 3
* ******,**********CO***(ICS DATA
71 mc
72 II NCH1
73 Cm 3 RESID IND/CO4 AGRIC
74 FR 18 700.0 600.0 550.0 450.0 350.0 250.0 150.0 90.0
75 FR 70.0 50.0 35.0 25.0 16.5 10.0 5.0 2.0 .5 .1
76 OF 400 490 530 640 000 1070 1480 1690
77 QF 1920 2170 2480 2850 3240 3640 4090 4900 5900 7100
78 Q0 15 400 600 730 960 1230 1530 1970 2500
79 Q0 3100 3490 3780 4290 5120 6020 7100
80 DG 1 3 0 1 2 3 5 7 28 49
81 DG 111 314 516 619 723 728 830
02 KK NCH2
83 Cm 3 RZSID IND/CO4 AGRIC
84 FR 12 95 81 60 45 25 11 5 2.5
85 FR 1 .5 .2 .1
86 Sr 843.6 844.8 845.6 846.0 046.6 847.3 847.9 848.4
87 SI 849.1 849.5 850.0 850.3
88 SD 8 845.0 845.5 847.0 847.6 848.3 849.0 849.6 851.0
89 00 1 1 0 720 1380 2710 5200 0000 10050 11250
90 DO 1 2 0 10.5 15.0 52.5 105.0 202.5 540 585
91 1z
8 100
V
V
16 200
27 - - LPTW
25 325
V
V
38 NCN1
43 300
52 300 ...........
57 . *--....
55 * 350
58 400 ...........
. .. V
.menmmmmtmmlmim~lolimimiillimJ~liiaa m
V
61 NCfl2
172
EXAMPLE PROLEM MD. 11
FLOOD DAMAGEANALYSIS
ROCKBED WARSHED
10 OUTPUT CNTROL VARIABLES
IPRWT 5 PRINT OROL
ZPLOT 0 PLOT CONTROL
QICAL 0.* rDROGRAPU PLOT SCALE
SG YES PRINT DIAGNOSTIC MESSAGE
IGLISH UNITS
DRAIKG AREA SWARE MILES
PRECIPITATION DEPTH INCHES
LU9GTE, n..VATION FEET
FLOW CUBIC FEET PER SECOND
STORAGE VOLUME ACRM-FEET
SURFACE A ACRS
TIPRATURE DEGRES FAHRENHEIT
jP NWLTI-PLAN OPTION
NPLAM 3 iUMER OF PLANS
JR I8LTI-RATIO OPTION
RATIOS OF RUNOFF
0.11 0.26 0.45 0.65 0.86 1.00 1.20 1.40 1.50
PEAK FLV4 ANDSTAGE (rn-P-PERIOD) SUMMARY FOR MULTIPLE PLAN-RATIO 0COP(C COMPTATIONS
FLows IN CUBIC FEE PER sscmS, ARm IN SQARE mILEs
TIME TO PEAK IN HYORS
RATIOS APPLIEDTrim
OPIRATION STATION AREA PLAN RATIO 1 RATIO 2 RATIO 3 RATIO 4 RATIO 5 RATIO 6 RATIO 7 RATIO I PATIO 9
0.11 0.26 0.45 0.65 0.06 1.00 1.20 1.40 1.50
HORGRAPR AT 100 35.10 1 FLOW 591. 1397. 2418. 3493. 4622. 5374. 6449. 7524. 6061.
TIME 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
2 FLOW 591. 1397. 2418. 3493. 4622. 5374. 6449. 7524. 6061.
TINS 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
3 FLOW 591. 1397. 2410. 3493. 4622. 5374. 6449. 7524. 8061.
TIME 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
ROUTD TO 200 35.10 1 FLOW 591. 1397. 2418. 3493. 4622. S374. 6449. 7524 061,
TIME 25.00 25.oc 25.00 25.00 25.00 25.00 25.00 25.06 25.00
2 F1W 367. 617. 864. 1052. 1206. 1317. 1467. 1573. 1627.
TIME 29.00 31.00 32.00 33.00 33.00 34.00 34.00 34.00 35.00
3 FLOW 591. 1397. 2418. 3493. 4622. 5374. 6449. 7524. 8061.
TIME 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
SPEAK STAGES IN FEET **
1 STAGE 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
TIME 0.00 0.00 0.00 00 0.00 0.00 0.00 0.00 0.00
2 STAGE 978.51 964.95 994.56 1003.99 1012.91 1020.02 1030.00 1039.32 1043.67
TIME 29.00 31.00 32.00 33.00 33.00 34.00 34.00 34.0 35.00
3 STAGE 0.00 0.00 0.00 0.00 0.00 0.00 000 0.0 0.00
TIME 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
DIVERSIOM ID FLOW 35.10 1 FLOW 0. 0. 0. 0. a
2
TIME
PLOW
1.00
0.
1.00
0.
1.00
0.
1.00
0.
I.06
0.
1.36
0. 0.
1.6
0.
I.FA
1166
0.
TIME 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
3 PLOW 360. 850. 1439. 1798. 2332. 2811. 3463. 4085. 4306.
TIME 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
f ,GRAPH
AT 325 35.10 1 PLOW 591. 1397. 2418. 3493. 4622. 5374. 6449. 7524. 4061.
TIME 25.00 25.00 25.00 25 00 25.00 25.00 25 0 25.0 10
2 FLOW 367. 617. 864. 102. 1206. 1307. 1447. 1573
TIME 29.00 31.00 32.00 33.00 33.00 34.00 34.00 34.00 35.0
3 PLOW 231. 547. 979. 1695. 2290. 2543. 2965. 3438. 3675.
TIME 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00 25.00
ROED TO NR1 35.10 1 FLOW 429. 976. 1742. 2680. 3660. 4313. 5232. 6156. 6701.
TIME 28.00 28.00 28.00 28.00 28.00 27.00 27.00 27.00 27.00
2 FLOW 305. 551. 784. 960. 1135. 1241. 1389. 1504. 1557.
TIME 34.00 38.00 39.00 41.00 41.0 41.00 42.00 3" 00 30
3 FLOW 199. 399. 675. 1129. 1 . 1868. 225. 266,. 2
TIME 27.00 28.00 28.00 28.00 20.00 28.00 20.00 28.00 20.00
173
HYNUOGRAPH AT 300 49 10 1 FLOW 790 1467. 222*J 4666 6174 7179 06*5. 10051 10768
TIME 25 00 25 00 25 00 25 00 25 00 25 00 25 00 25.00 25 00
2 FLOWI 790 1967 3231 4666 6174 7179 6is 10051 107d
TIME 25 00 25 00 25 00 23 00 25 00 25 00 25 00 2500 25.00
3 FLOW 790 1067. 3231 4666 6174 7179 0615 10051 107dm
TIME 25 00 25 00 250 0 25 00 25 00 25 00 25 00 250 0 25 00
2 COMBINED AT 300 04 20 1 FLOW 1162 26Q89 4697 6692 9339 10959 13250 t5529 16"3
TIlE 25 00 25 00 25 00 26 00 25 00 25 00 25 00 25 00 25 00
2 FLOW 979 2176 3649 591 &777 7833 9332 10025 11571
TIlE 25 00 25 00 25 00 25 00 25 00 25 00 25 00 25 00 25 00
3 FLOtd 974 2215 3005 5597 7500 8712 10420 *2175 13106
TIME 25 00 25 00 25 00 25 00 25 00 25 00 25 00 25 00 25 00
2 COMBINED AT 400 84 20 1 FLOW 1162 2689 4687 6892 9339 10959 13250 15529 1663
TIME 25 00 25.00 25.00 26 00 25.00 25.00 25"00 25.00 25 00
2 FLOW 979 217& 3649 5*6* 6777 7833 q322 lo05 S1571
TIME 25 00 25.00 25 00 25.00 25 00 25.00 25.00 25 00 25.00
3 FLOW 1333. 3065 5244. 7395 9832 11523. 13903 16261 17494
TIE 25 00 25.00 25 00 25.00 25 00 25 00 25 00 25,00 25,00
HYDROGRAPH AT RCH2 84 20 A LO 964. 1257. 1273 1291 1312. 1326. 1347 1369 1379.
TIME 28 00 33, 00 37 00 39.00 40 00 41.00 43 00 45 00 46 00
2 FLOW 802 1127. 1251 1252 1260 1265. 1274 1294 1290.
TIME 28 00 25.00 24.00 28.00 30 00 30.00 32 00 33.00 33.00
3 FLOW 935 1250. 1252 1263 *278. 1299. 1306 1323 1333
TIME 25 00 25.00 29 00 30 00 32 00 33 00 34 00 35 00 35 00
174
4DEAGB DATA VORPLAN 1 - ++
F32 FLow STAGE RNBID hDO/COM AMUC 70'AL
1 9 0 843.86 0.00 0.00 0.00 0.00
2 70.19 0. 845.27 387.69 5.65 0.00 393.34
3 48.51 0. 845.90 898.05 11.71 0.00 909.77
4 23.40 0. 846.65 1227.45 13.96 0.00 1241.41
5 8.77 0. 847.48 2451.80 45.22 0.00 2497.02
6 4.06 0. 848.05 4327.31 86.60 0.00 4413.91
7 1.36 0. 848.89 7559.55 187.16 0.00 7746.71
8 0.33 0. 049.74 9699.25 515.18 0.00 10414.43
9 0.13 0. 850.17 10422.49 553.97 0.00 10976.45
WEP ANNUAL DANAGI 1099.86 20.21 0.00 1120.06
++DA FAZ DATA PONPAY 2 -- 4+
175
12.12 Example Problem #12: Flood Control System Optimization
The data for the optimization model is shown in Table 12.12a. A number
of points should be noted about the data:
The input data in the appropriate HEC-l format and the output from the
model are shown in Table 12.12b. Note that the cost and optimization data
for the reservoir and pump are located in the stream network portion of the
input data deck, whereas, the local protection and degree of protection data
are located in the economic analysis portion of the data deck. The results
of the optimization analysis are shown in the output summaries at the end of
the computer output (Table 12.12b).
TABLE 12.12a
RECORD IDENTIFIERS
RESERVOIR DATA
PUMP DATA
CAPACITY-COST DATA
CAPACITY (cfs) 0, 5000, 5500, 7000, 8300, 9300, 12000, 15000, 21000 LC
COST ($ THOUS) 0, 103, 149, 122, 283, 340, 600, 1000, 3000 LD
UPPER PATTERN 0, 1, 2, 3, 5, 7, 28, 49, 111, 314, 516, 619, 723, 728, 830
LOWER PATTERN 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, .44, 3.5, 7.15
177
TABLE 12.12b
Example Problem #12: Input and Output
WIC1 IPUT PM 1
LIVE ID ....... 1 ....... 2 ....... 3 ....... 4 ....... 5 ....... 8 ....... 7 ....... 8 ....... 9 ...... 10
I ID XXUANIZ PF'ULWE 12
2 ID 8WOD COROL SYSTM OPIIZATION
3 ID N MR241HED
*DIAGRAM
4 IT 60 0 0 130
5 10 4
6 06 -15000 -8000 -17000
* ***M**************************MOO**** DIU PLAN AD RATIO DATA
7 jP 2
8 3R FIO6 .11 .26 .45 .65 .86 1.00 1.20 1.40 1.5
9 IK 100
10 IM POTENTIAL RESERVOIR INFU
11 DA 35.1
12 QI 24 24 24 26 33 s0 86 189 376 516
13 QI 594 657 710 760 801 839 910 1044 1287 1921
14 QI 2995 3953 4599 5077 5363 5374 5099 4603 3980 3325
15 QI 2719 2200 1844 154C 1251 994 777 605 471 365
16 QI 281 0 0 0 0 0 0 0 0 0
* PR POSED RBSVOIR DATA
17 KR 200
18 K PROPOSED RESERVOIR
19 RN
20 KP 2
21 is 1 STOR -1. 0.
22 SO 1 .023 .0504 29000 0
23 SL 975 35 .7 .5
24 S8 1105 35 2.8 1.5
25 SV 0 2500 4000 5200 6600 9000 11500 15500 21000 30000
26 SE 965 1000 1015 1030 1045 1060 1075 1090 1105 1120
27 SD 0 1500 2400 3000 3600 4350 4950 5550 6000 7200
28 R AC1
29 AG POTENTIAL CANNU IODIFICATION hCH
30 RS 1 1103 -1. 0.
31 SV 0. 50. 475. 940. 2135. 3080. 0. 0. 0. 0.
32 SQ 0. 200. 1020. 2050. 6100. 10250. 0. 0. 0. 0.
33 R 300
34 AG RINDOP FROM SI3DHASIN 300
35 BA 49.1
36 QI 32 32 32 35 44 67 114 252 501 688
37 QI 789 877 940 1013 1066 1119 1214 1392 1717 2561
38 QI 3993 4273 6139 6727 7163 7179 6789 6137 5308 4433
39 QI 3622 2930 2458 2053 1665 1325 1032 806 628 467
40 01 374
41 II 300
42 IN OMINEUD OPSTRRA INFWS
43 1C 2
44 n FIl2
45 AG PRPOSHI P0160C lART SITE
46 R 1 STOR -1. 0.
47 6V 0. 400. 30000.
48 SR 840 845 855
49 SQ 0 1250 1500
* e PLAN 2 POW DATA
50 9P 2
51 No 2 .023 .0504 100 10000
52 WP 843.5 3000
53 vc 0 250 500 1000 2000 6000 8000 10000
54 iD 0 670 1000 1600 2300 6000 7860 8670
a ******a~ta******ee..a.**aaaaaaa,.aaaa CS DATA
55 IC
56 R RCH1
57 ClH 3 RUmID IND/COM AGRIC
5 FR 18 700.0 600.0 550.0 450.0 350.0 250.0 150.0 90.0
59 FR 70.0 50.0 35.0 25.0 16.5 10.0 5.0 2.0 .5 .1
60 QC 400 490 530 640 800 1070 1480 1690
61 QP 1920 2170 2480 2850 3240 3640 4090 4900 5900 7100
62 QD 15 400 600 730 960 1230 1530 1970 2500
63 QD 3100 3490 3780 4290 5120 6020 7100
64 DG 1 3 0 1 2 3 5 7 28 49
65 OG 111 314 516 619 723 728 830
1 7R
66 LO 3 .023 .0504 21000 0
67 LC 0 5000 5500 7000 8300 9300 12000 15000 21000
68 LID 0 103 149 222 263 340 600 1000 3000
69 vu 3 0 1 2 3 S 7 28 49
70 DO 111 314 516 619 723 726 630
71 UL 3 0. 0. 0. 0. 0. 0. 0. 0.
72 DL 0. 0. 0. 0. 0.44 3.5 7.15
73 iU ItE2
74 CM 3 AMID MIO/O JIC
75 FI 12 95 81 60 45 25 11 S 2.5
76 rt 1 .5 .2 .1
77 SP 843.6 844.8 845.6 646.0 846.6 847.3 647.9 846.4
78 SP 849.1 849.5 850.0 850.3
79 SD 6 645.0 845.5 847.0 647.6 648.3 649.0 849.8 651.0
80 Da I 1 0 720 1380 2710 5200 8000 10050 11250
81 Dr 1 2 0 10.5 15.0 52.S 10S.0 202.5 540 585
82 Up 1 846.9
83 22
SCHEMATIC DIAGRAM OF STREAM NETWORK
INPUT
LINE (V) ROUTING (I'NG DIVERSIONOR PUMPFLOW
9 100
V
V
1? 200
V
V
28 RCHI
33 300
41 300 ............
V
V
52 .......
44 RCH2
EXAIPLE PNWBSST.4 12
FLOOD COOTRL 8135KM OPTIMIZATION
FOCMMMATBRtS8X
PTATION INTERVAL
COMP 1.00 SOORS
TOTAL TIME BASE 129.00 ROM
ENGLIzs UNITS
DRAINAGE ARM so"R MILES
PRECIPITATION DEPTE I
LaG5T, XLEVATIow FT
FLOW CUBIC I= PIU ANY r
STORAGE VOLWNE ACRE-IEET
SOWACE AREAACE
TUWRRATURN ONIn PAPNSRIT
JP NULTI-PUJ OPTION
NPLAl 2 NIUMER OF PLANS
JR ILTI-RATIO OPTION
RATIOS OF RUNS?
0.11 0.26 0.45 , 0.65 0.86 1.00 1.20 1.40 1.50
179
PEAK FLOI4 AND STAGE 4END-F-110)0 IUEARY FOR MU.LTIPLE PLAN-RATIO ECONOMIC COMUTATIN
FLMW IN CUNIC FEET PER SECOND. AREA IN SQUARE MILES
TIME TO PEAK IN HOURS
HYDROGRAPH AT t00 33 10 I FLOW 591 1397 2418 3493 4622 3374 6449 7524 8061
TIME 25 00 25 00 25 00 25 00 25 00 25.00 25 00 25 00 29 00
2 FLOW4 591 1397 2410 3493 4622 5374 6449 7524 8061
TIME 25 00 2500 25 00 25 00 25 00 29 00 25 00 25 00 25 00
ROUTED TO 200 35 10 I FLOWl 591 1397 2410 3493 4622 5374 6449 7524 0061
TIME 25 00 25 00 25 00 25 00 25 00 25 00 25.00 25 00 25 00
2 FLOW 367 617 066 1052 1206 1315 1467 1573 1627
TIME 29 00 31 00 32 00 33 00 33 00 34. 00 34 00 34.00 35 00
ROUTED TO RCHI 35.10 1 FLIW 429 978 1742 2690 3668 4313. 5232 6156 6701.
TIME 28 00 29 00 29.00 29 00 29. (O 27.00 27.00 27 00 27.00
2 FLO 303 551 795. 900 1135, 129. 1389 1504 1557
TIME 34.00 39 00 39 00 41.00 41 00 41 00 42 00 43 00 43 00
HYDROGRAPH AT 300 49.10 1 FLOW 790 1867 3231. 4666 6174 7179. 615 10051. 10768
TIME 25 00 25 00 25 00 25 00 25 00 25.00 25 00 25. 00 25 00
2 FLOW4 790 1867 3231 466 6174. 7179. 9615 10051 10768
TIME 25.00 25,00 23.00 25 00 25 00. 25 00 2500 25 00 2500
2 COMU INED AT 300 84.20 I FLOW 1162 2688. 4687 6992 9339 10959. 13250, 15529 16663
TIME 25 00 25 00 25 00 26 00 25 00 25.00 25.00 25.00 29.00
2 FLO 990 2176 3649 5191 6777 7933. 9332 10825, 11571
TIME 25 00 25 00 25 00 25 00 25 00 25 00 25.00 25 00 25 00
HYDROGRAPH AT RCH2 94.20 1 FLOW 964. 1257 1273 1291 1312 1326 1347, 1369 1379
TIME 29 00 33 00 37 00 39 00 40 00 41.00 43. 00 45 00 46 00
2 F lW 902. 1127 1251. 1252 1260 1265 1274. 1294 1290
TIME 28 00 25 00 24 00 28 00 30 00 30 00 32 00 33 00 33 00
* T 129.22 0.00
180
OF COMNENTr
8SNWMMT OO"
AVRMG SIBTU
ANNUAL NET BMUNITS * ** * -94.
181
IW73IMIATE VALUES (O OPTIKIZATION VARIABLES
OBJECTIVE
?rWCTI6o VAR I VAR 2 VAR 3
15000.000 8000.000 17000.000 TROT PNLTY MIN COST MUDM63O6JCTrFUCI
4964.224 1201.122 141.834 6.66813+06
14850.000* 8000.000 17000.000
LOCATION TARGETf COIP VAL DEVIATH PENALTY
RNc2 846.90 846.06 0.84 4964.24
TROT PRLTY AM COST AU DAM OBJCT FICTh
4964.237 1199.471 141.836 6.65993+06
-- "' ' m i | I l I I I~
MEA FLOW AND STAGE (D-Op-PRXOO) SWIBAT FOR mumLtIP PnaN-RATIO scmmc c~uoTATnaIs
FLOWS IN CUBIC FET PER 'I ARIA IN SQUAREIID
MILES
TIME TO PEAR IN BOORS
2 COMBINED AT 300 84.20 1 FLOW 1162. 2688. 4687. 6892. 9339. 10959. 13250. 15529. 16663.
TIM 25.00 25.00 25.00 26.00 25.00 25.00 25.00 25.00 25.00
2 FLOW 982. 2176. 3649. 5182. 6777. 7645. 9476. 11279. 12345.
TIME 25.00 25.00 25.00 25.00 25.00 25.00 25.00 26.00 27.00
RTED TO RRl2 84.20 1 FLOW 964. 1257. 1273. 1291. 1312. 1326. 1347. 1369. 1379.
TIME 25.00 25.00 25.00 26.00 25.00 25.00 25.00 25.00 25.00
2 FLOW 804. 1070. 1250. 1256. 1266. 1276. 1294. 1312. 1322.
TIME 28.00 26.00 26.00 29.00 33.00 35.00 37.00 38.00 39.00
56 * W111 *
FR PERCfhT EXCEEDANCE
700.0 600.0 550.0 450.0 350.0 250.0 150.0 90.0
70.0 50.0 35.0 25.0 16.5 10.0 5.0 2.0 0.5 0.1
QP 400. 490. 530. 640. 800. 1070. 1480. 1690.
1920. 2170. 2480. 2650. 3240. 3640. 4090 4900. 5900. 7100.
DANA" DATA
FLOW RESZD IND/COH AGRIC TOTAL
400.0 0.000 0.000 0.000 0.000
600.0 0.000 0.000 1.000 1.000
730.0 0.000 0.000 2.000 2.000
960.0 0.000 0.000 3.000 3.000
1230.0 0.000 0.000 5.000 5.000
1530.0 0.000 0.000 7.000 7.000
1970.0 0.000 0.000 28.000 28.000
2500.0 0.000 0.000 49.000 49.000
3100.0 0.000 0.000 111.000 111.000
3490.0 0.000 0.000 314.000 314.000
3780.0 0.000 0.000 516.000 516.000
4290.0 0.000 0.000 619.000 619.000
5120.0 0.000 0.000 723.000 723.000
6020.0 0.000 0.000 728.000 728.000
7100.0 0.000 0.000 830.000 630.000
183
rm Man 0
MID Tmm0 A=jI 10
97,. 0.00 0.00 0.00 .14 3.14
3 ,5.07 1742. 0.0 0.00 0.00 17.11 17.11
4 29. 20. 0.00 0.00 0.00 47.65 47.65
--. 36. 0.00 0. 0.00 436.01 436.01
1 3.77 4313. 0.00 0.00 0.00 621.4 621.66
1 1.30 5232. 0.00 0.00 0.00 723.62 723.62
0.33 6156. 0.00 0.00 0.00 740.82 740.82
9 0.11 6701. 0.00 0.00 0.00 792.26 792.26
UP AUmAL DAMIm 0.00 0.00 129.22 129.22
44-U UWDM
O K S NPS 2-4
16" FLOW. MUi
PAGE I=/CON AIC W5L
S6.0 307. 0.00 0.00 0.00 0.0 0
2 279.57 551. 0.00 0.00 0.00 0.00 0.00
3 65.07 767. 0.00 0.00 0.00 0.00 0.00
4 29.26 993. 0.00 0.00 0.00 0.00 0.00
5 9.60 1764. 0.00 0.00 0.00 0.00 0.00
6 3.77 2426. 0.00 0.00 0.00 0.00 0.00
7 1.30 3440. 0.00 0.00 0.00 0.00 0.00
6 0.33 4369. 0.00 0.00 0.0 0.00 0.00
9 0.11 4072. 0.00 0.00 0.00 0.00 0.00
UP NWEUL DAMAG 0.00 0.00 0.00 0.00
AVERAGE ANUAL 6MITS 0.00 0.00 129.22 129.22
73 i * W2 *
1 73
93. FLW.
0 ST
3..03
643,.66 MID
0.00 Mg. 00i AMC
000 Toy"
0.0
2 70.19 0. 645.27 367.69 5.65 0.00 393.34
03 46.51 0. 845.90 896.05 11.71 0.00 909.77
,4 23:46 0. 646.45 1227.45 13.96 0.00 1241.41
6.77 0. 647.46 2451.60 45.22 0.00 2497.02
6 4.06 0. 648.05 4327.31 66.60 0.00 4413.91
7 1.36 0. 646.69 7559.55 167.16 0.00 7746.71
6 0.33 0. 849.74 9699.25 515.16 0.00 10414.43
9 0.13 0. 650.17 10422.49 553.97 0.00 10976.45
UP AiWUAL DAIJP 1099.66 20.21 0.00 1120.06
4.+DAAGE DATA POR PLAN 2 -
FLAW. 117.0 MlID MrDM~ AGUIC TOTRL
0 43.21 0.0 .00 0.00 00
2 70619 0. 044.26 0.00 0.00 0.00 0.00
3 0. 645.00 6.75 0.10 0.00 6.65
4 2304 O 645.24 340.6 4.97 0.00 345.04
5 8.77 0. 645.66 769.61 10.97 0.00 600.56
6 4.08 0. 646.06 965.46 12.17 0.00 977.6
7 1.36 0. 646.75 1269.03 14.24 0.00 1263.26
6 0.33 0. 647.47 2423.43 44.42 0.00 2467.65
9 0.13 0. 647.87 3663.36 72.60 0.00 3735.96
UP AIUUAL DAMAGE 234.03 3.31 0.00 237.34
AVERAGE ANWAL SmlW1.s 665.63 16.90 0.00 82.73
184
NNaw AIML FLOOD -WCI 6UNARY
or COMPCT ONT
SUMMARY
OPTINIZATICI RSULTS
VAR 1 VAR 2 VAR 3
2616.11 21S0.86 4941.73
AVZAM3 ANNUAL
DANAM - EXISTING CODITIOS * * f 1249.
DANAGE -
AVERAGE ANNUAL OPTINIZKD SYSTEN *t * *ftf 237.
AvURMm ANIMAL DANMEz AncTICII (UwiTs) ***1012.
Note: The results of this test are dependent on the machine word size.
Results are likely to be within five percent of the answer shown
above.
19,5
Section 13
COMPUTER REQUIREMENTS
Figure 13.1 shows the sequence of operations for most jobs. HEC-1 uses up
to 16 I/O and scratch files. These can be stored on disk, tape, or whatever
medium is available. The program knows these files by their assigned unit
numbers. Table 13.1 shows the unit numbers used by HEC. These numbers can be
changed for a particular installation by changing their definition in BLOCK
DATA. All files are sequential.
Table 13.2 lists compile time and memory required for execution. Also,
execution times are given for the example problems described in Section 12 of
this manual.
RADIPT DATA;
WRIT TO WORING FILE
-9 18 1
186
TABLE 13.1
35 ISOP Scratch; data for flood control system optimization U 2400 real words
36 LSFIL Scratch; data for user-defined output tables U 301 real words
** Unit number is defined by user on KO or BI records (The unit numbers specified should not
conflict with other file definitions, for example, 21 and 22 are possible choices).
187
Table 13.2
I PC XT
(with 8087) CDC C r 175 Harris 500
3. Snownelt Runoff 80 .3 14
5. Routing Optimization 50 .1 7
9. Multiflood Analysis 70 .3 13
188
Section 14
REFERENCES
189
National Weather Service 1961. Rainfall Frequency Atlas of the United States.
Technical Paper No. 40, U.S. Department of Commerce, Washington, D.C.
1977. Five to 60-Minutes Precipitation Frequency
for the Eastern and Central United States. Technical Memo NWS HYDRO-35,
National Oceanographic and Atmospheric Atmospheric Administration, U.S.
Department of Commerce, Maryland.
1978. Probable Maximum Precipitation Estimates,
United States East of 105th Meridian, NOAA Hydrometeorological Report No.
51, Office of Hydrology, National Weather Service, Silver Spring, Maryland.
1979. DAMBRK: The NWS Dam-Break Flood Forecasting
Model. Technical Paper, Office of Hydrolgy, National Weather Service,
Silver Spring, Maryland.
1982. Application of Probable Maximum Precipitation
Estimates United States East of 105th Meridian, Office of Hydrolgy,
National Weather Service, Silver Spring, Maryland.
Snyder, F. F. 1938. Synthetic Unit Hydrographs. Transactions of the
American Geophysical Union, Vol. 19, Part 1, pp. 447-454.
Soil Conservation Service 1972. National Engineering Handbook, Section 4,
U.S. Department of Agriculture, Washington, D.C.
Soil Conservation Service 1975. Urban Hydrology for Small Watersheds.
Technical Release No. 55, U.S. Department of Agriculture, Washington, D.C.
Viessman, W. Jr., Knapp, J. W., Lewis, G. L. and Harbaugh, T. E. 1977.
Introduction to Hydrology, Dun-Donnelley Co., New York.
Woolhiser, D. A. 1975. Simulation of Unsteady Overland Flow. In, Unsteady
Flow in Open Channels. (K. Mahmood and V. Yevjevich, eds). Water
Resources Publications, Colorado.
190
Appendix A
Preface
This appendix contains a description of the input data for the HEC-1
computer program. It is only applicable to the 1981 version of the program.
Do not use this input format for any previous version (e.g. Dam Break) of the
program.
The yellow pages. A-1 through A-20, describe the general structure of the
input data, and data requirements and options for JOB DESCRIPTION and JOB
INITIALIZATION.
The blue pages, A-21 through A-92, describe the input data requirements
and options for HYDROGRAPH CALCULATIONS throughout a river basin. Record
types are arranged in alphabetical order.
The yellow pages. A-93 through A-111, describe the input data requirements
and options for ECONOMIC ANALYSES of flood damage, the required END-OF-JOB
record, and a SUMMARY of all input data records and variables.
A-i
HEC-1 INPUT DESCRIPTION
Table of Contents
1 Introduction
• Required records
A-iii
HEC-1 INPUT DESCRIPTION
12.1 PB and PI/PC Records - Storm Total and Distribution Option A-48
12.1.1 PB Record - Basin Average Precipitation A-48
12.1.2 PI aecord - Incremental Precipitation Time Series A-49
12.1.3 PC Record - Cumulative Precipitation Time Series A-49
12.2 PG Record - Storm Total Precipitation for a Station (Gage) A-50
12.3 PH Record - Hypothetical Storms A-51
12.4 PM Record - Probable Maximum Precipitation A-53
12.5 PS Record - Standard Project Precipitation (SPS) A-55
12.6 PR, PT and PW Records - Precipitation Gage Weighting A-56
12.6.1 PR Record - Recording Stations to be Weighted A-56
12.6.2 PT Record - Storm-Total Stations to be Weighted A-57
12.6.3 PW Record - Weightings for Precipitation Stations A-57
A-iv
HEC-1 INPUT DESCRIPTION
A-v
HEC-1 INPUT DESCRIPTION
18 Economic Data
A-vi
HEC-1 INPUT DESCRIPTION
INTRODUCTION
1 INTRODUCTION
This input description is ortanized into three major types of data: 1) job
description and initialization data, 2) hydrograph calculation data, and 3)
economic analysis data. This corresponds to the general sequence of data
necessary to build the digital model of a river basin as described in the next
subsection on Input Data Structure.
The first group (pages A-8 through A-20), JOB DESCRIPTION AND INITIALIZA-
TION DATA, begins with the I records and goes through the V records. The ID
and IT records are required and are described first. The other records are
optional and are described in a recommended input sequence, i.e., I, J, 0, V
records as desired.
The second group (pages A-21 through A-92), HYDROGRAPH CALCULATION DATA,
comprises all of the data necessary to simulate the various river basin
processes. The input data in this group are organized ALPHABETICALLY, beginning
with the B records and ending with the W records. The required and
recommended order to input these data are described in the next subsection,
Input Data Structure.
The third group (pages A-93 through A-106), ECONOMIC ANALYSIS DATA,
consists of data to be supplied after all of the hydrologic and hydraulic
calculations are completed. These data are optional and begin with the EC
record and are organized in the recommended sequence of input.
The last record described is the REQUIRED ZZ RECORD, pae A-107, to end
the job.
A-1
HEC-1 INPUT DESCRIPTION
INPUT DATA STRUCTURE
The input data set is divided into three sections - job description and
initialization data, hydrograph calculation data, and economic analysis data.
The third section begins with an EC record. All data following the EC
record are for calculation of expected annual damages.
Finally the job is) terminated by a ZZ record. Data for a new job
beginning with an ID record may follow immediately after the ZZ record.
The record sequence for a typical job is shown on the next page. A dash,
-, is used to indidate the second character of a record identification which
will be selected at the option of the user.
Continued
A-2
HEC-1 INPUT DESCRIPTION
INPUT DATA STRUCTURE
ID Job identification
IT Time specification
I-* Additional initialization data
0- * Optimization
ZZ End-of-job record
*Optional records
Continued
A-3
HEC-1 INPUT DESCRIPTION
INPUT DATA STRUCTURE
Data input for RUNOFF calculations will be retained and used for subsequent
runoff calculations until new data are read. Thus the data used in calculating
runoff need only be read once, unless they are to be changed for a new basin.
A typical record sequence for computing subbasin rainfall-runoff is:
( KK Etc. )
C )
*Optional records
Continued
A-4
HEC-l INPUT DESCRIPTION
INPUT DATA STRUCTURE
*Optional records
Continued
A-5
HEC-l INPUT DESCRIPTION
INPUT DATA STRUCTURE
When data are entered in FIXED FORMAT the record is divided into ten
fields of eight columns each, except field one. Variables occurring in field
one may only occupy columns 3-8 because columns 1 and 2 are reserved for the
record identification characters. Integer and alphanumeric values must be
right justified in their fields.
Data may also be entered in FREE FORMAT where fields are separated by a
comma or one or more spaces. Successive commas are used to indicate blank
fields. When entering time series data (flow, precipitation, etc.), more (or
less) than 10 values can be placed on a record.
A--6
HEC-1 INPUT DESCRIPTION
INPUT CONTROL RECORDS
The following records may be used to control the format and printing of
the input data. An input comment record is also described which may be
inserted anywhere in the input data stream.
NOTE - The asterisk (M) must be in column 1 and followed by the remainder of
the identification. If column 2 is blank, it is assumed to be a COMMENT
record.
A-7
ID HEC-l INPUT DESCRIPTION
JOB INITIALIZATION (I Records)
The ID and IT records are required to begin the job. The other records
(IN AND 10) are only used if those options are desired.
At least one ID record is required but any number may be used as desired
to title the output from this job. The title information is contained in
columns 3-80 inclusive and any characters or symbols may be used.
**REQUIRED
A-8
IT
HEC-1 INPUT DESCRIPTION
JOB INITIALIZATION (I Records)
The IT record is used to define time interval, starting date and time, and
*CAUTION: IDATE and ITIME are the time of initial flow conditions. No runoff
calculations are made from precipitation preceding this time.
Use 3-character code for month: JAN, FEB, MAR, APR, MAY, JUN, JUL, AUG, SEP,
OCT, NOV, DEC. Use of any other code for month means this is not a date, and
days will be numbered consecutively from the given day. Default is day = 1.
**REQUIRED
A-9
NHEC-1
Il INPUT DESCRIPTION
JOB INITIALIZATION (I Records)
The IN record is used to define time interval and starting time for time
series data which are read into the program on PC, PI, QO, QI, QS, MD, MS, MT
and MW records.
If an IN record is not used the time interval and starting time for all time
series will be the values specified on the IT record.
IN records may appear anywhere (exception: not after JD and before PI) in
the input stream. The same time interval and starting time will be used for
all time series data until these values are reset by reading new values on an
IN record.
When time series data are read from PC, PI, QO, QS, QP, MD, MS, MT, or MW
records, values to be used by the program are computed using linear interpola-
tion to match the tabulation interval specified on the IT record.
For times preceeding or following the given ordinates, the first or last
value is repeated as necessary to define NQ (IT-4) ordinates.
Data on PC, QI, QO, QP and QS records are instantaneous values. The first
value will occur at JXDATE and JXTIME.
Data on PI, MD, MS, MT and MW records are cumulative or average values
over a time interval. The first value on these records is for the time
interval beginning at JXDATE, JXTIME and ending at JXTIME + JIMIN.
A-1O
HEC-l INPUT DESCRIPTION 10
JOB INITIALIZATION (I Records) IM
The 10 record is used to control output for the entire job. The KO record
may be used to change output control for each hydrograph calculation.
A-11
JP
A-12
HEC-I INPUT DESCRIPTION
JR
JOB TYPE OPTION (J Records)
A-13
JD
HEC-l INPUT DESCRIPTION
JOB TYPE OPTION (J Records)
A-14
OuJ
HEC-l INPUT DESCRIPTION
OPTIMIZATION OPTION (0 Records) OR
Blank
A-15
OS HEC-1 INPUT DESCRIPTION
OPTIMIZATION OPTION (0 Records)
DO Diversion
SO Reservoir
WO Pump
LO Local protection projects and uniform degree of protection
A-16
OF
HEC-l INPUT DESCRIPTION
OPTIMIZATION OPTION (0 Records)
A-17
00 HEC-1 INPUT DESCRIPTION
OPTIMIZATION OPTION (0 Records)
A-18
vS
HEC-1 INPUT DESCRIPTION
USER-DEFINED OUTPUT TABLES (V Records)
VS and W records define tables which may be used to display selected time
series output. Each table may contain up to 10 columns of data as defined on
one pair of VS/TV records. Up to 5 tables may be output by using 5 successive
pairs of VS/W reecords.
A-19
VV
HEC-l INPUT DESCRIPTION
USER-DEFINED OUTPUT TABLES (V Records)
1. Observed flow
2. Calculated flow
3. Rainfall values
4. Rainfall loss values
5. Rainfall excess value
6. Storage values
7. Stage values
A-20
u - 1 1
BA
UUC-l XUPW DE3CRIPTION
BAINE 30V DATA (B ReC6946)
A-21
BASIN RNOF PAT& (B Records)
Base flow parameters (STRTQ, QRCSMN, and ITIOR) will be assumed equal to
zero unless this record is supplied. Once this record Is supplied, all
following subbasins will be assumed to have these values unless overriden by
another BF record.
NOTE - The definition of ETIOR has been changed from the old version of
HEC-1. The old value is QA/QI in the following equation:
A-22
OR
NEC-i IXPUJT DESCRIPTION
A-23
DIIOSTOD
Om (D leord)
A-24
DT
INC-i InPU DERIMION
DIVERSION DATA (D Records)
For multiplan simluatious UJP record), diverison data (DI and DQ records)
must be supplied for all plans. if no water is be diverted for a particular
plan, then the DQ record would contain only zero**. Diversicn hydrographo are
saved for all plans using the name in Field 1 of the 0T record.
A-25
NC-1 INPU DESCRIPTION
A-26
EEC-i INUTDESCRIPTIoD
DO
DIVERION DATA (D Records)
A-27
DD NEC-i INPUT DESCRIPTION
DIVERSION DATA (D Records)
A-28
HB
NEC-1 INPUT DESCRIPTION
HYDROGRAPH TRANSFORMATION (H Records)
A-29
* 9C
HL HYDROGRAPB TRANSFORMATION (H Records)
A-30
HO
EEC-3i INPUT DESCRLIPTION
Records)
H
myD&OGAPR TRASOETION (H
A-31
EHC-1 INPUT DSCRIPTION
KJOB STEP CONTROL (K Records)
**RBQUIRED
A-32
NEC-i INPUT DECRIPTION
KO
JOB STEP CONTROL (K Records)
Continued
A-33
4
OEI-i
INPUT DESCRIPTION
JOB STEP COWrToL (K Records)
A-34
!!i.
~K ElC-1 INPUT DESCRIPTION KF
JOB STEP CONTEOL (K Records)
Use this record to specify format for the hydrographs on unit 7.(See KO-4).
This format will be used until a new KF record is read. Default format is
(2DQ1,I6,918). KF record should not be used unless format is to be changed.
This file may be used to punch cards.
before writing.
A-35
KP EEC-i INPUT D19SCRIPTIOU
JOB STEP CONTRML (K Records)
KK
KP 1
KP 3
*Data for PLAN 2 is not provided and thus will be the same as
PLAN 1
A-36
r!ALU HEC-1 INPUT DESCRIPTION
LOSS RATE DATA (L Records)
One of four different rainfall loss rate procedures may be used for a
subbasin runoff computation. A different loss rate may be used for each
subbasin and/or plan. Snowmelt loss rate (N record) may be used in
conjunction with the exponential LE record) or uniform (LU record) loss rates.
A-37
LE . C-I INPUT DESCRIPTION
LOSS RATE DATA L Records)
A-38
WIIC'- IU DESCRIPTION
LM
LOSS RATE DATA (L Records)
A-39
IK-1 INPUT DESCRIPTION
LOSS RATE DATA (L Records)
*This factor should only be used for directly connected imperivous areas not
already accounted for in the curve number land use.
A-40
LH
HKC-l INPUT DESCRIPTION
LOSS RATS DATA L Records)
A-41
MA
HC-l INPUT DESCRIPTION
SNOWhELT DATA (N Records)
NOTE: Up to 10 records, one for each zone. Zones must be in equal elevation
increments corresponding to lapse rate coefficient TLAPS (M-1).
A-42
MC
HEC-1 INPUT DESCRIPTION
SNOWNELT DATA (H Records)
A-43
TEC-1 INPUT DESCRIPTION
These data are required for any snowfall/melt simulation. See IN record
description for discussion of time interval and number of values.
3 TEMPR(3) + Etc.
The MS, MD, and MW records are only used for the energy budget snowmelt
simulation. See IN record description for discussion of time interval and
number of values.
3 SOL(3) + Etc.
A-44
~MD
HEC-i INPUT DESCRIPTION
SNOWNELT DATA (N Records) MW
See HS record.
3 DEWPT(3) + Etc.
See MS record.
3 WIND(3) + Etc.
A-45
iiC-l INPUT DESCRIPTION
PRECIPITATION DATA (P Records)
ID
IT Etc., for job initialization
U-
KK Etc.
PG and PG + PI/PC record combinations can be included at any point in the data
set following the IT record. It is usually convenient to group them together
as a precipitation data bank before the first KK record. Different storms can
then be simulated by simply inserting different data banks, as long as the
gage identification and weightings are the sae.
Continued
A--46
HC-1 IIPUT DSCRIPTION
ID
IT
KK
KK
A-A7
EEC-i INPUT DESCRIPTION
_V PRECIPITATION DATA (P Records)
A-48
PC
2 PRCPR(2) + Etc.
4 PRCPR(4) + Etc.
A-49
PGNC-i INPUT DESCRIPTION
PRECIPITATION DATA (P Records)
All precipitation gages are total-storm stations. Some stations may also
have temporal distributions associated with the storm-total precipitation.
These stations are also called recording stations when referring to the
temporal pattern. The temporal distribution is defined on PI or PC records
imediately following a PG record.
PR, PT and PW records are used within each KK, BA, etc., record series to
specify weightings of precipitation station data to compute the subbasin-
average precipitation distribution.
A-50
I a
mnmm ma =nn
u I l~ll • nun M ml
PH
H-C-l INPUT DESCRIPTION
PRECIPITATION DATA (P Records)
The duration of the storm will be the duration for the last non-zero depth
which is specified. The first non-zero depth specified will be the most
intense portion of the storm. Depths must be specified for all durations
between these limits.
Continued
A-51
PH1
HEC-1 INPUT DESCRIPTION
" w PRECIPITATION DATA (P Records)
A-52
EEC-1 INPUT DESCRIPTION PM
PRECIPITATION DATA (P Records)
NOTE: Hydrometeorological Report No. 33 has been superseded by HMR No. 51 and
No. 52. Computer program HMR52 (NEC, 1984) may be used to calculate PES
hyetographs.
1000 10 .90
500 10 .90
200 10 .89
100 13 .87
50 15 .85
10 OR LESS 20 .80
Continued
A-53
PM IEC-i INPUT DESCRIPTION
PRECIPITATION DATA (P Records)
Duration of the computed PMS will correspond to the last non-zero percentage
entered. Minimum duration is 24 hours.
A-54
HEC-i INPUT DESCRIPTION
PS
PRECIPITATION DATA (P Records)
A-55
PR EC-1 INPUT DESCRIPTION
PRECIPITATION DATA (P Records)
These records are used to identify the gages and their relative weightings
for computing this subbasin's average precipitation.
A-56
HEC-1 INPUT DESCRIPTION PT
PRECIPITATION
DATA (P Records)
A-57
00 HEC-1 INPUT DESCRIPTION
HYDROGRAPH TIME-SERIES DATA (Q Records)
These records contain hydrograph time series data. The first value on the
record is at the starting time specified on the previous IN record.
Subsequent values are spaced at the time interval specified on the IN record.
The program reads all consecutive Q records and interpolates values for the
computation time interval and time period specified on the IT record. If the
computation time period extend! before or beyond the Q data supplied, the
first or last value will be repeated as necessary to produce a bydrograph for
the full time period.
2 QO(2) + Etc.
A-58
HEC-l INPUT DESCRIPTION QI
HYDROGRAPH TIIK-SERIES DATA (Q Records)
2 QI(2) + Etc.
These records are used to input a stage hydrograph for comparison with the
computed hydrograph. A rating table, on HQ and HE RECORDS, must also be
supplied. Comparison of hydrographs is a distinct operation which must be
separated from other operations with a KK RECORD.
2 QS(2) + Etc.
A-59
ofQPl NU DESCRIPTIONecrS
L HYDROGRAPH TIM-SZRIZ3 DATA (Q Records)
These records are used to input a pattern hydrograph for which local
2 QP(2) + Etc.
A-60
RN
HEC-1 INPUT DESCRIPTION
ROUTING DATA (R Records) R L
A-61
RM
HEC-1 INPUT DESCRIPTION
ROUTING DATA (R Records)
1 (AHSKK*60.) 1
A-62
SEC-i INPUT DESCRIPTION R S
ROUTING DATA (R Records)
Continued
A-63
tR
HRC-l INPUT DESCRIPTION
ROUTING DATA (R Records)
A-64
REC-i INPUT DESCRIPTION
RC
ROUTING DATA (R Records)
A-65
RX EEC-i INPUT DESCRIPTION
ROUTING DATA (R Records)
Left bank and right bank of channel are assumed to be located at points 3
and 6, respectively, of the cross section.
A-66
HEC-l INPUT DESCRIPTION RY
ROUTING DATA (2 Records)
Left bank and right bank of channel are assumed to be located at points 3
and 6, respectively, of the cross section.
&-67
RK HC-l INPUT DESCRIPTION
ROUTING DATA (R Records)
3 N + Channel roughness.
A-68
HEC-l INPUT DESCRIPTION
RT
ROUTING DATA (R Records)
NOTE - The variables used for this routing method are dependent on the
computation time interval. The user should make proper adjustments when using
different time intervals.
0 LAG method.
0 Tatum
A-69
HEC-i INPUT DESCRIPTION
STORAGE ROUTING DATA (S Records)
S records are used to provide storage and outflow data for storage routing.
1. Discharge on SQ records
2. Weir and orifice data on SS and SL records
3. Ogee spillway data on SL, SS, SG, SQ, and SE records
When spillway data (weir or ogee) are provided, the program computes a
steady flow rating curve, then interpolates from that rating curve during the
routing calculation. Elevation data may be input for storage or outflow by
following SV or SQ records with SE records.
A--70
SV
HEC-I INPUT DESCRIPTION
SA
STORAGE ROUTING DATA (S Records)
These records are to be used if the reservoir volumes are known. Do not
use if SA records are supplied.
These records are used if the reservoir volumes (SV record) are not
known. Do not use if SV records are supplied.
A- 1
SE
HEC-1 INPUT DESCRIPTION
s Q STORAGE ROUTING DATA (S Records)
The SQ record gives outflow data for storage routing. Values should
correspond to storage data, or if elevation data are provided for both storage
and outflow, the program will interpolate discharges for the given storages.
The SQ and SE records are also used to specify tailwater data for the ogee
spillway option.
A-72
HEC-1 INPUT DESCRIPTION S L
STORAGE ROUTING DATA (S Records)
A-73
HEC-1 INPUT DESCRIPTION
$S STORAGE ROUTING DATA (S Records)
This record is used to compute flow for weir or ogee spillways. If the
dam overtopping summary is requested (ST record), the spillway crest elevation
should be provided on this record.
A-74
.EC-1 INPUT DESCRIPTION ST
STORAGE ROUTING DATA (S Records)
This record is used to compute flow over the top of a dam. Flow computed
using the weir coefficients specified on this record is added to outflow
computed from the spillway (SQ, SS, SL, or SG records). Use of this record
calls for the dam overtopping sumary (spillway crest elevation should be
provided on SS record). This record is required if the non-level top-of-dam
option (SW/SE records) to used. The discharge over the top of dam is added to
the discharge elevation relationship generated by the program (SL, SS, SG
options) or specified by the user (SQ, SE option).
A-75
SW NEC-1 INPUT DESCRIPTION
SE STORA E ROUTING DATA (S Records)
A-76
HEC-l INPUT DESCRIPTION
STORAGE ROUTING DATA (S Records)
SG
Continued
A-77
so HiEC-1 INPUT DESCRIPTION
NOTE- SQ and SE records to define the tailwater must follow this SC record.
If a low-level outlet is specified, it should precede the SC record to
prevent error message.
'A- 78
\
SB
HEC-1 INPUT DESCRIPTION
STORAGE ROUTING DATA (S Records)
NOTE - Tables and plots of dam-breach hydrographs for each plan are generated
automatically when IPRNT (10-1 or KO-1) is less than 4. Those tables and
plots show how well the breach hydrograph is represented by the normal time
interval specified on the IT record.
A-79
so KKIIC-l INPUT DESCRIPTION
STORAGE ROUTING DATA (S Records)
A-80
USC-1 INPUT DESCRIPTION
STORAGE ROUTING DATA (S Records)
SD
A-8i
U I HEC-1 INPUT DESCRIPTION
UNIT GRAPH/KINEMATIC DATA (U Records)
The given unit hydrograph must have been derived for the time interval on
the IT record (IT-l, IT-2). For example, if the time interval is 15 minutes,
then a 15-minute unit hydrograph must be used.
A--82
I l miIm II iI m n
HEC-1 INPUT DESCRIPTION
UC
UNIT GRAPH/KINEMATIC DATA (U Records)
A-83
US HEC-1 INPUT DESCRIPTION
UNIT GRAPH/KINEKATIC DATA (U Records)
A-84
II I I I I I
EEC-1 INPUT DESCRIPTION UA
UNIT GtUM/KIN&TIC DATA (U Records)
This time-area data may be used with either the Clark or Snyder methods.
This data may be in any units, since area is scaled to the subbasin area and
time is scaled to time of concentration. The areas contribute to runoff at
the basin outlet at equally spaced time intervals. A synthetic time-area
curve will be used if the UA record Is not supplied.
A-ms
HXC-i INPUT DESCRIPTION
4 . UNIT OWNP/KINCHATIC DATA (U Records)
A-86
7, RU-i INPUT DESCRIPTION
UNIT GRAPH/gINSHATIC DATA (U Records)
UK
The first and second loss rates specified on a previous L record will be
used for the first and second UK subcatchments, respectively.
A-8 7
UNIT GAHKNKTCDT URcrs
A-88
WP
MC-1 inPU DgSCR1PTIOM
PUtRI DATA (W Records)
The program checks the elevation at the end of the previous time interval
to see if a pump should be turned on or off. The use of the WP record with
the multiplan capability requires some special conventions. A single WP
record with a non-zero (can be set very small) pump flow is required (PIJPQ,
field 2) for Plan 1. All other plans (Plan 2, 3, etc.) must specify first a
HP record with zero PUHPQ and them a second WP record with the desired pumping
rate; for example:
Feld 1 2 3 4
HP 0
WP 843.5 3000 PNPQ1
A-89
WW ~ EEC-i ImPU DES8CRIPTION
A-90
N'C-i INPUT DISCIIPTIOU WO
PUMhP DATA (V Records)
A-91
W C
NEC-i INPUT DESCRIPTION
WD PUMP DATA (W Records)
A-92
HEC-l INPUT DESCRIPTION
ECONOMIC DATA
18 ECONOMIC DATA
Data for economic evaluation of flood damage is placed in the data set
following the last hydrograph calculation and before the ZZ record. The first
record in the economic data is an EC record, and all records between the EC
and ZZ records are economic-data records.
H
PN * Plan names
FR Frequency data
DC Damage data
*Optional records
A-93
E C HEC-1 INPUT DESCRIPTION
ECONOMIC DATA
CN
A-94
HEC-1 INPUT DESCRIPTION
PN
ECONOMIC DATA
This record is used for description of the plans. One record is used for
vach plan. A maximum of 5 plans (PN records) may be used.
A-95
WN
TN
HEC-l INPUT DESCRIPTION
ECONOMIC DATA
See WN record.
A-96
KK
HEC-1 INPUT DESCRIPTION
JOB STEP CONTROL (K Records)
A-97
WT
FR HEC-1 INPUT DESCRIPTION
ECONOMIC DATA
This record is used to identify the watershed and township for the stream
station given on the KK record. Watershed and township designations will be
the same for all stations until a new WT record is read.
This record is required for the first station. These frequency values
will be used until changed by a new FR record.
** REQUIRED
A-98
QF
HEC-1 INPUT DESCRIPTION SF
ECONONIC DATA
1 Not used.
2 Not used.
This record should be used only if peak stage have been calculated in the
hydrologic portion of HEC-1. This record is required for each station if QF
record is not provided.
1 Not used.
2 Not used.
A-99
SQ HEC-1 INPUT DESCRIPTION
QS ECONOMIC DATA
1 Not used.
1 Not used.
2 Not used.
A-100
.EC-i INPUT DESCRIPTION SD
ECONOMIC DATA OD
Do not use this record if flow-damage data are to be used. Provide one SD
record for each station. If stage-damage data change for each plan, a new SD
record must be provided for each plan.
1 Not used.
1 Not used.
A-101
HEC-1 INPUT DESCRIPTION
DG ECONOMIC DATA
Damage data must be provided for each station. One (two if NDMG is
greater than 8) record is required for each damage category.
1 Not used.
**REQUIRED
A-102
EP
HEC-l INPUT DESCRIPTION
ECONOMIC DATA
This record is required to indicate the end of data for a plan. The
current plan will be evaluated and new data will be read for the next plan.
If there are no additional data, the last data set read will be used to
compute expected annual damages for any plan which has not been evaluated.
- Labels such as Plan Name (PN) and Damage Category Name (CN) can be
changed. Plan Names could be specified for all plans in the first group
of data (for the first plan).
A-103
LO HEC-1 INPUT DESCRIPTION
ECONOMIC DATA
A-104
LC
HEC-l INPUT DESCRIPTION
ECONOMIC DATA LD
A-105
D U HEC-l INPUT DESCRIPTION
ECONOMIC DATA
DL
Pattern damage table for minimum design level (XLPKN) for local protection
project.
I Not used.
Pattern damage table for maximum design level (XLPMX) for local protection
project.
1 Not used.
A-106
HEC-1 INPUT DESCRIPTION
ECONOMIC DATA
D P
A-107
Z Z HEC-1 INPUT DESCRIPTION
END-OF-JOB RECORD (** ZZ RECORD)
This record identifies the end of an HEC-1 job and causes summary
computations and printout to occur. Another job may be started with another
ID, IT, etc., record series if desired. If another job does not follow, the
control is passed back to the computer operating system.
**REQUIRED
A-108
HEC-1 INPUT DESCRIPTION
HEC-1 INPUT RECORD SUMMARY
*LIST A-7
*NOLIST A-7
*FIX A-7
*FREE A-7
* (comment beginning in column 3 A-7
*DIAGRAM A-7
ID (TITLE) A-8
IT NMIN IDAIK ITIME NQ NDDATE NDTIME A-9
IN JXMIN JXDATE JITIME A-10
In A-11
10 IPRT IPLT QSCAL A-I1
JP NPLAN A-12
JR IRTIO RTIO A-13
JD STRM TRDA A-14
VS ISTA A-19
VV SMVAR A-20
DR ISTAD A-24
DT ISTAD DSTRMX A-25
DI DINFLO . A-26
DQ DIVFLO A-26
DO IOPTD DANCST DDSCNT DVRMX DVRMN A-27
DC DCAP A-28
DD DCST A-28
Continued
A-109
HEC-1 INPUT DESCRIPTION
HEC-1 INPUT RECORD SUMMARY
PB STORM A-48
PI PRCPR A-49
PC PRCPR A-49
PG ISTAN PRCPN ANAPN ISTANX A-50
PH PFREQ TRSDA PNHR .A-51
PM PMS TRSPC TRSDA SWD R6 R12 R24 R72 R96 A-53
PS SPFE TRSPC TRSDA SWD A-55
PR ISTR A-56
PT ISTN A-57
PW WTR A-57
QO QO A-58
QI Q A-59
Qs Qs A-59
QP QP A-60
RN A-61
RL QLOSS CLOSS A-61
RM NSTPS AMSKK X A-62
RS NSTPS ITYP RSVRIC x A-63
RC ANL ANCH ANR RLNTH SEL ELMAX A-65
RX x A-66
RY Y A-67
RK L S N --- SHAPE WD Z A-68
RT NSTPS NSTDL LAG A-69
Continued
A-110
HEC-l INPUT DESCRIPTION
HEC-l INPUT RECORD SUMMARY
SE ELEV ... A- 72
SQ DISQ ... A-72
SL A-73
SS CREL SPWID COQW EIPW A-74
ST TOPEL DAIIWID COQD EXPD A- 75
SW WIDTH A- 76
SE ELVW ... A-76
SO IABCOA ISPITW ISPCTW lIGATES SS DESHD APEL APWID APLOSS PDPTH A-77
SB ELBM BRWID Z TFAIL FAIL.EL A- 79
SO IOPTR RANCST RDSCNT CAPHI CAPHN A-80
SD RCST ... A-81
UI QUNGE . A-82
UC TC R A-83
us TP CP A-84
UA QCLK .. . A-85
UD TP A-86
UK L S N A A-87
RK L S N CA SHAPE WD Z UPSTQ A-88
aC A-94
CN NCAT NECAT ... A-94
PN IPLN NMPLN A-95
UN NUAT WID .
- A-96
TN NWN TID . . .A-96
KK ISTAQ NAME .. . A-9 7
WT IWAT ITWN A-98
FR NFRQ PFREQ . . . A-98
QF --- --- QFRQ A- 99
SF - --- SFRQ ... A-99
SQ -- NSTG STGQ ... A-100
QS QSTG . . .A-100
SD -- NDNO SDMG . . .A-101
QD --- NDNG QDMG A-101
DG --- IPLN DAMG DANG A-102
EP k-103
LO IDPTLP XANCST XDSCNT LPMX ILPEN A-104
LC ILCAP A-105
LD XLCST A-105
DU ICAT TUDAMO A-106
DL ICAT TLDAMG A-106
DP DOPRT TRGT A-107
ZZ A-1O8
A-1l1
Appendix B
KC-i USAGR WITH MC DATA STORAGE SUT¥
3.1 Introduction
The HEC Data Storage System (DSS) (HKC, 1983) has been developed to allow
transfer of data between NC programs. The data are identified by unique
labels called PATiUANEKs which are specified when the data are created or
retrieved. Thus, a hydrograph computed by HEC-1 can be Labeled and stored in
DeS for later retrieval as input data to HIC-5, for instance. The D88 has
several utility programs for manipulating data. These programs enable editing
of information, changing pathnames, purging unwanted data sets and insertion of
other data sets. Graphic and tabular portrayal of DSS data are also available.
The pathname is separated into six different parts by a slash "/" delimiter
so that each part refers to a specific, unique identifier. One convention
that has been developed to simplify definition of psthname parts for typical
hydrologic data is shown below: .
In general, DSS software finds the data associated with a pathname by using
each of the six parts to search the DSS file structure, which is hierarchical,
or "tree-like." An example of a pathname for a time-series data record is:
/WISSISSIPPI/CAIRO/STAGK/01JAN85/HOUR/OBSERVED/
HZC-i can interact with DSS as follows: retrieve runoff parameters stored
in DSS by program HYDPAR (Corps of Engineers, 1978); retrieve and/or store
time-series data; and store flow-frequency curves. The access to this data is
B-1
...
accomplished using the BZ, ZR and ZW records in the HzC-1 input data set.
Table B.1
ID
IT
I0
JP (required for multiplan simulation)
JR (required for multiratio simulation)
KK
KP (only required if multiplan simulation)
ZR
BZ
L (only required if Snyder unit graph is used)
KP (only required if multiplan simulation)
ZR
BZ
ZZ
The BZ and ZR records can be used in either fixed or free format modes
independent o. the input mode for the rest of the data. As an example of the
BZ and ZR record formats, consider the pathname,
A A a [
/IISSISSIPPI/CAIRO///1985/PLAN A/
5-2
Bz
the 32 and ZRrcord would then have the following fixed form
1 13TAuCAIO (Part 3)
DZ D-CAIO
ZR AmUISSISSIPPI 3-1985 F-PLAN A
Note, that only parts A, K and F are entered on the ZR record. Fattines part 9
is entered on the DZ record and parts C and D are left blank.
The format and content of the Il and ZR records for UYDPAR parameter
retrieval in fixed format are as follows:
* A-Alphanumeric data
3-3
-ad the iaput for the free format is as follow:
0 ID 8Z Record identification.
0 - ZR Record identification.
The time-series data that can be retrieved with the ZR record are cumula-
tive or incremental precipitation end discharge hydrographa, corresponding to
data which can be specified on PC, PI, Q1 or QO records. The record sequence
needed to perform this operation is shown in Table 5.2. This option is useful
in either stream network or multiplan-iultiratio simulations.
Table 7.2
ID
IT
I0
jP (required for mwltiplan simulation)
JR (required for multiratio simulation)
3-4
ZR'
YPtbame parts D and I are not used. The program uses information on the
3? record in the place of information normally specified with parts D and 3.
A an eample application, consider the pathnme needed to retrieve an
observed
+ bydrograph: A I t
/MISSISSIPPI/CAIROFLOW/ /1OS8/
where all the pathname part descriptors and the type of time-series data is
specified by the "*Qo". Note that the additional parameter "-aa", must set
the value "as" equal to PC, PI, Q1 or QO to indicate the type of time-series
data.
0 ID ZR Record identification.
8-5
S 3ahl Tim interval of 053 data (e.g., LSMIE)
computation interval spec ifie on IT record. Must
be a standard DII tims interval; SmiN, la1mi,
Flow, storage or stage time-series data may be stored in D83 using the ZV
record. Th ZN cowvention is similar to the use of the ZR record (see Table
5.2:). Wing the previous eameple for the ZR record, the ZM record specifies
tbo petbom as:
Zu DUCAIRO cgs=O
ZN B=CAI30 C3AM
Mote that parts A and F need not be repeated. If part B were not use&, then
the station name oft the IKK record would be used for location now..
As in the case of the ZR record, the ZU data my be used in the free field
format mode. However, the application of the ZW record differs slightly in
that f or each plani all ratios of the computed: time series awe saved (as
opposed to a single time-series trace for the ZR record). The pathname part F
need not be repeated for each plan, as the program, automatically assumes the
description given for plan 1. As in the case of the ZR record, a NP record
must be used with each Z% raerd' for each plan.
B-6
,Z W
Table 8.3
ILD
IT
10
JRP
9C
UZ
3-7
. M A - Flow Frequency Curves to DSS (Fixed Format Option)
-Writing
YW-
aiale V*W DescritnO
0 - ZW Record identification.
0 ZW Record identification.
This record must always follow a PU record in the economic data. The
conventions for specifying this record are analogous to the reading of HYDPAR
data.
8-8
ame fle, 29, A-22 exmle, 169
butu.-mvsep precipitation, S. Be. at, ASmodification. 64
Oifuraited cisels (in Divarsions)
mmulInfltrtion
31 ~lGage precipitation, 9, 80, A-50
Channel storage (no routine) Holtan infiltration, 19, A-41
Clark unit grph. 20, A40 Hyrograhs, 6
g exe 90. 120 analysis, 43
Cmnfluences (cmining),* 6, A-30 balanced, 7, A-29
muter-rquarumts, 3. M8 input. al. A-"9
Costs for projects, 69, A-IS, A-SI, A-92, A-105 ratios, A-21 (sm Ikultiflood too)
exampe. 116 recession of, 30, A-Il
Cross sections, 33, A-66 trnfraion, 1
Cum nuer (SCS) (sm Intercept ion/i nfi ltration) Hyeto
user input. a1. A-*S
DOw break, 56, A-i9 synthetic, 10, 81, A-51
exmle, 138
Do safety, 53, 88, A-74 Initial and uniform loss, 16, A-37
breach, 56, A-79 Initial flow, 29, A-42
examles, 133, 138 ample, 90
non-level crest, 54, A-76 Input data overview. 72, A-2
overtopping, 53, A-75 Intercenption/infilitration, 15
simulation limitations, 51 exampes, 90
Dowe (see economics) WEC eponential, 16, A-US
Damage reach, 63 Holtan, 19. A-41
Data repetition, 76 initial and uniform, 16, A-Si
Depth-4uration precipitation, 12 parameter estimation (optimization), 43, A-15
Depth-area precipitation, 55, A-14 SCS curve mber, 18, A-40
exmle, 129
Design storm, 10, A-S1 Kinematic wove, 23
exmle, 136 exoele, 106
Diagrami (stream network), 65, A-i land-surface rnoff, 23, A-Si
exmle, 96, 155, 162,172, 179 routing, 41, A-IS
Direct runoff (see Land-surface rnoff)
Diversion, 7, 42, A-24 Land-surface runoff, 5, 51
exmle, 90, 158, 159 kinemtic wave, 23, A-Si
optimization, 68 uni t graph, 20. A-82
Limitations of propra, *2, 57
Economics (flood daeig), 62, 83, A-102 Local protection, 69, A-104
exampe, 169 Loss rate (we Intercuption/infiltration)
Error massages, 86
Evapotranspiration, 15 Mlanning's equation, 33
Exceaence frequency. 64, A-98 Malt (me Snowmalt)
Excess precipitation, 16 Hodif led Puls, 32, A-03
Expected annual flood dewmage, 66 Muskingum, 32, A-62
exmle, 169 Multiflood, S1, A-13
eample, 153
Flood damage, 62, 65, 83. A-101 flultiratio (see Nultiflood)
Floodplain storage (see Routing) Multiplan, SI, A-12
Flood Routing (we Routing) ex&Vle, 158
Flow target optimization, 11, 116, A-107
Frequency cves, 64, A-98 Norml-dpth rauting, 33, A-IS
example, 13
Snowfall 13
)
am reet"Snommlt. 13, A-42
pumtes 4, 42 A-15 dos dy 14, A4
error anlysis, 49 exmle, 114
Gamin., 120, 126 energ budget, 14, A-44
fnwd cotrol systu size, 6S, 69, 84, X-16 Snyder unit graph, 21, A4W
SUN*% 176 ezampe, 90
Otput. as Soil Conservation Service
error Messages, IS exampe, 90
doovertopping iry, 137, A-74 curve numer, 18, A-40
dam-breech anlysis, 151, A-79 unit graph, 22, A-SB
user-desipned tabl, 152, A-19 Soil moisture (see interceptian/infiltratian)
Ovrland neeo (soe Land-surface runoff) Spatial data (SM) usage, A-23
Spi I)ways
Parmeter calibration, 43 weir, 37, 54, X-94
routing, 49, A-15 exmle, 133
selection process, 47 traezoidul. 37. A-77
unit graph and loss rate, 43, A-15 agpt, 37, A-17
Percolation, 15. 19 Stage hydrograph usage. A-59
Precipitation, 8, 81 Standard Project Storm, 10, AS55
gage wightings, 9, 80, A-SS Storage-outflow routing (see Routing)
historical staorm, 9,80 Storage of reservir, 37, A-?Il
depth-area, 59. A-14 Staous (see Precipitation)
examle, 129 Straddle-stagger (see Routing)
depth-duration. 13, A-Si Strea network model
input data, 81. A-*6 definition, 4, A-32
partial duration, 13 exmle data set, 79
synthetic, 10. A-51, A-S3, A-55 exampe, 90
Prdibi1e Maximm Storm, 12, A-S3 Clark, 20, A-03
Pumps. 7, 42, A-0 parmter estimetion (optimization)
optimization. 68 43, 82, 120, A-15
mamles, 156, 169, 176 SCS, 22, A-US
Snyder, 21. A-4
Rainfall (see Precipitation)
Rating curve use. A-31 Tatum (see Routing)
Recession flow, 29, A-22 Tamerature (see SnwAmelt)
Rleservoi r, 6 Thiessen weighting (se Precipitation-gage weighting)
low-level outlet, 37, A-73 Time-series input data, 76, A-10
optimization. UI Time-area diagram, 21, A-US
example, 176 Transposition coefficient, 12
routing, 32, 35, A-7l
spillway, 37, A-74 Uniform level of protection, 69, A-l07
storage, 36. A-7l Unit hydrogrqahs, 20
Routing, 6, 32, 92 example, 90
average-ag, 41, 90, A-69 Clark, 21, A4U
given storage vs. outflow, 33, 95, A-63 parmter estimetian (optimiization)
kinematic wave, 41, 106, A-68 43. 62, 120, A-IS
level-pool reservoir, 35, 1W3 A-71 SCS. 22, A-US
Nodified Puls. 32. 95, A-63 Snyder, 21, A-94
Muskingum, 32, 126, A-62 User-designed output, 152. A-19
not Il depth, 33, 136, A-ES
parmeter estimtion examle, 126 Working RIO (see Rtutig
resevoIrs, 32, 35
streddle-stagge, 41, A-69 Zones for snumlt, i3, A-42
Tat, 41, A-49 examle, 114
working RO, 35, A-63