SPOKANE RIVER BASIN MODEL PROJECT
Volume IV - User's Manual for Steady-State Stream Model
by
John L. Shepherd
E. John Finnemore, Ph.D.
Systems Control, Inc., Palo Alto, California
for the
ENVIRONMENTAL PROTECTION AGENCY
Contract No. 68-01-0756
October 1974
-------
EPA Review Notice
This report has been reviewed by the Environmental Protection Agency
and approved for publication. Approval does not signify that the
contents necessarily reflect the views and policies of the Environ-
mental Protection Agency, nor does mention of trade names or commer-
cial products constitute endorsement or recommendation for use.
ii
-------
ABSTRACT
Three existing mathematical models, capable of representing water quality
in rivers and lakes, have been modified and adapted to the Spokane River
Basin in Washington and Idaho. The resulting models were named the Steady-
state Stream Model, the Dynamic Stream Model, and the Stratified Reservoir
Model. They are capable of predicting water quality levels resulting from
alternative basinwide wastewater management schemes, and are designed to
assist EPA, State, and local planning organizations to evaluate water qual-
ity management strategies and to establish priorities and schedules for
investments in abatement facilities in the basin.
Physical data and historical hydrologic, water quality and meteorologic
data were collected, assessed and used for the model calibrations and
verifications.
The modified models are all capable of simulating the behavior of various
subsets of up to sixteen different water quality constituents. Sensitivity
analyses were conducted with all three models to determine the relative
importance of a number of individual model parameters.
The models were provided to the EPA as computer source card decks in
FORTRAN IV language, with accompanying data decks. All development work
on, and applications made with, these models were fully documented so as
to permit their easy utilization and duplication of historical simulations
by other potential users. A user's manual with a complete program listing
was prepared for each model.
This report was submitted in fulfillment of Contract No. 68-01-0756 under
the sponsorship of the Environmental Protection Agency.
The titles and identifying numbers of the final report volumes are:
Title EPA Report No.
SPOKANE RIVER BASIN MODEL PROJECT DOC /74
Volume I - Final Report
SPOKANE RIVER BASIN MODEL PROJECT DOC /74
Volume II - Data Report
SPOKANE RIVER BASIN MODEL PROJECT DOC /74
Volume III - Verification Report
SPOKANE RIVER BASIN MODEL PROJECT DOC /74
Volume IV - User's Manual for Steady-state Stream Model
SPOKANE RIVER BASIN MODEL PROJECT DOC /74
Volume V - User's Manual for Dynamic Stream Model
SPOKANE RIVER BASIN MODEL PROJECT DOC /74
Volume VI - User's Manual for Stratified Reservoir Model
111
-------
CONTENTS
SECTION PAGE
I. INTRODUCTION 1
II. MODE OF OPERATION 3
Purpose 3
Restrictions 3
Solution Technique 4
Model Usage 4
Model Verification and Sensitivity 11
III. PROGRAM DESCRIPTION 13
Main Program 13
Subroutine BLEND 13
Subroutine CMIN 13
Subroutine DOEQU 13
Subroutine DOSAT 15
Subroutine GETCON 15
Subroutine K2CAL 15
Subroutine NEWIN 15
Subroutine PLOT 15
Subroutine PLTSET 15
Subroutine REINI 20
Subroutine REPRE 22
Subroutine RMATC 22
Subroutine RUNON 22
Subroutine SCAN 22
Subroutine SETOPT 23
Subroutine TRIED 23
IV. INPUT REQUIREMENTS 25
V. OUTPUT DESCRIPTION 57
VI. DEFINITION OF PROGRAM VARIABLES 59
VII. SAMPLE INPUT DECK 75
VIII. SAMPLE OUTPUT 81
IX. PROGRAM LISTING 105
X. REFERENCES 141
XI. ABBREVIATIONS 143
v
-------
FIGURES
NO. PAGE
1. Diagram of a Typical River System 5
2. General Program Flow Diagram 14
3. Flowchart - Subroutine GETCON 16
4. Flowchart - Subroutine NEWIN 21
5. Flowchart - Subroutine SETOPT 24
6. Input Data Deck Organization 26
vi
-------
TABLES
NO. PAGE
1. Coefficients for Computation of the Reaeration
Coefficient 8
2. DOSCI Input 27
3. Definition of Constituent Selection Option, ICOMB ... 56
4. Definition of Common Variables 59
5. Definition of Local Variables 70
6. Example Input Data Deck 76
7. Example Output 82
vn
-------
SECTION I
INTRODUCTION
The Steady-State Stream Model (DOSCI), documented herein, is an exten-
sive modification of the program DOSAG-I [Reference 8], and is one of
three water quality simulation programs developed by Systems Control,
Inc. (SCI) for use on the Spokane River Basin in the states of
Washington and Idaho. Portions of DOSAG-I, as developed by the Texas
Water Development Board in 1970, remain intact in DOSCI. Subroutines
which are unchanged or only slightly modified include BLEND, CMIN,
DOSAT, K2CAL, REINE, RMATC, SCAN and TRIED. DOEQU, REPRE and RUNON
have been extensively modified. Five newly added subroutines are
named GETCON, NEWIN, PLOT, PLTSET, and SETOPT.
These modifications have given DOSCI the capability of modeling as many
as sixteen quality constituents simultaneously. The algorithms used to
model the quality constituents are described in detail in Volume I (Part
III) of this report.
The format of the additional input required to execute DOSCI has been
made consistent with the format of the input required to run DOSAG-I.
Many input data cards required by DOSAG-I have been made optional by
the incorporation of default values for the variables involved. Por-
tions of this manual are taken from Reference [8], and the format of
this manual has been made generally consistent with the format of
Reference [8], subject to the constraints imposed by EPA documentation
specifications.
DOSCI has been extensively tested on a wide variety of stream condi-
tions on both the UNIVAC 1108 and the IBM 370/155. Run time is
minimal and the chief use of the program is to give a rapid evaluation
of a given set of stream conditions.
-------
SECTION II
MODE OF OPERATION
PURPOSE
The purpose of this model is to calculate the concentrations of as many
as sixteen water quality constituents in a particular stream system,
including the biochemical oxygen demand and the dissolved oxygen con-
centration. If desired, the minimum dissolved oxygen concentration in
the stream system may be checked against a pre-specified target level
dissolved oxygen concentration. If the minimum dissolved oxygen level
is below the target dissolved oxygen level, the program will compute
the required amount of flow augmentation to bring the dissolved oxygen
level up to the target level in the entire system. The user of the
program specifys the locations within the stream system at which dilu-
tion water is available for flow augmentation. The program is designed
to be run for varying climatic and hydrologic conditions during a twelve
month period. Thus, it is possible to enter up to twelve different
temperatures and corresponding discharges to each of the headwaters
within the stream system being modeled.
The Steady-State Stream Model (DOSCI) is quite flexible and can be
adapted to any stream system. One major restriction is that large
impoundments such as reservoirs cannot be considered by this program.
The output from a single run of the DOSCI dissolved oxygen simulation
model will provide the user with a complete description of the dissolved
oxygen resources of the stream system investigated, and the required
dilution water needed to bring the system up to the target level
dissolved oxygen concentration. An additional user option available
with this program is the ability to find the dissolved oxygen distribu-
tions for varying levels of treatment (waste treatment plants) in the
simulated river basin.
RESTRICTIONS
1. Maximum number of headwater stretches = 10, and headwater
stretch // must be <_ 10
2. Maximum number of junctions = 20; jets should be numbered in in-
creasing order downstream to decrease program run time and output
3. Maximum total number of reaches = 100; maximum number of reaches
per stretch = 20
4. Maximum number of stretches = 20
5. Maximum of twelve months of routing for temperature and headwater
flows; a minimum of one month must be used.
6. Maximum number of dissolved oxygen targets = 4. A minimum of one
dissolved oxygen target must be specified. This dissolved oxygen
target may be entered as a negative number if no flow augmentation
is desired.
7. Maximum of four degrees of waste treatment.
-------
SOLUTION TECHNIQUE
The algorithms used by DOSCI in the calculations of the various constit-
uent concentrations in a stream system are described in detail in Part
III of Volume I of this report. Briefly, given a stream reach, the
procedure is as follows. All inflows to the head of the reach are mixed
together and the resulting concentration for each constituent being
modeled is calculated. In the simplest case the inflow to a reach is the
outflow from the upstream reach. A reach may have as many as three in-
flows, however, as is the case for a reach located immediately below a
junction and with a tributary inflow. In any case the resulting mixture
is allowed to react for a time At according to the quality equations
where At is the travel time through the reach. The process is then
repeated for the next downstream reach. This process continues until the
entire stream system has been modeled.
MODEL USAGE
In order to use the DOSCI program, the user must take the stream system
which he proposes to simulate and break it down into the elements which
are used as input to the program. Figure 1 shows a schematic diagram of
a typical river system which has been decomposed into the elements
required to model it using DOSCI. There are essentially four major
elements into which a system must be decomposed so that it can be modeled
using this program. These elements are:
1. junctions - the confluence between two streams within the river
basin being modeled,
2. stretches - the length of a river between junctions,
3. headwater stretches - the length of a river from its headwater to
its first junction with another stream, and
'4. reaches - the subunits which comprise a stretch (headwater or
normal).
A new reach is designated at any point in the stretch where there is a
significant change in the hydraulic, biologic, or physical characteris-
tics of the channel, including the addition of a waste load, or the
withdrawal of water from the stream.
After a stream has been represented schematically, it is necessary to
specify the hydraulic and physical characteristics of each reach in the
stream system. This step involves reading into the program various
coefficients to describe all of the factors which are involved in the
various reactions occurring among the quality constituents being modeled
in the reach. It should be noted that this is some of the most important
information required for the simulation process. The results obtained
from the simulation of the dissolved oxygen resources within a river
system, using DOSCI as a modeling medium, are only as accurate as the
input data provided for the modeling process. It thus behooves the user
to take great care in specifying the coefficients to be used in the
program for simulating the stream system.
-------
Junction ....
Stretch 1,11..
Retch 1,2, .
Headwater . . .
Discharge or
Withdrawal . .
End of System .
Waste Discharge
Withdrawal
FIGURE 1. SCHEMATIC DIAGRAM OF A TYPICAL RIVER
SYSTEM
-------
Two equations are used to describe the hydraulic characteristics of
each reach in the river system. The first equation represents the
relationship between discharge and depth. It is assumed that both of
these relationships can be represented by the following exponential
equations:
V = A^5! (1)
D = A2QB2 (2)
where
V = mean velocity in a reach (fps)
Q = mean discharge in a reach (cfs)
D = mean depth (feet)
A ,A ,B ,B = coefficients
The above regression coefficients are used as input data into the program.
These coefficients must be developed from data obtained from the actual
stream system. The most readily available data of this type are those
collected by the U.S. Geological Survey at each of its streamflow gaging
station sites. Data from other sources, but of the same type, in a given
river basin, can also be very useful to the modeling process. If the type
of data necessary to develop these relationships are not available for a
given stream, it will be necessary to estimate these coefficients based
on the known topography and physical characteristics of the stream. How-
ever, this method is subject to serious error and is not recommended
unless absolutely necessary.
An extremely important factor in the dissolved oxygen modeling of a
stream system is the reaeration coefficient, K_, which is used in the
calculation of the rate of diffusion of dissolved oxygen into the stream.
There are five options available in this program for specifying the
reaeration rate coefficient. One option is to read it in for each of
the reaches in the stream system. The read-in values might be used if,
for example, field surveys of the stream to be modeled have been made
and values of the reaeration coefficient have been computed from the
results of these. However, the reaeration coefficients determined from
surveys of this type are only useful for the discharges measured during
the survey period. A change in discharge in the stream would probably
result in greatly different values for this coefficient. The program
user may also choose to estimate K2 values for each reach, based on the
known physical and hydraulic characteristics of the stream being modeled.
Obviously this method is very subjective and may involve large errors.
Several investigators have found that the reaeration coefficient, K ,
can be represented by a relationship as shown in Equation 3.
(3)
-------
where
A_,B»,C_ = coefficients
This relationship postulates that the reaeration rate coefficient is
directly proportional to the mean stream velocity and inversely propor-
tional to the mean depth. It is based on two observed phenomena:
• increasing velocity and turbulence increases the surface renewal
rate of dissolved oxygen and promotes mixing and dispersion of the
oxygen throughout the depth of the stream, and
• increased depth decreases the rate of dispersion of dissolved
oxygen throughout the water mass, thus resulting in lower quanti-
ties of oxygen being transferred from the atmosphere.
Several investigators have presented the necessary coefficients for use
in Equation 3 to calculate the reaeration coefficient for any stream in
which the mean velocity and mean depth are known. Table 1 presents these
coefficients as determined by four investigators in both field and
laboratory tests. The coefficients developed by Churchill in 1962 and
by Langbein and Durum in 1967 are probably the best known for the
computation of the reaeration rate coefficient for general model use.
If Equation 3 is preferred for the calculation of the reaeration rate
coefficient, the appropriate coefficients in this equation are read into
the program for every reach in which it is desired to calculate K in
this manner. A note of caution should be observed when using an equation
of this form. Mean stream depths of less than one foot cause the reaera-
tion coefficients predicted by the equation to be higher than are normally
observed under actual field conditions. If the stream being modeled has
significant areas in which the mean depth of flow is less than one foot,
the user is advised to employ alternative methods for computing the re-
aeration coefficient in these areas.
Another technique for computing the reaeration coefficient, also avail-
able to the user of this program, is a direct proportionality between
the reaeration coefficient and the stream discharge. This relationship
is shown in Equation 4.
B,
K2 = A4Q * (4)
A relationship of this type may be developed from data obtained from a
field survey of a stream in which mean velocity and depth were not
determined but discharge was known. The coefficients used in the
discharge-reaeration coefficient equation must be computed from measured
field data. Equation 4 is not generally applicable to most river systems,
because many of the factors which effect the reaeration coefficient are
not adequately described by a simple discharge-reaeration coefficient
relationship.
-------
TABLE 1.
COEFFICIENTS FOR COMPUTATION OF THE REAERATION COEFFICIENT
Investigator
Churchill, et al. (1962)
Langbein and Durum (1967)
O'Connor and Dobbins (1958)
Owens and Gibbs (1969)
A3
5.026
3.3
12.90
9.5
B3
0.969
0.50
0.50
0.67
C3
1.673
1.33
1.50
1.85
-------
The fourth technique available for computing the reaeration coefficient
for each reach is based on the investigation by Thackston and Krenkle.
This technique was developed experimentally by determining the reaera-
tion coefficient in laboratory channels, using as parameters the mean
velocity, channel slope, and mean depth. Equation 5 shows the solution
for the reaeration coefficient as developed by Thackston and Krenkle,
and as is used in this program:
= . 00125 l (5)
where
S = mean channel slope (feet/feet)
Use of this option requires that the program user specify a mean channel
slope for each reach in the river system being modeled. Equation 5
indicates that the reaeration coefficient is proportional to the shear
velocity developed within the stream. Thackston and Krenkle applied
this equation to their laboratory data and showed that it gave a reason-
ably good description of the reaeration rate in the channel. However,
only limited data were available to verify the predictive capability of
this equation in an actual stream system. Studies at the Texas Water
Development Board have indicated that Equation 5 may tend to give higher
values for reaeration coefficients than are actually measured during the
field evaluation of Texas streams.
The fifth technique available calculates K. by Equation 6.
QfiQ
K = 11.61 V A
D1.673 (6)
The program user may specify any of the five methods described above for
the prediction of the reaeration rate coefficient for a given reach.
The use of Equations 4, 5 and 6 for reaeration coefficient computation
require the user to specify the appropriate coefficients for the equa-
tions. Equation 6 requires the program user to specify the mean channel
slope for each reach for which this equation is to be used. The user
may elect to use the same technique for calculating the reaeration co-
efficient for all reaches in the stream system or he may use a different
method for each of the reaches in the system, depending upon the degree
of knowledge obtained of the physical and hydraulic characteristics of
each reach.
Waste discharges are entered into the system by specifying a separate
reach at each location at which a discharge takes place. The reach
specified should be of zero length and should be located at the site of
the actual waste discharge in the prototype system. The user specifies
the waste discharge volume in cubic feet per second, and the concentra-
tions of all constituents in milligrams per liter. A provision is
available in this program to reduce the treatment factor which is read
into the program in percent. If the treatment factor is to be used, it
-------
is assumed that the concentrations of the waste, as read into the pro-
gram, are actually the concentrations present in the raw wastewater
prior to waste treatment. If the user desires to suppress this option,
he inserts a treatment factor of 0.00. In effect, this means that the
values specified for each waste treatment plant will not be changed in
any manner by the model. (For a further discussion of the inflow
treatment factor see the FILE F2 description in Table 2.
The model has provisions for withdrawing water at any location within
the stream system. The water is withdrawn from the stream at the quality
existing at the location of the withdrawal as determined by the model.
The withdrawal is specified in the same manner as a waste input to the
stream system. A separate reach of zero length is set up for each with-
.drawal in the system. A negative flow is specified on the waste and
withdrawal input cards to indicate that water is being withdrawn from
the system. The treatment factor values are not taken into considera-
tion by the model for withdrawal.
The user of the DOSCI program has several options available which will
enable him to simulate several different stream conditions in one
computer run, without reading in additional data. The user of the
program may read in up to four waste treatment factors. The program will
calculate a new dissolved oxygen profile based on the organic load re-
leased from each plant after the treatment factor has been applied.
This process is repeated for each of the treatment factors entered in
the program. The user may also specify up to four dissolved oxygen
target levels, which are the minimum permissable dissolved oxygen
concentrations in the stream system. By specifying a positive dis-
solved oxygen target level, the user also indicates that he wishes the
program to calculate the flow augmentation requirements, if any, needed
to meet this target level. The other option available to the user is
that up to twelve different temperatures and corresponding headwater
flows may be specified, and the program will completely model the stream
for each value. It should be noted that the input headwater temperature
(FILE H) is used as the temperature in a given reach only if no
temperature is specified for that reach (FILE J).
The above three options enable the program user, in a single run, to
perform a large number of simulations of the stream system to determine
the effects of various waste loadings, temperatures, and dissolved
oxygen target levels on the dissolved oxygen concentration within the
system.
The procedure used by DOSCI to determine augmentation requirements
should be briefly mentioned. The model begins routing organic wastes
and dissolved oxygen from the uppermost point in the stream system and
proceeds downstream. As the simulation progresses downstream, reach by
reach, the calculated dissolved oxygen concentration is checked against
the target dissolved oxygen level specified by the user. When the model
discovers a dissolved oxygen concentration below the target, it stops
at this reach. The model then searches all of the upstream headwaters
10
-------
to see which, headwaters have water available for flow augmentation
purposes. The model then estimates, using a parabolic relationship be-
tween dissolved oxygen, deficit and the target dissolved oxygen level,
the quantity of water required for flow augmentation to increase the
minimum dissolved oxygen, concentration to the target level. The volume
of water required is then divided equally among all of the headwaters
from which water is available for augmentation. These new flows are then
re-routed through the stream system. If the amount of augmentation was
insufficient to raise the minimum dissolved oxygen concentration above
the target level at the reach being investigated, the process is repeated
until the target dissolved oxygen concentration is attained. After the
target dissolved oxygen concentration has been satisfied, the program
proceeds downstream until it comes to another dissolved oxygen concentra-
tion level below the target and then the process is repeated as before.
The program is designed so that it can only augment from a headwater
stretch, and the augmentation requirements are divided equally among all
the headwater stretches which have augmentation availability. The flow
augmentation option is suppressed by specifying a negative value for the
target dissolved oxygen concentration.
MODEL VERIFICATION AND SENSITIVITY
The DOSCI model was verified on five different stream regions in the
Spokane Basin in the states of Washington and Idaho. These regions
varied from deep, slow flowing streams such as the Spokane River above
FDR Lake to shallow, steep streams such as the upper reaches of the South
Fork of the Coeur d'Alene River. Two different simulation periods were
chosen for each of the five regions. Quality constituents modeled in-
cluded coliforms, zinc, chlorides, BOD, NH -N, NO—N, NO -N, PO.-P and
DO. Waste discharges modeled included urban sewers and mine processing
outfalls. In general the simulated values of the constituent concentra-
tions throughout the system matched the observed values very well. A
sensitivity analysis was also performed, using one of the five regions
as the study area. The results of the verification runs and the
sensitivity analysis may be found in Volume I (Parts IV and V) of this
report.
11
-------
SECTION III
PROGRAM DESCRIPTION
The subroutines employed by DOSCI are described in this section. Flow-
charts are provided for the new subroutines with the exception of the
plotting routines. The nontrivial equations occur in subroutines K2CAL
and GETCON. The equations used by K2CAL are described in Table 2 (File
F-4). The equations used by GETCON are described in Volume I, Part III.
The important local variables are described in Table 5. All subroutine
interlinkages and simplified flowcharts for all subroutines are presented
in Figure 2, which summarizes the modifications to the original flowchart
in Reference [8].
Main Program
The main program, with the help of subroutines NEWIN and PLTSET, reads
and echo prints all input data. It sets the basin conditions to be
considered, and calls the necessary subroutines for the modeling of the
stream system. The main program repeats the modeling process for each
of the options specified for program operation. Figure 2, which is
based on Figure 4 of Reference [8], shows the general program subroutine
interlinkages.
Subroutine BLEND
This subroutine computes the discharges and the quality constituent
concentrations entering a stretch immediately downstream from a junction.
Discharges are computed by simple addition; constituent concentrations
are computed by a mass balance at the junction.
Subroutine CMIN
For each reach, this routine finds the sub-reach having the lowest
dissolved oxygen concentration, and the river mile at which this concen-
tration occurs.
Subroutine DOEQU
When DOEQU is called to process a given reach the following procedure is
followed. If the reach has length zero, the end conditions are set equal
to the start conditions, the appropriate quantities are stored for future
reference, and DOEQU returns control to the calling routine. If the
reach length is not zero and only conservatives are being modeled, the
concentrations of the conservatives are calculated at the end of the
reach by a simple mass balance, the appropriate quantities are stored for
future flow augmentation computations and printout, and control is returned
to the calling routine. If the reach length is not zero and nonconserva-
tives are being modeled, the reach is divided into ten equal subreaches
and each subreach is then considered as follows. The velocity, depth,
13
-------
1. Find* lh« hcid vat«i
•trttcK«i th»c «cu*
2, Cilia BLENT).
3. 5«(a Bachlm
Idcntldi-r* to
lrdlc*ti th* nuctir
nntaclon upatrei*.
4. C.ll. TRIED.
S. Ctlli SCAS.
6. S«ti final condl-
2. Calif DOSAT.
4. C.lti KEIKI.
3. C<1U RMATC.
7. Write* final •<
t. C«lli FLTSET.
Siti Initial con-
dition, of it fitch
• qual to final
e««dlo| itnccb.
Calli TKIBD.
Cilia tUS.
Sita final con-
dition! to cal-
culated condition*.
Kipiata fro* (1)
•tr«tch.
»•• THBO •*•
Chick* to Ma If thtrt
an any dovnatriu
nachci 10 b« con*Id*I
J. Coaaldiri o«*t dovn-
3. Calli KZFRZ.
4. Calli DOEQU.
dlaaolwd oxy(*o of
nach.
1. Quckj for riich*i ID
• itt.tcb.
3. Calculac*B DO deficit
for r*ach.
4. Cilculacei imunt of
dilution uattr CMdid
dtflclt.
5. Check* to
«( •ugceataclon
•vallibli.
t. Divldt. dllutloi
1 or 10 subreacl.es
and considers each
iLlvar •!!• at
b«(innlo( of «.
aubrcach la
calculated.
tin (low.
S. Call C2CJU..
6. Call GETCOB.
7. Call CMIK.
9. Scon r*ach 10^
for prloilni Ui
final .a —ry.
FIGURE 2. GENERAL PROGRAM FLOW DIAGRAM
-------
and travel time for the subreach are computed, subroutine K2CAL is
called to calculate K?, and subroutine GETCON is called to compute the
concentrations of the nonconservatives at the end of the subreach, given
the concentrations at the beginning of the subreach. Any incremental
flow is added, the concentrations are adjusted accordingly by a mass
balance, and the procedure is repeated on the next subreach. When the
last subreach is completed appropriate quantities are stored, and control
is returned to the calling routine.
Subroutine DOSAT
This subroutine calculates the dissolved oxygen saturation concentration
for each reach based on the temperature and elevation in that reach.
Subroutine GETCON
Subroutine GETCON is called by subroutine DOEQU to calculate the changes
in concentrations which occur in a stream reach over a time interval At,
where At is the travel time through the reach. The algorithms used are
explained in detail in Part III of Volume I. The inputs to GETCON in-
clude stream depth and velocity, average light intensity at the surface
during At, and the appropriate system constants and reach variables, as
well as the concentrations of all constituents at the beginning of At.
A flow chart of subroutine GETCON is given in Figure 3. Equations
referenced as A.NN may be found in Volume I, Part III.
Subroutine K2CAL
This subroutine calculates the reaeration rate coefficient for each reach
in the river system being modeled. The technique used for calculating
the reaeration coefficient is specified by the user in the input data
(see Table 2, File F-4).
Subroutine NEWIN
This subroutine reads and echo prints the data input to FILE J and FILE
K. Subroutine SETOPT is called to set various option flags based on the
inputs to FILE K. A flowchart of subroutine NEWIN is given in Figure 4.
Subroutine PLOT
This subroutine is called to plot various quantities on the line printer.
See the description of input FILE L for a list of the quantities which
may be plotted.
Subroutine PLTSET
This subroutine has two uses. It is called to read and echo write FILE
L input data. It is also the driver for subroutine PLOT and is used in
this respect for plotting desired quantities.
15
-------
IS
THIS THE
FIRST TIME
IN
INITIALIZE
Zero Out Deltas
TCOR = (T-20)
A.82
ARE
COLI FORMS
BEING
MODELED ?
COLIFORM DECAY
A. 8 A.
10 A. 7
BOD DECAY
A.15 A.14
A.12 A.16
A.13 A.17
FIGURE 3. FLOWCHART - SUBROUTINE GETCON (Page 1)
16
-------
IS
PO..-P \YES
BEING MODELED
IS
NH3-N \YES
BEING MODELED
FIGURE 3.
FLOWCHART-SUBROUTINE GETCON (cont'd)
17
-------
IS
N02-N
BEING MODELED
N02-N DECAY
A.32 A.33
A.31
IS
N03-N
BEING MODELED
7
IS
PCU-P
BEING MODELED
i
IS
NH,-N
'/OLITIZATIO!
OCCURING ?
f!H,-N VOLITIZATION
ARE
HYTOPLANKTO
EING MODELE
IS
IT DARK
DURING TIME
STEP ?
ALGAL GRO/JTH
A.4"5"
A.46
A.47
A.48
A.49
A.50
A.51
A.52
A.53
A.54
A.55
A.44
A.60
FIGURE 3.
FLOWCHART-SUBROUTINE GETCON (cont'd)
18
-------
IS
BOD \ YES
BEING MODELED > '
^v YES
ODELED^ — *•
ALGAL RESPIRATION
A. 66
YES
ALGAL RESPIRATION
A.65
ALGAL
RESPIRATION.
A.93 A.29 A.34
FIGURE 3.
FLOWCHART-SUBROUTINE GETCON (cont'd)
19
-------
ARE
HEAVY \ YES
METALS BEING
MODELED
FIGURE 3.
FLOWCHART-SUBROUTINE GETCON (cont'd)
20
-------
CENTRYj
READ 7
REACH /
VARIABLES/
WRITE
REACR
VARIABLES
FODMISUELV
LANEOUS
VVARIABLES
CALL
SETOPT
WRITE SUMMARY
OF PARAMETERS
TO BE USED
fRETURN^\
FIGURE 4. SUBROUTINE NEWIN
21
-------
Subroutine REINI
This subroutine flags each headwater stretch to indicate whether water
is available for flow augmentation. It also sets the initial conditions
of a stretch equal to the final conditions of the preceeding stretch.
After the calculations for each stretch, are completed, including flow
augmentation requirements (if any), this subroutine sets the final condi-
tions for each stretch equal to the calculated conditions prior to print-
out. This procedure is repeated for each of the headwater stretches.
Subroutine REPRE
This subroutine calculates the concentrations at the head of a reach from
the end conditions of the upstream reach and the concentrations of any
tributary inflow using the mass balance technique. The treatment factor
is applied to the inflow if required, i.e., if the DO level of the
tributary inflow is input as a negative number (see footnote to FILE F-2
input data).
Subroutine RMATC
Subroutine RMATC identifies the headwater stretches as they enter a junc-
tion. As computations in the program progress downstream, RMATC stores
and flags the headwater sources which are available for flow augmentation.
This flagging operation is performed each time a junction within the river
basin system is encountered.
Subroutine RUNON
This subroutine initializes the concentrations in each headwater of the
system and calls DOSAT to calculate the DO saturation level in each
reach. This subroutine also prints the headings for the intermediate
and the final summaries, and prints out all of the information in the
final summary after all calculations have been completed. Subroutine
PLTSET is called to plot any quantities requested in FILE L.
Subroutine SCAN
Subroutine SCAN, operating on one stretch at a time, computes the flow
augmentation required (if any) to meet the desired target dissolved
oxygen concentration at all points in the stretch. It first checks for
all available augmentation upstream of the stretch under consideration.
The subroutine then determines the minimum dissolved oxygen level in
each of the reaches in the stretch. It calculates and stores the dis-
solved oxygen concentration. It then estimates the quantity of dilution
water required to satisfy the greatest deficit within the stretch. The
required amount of dilution water is divided equally among all sources
available for flow augmentation. SCAN then adds the dilution water to
each of these upstream sources.
22
-------
Subroutine SETOPT
This subroutiaa is called by subroutine NEWIN to set various option
flags as required by the FILE K Inputs. A listing of the constituents
which are to be modeled is output. A flowcliart of SETOPT is presented
in Figure 5.
Subroutine TRIED
This subroutine processes each reach in a stretch by calling REPRE to
calculate conditions at the head of the reach and DOEQU to calculate the
changes occurring throughout the reach. The minimum dissolved oxygen
concentration in the reach is stored for use by flow augmentation calcu-
lations.
23
-------
CENTRYj
SET OPTIONS
BASED ON IN-
PUT PARA-"
METERS
WRITE
OPTION
SUMMARY
RETURN
)
FIGURE 5. SUBROUTINE SETOPT
24
-------
SECTION IV
INPUT REQUIREMENTS
The inputs required by DOSCI are described in this section in Table 2
and Table 3. Any variable which does not have a default value must be
input. If a value of zero is desired for an input variable which has a
non-zero default value, a small positive number (e.g., 0.00001) must be
input. Restrictions on input data are discussed in Section II. Equa-
tions referenced as (A.NN) in Table 2 may be found in Part III of Volume
I.
Figure 6 illustrates the physical arrangement of the DOSCI input deck.
A typical DOSCI input deck is presented in Table 6.
25
-------
FILE L Plotting Variables
FILE K Miscellaneous Variables
FILE J Reach Variables
FILL I Mean Monthly Headwater Flows
FILE M Mean Monthly Temperatures
Treatr'jnt Factors and
DO Tarrj?t Levels
FILE F-4 Data 1,21 through Data 1,23)
FILE F-3 Data (1,15) throunh Data (1,20)
Data (1,7) through Data (I.
Data (!,2-'>) throunh n = U (1,32)
FILE F-l Data (1,1) throuTh 33ta (1,6)
FILE E Headwaters
FILE D Junctions
FILE C Reach Order
FILE B Physical Description
FILE A Basin Title
POSCI PROGRAM 3ECK
FIGURE 6. INPUT DATA ORGANIZATION
26
-------
CARD
#
1
2
1
2
CARD
COLUMN
1-8
9-72
1-8
1-8
14-15
25-26
35-36
45-46
55-56
75-80
1-8
TABLE 2 .
DOSCI INPUT
VARIABLE DEFAULT
FORMAT NAME UNITS VALUE
File A - Run Title
2A4 DUM
18A4 TITLE
2A4 DUM
File B - Physical Description
2A4 DUM
12 NINIT
12 NJUNC 0
12 NREA
12 NTRIB
12 ICK 0
F6.1 ELEV (FEET) 0.
2A4 DUM
DESCRIPTION
FILE A
Title for run
ENDFILE
FILE B
Number of headwater reaches
(
-------
TABLE 2. (Continued)
ro
oo
CARD CARD
# COLUMN
1 1-8
12-13
21-80
2-NTRIB Each
reaches.
NTRIB+1 1-8
1 1-8
22-23
30-31
40-41
50-51
FORMAT
2A4
12
20(13)
of the
2A4
2A4
12
12
12
12
VARIABLE
NAME UNITS
File C - Reach Ordering
DUM *
I
(IORD(I,J),
J=l,20)
NTRIB stretches requires a card
DUM
File D - Junctions
DUM
I
JUNG (1,1)
JUNC(I,2)
JUNG (I, 3)
DEFAULT
VALUE DESCRIPTION
FILE C
Stretch number (<20)
Reaches in stretch I (upstream
to downstream, no more than 20
reaches in a stretch)
(same format as card #1) defining its
ENDFILE
FILE D
Junction number (<20; junctions
should be numbered in increasing
order downstream)
Upstream stretch for junction I
Upstream stretch for junction I
Downstream stretch for junction
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
File D - Junctions
2-NJUNC Each of the NJUNC junctions requires a card (same format as card #1) defining its
stretches.
NJUNC+1 1-8 2A4 DUM
If there are no junctions, only the ENDFILE card is required.
ENDFILE
File E - Headwater Data
NJ
1-8
9-13
14-18
1-10
11-20
21-30
2A4
15
15
F10.4
F10.4
F10.4
DUM
I
LAUG(I)
CONDZ(I,1) (MG/L)
CONDZ(I,2) (MG/L)
CONDZ(I,3) (MG/L)
0
0
0
0
31-40 F10.4 CONDZ(I,4)
(MG/L)
FILE E
Headwater stretch number (must
be £10)
0 = no augumentation
1 = augment flow if necessary
DO concentration in headwater I
BOD concentration in headwater I
NH_-N concentration in headwater
I
NO--N concentration in headwater
I
-------
TABLE 2. (Continued)
U)
o
CARD CARD
// COLUMN
2 41-50
51-60
61-70
71-80
3 1-10
11-20
21-30
31-40
41-50
51-60
FORMAT
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
VARIABLE
NAME
File E -
CONDZ(I,5)
CONDZ(I,6)
CONDZ(I,7)
CONDZ(I,8)
CONDZ(I,9)
CONDZ(I,10)
CONDZ(I,11)
CONDZ(I,12)
CONDZ(I,13)
CONDZ(I,14)
UNITS
Headwater Data
(MG/L)
(MG/L)
(MG/L)
(MPN/100)
(MG/L)
(MG/L)
(MG/L)
(MG/L)
(MG/L)
(MG/L)
DEFAULT
VALUE DESCRIPTION
0. NO»-N concentration in headwater
I
0. P0/,~p concentration in headwater
0. CHLOR-A (phytoplankton) concen-
tration in headwater I
0. COLIFORM concentration in head-
water I
0. HEAVY METAL 1 concentration in
headwater I
0. HEAVY METAL 2 concentration in
headwater I
0. HEAVY METAL 3 concentration in
headwater I
0. TOTAL NITROGEN concentration in
headwater I
0. CHLORIDE concentration in head-
water I
0. HEAVY METAL 1 ION concentration
headwater I (not used if PIHM1
0, FILE K card 9)
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
File E - Headwater Data
61-70
71-80
F10.4
F10.4
CONDZ(I,15)
CONDZ(I,16)
(MG/L)
(MG/L)
0.
0.
HEAVY METAL 2 ION concentration
in headwater I (not used if
PIHM2 / 0, FILE K card 10)
HEAVY METAL 3 ION concentration
in headwater (not used if
PIHM3 ± 0, FILE K, card 10)
4-3*NINIT
Each of the NINIT headwaters requires 3 cards (same formats as cards #1, #2, #3) defining the
concentrations of its constituents. Although no concentrations need be specified on the cards
(in which case the default values apply), the cards must be included.
3*NINT+1 1-8
2A4
BUM
File F-l Reach Data
1-8
17-18
23-30
33-40
43-50
53-60
2A4
12
F8.0
F8.0
F8.0
F8.0
DUM
I
DATA (1,1)
DATA (I, 2)
DATA (I, 3)
DATA(I,4)
(MILES)
(MILES)
(HOUR"1) . 008
(HOUR"1) 0.
ENDFILE
FILE F-l
Reach number (<99)
Reach length
River mile to head of reach
BOD decay coefficient
BOD settling coefficient
-------
TABLE 2. (Continued)
CARD CARD VARIABLE DEFAULT
// COLUMN FORMAT NAME UNITS VALUE
DESCRIPTION
File F-l Reach Data
63-70
73-80
F8.0
F8.0
DATA(I,5) .255 Coefficient on flow for velocity
DATA(I,6) .414 Exponent on flow for velocity
2-NREA
Each of the NREA reaches requires a card (same format as card #1) defining the above data.
u>
NJ
NREA+1 1-8
1-J
2A4
2A4
12-20 19
DUM
File F-2 Reach Data
DUM
I
21-30 F10.4 DATA(I,7) (CFS)
1-10 F10.4 DATA(I,8) (MG/L)
11-20 F10.4 DATA(1,9) (MG/L)
ENDFILE
FILE F-2
Reach number (£99)
Flow entering reach I (see F-2
note below)
DO concentration in flow for
reach I
BOD concentration in flow for
reach I
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
2 21-30
31-40
41-50
51-60
61-70
71-80
3 1-10
11-20
21-30
31-40
FORMAT
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
VARIABLE
NAME
File F-2
DATA (I, 10)
DATA(I,11)
DATA (I ,12)
DATA (I, 13)
DATA (I, 14)
DATA (I, 24)
DATA (I, 25)
DATA(I,26)
DATA(I,27)
DATA(I,28)
UNITS
Reach Data
(MG/L)
(MG/L)
(MG/L)
(MG/L)
(MG/L)
(MPN/100)
(MG/L)
(MG/L)
(MG/L)
(MG/L)
DEFAULT
VALUE DESCRIPTION
0. NH--N concentration in flow
reach I
0. NO -N concentration in flow
reach I
0. NO_-N concentration in flow
reach I
0. PO/.~^ concentrati°n i-n flow
reach I
for
for
for
for
0. CHLOR-A(phytoplankton) concen-
tration flow for reach I
0. COLIFORM concentration in flow
for reach I
0. HEAVY METAL 1 concentration
flow for reach I
0. HEAVY METAL 2 concentration
flow for reach I
0. HEAVY METAL 3 concentration
flow for reach I
0. TOTAL NITROGEN concentration
in
in
in
in
flow for reach I
-------
TABLE 2. (Continued)
CARD CARD
// COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
U)
-C-
41-50
51-60
61-70
71-80
4-N
File F-2 Reach Data
F10.4
F10.4
F10.4
F10.4
DATA(I,29)
DATA(I,30)
DATA(I,31)
DATA(I,32)
(MG/L)
(MG/L)
(MG/L)
(MG/L)
0.
0.
0.
0.
CHLORIDE concentration in flow
for reach I
HEAVY METAL 1 ION concentration
in flow for reach I (not used
if PIHMI ^ 0, FILE K card 9)
HEAVY METAL 2 ION concentration
in flow for reach I (not used
if PIHM2 ± 0, FILE K, card 10)
HEAVY METAL 3 ION concentration
in flow for reach I (not used
if PIHM3 ± 0, FILE K, card 10)
Each reach with a nonzero input or withdrawal requires 3 cards (same formats as cards #1, #2,
//3) defining the above quantities. Although no concentrations need be specified on the cards
(in which case the default values apply), the cards must be included.
The last card must be followed by the following:
N+l 1-8 2A4 DUM ENDFILE
If no reach has a nonzero inflow, only the ENDFILE card is required.
-------
TABLE 2. (Continued)
F-2 Note:
If the length of reach I is zero (card #1, File F-l) and the flow is positive, the reach is
considered to be a point source. If the DO concentration is input as a negative number, the
treatment factor (card #1, File G) is applied to all constituents being modeled and the DO
concentration is treated as a positive quantity. If the DO is not input as a negative number,
the treatment factor is not applied. If the length is zero and the flow is negative, the reach
is considered to be a point withdrawal. The constituents are withdrawn from the system at the
concentrations determined by the model.
un If the length of the reach is positive and the flow is positive, the flow is added uniformly
to the system over the reach length. This is accomplished by considering the flow to be divided
evenly amongst ten equally spaced' point sources along the reach length. Each point source is
then treated as described above.
Similarly, if the length is positive and the flow is negative, the flow is withdrawn uniformly
over the reach length by considering ten equally spaced point withdrawals and treating each
point withdrawal as described above.
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
1 1-8
17-18
23-30
33-40
43-50
53-60
63-70
73-80
FORMAT
2A4
12
F8.0
F8.0
F8.0
F8.0
F8.0
F8.0
VARIABLE
NAME UNITS
File F-3 Reach Data
DUM
I
DAT A (I, 15) (HOUR"1)
DATA (I, 16)
DATA (I, 17)
DATA (I, 18)
DATA (I, 19)
DATA(I,20)
DEFAULT
VALUE DESCRIPTION
FILE F-3
Reach number (<99)
0. K- (reaeration) value for
option #1 (see File F-4 for a
definition of option values)
3.3 Velocity coefficient for K^
equation for option #2
1.0 Velocity exponent for K9
equation for option #2
1.33 Depth exponent for K? equation
for option #2
0. Flow coefficient for K«
equation for option #3
0. Flow exponent for K equation
for option #3
2-N If the default values are not satisfactory for a given reach, a card (same format as card #1)
must be included for that reach.
-------
TABLE 2. (Continued)
CARD CARD
// COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
File F-3 Reach Data
If all default values are satisfactory, only the ENDFILE card is required.
N+l
1-8
2A4
DUM
ENDFILE
File F-4 Reach Data
1-8 2A4
17-18 12
25-26 12
31-50 5A4
53-60 F8.0
63-70 F8.0
73-80 F8.0
DUM
I
K20PT(I)
RIDENT
DATA(I,21)
DATA(I,22)
DATA(I,23)
.259
.414
FILE F-4
Reach number (£99)
Option for calculating reaera-
tion constant K (see F-4 note
below)
Reach identification
Flow coefficient for depth
equation
Flow exponent for depth
equation
Channel slope for K equation
for option 4
-------
TABLE 2. (Continued)
CARD CARD VARIABLE DEFAULT
# COLUMN FORMAT NAME UNITS VALUE DESCRIPTION
File F-4 Reach Data
2-NREA Each of the NREA reaches requires a card (same format as card #1) defining the above data.
NREA+1 1-8 2A4 DUM ENDFILE
F-4 Note:
The reaeration constant K is calculated as follows for reach I
£ Option #1 K2 = DATA(I,15)
Option #2 K0 = ^- 4^
where A = DATA(I,16)
V = velocity (ft/sec)
B = DATA (I, 17)
D = (ft)
C = DATA(I,18)
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
Option #3
= (A.Q) (
where A = DATA(I,19)
Q = flow (CFS)
B = DATA(I,20)
Option #4
where V = velocity (ft/sec)
2
g = gravity (ft/sec )
D = depth (ft)
S = DATA(I,23)
-------
TABLE 2. (Continued)
CARD
#
CARD
COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
Option #5
T,
K9 =
2
n fii
,11. ol V
969
where V = velocity (ft/sec)
D = depth (ft)
A = DOK2 (I) (see File J)
o
File G - Treatment Factors and DO Levels
1-8
11-15
2A4
F5.0
DUM
DOL(l)
(MG/L) 0.
16-30
31-35
3F5.0 (DOL(I),
1=2,4
F5.0
TRFAC(l)
(MG/L)
FILE G
DO concentration level for 1st
run (if the DO concentration
falls below this level in any
reach, augmentation will occur
if allowed)
DO levels for 2nd, 3rd, and 4th
runs (if DOL(l)=0, for 1=2,3,4
run I will not be made)
Fraction of constituents removed
at treatment stations (see foot-
note to File F-2) for 1st case
-------
TABLE 2. (Continued)
CARD
#
CARD
COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
36-50
l-l
File G - Treatment Factors and DO Levels
3F5.0
2A4
(TRFAC(I)
1=2,4)
DUM
0.
Fractions of constituents removed
for 2,3,4 (if TRFAC(I)=0, for
1=2,3,4, case I will not be run)
ENDFILE
If the default values are satisfactory, only the ENDFILE is required.
Run 1-Case 1, i.e., the DOL(l) - TRFAC(l) combination is always run. In addition, for any non-
zero TRFAC(K),K=2,3,4, the DOL(l) - TRFAC(K) case is run. For 1=2,3,4 if DOL(I)^0, the DOL(I)
- TRFAC(l) case is run plus the DOL(I) - TRFAC(K) cases for TRFAC(K)^0, K=2,3,4. Hence a maximum
of 16 cases may be run. The default values result in only the DOL(l) - TRFAC(l) case, i.e., no
augmentation and no treatment.
-------
TABLE 2. (Continued)
CARD
//
1
CARD
COLUMN
1-8
21-80
VARIABLE
FORMAT NAME
File H - Mean
2A4 DUM
12F5.0 (TEMMO(I),
1=1,12)
UNITS
Monthly
(CENT0
DEFAULT
VALUE
Temperatures
) o.
DESCRIPTION
FILE H
Mean monthly temperatures (Oct.-
Sept.) for the system being
I-t
2A4
DUM
modeled (if no temperature is
assigned to a given reach in
File J, the average is assigned
to that reach
ENDFILE
For each case executed (as defined in File G), a subcase is run for each TEMMO(I)^0. Hence at
least one TEMMO(I) must be nonzero.
File I - Headwater Flows
1-8 2A4 DUM
12-13 12 I
21-80 12F5.0 (HWFLOW(I.J), (CFS)
J=l,12)
FILE I
Headwater stretch number (
-------
TABLE 2. (Continued)
CARD
//
CARD
COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
File I - Headwater Flows
2-NINIT
Each headwater requires a card (same format as Card #1) defining the flows.
NINIT+1 1-8
2A4
DUM
ENDFILE
File J - Reach Variables
1-8
10-13
14-19
20-25
26-31
32-37
38-43
2A4
13
F6.0
F6.0
F6.0
F6.0
F6.0
DUM
I
COLK(I)
NH3K(I)
N02K(I)
N03K(I)
P04K(I)
(HOUR 1)
_i
(HOUR )
_i
(HOUR )
(HOUR" )
_i
(HOUR )
FILE J
Reach number (£99)
.004 Coliform reaction coefficient
for reach I (A.7)
.004 NH_ reaction coefficient for
reach I (A.26)
.015 NO reaction coefficient for
reach I (A.31)
.0014 NO reaction coefficient for
reach I (A.37)
.0009 PO, settling coefficient for
reach I (A.75)
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
44-49
50-55
56-61
62-67
68-73
74-79
VARIABLE
FORMAT NAME
File J
F6.0 EXTK(I)
F6.0 DOK2(I)
F6.0 HMIK(I)
F6.0 HM3K(I)
F6.0 HM3K(I)
F6.0 TEMREA(I)
DEFAULT
UNITS VALUE
- Reach Variables
(HOUR"1) .04
(HOUR"1) 1.
(HOUR"1) .004
(HOUR"1) 0.
(HOUR"1) 0.
(CENT0) File H value
DESCRIPTION
Extinction coefficient for reach
I (A. 51)
K coefficient for option 5 (see
footnote to File F-4)
Settling coefficient for HEAVY
METAL 1 for reach I (A. 2)
Settling coefficient for HEAVY
METAL 2 for reach I (A. 2)
Settling coefficient for HEAVY
METAL 3 for reach I (A. 2)
Temperature in reach
2-N If the default values are not satisfactory for a given reach, then a card (same format as card
+1) defining the above variables must be input for that reach.
If all default values are satisfactory, only the ENDFILE card is required.
N+l 1-8 2A4 DUM ENDFILE
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
1-8
10-20
1-8
11-20
21-30
31-40
1-8
11-20
2A4
F10.4
2A4
110
110
110
2A4
110
File K - Miscellaneous Variables
DUM
10
DUM
IFN
(LANGLEYS/MIN)
ICOL
ICOMB
DUM
IHEAVY
FILE K
Average light intensity (A.53)
FILE K
PHYTOPLANKTON growth function
option
0 = growth limited by NO.,-N
concentration
1 = growth limited by NH--N
concentration
2 = growth limited by maximum
of NH -N and NO -N
COLIFORM option
0 = don't model COLIFORMS
1 = model COLIFORMS
Constituent selection option
(see Table 3)
FILE K
Heavy metal option
0 = model no heavy metals or
ions
N = model N heavy metals and
their associated ions
(N=l,2,3)
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
3 21-30
31-40
4 1-8
11-20
VARIABLE DEFAULT
FORMAT NAME UNITS VALUE
File K - Miscellaneous Variables
110 ITOTN
110 ICHLOR
2A4 DUM
110 INH 1
DESCRIPTION
TOTAL NITROGEN option
0 = don't model TOTAL NITROGEN
1 = model TOTAL NITROGEN
CHLORIDE option
0 = don't model CHLORIDES
1 = model CHLORIDES
FILE K
NH reaction order
2 = 2nd order
21-30
110
IN 2
NO reaction order
1 = 1st order
2 = 2nd order
31-40
110
IN 3
NO reaction order
1 = 1st order
2 = 2nd order
41-50
110
IP
PO. reaction order
4
1 = 1st order
2 = 2nd order
-------
TABLE 2. (Continued)
CARD
CARD
COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
File K - Miscellaneous Variables
1-8 2A4 DUM
11-20 F10.4 THKCOL
21-30 F10.4 ABOD
31-40 F10.4 AHM
41-50 F10.4 CHMOC
51-60 F10.4 THKNH3
61-70 F10.4 VOLITK
71-80 F10.4 THVOLK
(MG/L)
1.07
0.
0.
20.
1.10
.01
.17
FILE K
Temperature correction constant
for (A.8) coliform (COL)
reaction coefficient
Coefficient on BOD in COL
calculation (A.10)
Coefficient on HEAVY METAL 1
(HM1) in (A.10) COL calculation
HM1 concentration limit in COL
calculation (A.10)
Temperature correction constant
for (A.27) NH -N decay
coefficient
Exponent for NH--N volitization
(A.35)
Temperature correction constant
for (A.35) NH -N volitization
-------
TABLE 2. (Continued)
CO
CARD CARD
# COLUMN
6 1-8
11-20
21-30
31-40
41-50
51-60
61-70
FORMAT
2A4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
VARIABLE
NAME
File K -
DIM
BODC
BOON
BODPC
BODOQ
NOREFR
GRMAX
DEFAULT
UNITS VALUE
Miscellaneous Variables
106.
16.
.5
(MG O./MG BOD) 1.5
.5
(HOUR"1) .1
DESCRIPTION
FILE K
Carbon to phosphorus ratio in
BOD (A. 18)
Nitrogen to phosphorus ratio in
BOD (A. 19)
Dry weight fraction of carbon
in BOD (A. 18)
BOD - oxygen quotient (A. 16)
Non-refractory part of BOD (A. 17
Maximum fractional growth rate
71-80
1-8
11-20
F10.4
THGRMX
1.07
2A4 DUM
F10.4 CHMOA (MG/L)
20.
for phytoplankton at 20° centi-
grade (A.45)
Temperature correction constant
for GRMAX (A.45)
FILE K
HM1 limit for phytoplankton
growth (A.46)
-------
TABLE 2. (Continued)
CARD CARD VARIABLE
# COLUMN FORMAT NAME UNITS
File K -
7 21-30 F10.4 HMKA
31-40 F10.4 MP04
41-50 F10.4 MIN03
51-60 F10.4 M2N03
61-70 F10.4 MNH3
71-80 F10.4 ML
DEFAULT
VALUE DESCRIPTION
- Miscellanous Variables
(MG P04-P/L)
(MG N/L)
(MG N03-N/L)
(MG NH -N/L)
(LANGLEYS/MIN)
.01 HM1 coefficient for phytoplankto
(A. 46) growth calculation
t
.03 Michaelis-Menton constant (A. 47)
.028 Michaelis-Menton constant (A. 47)
.045 Michaelis-Menton constant (A. 48)
.045 Michaelis-Menton constant (A. 49)
.03 Light intensity calculation
1-8
11-20
21-30
31-40
41-50
2A4 DUM
F10.4 APR
F10.4
F10.4
F10.4
NR
ASR
AND
(HOUR DEG.C
(FT/HOUR)
(HOUR" )
.0001
.05
.001
factor (A.50)
FILE K
Chlorophyll-A to phosphorus
ratio in phytoplankton (A.57)
Phytoplankton respiration factor
(A.63)
Phytoplankton sinking rate (A.68)
Fractional death for phytoplank-
ton (A.72)
-------
TABLE 2. (Continued)
Ul
o
CARD CARD
# COLUMN
8 51-60
61-70
71-80
9 1-8
11-20
21-30
31-40
41-50
51-60
61-70
71-80
FORMAT
F10.4
F10.4
F10.4
2A4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
VARIABLE DEFAULT
NAME UNITS VALUE
File K - Miscellaneous Variables
ATD . 001
BRRBOD (MG/M2-HR) 61.
2
BRRP04 (MG/M -HR) .125
DUM
2
BRRNH3 (MG/M -HR) .108
BENOD (MG/M2-HR) 15.
AHM2 0.
AHM3 0.
ATD2 0.
ATD3 0.
PIHM1 0.
DESCRIPTION
Phytoplankton toxic death
coefficient for HM1 (A. 73)
BOD benthal release rate (A. 79)
PO.-P benthal release rate (A. 78)
4
FILE K
Nitrogen benthal release rate
(A. 77)
Benthal oxygen demand (A. 81)
Coefficient on HEAVY METAL 2
(HM2) in COL calculation (A. 10)
Coefficient on HEAVY METAL 3
(HM3) in COL calculation (A. 10)
Phytoplankton toxic death
coefficient for HM2 (A. 73)
Phytoplankton toxic death
coefficient for HM3 (A. 73)
Fraction of HM1 in ion form (A. 3)
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
10 1-8
11-20
21-30
31-40
41-50
51-60
61-70
71-80
FORMAT
2A4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
F10.4
VARIABLE DEFAULT
NAME UNITS VALUE
File K - Miscellaneous Variables
DUM
PIHM2 0.
PIHM3 0.
CHM02C (MG/L) 0.
CHM03C (MG/L) 0.
CHMOA2 (MG/L) 0.
CHMOA3 (MG/L) 0.
HMKA2 0.
DESCRIPTION
FILE K
Fraction of HM2 in ion form
(A. 3)
Fraction of HM3 in ion form
(A. 3)
HM2 concentration limit in COL
calculation (A. 10)
HM3 concentration limit in COL
calculation (A. 10)
HM2 limit for phytoplankton
growth (A. 46)
HM3 limit for phytoplankton
growth (A. 46)
HM2 coefficient for phytoplank-
11 1-8 2A4 DUM
11-20 F10.4 HMKA3
ton growth calculation (A.46)
FILE K
HM3 coefficient for phytoplank-
ton growth calculation (A.46)
-------
TABLE 2. (Continued)
CARD CARD
# COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
File K - Miscellaneous Variables
11 21-30 F10.4 THN03K
31-40 F10.4 THP04K
1.12
1.084
1-8
2A4
DUM
Temperature correction constant
for NO -N decay coefficient
(A.38r
Temperature correction constant
for PO.-P settling coefficient
(A.76.1)
ENDFILE
Although no values need be specified (in which case the default values apply), all 11 cards
must be included.
File L - Plot Variables
1 1-8 2A4 DUM
11-15 15 NR
16-20 15
NIND
FILE L
Number of reaches for which
quantities are to be plotted
(£99)
Number of quantities to be
plotted per reach (<43)
-------
TABLE 2. (continued)
CARD CARD
# COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
VALUE
DESCRIPTION
Ln
2-N
N-M
File L - Plot Variables
1-8 2A4
11-80 1415
DUM
NR)
14 reaches per card are input
1-8 2A4 DUM
11-80 1415
NIND)
14 quantities per card are input
FILE L
Reach number for which plots
are wanted
FILE L
Quantities to be plotted
4 = flow at end of reach
minimum DO in reach
DO at reach end
BOD at reach end
DO at reach start
BOD at reach start
5
7
11
14
15
16
17
18
19
20
NH -N at reach start
NH -N at reach end
NO -N at reach start
NO -N at reach end
NO_-N at reach start
-------
TABLE 2. (Continued)
CARD
#
CARD
COLUMN
FORMAT
VARIABLE
NAME
UNITS
DEFAULT
NAME
DESCRIPTION
File L - Plot Variables
N-M
21 = NO -N at reach end
22 = PO.-P at reach start
4
23 = PO.-P at reach end
4
24 = phytoplankton at reach
start
25 = phytoplankton at reach
end
26 = coliform at reach start
27 = coliform at reach end
28 = HM1 at reach start
29 = HM1 at reach end
30 = HM2 at reach start
31 = HM2 at reach end
32 = HM3 at reach start
33 = HM3 at reach end
34 = total nitrogen at reach
start
35 = total nitrogen at reach
end
-------
TABLE 2. (Continued)
CARD CARD VARIABLE DEFAULT
# COLUMN FORMAT NAME UNITS NAME
File L - Plot Variables
N-M 36
37
38
39
40
41
42
43
DESCRIPTION
= chlorides at reach start
= chlorides at reach end
= HM1 ions at reach start
= HM1 ions at reach end
= HM2 ions at reach start
= HM2 ions at reach end
= HM3 ions at reach start
= HM3 ions at reach end
Ln
Ui
M+l
2A4
DUM
ENDFILE
If no plots are required, only the ENDFILE card is needed,
-------
TABLE 3.
DEFINITION OF CONSTITUENT SELECTION OPTION, ICOMB (FOR CARS 2 OF TABLE 2, FILE K)
ICCMB 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
CO XXXXXXXXX XXX XXXXX XXXXX
BCD XXXXXXXXX XX
N-I.J-N xxx xxx xxxx x x
NO.-N X X
CTv
No3-;; xxxxxxxxx xxxxxx xxx
POA-P xxx x xx x xxxxxx
PHVTO-
?L/C;KTC:I xx x xxx
X indicates that the constituent will be nodeled iinder the indicated ICOMB option.
-------
SECTION V
OUTPUT DESCRIPTION
The output generated By DOSCI is printer output and consists of the
following:
1. All of the data input in File A through File F (see Table 2),
appropriately labeled
2. A summary of the computed conditions in each reach including
(a) reach length (mi)
(b) river mile at head of reach
(c) flow in reach (cfs)
(d) minimum DO level in reach (mg/L)
(e) river mile where minimum DO occurs
(f) temperature of reach (°C)
(g) concentrations at head and end of reach of all constituents
being modeled _.
(h) reaeration value for reach (day )
(i) travel time through reach (day)
(j) velocity (fps), width (ft), depth (ft), and slope of reach
3. augmentation required
4. plots of requested concentrations versus river mile. Sample
output is shown in Table 7 -
57
-------
SECTION VI
DEFINITION OF PROGRAM VARIABLES
TABLE 4.
DESCRIPTION OF COMMON VARIABLES
FORTRAN
NAME
DEFINITION
UNITS
COMMON/BLOCKl/
CRMIN(IA)
DATA(IA,J)
J=l
J=2
J=3
J=4
J=5
J=6
J=7
J=8
J=9
J=10
J=ll
J=12
J=13
J=14
Minimum DO concentration for reach MG/L
IA
-1
-1
Length of reach IA MILES
River mile at head of reach
Decay coefficient for BOD in reach HOUR
Settling coefficient for BOD in HOUR
reach
Flow coefficient to calculate
velocity in reach
Flow exponent to calculate velocity
in reach
Flow (incremental or point source) CFS
into reach
DO concentration in inflow MG/L
BOD concentration in inflow MG/L
NH -N concentration in inflow MG/L
NO -N concentration in inflow MG/L
NO.,-N concentration in inflow MG/L
PO.-P concentration in inflow MG/L
4
Algae concentration in inflow MG/L
59
-------
TABLE 4. (Continued)
FORTRAN
NAME
DEFINITION
UNITS
J=15 Value for the reaeration coefficient HOUR
K , if option 1 for calculating K
is used (base e)
J=16 Coefficient of velocity if option 2
is used for calculating K
J=17 Exponent of velocity if option 2 is
used for calculating K
J=18 Exponent of depth if option 2 is used
to calculate K?
J=19 Coefficient of discharge if option 3
is used to calculate K
J=20 Exponent of discharge if option 3 is
used to calculate K_
J=21 Coefficient of discharge to calculate
the average depth of water for reach I
J=22 Exponent of discharge to calculate the
average depth of water for reach I
J=23 Slope of channel if option 4 is used
to calculate K
J=24 Coliform concentration in inflow MPN/100 ML
J=25 HM1 concentration in inflow MG/L
J=26 HM2 concentration in inflow MG/L
J=27 HM3 concentration in inflow MG/L
J=28 N concentration in inflow MG/L
J=29 Cl concentration in inflow MG/L
J=30 HM1 ion concentration in inflow MG/L
J=31 HM2 ion concentration in inflow MG/L
J=32 HM3 ion conentration in inflow MG/L
60
-------
TABLE 4. (Continued)
FORTRAN
NAME
FINIS (I, J)
J=l
J=2
J=3
J=4
J=5
J=6
J=7
J=8
J=9
J=10
J=ll
J=12
J=13
J=14
J=15
J=16
J=17
J=18
J=19
DEFINITION
Reach I number
River mile to head of reach I
Length of reach I
Total discharge leaving reach I
Minimum dissolved oxygen concentra-
tion in reach I
River mile where minimum dissolved
oxygen concentration occurs in reach I
Dissolved oxygen concentration at
the end of reach I
Reaeration coefficient, K , value
for reach I (base e)
Time of travel for reach I
Mean velocity in reach I
BOD at the end of reach I
Not used
Depth of water in reach I
Dissolved oxygen concentration at
start of reach
BOD concentration at start of reach
NH_-N concentration at start of
reach
NH--N concentration at end of reach
NO -N concentration at start of reach
NO_-N concentration at end of reach
UNITS
MILES
CFS
MG/L
MG/L
HOUR'1
DAYS
FPS
MG/L
FEET
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
61
-------
TABLE 4. (Continued)
FORTRAN
NAME
FINIS (I, J)
J=20
J=21
J=22
J=23
J=24
J=25
J=26
J=27
J=28
J=29
J=30
J=31
J=32
J=33
J=34
J=35
J=36
J=37
J=38
DEFINITION
(Continued)
NO -N concentration at start of reach
NO«-N concentration at end of reach
PO.-P concentration at start of reach
4
PO.-P concentration at end of reach
-4
Algae concentration at start of reach
Algae concentration at end of reach
Coliform concentration at start of
reach
Coliform concentration at end of
reach
HMl concentration at start of reach
HMl concentration at end of reach
HM2 concentration at start of reach
HM2 concentration at end of reach
HM3 concentration at start of reach
HM3 concentration at end of reach
N concentration at start of reach
N concentration at end of reach
Cl? concentration at start of reach
C19 concentration at end of reach
HMl ion concentration at start of
UNITS
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MPN/100ML
MPN/100ML
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
reach
J=39
HMl ion concentration at end of reach MG/L
62
-------
TABLE 4. (Continued)
FORTRAN
NAME
DEFINITION
J=41 HM2 ion concentration at end of
reach
J=42 HM3 ion concentration at start of
reach
J=43 HM3 ion concentration at end of
reach
J=44 Internal flag used if total nitrogen
concentration is greater than sum of
components.
RIDENT(IA,J)
J=l-5 Reach identification for reach IA
K20PT(IA)
COMMON /BLOCK2/
JNIT(IA)
F(IA)
C(D
IONE(K)
TITLE (I)
IORD(I,J)
J=l-20
JUNG (I, J)
J=l-3
INIT(I)
Option used for reaeration calcula-
tion for reach IA
Headwater identifier for flow
augmentation calculations
Not used
DO concentration in subreach I
Augmentation logic array
Not used
Run title
Reaches in stretch I
Stretches defining junction I
Headwater stretch identifier
UNITS
FINIS(I, J) (continued)
J=40 HM2 ion concentration at start of
reach.
MG/L
MG/L
MG/L
MG/L
MG/L
63
-------
TABLE 4. (Continued)
FORTRAN
NAME
IAUG(I)
DOL(I)
TRFAC(I)
COMMON/BLOCKS/
CONDZ(I,J)
J=l
J=2
J=3
J=4
J=5
J=6
J=7
J=8
J=9
J=10
J=ll
J=12
J=13
J=14
J=15
J=16
CONDI (I, J)
J=l-16
DEFINITION
Augmentation logic array
Minimum permis sable DO level
Treatment factors
DO concentration in headwater I
BOD concentration in headwater I
NH -N concentration in headwater I
NO -N concentration in headwater I
NO -N concentration in headwater I
PO.-P concentration in headwater I
4
Algae concentration in headwater I
Coliform concentration in headwater I
HM1 concentration in headwater I
HM2 concentration in headwater I
HM3 concentration in headwater I
N concentration in headwater I
Cl« concentration in headwater I
HMl ion concentration in headwater I
HM2 ion concentration in headwater I
HM3 ion concentration in headwater I
Concentrations at start of reach I,
UNITS
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MPN/100ML
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
same constituents as CONDZ
64
-------
TABLE 4. (Continued)
FORTRAN
NAME
DEFINITION
UNITS
CONDE(I.J)
J=1-16
COMMON/BLOCK4/
TEMMO(I)
IDMCH(I.J)
HWFLOW(I.J)
J=l-12
TRFACN(I)
SEASON(I)
COMMON/BLOCKS/
JJ
KK
II
B
RMLOW
MAX
CLOW
NTRIB
NREA
NINIT
NJUNC
ELEV
Concentrations at end of reach I, MG/L
same constituents as CONDZ
Monthly mean stream temperature °C
Augmentation logic array
Monthly flows in headwater I CFS
Not used
1^0 means month I will be run
Junction identifier for flow
augmentation
Junction identifier for flow
augmentation
Junction identifier for flow
augmentation
Not used
Location of minimum DO in a reach RIVER MILE
Number of subreaches in a reach
Minimum DO in a reach MG/L
Number of stretches in the system
Number of reaches in the system
Number of headwaters in the system
Number of junctions in the system
Basin elevation FEET
65
-------
TABLE 4. (Continued)
FORTRAN
NAME
DOLEV
TF
TEMP
CSA
M
DEFINITION
DO target level
Treatment factor
Average temperature
Saturation concentration for DO
Flag to indicate end of computations
UNITS
MG/L
°C
MG/L
QUP
FINL
JA
IA
ICK
FINLN
DELQ
NI
NJ
K2
VEL
NSEAS
NRUN
TFN
for a stretch
Total upstream flow entering a reach
Not used
Stretch subscript
Reach subscript
Intermediate printout flag
Not used
Incremental flow used in subreach
calculations
Input unit
Output unit
Reaeration coefficient
Reach velocity
Month of run
Run number
Not used
CFS
CFS
HOUR 1
FT/SEC
66
-------
TABLE 4. (Continued)
FORTRAN
NAME
COMMON /CONBEG
DO
BOD
NH3
N02
NO 3
P04
ALG
COL
HM1
HM2
HM3
HM
TOTN
COMMON /CONEND
DOE
BODE
NH3E
N02E
N03E
P04E
DEFINITION
Initial DO concentration for GETCON
Initial BOD concentration for GETCON
Initial NH -N concentration for
GETCON
Initial NO--N concentration for
GETCON
Initial NO -N concentration GETCON
Initial PO.-P concentration GETCON
4
Initial algae concentration for
GETCON
Initial coliform concentration for
GETCON
Initial HMl concentration for GETCON
Initial HM2 concentration for GETCON
Initial HM3 concentration for GETCON
Initial HM concentration for GETCON
Initial N concentration for GETCON
Final DO concentration from GETCON
Final BOD concentration from GETCON
Final NH--N concentration from GETCON
Final N09-N concentration from GETCON
Final NO--N concentration from GETCON
Final PO.-P concentration from GETCON
UNITS
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MPN/100ML
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
MG/L
67
-------
TABLE 4. (Continued)
VARIABLE DEFINITION UNITS
ALGE Final algae concentration from GETCON MG/L
COLE Final coliform concentration from MPN/100ML
GETCON
HM1E Final HM1 concentration from GETCON MG/L
HM2E Final HM2 concentration from GETCON MG/L
HM3E Final HM3 concentration from GETCON MG/L
HME Final HM concentration from GETCON MG/L
TOTNE Final N concentration from GETCON MG/L
COMMON/CONST/
Variables in CONST are defined in FILE K input (see Table 2).
COMMON/RCHVAR/
Variables in RCHVAR are defined in FILE J input (see Table 2).
COMMON/OPTION/
IFN Algae (phytoplankton) growth function
option
0 = growth limited by NO -N
concentration
1 = growth limited by NH -N
concentration
2 = growth limited by maximum of NH -N
and NO -N
IK2 K selection option
ICOL 0 = don't model coliforms
1 = model coliforms
ICOMB Constituent selection option (see
Table 3)
INH3 0 = don't model NH -N;
1 = model NH3-N
68
-------
TABLE 4. (Continued)
VARIABLE
DEFINITION
UNITS
IPO.
4
IALG
IFIRST
COMMON/OPT2
IHEAVY
ITOTN
ICHLOR
COMMON/OPT3
IP
INH
IN2
INS
0 = don't model NO -N;
1 = model N02-N
0 = don't model NO -N;
1 = model N03-N
0 = don't model PO.-P;
1 = model PO.-P
4
0 = don't model algae;
1 = model algae
Logic flag for GETCON
=1 means model I heavy metals
0 = don't model total nitrogen;
1 = model N
0 = don't model Cl ;
1 = model C12
=1 means model PO.-P with I'th order
4
reactzon
(I = 1 or 2)
=1 means model NH -N with I'th order
reaction
(I = 1 or 2)
=1 means model N02~N with I'th order
reaction
(I = 1 or 2)
=1 means model NO -N with I'th order
reaction
(I = 1 or 2)
69
-------
Local variables employed by DOSCI are described in Table 5 under their
respective subroutine.
TABLE 5.
DESCRIPTION OF LOCAL VARIABLES
VARIABLE
SUBROUTINE
AVK2
DELQ
DELQ1
DEP
HSUM
TOTFLO
TRAV
TSUM
VEL
VSUM
SUBROUTINE
ARR
BODMC
BODMTL
BODNWR
BODWT
DALND
DEFINITION
DOEQU
Average K_ value for a reach
Flow in subreach
Incremental flow into subreach
Depth of subreach
Sum of subreach depths
Total flow at end of reach
Travel time through subreach
Travel time through reach
Velocity in subreach
Sum of subreach velocities
GETCON
Algal respiration rate
BOD convertible to inorganic forms
BOD material in BOD decay
BOD nitrogen weight ratio
BOD weight
Algae change due to natural death
UNITS
DAY-1
CFS
CFS
FEET
FEET
CFS
DAYS
DAYS
FPS
FPS
HOUR"1
MG/L
70
-------
TABLE 5. (Continued)
VARIABLE
SUBROUTINE
DALTOX
DN
DOBEN
DOD
FACHM1
DEFINITION
GETCON (Continued)
Algae change due to toxicity
Nitrogen demand due to algal growth
Benthal DO demand
DO reaeration change
Heavy metal factor on coliforms and
UNITS
MG/L
MG/L
MG/L
MG/L
FACHM2
FACHM3
FL
FLIM
FN
FNH3
FN03
FP
GRLIM
IBAR
NOUT
OS AT
RAT
algae reactions
Heavy metal factor on coliforms and
algae reactions
Heavy metal factor on coliforms and
algae reactions
Algal growth limitation function due
to light
Minimum of FL, FN, FP
Algal growth limitation function due
to nitrogen
Algal growth limitation function due
to NH3-N
Algal growth limitation function due
to N03-N
Algal growth limitation function due
to P04-P
Total algal growth limiting function
Maximum light intensity
Output unit
DO saturation level
Ratio used in settling calculations
(LANGLEYS
/MIN)
71
-------
TABLE 5. (Continued)
VARIABLE
DEFINITION
UNITS
SUBROUTINE
TCOR
TT
SUBROUTINE
TEMSET
XCK
XDOK
XEX
XNH
XN2
XN3
XP
SUBROUTINE
A
DX
FMAX
FMIN
IPLOT
LAB
N
SMIN
SMAX
T
GETCON (Continued)
Temperature correction term
Temperature
NEWIN
Default temperature value
Default coliform reaction coefficient
Default reaeration factor
Default extinction coefficient
Default NH,-N reaction coefficient
Default NO--N reaction coefficient
Default NO -N reaction coefficient
°C
°C
HOUR
-1
FEET
-1
HOUR
-1
HOUR
-1
HOUR
-1
Default PO.-P reaction coefficient
4
HOUR
-1
PLOT
Dependent variable array to be plotted
Plot increment
Maximum value to be plotted
Minimum value to be plotted
Print buffer for plotting
Plot label
Number of values to be plotted
Minimum grid limit
Maximum grid limit
Independent variable array
72
-------
TABLE 5. (Continued)
VARIABLE
DEFINITION
UNITS
SUBROUTINE PLTSET
IFLAG
INDS
INR
IOUT
NIND
NR
TITLE
XPLT
YPLT
0 = plot;
1 = read plot inputs
Array of indices of quantities to be
plotted
Array of reach numbers where plots
are requested
0 = plot;
1 = no more plots
Number of quantities to be plotted
Number of reaches where plots are
requested
Array of plot titles
Independent variable array passed to
subroutine PLOT
Dependent variable array passed to
subroutine PLOT
SUBROUTINE REPRE
FAC
FACFLO
TOTFLO
VAR
XMAX
Treatment factor
Incremental flow for subreaches
Total flow for reach
Concentration of flow into subreach
Number of subreaches in reach
SUBROUTINE RUN ON
CHK Velocity check parameter
IR Reach number
CFS
CFS
MG/L
FT/SEC'
73
-------
TABLE 5. (Continued)
VARIABLE DEFINITION UNITS
SUBROUTINE RUNON (Continued)
QAUG Augmentation required CFS
WIDM Reach width FEET
XSLOP Reach slope
SUBROUTINE SCAN
LL Number of headwaters available for
augmentation
QADD Total augmentation required CFS
QPLUS Augmentation required per headwater CFS
Z DO deficit' ' MG/L
74
-------
SECTION VII
SAMPLE INPUT DECK
A listing of a sample input deck is provided in Table 6. The deck is for
a DOSCI run simulating conditions in River Region 2 (Main Stem and South
Fork of Coeur d'Alene River and Hangman Creek) during August, 1969. The
region is divided into fifty two reaches. Numerous point source and in-
filtration flows are modeled. Three municipal outfalls are modeled.
Many of the inflows are polluted due to the extensive mining activities
in the area. Reaeration values are calculated by Option 2 (see Table 2
File F-4). Nominal values of all reaction coefficients are used with the
exception of NH -N volitization, BOD benthal release, and benthal oxygen
demand. Plots are produced for nine quantities in twenty two reaches.
The resulting output is presented in Table 7-
75
-------
TABLE 6. SAMPLE INPUT DECK
FILE A
EN'DFILE A
FILE B
ENOFILF. B
FIL? C
FILE C
FILE C
FILE C
FILE C
ENDFILE c
FILE D
FIL& D
PROFILE D
FILF. F
n,2
0.
FILE F
8,2
0.
FILE F
8,2
4.
ENDFILE E
FILE' FM
FILE FM
.FILE FM
FILE FM
FILE FM
FILE FM
FILF. FM
FILE FM
FILF. F«i
FILE F-i
FILE FM
FILE FM
FILE FM
FILE FM
FILE FM
FILt FM
FILE FM
FILE FM
FILE FM
FILF. FM
FILt FM
FILF FM
FILE FM
FILt FM
FILE FM
FILE FM
FILE FM
FILE FM
FILE FM
FILE FM
FILE FM
FILE FM
FILE FM
FILE F.I
FILE F-J
FILF FM
FILE FM
SPUX ANF
3
1
?.
3
a
5
1 0
.ft
2 0
,8
3 0
,8
1
2
3
a
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
20
2b
2b
• 27
28
29
30
31
32
33
30
35
36
37
MIVF.R RT
2
1 2 3
16 17 IB
23 20 25
? 6 27 ?. fl
03 an 05
2
1
,03
.03
,03
6.7
.2
0.0
2.6
0.0
7.6
0.0
o . o
o. n
2.5
0.0
„.. 2 , 0
0,0
3,5
1,1
4,1
0,0
3,6
0,0
1.9
0,0
5,5
2,4
0,0
4.0
1.3
.7
0.0
2.7
0.0
1,3
0,0
2.0
0.0
3. a
0.0
.9
.GIOM Mm-"'
52
a •, 6
19 r?0 21
29 30 31
Uh
-------
TABLE 6. (cont'd)
FILE-. F-*l
FTLE- FM
FILE F»\
FILE F*l
FILE F«l
FILt F«l
FILt F«i
FILE-. F-l
FILE- F'M
FILE F«l
FILt. f-«l
FILE FM
FILE F«l
FILL F-l
FILE F"l
ENOFILE F-l
FILE F«2
8,0
0,0
FILE P*2
8,0
0,0
FILE F«2
6,0
0,0
FILE F*2
8,0
0,0
FILE F-2
8,0
0,0
FILt K»2
8,0
0,0
FILE F"2
8,0
0,0
FILE F*2
8,0
0,0
FILE F-2
8,0
0,0
FILE- F-2
8,0
0.0
FILF. F-2
8,0
FILE F»2
8.0
FILE F-2
3.0
FILE F-2
8,0
FILF. F*2
8,0
38
V)
00
«1
'12
03
uu
05
06
07
08
09
50
51
52
1.0
0,0
1.0
0.0
1.0
0,0
l.o
0.0
1.0
o.o
1.0
0.0
1.0
0,0
1.0
0.0
1.0
0.0
1,0
0.0
1.0
1.0
1.0
1.0
1.0
2
3
5
6
7
8
9
11
12
13
10
16
17
IB
10
0.0
1 .9
2,7
0.0
2.J
3.6
7,9
0,0
3,6
0,0
6,1
0,0
7,2
0.0
8,0
0.
.03
0,0
36,
0,0
19.
0,0
12.
0.0
o.
,03
0.0
o.
.03
0,0
36.
0.0
18.
.03
0,0
0.
.03
0,0
60,
0,0
0,
.03
6.
.03
6.
.03
6.
.03
6.
.03
6.9
6,9
'3,0
2.3
2.3
167.3
160,2
156.3
1-5J, .3
I'j2,7
l'S2,7
l'»6.6
106,6
139.0
1.59,0
.006
0,0
0,0
0,0
0,0
,006
0,0
.006
0,0
0.0
J006
0.0
.006
0.0
0.0
.006
,006
,006
,006
,006
,02
0.0
.02
0.0
.02
0,0
,02
0,0
.02
0,0
.02
0,0
,02
0.0
,02
0,0
.02
0,0
.02
0.0
,02
0,0
,02
0,0
,02
0,0
.02
0.0
,0?
,03
0,0
,03
0,0
,03
0,0
.03
0.0
.03
0.0
.03
0.0
,03
0,0
.03
0.0
.03
0,0
.03
0,0
,03
,03
,03
,03
.03
.31
,21
.21
,21
,21
,1
.1
.071
,067
,071
,067
.071
,067
,071
.067
0.
0,0
0,
0,0
0.
0,0
o.
0,0
0,
0,0
0,
0,0
o.
0,0
0.
0,0
0,
0.0
0.
0.0
0,
0.
0,
0,
0.
200,
0,0
200,
0,0
200.
0,0
200,
0,0
200,
0,0
200,
0,0
200,
0.0
200,
0,0
200.
0,0
600.
0,0
200,
200,
200.
200.
200.
77
-------
TABLE 6. (cont'd)
0.0
FILE F"2
8.0
FILE F-2
8,0
FILE F*2
e.o
FILE F-2
a.o
3,
FILE K«2
8,0
3,
FILE F-2
8,0
FILE F-2
8,0
1,0
FILE F-2
6.0
,8
FILE F«2
1,
,63
FILE F»2
1.
.6
FILE F«2
6,0
.3
FILE F«2
8,0
.3
FILE F-2
8.0
.3
FILE f-2
1,
100,
FILE F-2
8,0
100,
FILE F-2
8,0
5,1
FILE F-2
8,0
n w
** I 7
FILE F-2
8,0
1.6
FILE F-2
8,0
1,8
FILE F-2
8,0
20
1.0
21
1.0
22
1.0
23
1,0
24
1,0
25
1.0
26
1.0
28
1.0
30
200.
32
200,
34
1,0
35
1,0
36
1.0
38
200.
39
1.0
41
1.5
45
20.
47
20.
49
20.
51
20.
4,
.03
4.
.03
6.
.03
4.
.03
4.
.03
4.
.03
4.
.03
6.
.03
.3
.03
.3
,03
33.
3.
.03
4.
.03
.4
.03
7.
.03
6.
.03
4.
.03
4.
.03
4.
.03
4.
.03
,006
.006
,006
,006
,006
.006
.006
.006
,006
.006
.006
.006
.006
.006
,006
.006
.006
.006
.006
.02
0.0
,02
0,0
,02
0,0
.02
0.0
,02
0.0
.02
0,0
.02
0.0
.02
30,
,02
0,0
.02
0.0
,02
,02
.02
0.0
.02
.02
.02
.02
.02
0,0
.02
0.0
.02
.03 0.
.03 0.
»C3 0,
.03 0,
.03 0.
,03 0.
,03 0,
.03 0.
.03 0,
.03 0.
.03 0.
.03 0,
,03 0.
,03 0.
.03 n.
,03 0.
.03 0.
.03 0,
,03 0,
,03 0.
200,
200,
20rt,
200,
200,
200.
200,
200.
15000.
15000,
200,
200,
200,
15000,
2on,
200,
200,
200,
200,
200,
78
-------
TABLE 6. (cont'd)
2.6
ENDFILE F-R
FILE F-3
FILL F«3
FILE F-3
EHOF1LE F-3
FILE F-4
FILE F»4
FILF. Fi-4
FILE F»4
FILE F-4
FILE F»4
FILE F«4
FiLt: r«a
FILE F-4
FILE F«4
FILE F-4
FILE F«U
FILE F-U
FILE F*4
FILE F-4
FILb F-4
FILE F-4
FILE F-4
FILfc Fn4
FILE F»4
FILE F-4
FILE F»4
FILF F«4
FILE F-4
FILE F-U
FILE F-4
FILE F-4
FILE F*4
FILF F-4
FILE F«4
FILE F-4
FILF F-4
FILE F-4
FILF F-4
FILE F-4
FILE F-4
FILE F-4
FILE F-4
FILE F-4
FILE F-4
FILE F-4
FILE F-4
FILE F«4
FILE F-4
FILE F*4
FILE F«4
FILE F«-4
FILE F-4
FILE F»4
FILE F-4
FILE F»4
FILE F-U
ENOFILE F.4
ENOFILE G
4fl
50
52
1
2
3
U
5
6
7
8
9
10
11
1?
13
14
15
16
17
18
19
20
21
22
?3
24
?5>
26
27
28
29
30
31
32
33
3a
35
3fr
37
3B
39
40
41
42
43
44
45
46
47
46
49
50
51
52
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
?
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
?
2
2
2
2
2
2
2
2
2
2
0.0
2.5
2.5
2.5
G4110 TO G4U5
0411-5 TO PRITCHARD C
PKITCHAHD CKtfcK
PRITCHARD C TO BEAVR
BEAVtR CKEEK
RMVER TO GHAHM CRK
GRAHM TO STf-AMBOAT
STEAMHOAT CRFEK
STEAMBOAT TO S-F-C
8MITH-F ALL-COUGAR C
S-F-C TO N FORK
N FORK COEUR D'ALENE
N FORK TO G4130
G4130 TO S FRK COAL
M34 TO DEADPAN
OtAOMAN GULCH
DEAD TO HUNTER'S G
GOLD HUNTERS GULCH
GOLD TO SLAUGHTER
SLAUGHTERHOUSE CREEK
SLTHOUS TO CANYON C
Mfe.4 TO GfcM OUTFALL
GEM OUTFALL
GEM OUTFALL TO 3 FRK
CAMYOM C TO r,4li1,S
G4131.5 TO LAKE CRK
LAKE CREEK
LAKE C TO SILVEHTON
SILVERTON OUTFALL
OUTFALL TO QSHUKN 0
OSBUR'J OUTFALL
OUTFALL TP BIG CREEK
RIG CRFf-K
BIG CMFF.K TO MILD C
TO KF.LLOGG OUT
KELLOGG OUTFALL
OUT TO G4133
C4133 TO PlMt CREEK
PIN'E CREEK
PIk.'E C TO CDALN RIV
S FRK TO G4135
G4135 TO 4JULY CRK
4TH OF JULY CREEK
4JULY C TO ROSE CPK
ROSE CREEK
ROSE C TO KILLARNtv
KILLARNEY OUTLET
KILL 0 TO BLACK LAKE
BL^CK LAKE OUTLf-T
OUT TO CDALENE LAKE
,U7
,13/
.n;
,137
.137
.137
,137
.147
.155
,131
.162
,16fl
,112
,1 12
,140
,140
,154
,140
,154
.21
,140
.21
jl^l
,19J
,191
,191
,191
.I'M
.191
,191
,191
,194
,200
.200
.200
,215
.222
.270
,236
,1S3
,146
,600
,215
.600
.215
,600
,215
,600
.215
79
-------
TABLE 6. (cont'd)
f-tl.E H
ENOFILE H
FILE i
FILE I
FILE I
F.NOFILC I
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE 'J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
FILE J
CTI C t
FILE J
FILE J
FILE J
FILE J
FILE J
ENOFILE J
FILE K
FILE K
FILE K
FILE K
FILE K
FILE K
FILE K
FILE K
FILE K
FILE K
FILE K
ENOFILE K
FILE L
FILE L
FILE L
FILE L
ENOFILE L
M
i
2
3
1
2
U
b
a
10
12
1<(
15
16
18
20
22
23
25
26
27
29
31
33
35
37
39
ltd
HZ
«5
«a
tib
as
50
52
*
2 1 ?
1 0 1
1112
,01 ,00001
22 9
52 50 «8 Hb UU U3 15 ia 10
39 37 35 33 29 27 26 22
7 11 17 19 21 23 27 29 37
19,
117,
1.
13,
19.7
19.
19.
20,
17.5
18.5
IB, 5
19.
19,
16. U
16. 4
16. U
15.3
17,
17.
18.3
17,5
16.
16.
19,
16.7
17.
16,5
17.3
SI.
17.3
16.7
19, 5
20.
.1
.onooi
6 (a 1 U2 00
80
-------
SECTION VIII
SAMPLE OUTPUT
A listing of a sample DOSCI output is provided in Table 7. The output
results from the input deck listed in Table 6 and is for the most part
self explanatory. All input conditions are printed out. The three out-
falls are point sources 30, 32, and 38. No augmentation is requested
and no treatment factors are applied. DO, coliforms, BOD, NH -N, NO -N,
NO -N, PO.-P, zinc, and chlorides are modeled.
81
-------
TABLE 7. DOSCI OUTPUT
(• T L F * - » A ' I 'J T I T i. F * » *
F7LF 1
CtKP
TYPE
FILE B
Hi VFH d'SlN
* • * ULE n - PHYSICAL
NO. OF NP. OP NO. UF NP. uc
MFinmlFKS ,IU«J(,TInNS hF4Cr>FS STPtTcMt! Fll1* FJNAl.
M*K ne 10 MA». of in "4> Pf iftg «4x UF ?o SUMMARY ONLY
«t»N
FLFV.
(FT)
2?bO.O
oo
r-o
C4HO
TYPE
FILE
FILE
FILE
^ ILF
'.'0.
» * * MLE C - ut.ALH Ph(^tR * » *
0"OKW Uf 4LL "LACHES IN
STDtTC"
1 ? i (i b «. 7 P v in i i i? li i-'Pv/10") lMli/L5 CKG/L) ("li/L) (rC/L)
C-3/L)
CO
fie
."i^O
,n3^o
. 0 ^ *in
, ii n t n
I.F *• (1 )
."•"lO
."^11)
.0^0"
.i^i'n
- .no^o
- . " 0 ^ ()
?0^
2(jn
J g(l
• ^
•0
.01
.on
* *
-.10
. c •> -.rc - . ;• ^ - . ^ o
TYPf
CiF
»ticw
SEfTtlNG PfJ u (-PH
-------
TABLE 7. DOSCI OUTPUT
t < » |.T|.F » - HA?1N TITLE * « »
CAHR
TYPF
FJLF
NAHf OF
KIVFK riiSlN
KIVER RFI.TO*' NU>'BFK J AuT.UST . J Obi
• * hlLE B - PHV5IC«L
00
to
CARP
TYKE
FILE B
EN^FILE
CAKO
TYPE
MLt C
FILE C
FILE C
FILE C
I- ILF C
FN^FILE
NO. OF NO, OF NO. OF Nfl. UF I'JiEHT 1
FAnhJTFK^ .lU^cH^S HFAf^FS STDtTcwtf Fl!1' FJNA1.
"A< Oh 10 ^A< [^ 10 "A* Of 100 MSX UF ct) SUM«A»T ONLY
».(0. 0^
ST"tTCH
* » » I-TLE C - "t-ALH ONntR * » »
0 ° 01~ k UF ALL
M..B IN tACH S7
'-TOO
FLFV.
(FT)
1 ? 3 3b 37 38 3° 00 ill 12 -0 -0 - 0
U3 uii as 14•> u/ us u^i so si s? -0 -0 -0 -0 -0 -n -0 -0 -0 -0
***FHtD- JUfiCTIOva
C4HO
TYPE
FILE 0
HLE 0
ENOI- ILE
MJ. Op
U I' S T K c A "
S7utTCW
1
2
* * • clLt h - HtAP*«TtKS * « »
CAWD STCH AIH;M oo boo
TYPE NO OPT (Mf./L) C"(./
n: M •* - N MUP-N f,Oj-\
("6/L3 (M^/I. )
B.2D
6.£|0
8.20
pf ILE
P PMV70
.81 .1300 ,00*>1 .1^10 .0^0n -»1C10
.80 .1300 . 0n t> 0 .1210 .0^00 -.1Q10
?00."
2 I) 0 . 0
?00.0
MM] h"<;
(HUVL) C^r./L)
TtT v r,
(r r, / L ) C •• i
.00 -.10 -.00 -.10
.00 -."O -.01 -.10
-.00 -.00 -.00 -.10
/L) Clf./L)
.01 -.10
"U/u) C-S/L)
-.0" -.00
-.ui -.00
-.Li - .00
Ml) - IH]A(IM) 7"*1! n A T A C I t 6 ) • » *
C 4HO
TYPF
NT. LtN(,TMUF»lvtP"lLt bOy nOU COtF.
PF St«CH TP "L^u1 HFArtIU« SETTLING Ofj U fO
txp,
ON „• F
-------
oo
OJ
F-l
F-l
F.I
FILF F'-i
FILT F.|
FILF F-l
FILF F-l
KILE F-l
FILE
FILF
F ILF
FILE
FILF
FILF
FILF
FILF
FILF
FILE
FILF
FILF
FILE
FILF
FILE
FILF
FILF
FILF
FTLF
FILE
FHF
FILE
FILF
FILF
F ILE
FTLF
FILE
F ILF
FILE
F ILF
FILF
^ ILF
MLF
F-l
F-l
F.|
P.]
c-l
F.I
FILF
FILF F-
F'LF F.
FILF
FILF
F ILE
FILF
FILE
FILF
FILF
F TLF
F-t
ILF
tSCH FLOW 00 "PO
MO tCFbl (MG/L3 CMC,/L) CMG/L)
ACH (fILES) OF KF4CH
1
2
3
1
5
h
7
fl
t)
10
11
12
13
10
15
16
17
1"
1 9
20
2 1
2?
23
20
25
2fr
27
29
2Q
31
31
i?
43
30
35
js
37
39
39
uO
01
u?
«3
UO
US
U 6
07
Oft
t°
bl
M
b2
* F IL
"3-N
G/L) !
*. 7 201.1
• 2 143.3
.0 I V3. I
2 , ft 1^3,1
.1 19ft. 5
7 , ft 1 -.3
3. ft lbft.3
.0 1 b?.7
>•.! lb?.7
. n 1 jf..<>
7,2 luk.^
,n 1 i".o
«.n ii°. o
F FC2) - DiUCI-7)
NO,,-M Nn3-N ^na.^
"O/L) (ir,/L) l"t./L!
CUFF'.
.no"
.no"
.109
. no "
. 009
. n 08
, no"
.no"
.00"
, nQ"
.109
. 00"
.nu9
.10"
.008
. noe
.no"
.009
,0u"
.00"
.00"
.no"
.no"
.10"
. no9
. no"
.009
. no"
. no"
.no"
.no"
.no"
.no"
.no"
.no"
.no"
.no"
.no"
. no"
.no"
.no"
.00"
.no"
.10"
.10"
.10"
, no"
.no"
. no"
.no"
. n o B
.OOP
THKII D4'
PhYin
) (Mf,/L)
COEF.
- .oooioo
- , 000001
-. ooooon
-.000000
-.000000
- . onoooo
- . OOoOOO
- .o&oioo
- . OOcOOO
-.QlCOOO
-.OOoOOO
- . oioioo
— .QlfllOl
— ,000000
-.OOOOOO
- , QOOOuO
- . o^oooo
- . OOoOOO
- . r, r o o o o
- , OOflOUO
- , o o o 1 01
- . 00onoo
- . OOoOOO
- , OOOOOO
- . OnpOOO
- , o 10 1oO
- , OOflOOl
- , 0 1 0 001
- . 010001
- .oioooo
- . o 1 0 0 U 1
-.OlQOOl
- . 00oni)1
-.oocioo
- . of o no i
- . oio°oi
- . ononoi
.. 0 0 0 1 01
.. 0 0 0 0 0 0
..OOoOOO
- .OOflOOO
- .OOoOOO
- .00o"00
- . OOftOOl
- . o io o 01
- . 0 1 0 0 t; 0
- .oopoon
- . ooo°oo
- . 0 n C 1 0 1
- . onoooi
- . onoooo
- . oooooi
Titl.,").
COLIFOH"
C"PN/ 1 CO
VELOCITY vi
,1«2
.102
.102
.102
. 1 «2
.102
.102
.131
.uo
.123
. 1 06
.121
. IbS
, 101
, 101
. *50
,500
.500
, 5 0 1
025
.5Q1
, 025
,521
,001
,520
.311
.310
.310
.311
.310
,311
.MO
.311
.310
* 2o1
.310
, 2f.1
.310
,?] 1
.210
.211
.210
. 100
. 1 01
,171
, 0 h7
,0 M
,0n7
,171
,0tj7
.071
,167
OM4CTt2allHKIJ
s WMI H*2
) ("G/L1 (»G/L
.010
.111
. °11
.010
.'MO
* o i o
. " 10
,010
.oil
.111
, a i (i
.011
.010
.010
.011
.oj o
.010
.Ho'
.(110
.010
.010
.010
.111
,110
,111
,011
I 0 11
."1"
,oio
. o 10
. II 1 U
.111
,010
,010
.'MO
,010
. 01'I
,010
,010
.'MO
.ojo
.011
HM3
TfT '!
OC/L)
HUJ1
(ir./LJ
-------
00
^-2
F- 2
F-a
F-a
F-?
F-a
F-a
F-?
F-?
F-?
F-a
F-a
F-a
F-a
F-?
F-a
F-?
F-a
F~a
F-a
F-a
F-a
F-a
F-a
F-a
F-a
F-a
F-a
F-a
F-2
F-a
2
3
5
6
7
S
<)
1 1
12
13
1U
16
17
1 B
1 9
?0
aa
?3
?u
a 6
?8
30
32
3 U
36
18
39
"1
• Ub
(17
(J 9
^1
,
".0 8.00
?6.o e.oo
i b.o fl.oo
12.0 ".00
u.o a. no
u.o fi . n o
36.0 6.00
13.0 6, no
u.O fl.OO
60.0 B.OO
U.O 6.00
6.0 8.00
6.0 6.00
6.0 6 , n o
6.0 6,00
U.O 6.10
U.O 6.0o
6.0 6.00
u.O 6.00
u.o 6. no
u.o ».oo
u « 0 ^ 1 0
0.0 fl » 0 0
.3 1,00
.3 1.00
36*0 ^ » 0 0
3.0 6.00
U.O 6.00
.u i , no
! • 0 6,0o
0.0 6.00
u.o 6 . n0
u.o 6.00
U.O 6.0Q
u.O R.10
FNP.MLE
C4KO
T'PF
HIF F-3
FILF F-3
FIL' F-3
F'LF F_3
FILK F-3
FIL1 c-3
ML* F-3
FIUF F-3
F ILF F-3
FILF F _ }
F'LF F.-J
F ILF F_J
FIL' F."!
h ILF F-3
1 .00
l .on
1 .00
1 .00
l.on
1 .01
1 .00
i .on
1 .00
1.01
1 .00
1 .00
l.on
1.00
1.01
1.01
I .00
I .00
I .00
l.on
I .00
i • u o
l.on
210.00
ano.oo
i A ft
1 • U "
1 .01
1 .00
ano.oo
1 .on
1 .bO
a o • o o
ao.on
a o • o n
20 .00
*
.0300
-.nono -
-,nooo -
- . oono -
.0300
.1310
-. oono -
.1300
. n30o
- . n 0 i 0 -
.1310
. n3no
. n 3 1 o
. n 3 n o
. n 3 n o
,0300
.0300
.1310
.1310
. 0 3 n o
0300
,0300
.0310
.0300
.1300
, njno
. n 3 n o
.nino
.1310
.1310
.0310
.lino
,03no
* * FILE
NO.
OF ^
Rt*CH o
1
?
u
h
«
10
1?
1"
I S
16
1 n
ao
aa
23
25
.0160
. 0 1 0 0
. 0 1 0 1
.0000
.OloO
,0n60
. o n o n
, 0 n f>n
. OnbO
* 0 n o n
. C0f>n
. on^i
. 0 0 M
. o n 6 n
.0160
. 0161
. C i h n
, 0 n n n
. 0 i h i
ioien
, o n f n
. o o hi
!o'n»!n
. o n h n
. o 0 6 n
. ont,n
. 0 n hi
. C ">i1
, o n h n
. o n i, i
. C n ei n
F(3> •
VM Vf
0« K2
PTTiiK'-
.non
.OQl
• n o n.
.non
.non
.non
.non
.101
.001
.Oon
.non
.000
.101
,ion
.non
.oano .0300
.naiO ,0300
.nano .o3oo
.nano .O3oo
.nano .0300
,1210 ,03 00
.1210 ,0300
.0210 .0300
.1210 .0301
,naiO . 03on
,1210 , 0 3 0 C
,nano .0300
,nano .0301
.nano • 0 1 0 f
nano o 3 o c
• n a n o .o3oc
.lano ,o3oO
.nano .0100
. n a 1 o . . u 3 o n
,nar|0 .oioo
. naoo ,o4oi
.nano ,o3oo
•121Q .030C
.nano .olO1"
. n <> f u .o3ijl
.nano .0301
. oano .o3oo
,0210 ,0301
. i a " o ,o3oi
.1210 .0301
.1210 .03oC
. nano ,o3oo
DiTAU,^)
COFF. ON
V F0« "2 V
npTION-2 0
1 .ion
1.301
i.ion
3.301
3.100
3 . 3 0 n
3 . 3 0 0
3 . 3 0 n
i. 5on
1.300
3. ion
3.301
3.101
3 . 1 U 0
1.300
.nooo
,0000
.1010
.0010
.1010
, nooo
,nono
.nono
.0010
.1010
.1000
.1010
.1010
!n&no
p i) n fi
.nooo
.1000
.1010
.1010
. n o n o
.1010
. 1 0 1 0
,1000
.1010
.1010
. i 0 1 0
.1010
.1010
.1010
.nono
.1010
.1000
1 H K 1.1 0 4 T
F.XP. ON
F03 K2
pTInN-a
1.100
1 .000
1.101
I . non
1 . n c n
1 ,ncn
1 ."on
1 .101
1 . 101
l .non
1 .100
1 .101
1 .non
1 .noo
1 .100
200.0
ano.o
aoo.o
200,0
aoo.o
aoo.o
aoo.o
200.0
BOO .0
210.0
ano.o
aoo.o
P 0 U 0
ano.o
t* 0 0 , 0
200 .0
210.0
210-0
2°0 . 0
ano.o
3 fi M f\
e lp v . U
aoo.o
J^oOo .0
15010.0
200. 0
ano.o
1 5 G 0 U . 0
ano.o
aoo.o
aoo.o
210.0
aio.o
aio.o
I (1,20) »
E"P. ON
D FO" *a
OpTlOi-i-2
1 .331
1 . 131
I .3iO
I .3^0
I .3^1
1 . 3 j n
I .3iO
1 . ?^n
l . ^in
1 .330
1 .330
1 ,33n
1 . 3 i 0
1 .3in
1 ,3iO
.00
.00
.00
.00
.00
.on
.00
.01
.00
.00
-.00
-.00
-.00
- n n -
. 0 •
-.00
• • o o -
-.00
-.00
3.01
i.OO
-.01
1 fi 0 ~
1 • \! '
.BO
.65
.60
\ rt _
* J 1
,3n
110.01
1 fl 0 « 0 0 *
i. in
u.70
1.61
1.60
2.61
* *
CUFF. ON
0 F0« K2
OPT10N-3
,100
, 100
,101
,001
.100
,noo
,000
.OGI
.00,1
. non
.101
.ncn
. n o i
,001
.000
.00 .00
.10 .01
.no .00
.no .00
.00 .00
.no .01
.10 .00
.no .00
.00 .01
.10 .01
.10 - . o 1
.no -.on
.no - . co
,ny - . oO
.10 -.01
.10 -.00
, fl u - . i_i n
.00 -.ijO
. * 0 - . 0 0
. n o ^ . 0 0
.10 -.00
.10 - . on
.no - .01
,nu . .on
.10 -.on
.10 -.00
.no -.on
.no -.on
.00 -.on
.no - .ol
.10 - .00
.10 - .00
EXP. OH
C FU^ K2
npiiDN-i
.101
,ncn
.101
. on
. 01
u 1
. o I
. on
.101
.101
.100
.101
. non
.101
,000
.00
.10
.00
.10
.00
,10
.10
.10
.10
.10
-.no
-.no
-.10
-.10
-.no
-.no
-.no
-.10
-.10
-.10
-.10
-.no
- . 10
- . 1j
-.10
-.no
-.no
-.no
-.no
-.10
-.no
.01
.01
.00
.01
.01
,01
.on
.01
.01
,01
.01
.01
.01
.on
.01
.01
.01
.on
30.01
.01
.01
-.01
.;>)
-.C-n
-.on
-.on
-.1.1
.01
.on
.00
.10
,10
.CO
.no
.10
.10
.10
.10
.oc
.10
-.no
-.10
-.1C
-.CO
-.10
-.10
-.10
-.10
-.10
-.00
-,10
-.10
- . 10
-.no
-.10
-,no
- , 10
-.10
-.10
-.no
.oc
.on
.on
.01
.01
.en
.00
. v n
.01
.01
-.on
-.00
-.01
-.un
-.01
-.01
-.01
- . vl
-.on
-.01
-.on
-.01
-.Ci
- . un
- . v 1
-.on
- . c n
- . (. 1
-.un
-.01
-•""
.00
.10
.10
.10
.10
.00
.CO
.10
.00
.10
-.00
-.no
-.10
-.no
-.00
-.00
-.10
-.50
- . '"' 0
-.10
-.00
-.00
-.no
-.no
-.10
-.no
-.CO
-.10
-.10
-.10
-.10
-.00
-------
FILF
F I LF
FILE
FILF
FILF
FILE
FTLF
F TLr
FILC
FILC
F ILC
FILE
FILF.
1- ILE
F ILr
FILF
F-
F-
F.
F-
F-
F.
F-
F.
F.
F.
F-
F-
F.
F.
F.
F-
3
3
3
3
3
3
3
3
3
i
3
3
3
3
3
3
2f>
27
£9
31
33
3S
37
3"
UO
U?
U3
ad
dA
(if*
50
5?
.noo
.000
.000
.000
.000
.001
. 0 (i 0
.000
• Ou 0
.000
.000
.000
.000
-.'ooo
-.000
-.000
3
3
3
T;
3
3
3
3
\
<
3
3
3
?
?
2
.300
.300
. 3u 0
.300
.300
.300
.300
. 3 0 0
. 3 j n
.3l.1
. ?00
. 3on
.300 ,
.son
.son
.500
1 .000
.000
.100
.000
. noo
.noo
.nun
.000
. o c o
.000
. OQO
.noo
.OUO
.000
1 .000
1.000
1 .331
1 .330
1 .330
1 .330
1 . 330
1 . 330
1 .330
1 .3.50
I ,3in
1 .330
1,330
1 .330
1 .330
1 .330
1 .330
1 .330
.0-00
.000
,000
.000
.000
.000
.000
.OOP
.000
.000
.000
.000
.000
-.oon
-.010
-.000
.00"
.OQO
.000
, OQ1
,000
,000
, o 0 0
.ooC
. n o n
.100
. ooo
.000
.000
-.000
-.000
-.000
ILE
* « » FIU
THKU
CO
Ln
C JRO
TYPE
F T I F
r I L k
F ILr
FTlF
FILF
F T Lp
FTLF
F T L c
F ILU
FILC
FTLF
FILF
t | L c
r I L
F TLr
FTLF
F 'L1-
F ! Li-
F :LE-
t- T Lr
F T Lr
F 'Lc
FTLF
FTLF
F IL1
F I L >•
F TLC
F 'L-
F T Lc
FILF
FTLr
FILF
FTLC
KlLi-
F ILF
F ILF
F-(l
F n
r — d
F-d
F.U
F-d
f -d
F-d
r-U
F - a
F-d
F - U
F-d
i- - a
t- - a
- - d
F- J
F-d
F — (J
f . u
F-d
F.U
F-d
F.U
F-d
c - t
r- a
' - tl
F-d
F-d
i- - d
F-d
F. d
F-U
F.d
i.n.
OF
1
•5
d
5
h
7
R
Q
10
1 1
i -i
1 '
1"
1 1;
; h
17
1"
1"
i
,;?
^3
jjd
2<,
^ ^
f 7
(f q
d Q
30
3!
3?
33
3"
3S
npnoN
FLH
? Ud
2 HP
? H °
f BF
? bF
Ni"E
Of
110 in r,i4i is
ITCUA°U I 10 Ht»vK
AVt° fhF^K
4 V t u' T L1 i.. >•• 4 n K C " K
? b»iHM C«tFH
? f,t:
? ST
? ST
? b "
T t
C 1'
? I-
? L' d
? Ill
? L>1
? UC
? Li"1
? l'n
? bl
; bi
? M h
? bf
? UF
? t-i
? Ij •'
? L4
? L ft
2 bl
? 01
? C?
? 0'
i'J1J 1^ ^]FAMfnAT
f 4 •'«'-' A I L u' I E A
Lt"nu.'. 1 Tl~. b-F.c
1^-|. 4LI -rilMkl-i C
F - 1 T" IV rLc-\
hP«^ Tfl 1,01^0
1 3 \.. T r. <; 1 1, K r u » L
u TP ^tir"4v
A "••>.•, (,1'Lff
jn in f, 1 1 , . T L •< 1 S T,
L^ wjriF^s r, u L L H
1 n TL bl. 4"buT r-i
Aii[, ^lik-^u'^F f^FtK
Tuii.'S 1" rA'-Y'V. C
. J T j ^F.i o')TF4LL
M C" Ic At L
^• u''T F A' L T" 5 FnK
1 . L ' T T^ i, t 1 < 1 . i
I 1 . s T - 1 A >• £ C K ^
„ r >* f- ^ *
K c TL ^TLytyTrN
L rJT"iv CI'TFiLL
T 41. L ln rSsJ>", U
,ii ^ ' r j T ; A L 1
T F i I. L TO HiR CrtEtK
? uTljC0t_Fn
? hT
I, CwtFK 10 "ILU C
rcEi- .
PN t FPK
. 137
,137
! 137
. 1 37
. 1 37
. 1 ^7
1 N 7
. 1 "7
. 15S
.131
. 16?
1 b ^
.n?
.11?
. 1 -o
.1-0
. 1 '"
, 1 1*0
. 1'j 'J
, 1 d 0
. ; 5d
.211
.140
,?1 0
. : v i
. 1 v\
. 1 V 1
.IV!
. 1 VI
.1 VI
,19!
. 1 VI
. 1V1
, 1 9U
FfP.
ON i: FOR
.did
.did
. d i u
.aja
u * d
, (1 1 U
did
. 11 1 n
.t id
.did
.iild
.did
!"'"
.111
.did
.tit
.tit
.did
.tit
. t Id
, 1 1 1
. 'j 1 1
. t i t
.did
.tit
, a 1 1
.tit
.did
.d) d
.did
.did
, d i d
.did
.did
0^,'!F.L
rpTin;;-
-.000
~ . n Q n
- .0 jO
-.000
- . r. 0 0
- .000
-.001
- . 0 0 0
- . "00
- . •> (, o
-.000
-.000
-.ooo
-.000
- . " u 0
- . " 0 "
-.000
-.000
-."00
- . "oo
-."01
-."on
-.0(iO
-.000
- . n J 0
- . " ~i '"
- . i u 0
-.00"
- , " u 0
- . 0 '.I 0
- . " 0 'I
- . ' o n
- . " u o
- . o o o
-.000
-------
OO
LARD
TYPE
FILE
F TLF
f IuE
HLF
ULE
FILF F-U
FILF F-U
FILE F-U
FILF F.U
FILE F-U
FILE F.a
FILF F.U
FTLF F.U
FILF F-U
FILE F.U
FILE F.U
F TLE F.U
FILF F-U
FTLF F.U
FILE F.U
ENDFILE
CARn
TYF>F
FTLF G
EMPFILF-
TYPF
FILF H
ENTF ILE
CIHD
TYPE
FILE I
FTLF. I
FILF I
ENTFILE
"t AC^ CfLTF OH"
NUMBER HEiCTio"
J 1 .OOUO
J ? .OOUO
J u .OOUO
J h .00^0
J fl .OOUO
36 ? nILO CnFL" .200 .UlU
37 2 MLn TO i-F'LLOr.G OUT ,?00 .uin
3* 2 KFLUir.i, OUTFALL .300 .aia
al 2 PINE ChFtK J270 .uiu
a 2 ? P T *^ F. rTUCnALNHIV . ? 3 fa .Uiu
a •? 2 s F P K T n r. u i j 5 .153 . a i a
uU 2 C.ui3b 10 UJULY TKK .!<(& .oiu
u5 2 UT^ Oc Jl'L> CUEEK ,6oO .uju
UA ? tjjoLY c TO BuSt CHK .215 .uju
U7 ? RHSF rR^^k ,600 ."ia
U9 ? KILIA°I»FY uuTLtT .hoo ,'jju
50 2 KTLI o H, ntio LAKE . ? 1 ^ . a i a
51 2 biATkLiKF: OUTLET .600 .uiu
52 2 OUT TO ruiLENt LAKt .215 .aia
» • * FILE t - H»SIN CHARACTERISTICS * * *
MINIMUM iLLOi-A<>LE FRACTION REMOVED BY TREATMENT
0.0, LFvEL (MO/L5
-.0 -.0 -.0 -.0 -.nOO -.000 -.OQO
* « • FH.t M - M'-A" >-L'NTHLY Tt"PE-(ATUF.'F.S * » *
OCT "QV DFC ,TA" FFB "AR APR MAY JUN ,TUL AUC
-.0 -.0- -.0 -.0 -.0 -.0 -.0 -.0 -.0 -.0 19.0
* * * FILF I - "fJAN MONTHLY HFAnoATTR FLOWS * * »
»IQ. OF HFADwATER fLU^'S I'J rFS
ST"tTcH UCT MuV DEC .'AM FEb MAR APR MAY JUN JUL AUC,
1 -0. -0. -0. -0. -0. -0. -0. -0. -0. -0. 117,
i .0. -o. -0. -0. -o. -0. -0. -0. -0, -o. 1.
j .0. -(!. -n. .j. -0. -0. -0, -0, -0, -0. 13.
Rt«C^10N HFiClTON " t 4 L T I n N SFTTi. TN
COFF COEF CLFF C^tF
.OQUQ .0150 .001U .OnOQ
.00^0 .0150 .00 1 a .OOJT
.OQUO .015-0 .OOlu .000''
n LHLF FU"
OPTION 5
.ouoo l.oo
,OyOO 1.00
. 0 COQ 1.00
.0^00 1.00
.Our>Q 1.00
-.000
-.000
-.000
-.000
-.000
-.000
-.000
-.000
- .000
-.000
-. 000
- .000
-.000
-.000
-.000
-.000
-.000
-.000
-.0
SEp
-0.
-0.
-0.
HEAVY HrAVY HEAVY
KF.T t «FT Z NET 3
COE.F CrtF t^EF
.00^0 -.0000 -.0000
,00"0 -.00 0 o -.0000
.OOUO -.0000 -.OCOQ
.OOUu -.0000 -."000
CCtAI
UEG)
1".70
19,00
lc."o
17.50
-------
HLF
ULE
FILF
ULE
HLF
t ILE
^ ILF
FILE
HLF
FILF
ULE
ULE
FILE
ULE
HLF
FILE
ULE
FILE
FILE
FILE
ULE
HLF
HLF
HLF
ULE
ULE
j
j
j
j
j
j
j
j
j
j
j
j
J
J
j
J
j
J
j
j
j
j
j
J
j
J
II i
10
1?
1"
is
16
IP
2(1
2?
23
25
i*
t>7
2"
31
3?
3S
37
3«
UO
u?
ui
UU
Uh
08
bO
b?
.OOUO
.OOUO
. 0 0 u 0
.oeuo
.OOuO
.oouo
.OOUO
.OOuO
.OOuO
.oouo
.OOuO
.OOuO
.OOuO
.OOuO
.OOUO
.oouo
.oouo
.OOuO
.OOUO
.oouo
.OOuO
.0"UO
.oouo
.OOUO
.OOUO
.oouo
.Oouo
.OOUO
.oouo
.OOUO
.oo"o
.OOUO
.OOUO
. 0 0 u o
.00"o
.OO^O
.OOUO
, 0 0 u 0
.COUO
.0010
.0010
.OOUO
.00"0
,00uo
.oouo
.0010
.00"0
. 00"0
.OOUO
.oouo
.0010
.OOUO
0150
OlbO
OlbO
0150
0150
0150
0150
OlbO
01 bO
OlbO
OlbO
OlbO
OlbO
OlbO
OlbO
0150
0150
0150
0150
01 50
OlbO
OlbO
0150
0150
0151
0150
.OOIU
.OOIU
.OOIU
.ooiu
.ooiu
.ooiu
.ooiu
.oolu
.ooiu
.ooiu
.OOIU
.OOIU
.ooiu
.ooiu
.ooiu
.ooiu
.0014
.ooiu
.ooiu
.OOIU
.00! u.
.0011
.ooiu
.OOIU
.001U.
.ooiu
.0009
.0009
.0009
.0009
.0009
.0009
.0009
.00(59
. 0 0 0 9
.0009
. 0 0 u 9
.0009
.0009
. 0 0 O 9
.0009
,00u9
,0009
.0009
.0009
.OOC9
.0009
. 0009
.0009
.0009
.0009
.0009
.OUOO
.OuOO
.OuOO
.OuOO
.OuOO
.OuCO
.OuOO
. Ou 0 Q
.OuOO
.OuCO
.OuOO
.OUOO
.OuOO
.OuOO
.OuOO
.OUOO
.OuPO
.OuOO
.OUOO
.OuOO
.OUOO
.OUOO
.OUOO
.ouoo
.OuOO
.OuOO
1 .00
1 .00
1 .00
1.00
1 .00
1 .00
1 .00
1 .00
1 .00
1.00
1 .00
l.oo
1 .00
1.1)0
1 .00
1 .00
1 .00
1 .00
1.00
1 .00
1 , 1)0
1.00
1.00
1 .00
1 .00
1 .00
.oouo
.OOUO
. nouo
.oouo
.OOUO
.0010
. 0 0 UQ
.OOUQ
.COUO
. o D u o
.oouo
.oouo
.OOUO
. 0010
.oouo
.0010
.0010
.0010
.OOUO
. OQIO
.COUO
.OOUO
. 0 0 << 0
,00U0
.0010
.0010
-.OQOQ
-.OOOQ
-. "000
-.0000
- . o o n o
-.1000
-.nooo
-.0000
-."000
- . 0 0 0 o
-."COO
-,"OOQ
- . 0 0 00
-."000
- . o o o u
-. OQOO
- . 0 ij 0 0
- . "000
-.0000
- . o 0 o o
-. "0"0
-.OQOO
-.0300
-.0000
- .0000
-.0000
— .
- .
-,
- .
- .
- .
- ,
- .
- .
- .
-.
- .
- .
- .
- .
- ,
- .
- .
- .
- .
- ,
- .
0000
0000
ocoo
" n 0 o
nooo
0000
onoo
0000
1000
0000
ocoo
0000
1500
oooo
"000
"COO
00 00
OQPO
0000
"COO
0 0 " 0
ojno
0000
0000
0000
0000
18
Ifl
19
19
1 ft
1ft
1 ft
15
17
17
IP
17
1ft
1 ft
1 <•
1 ft
17
1ft
17
21
17
1ft
19
1°
19
2"
.^0
.50
.00
.00
. 1 0
.10
. U 0
.30
.00
. Ol;
.30
.50
,00
.00
. oo
.70
.00
.bO
.30
.00
.50
.70
.20
.?0
.^c
.00
00
-------
» * * FILE * • MISCALLAVEHUS VARIABLES » * *
THE FOLLQWIIS COMSTITU£MT5 ARE RElNG "u"tLtri DISSOLVED OXYUFN
CULIFORMS
nyn
NHl-N "OOElEn HY 1ST ORPFR KFACTION
MQP-N MODLLEO RY IST OPQFH
PO"-D MOOtLtD BY ?ND ORpEH I
1 "tAVY MtTAL (AND ITS ASSOCIATED IUN)
CMLURIOtS
OO
CD
^.OtF^ lCTt«T Okl BOD IN CntlFnf'" CAI.CUL'TIO'I =
l-T£Fh ICTL^T 0N ^L'vY METAL I IN c^LUIW C Al (. I ]L » T I ON c
nFAVy KET»L 1 CTNTL' T- UHy kf-ifTIO" COEFFICIENT I
bCj CVYn^r n.j^TIt"! =
NDS-hFfPAC1nKY FWiCTION OF PKT. A^IC uATt»14L =
c.ihRo>j TO ^«u5P«cDus RATIO j K unu =
MTBor.tN To pwo^PHO'uS osTjo IN bnD =
U«v kEIGhT r^JCTI?'. OF CABRl^, IN "on =
lFM"r.BATuOt (.nh-FtfTTOM CCVaTfJT F^K f.Hj DF(;AY CO t FF I C 1 1 '•'T =
'i =
F^-< N"3 VOi, I T T z t T I Qv
bid nFNThiL H'LFAft Ktlf (r'i: / <; rjl) A " t MFTf^'-h".') =
r -.pr,fH"uOU" Up-'ITHJI 5tLt«SF - l T l ( "[,/OIIHnF "tTtB
'•ITcur,tM PE""THAL Dtl£ASF pATt ('t!;/a(JU4^F
D E ', T h 1 L CXYr'tN nt^ANy (H^/il.'ut^E ftH^-W
l-aACTION OF MEAVy PfT»L 1 I>J ION CO^M =
C.NDF ILfc
1 .
.0^000
2 0 . 0 n 0 c 0
1 , U 7 Q ^ 0
1 . •> c 0 r 0
. b 1 0 f> 0
1 c . u '1 0 n 0
1 . 1 n 0 f- 0
!;7o^o
i. i v o
. i ? 3 <• c
i ,•-•--r c
. * ten
.00o'1l
- , 0 0 0 A 0
-------
» * » FILE K - MJSCFLL A«.'EPUS VARIABLES * * *
THE FOLLOWING CONSTITUENTS ARE BEING "u"tLEf ni^S^LVEn OxYOFN
itp BY 1ST OP0ER Hf1
'-N MOOLLED BY IST UBDFH
u-o MJPtLtP BY tHO OPfEK RFACTIO^
Jt»VY MtTAL tA«D ITS ASSOCIATED ION)
C9LFI-
CO
00
OM POD IN COLIFOPM. CA[ :HI.ATIO'I -
o" HEIVY METAL i IN cnLunw" CSLL
«ET'L i cr\rtMi"tTin', LTMIT (^I,/LI IN C^LIFOH- CALCULATION
OF
"QN TO HHUSr'"C!Do5 ^ATJO IN onu i
Bortw To PMOSPHOSUS BATIO IN I,PD =
.FIGHT Tn|acTIO\ OF CA"nrS I M "I* =
C I tk'T =
v PnOcFSS
T 7 U'J r.U'.'STA"T K?* ^nl 'JfC'Y C°t r >• ' C T t^T =
t KIT*7 t '' r; / ^ U '.! A D t "F T >•" ^-h1-' ) =
~"tLt .SSF '-'iTt ( "|, /iOll (i ^F "t. T t "-K1 1 =
t>FNTH»L
hBACTIO'4 OF hEAVy HFJSL I 1'' IOM FO°M =
L.NOF ILL
2 0 . 0 c 0 c 0
1 . u 7 o " 0
1 .i^O^O
. b ^ 0 n 0
106.o"0"0
I c. . u 0 0 n 0
. b n 0 " 0
i . ; * (.• t o
. 0 n u o 1
- t 0 fl 0 A u
-------
FILL I --- PLUT VA«l«bLES
""ILL BE pLOTTtD FOR THF FOLLOhlNC. "E'CWtS
52 SO ufl ut DO Ui 15 !« 10 B 6 1 «Z 10
39 37 35 33 29 27 2* ?£
THt FOLLOWING QLlAMT!TT£S "ILL Bf
TO AT Bt6CH END [wr./
°00 4T rfFsCH LVD (
Nn3 «T HFTH fjp {^li/Ll
NO? AT HFiCM E'JD (Mb/L)
••'OS AT HfiCM LNO (f'G/L)
POU 'T hfiCM E.VU (UG/L)
COL AT HEACH t*O fu/L)
HMJ AT HEATH t^D
CL- AT KEiCM fO
OO
-------
KlvFH REGION N'U"BEX 2 AfGUS I t
NO.
Of
NUM^E" OF RUN = fl
SFASO^l OF yP. = »UF
IDE*ITIHC>TION HIVER MILE
OF TO HFAP
REACH PEACW
i
2
U
6
8
10
1?
I"
15
Ifr
18
20
2?
23
25
2*
27
29
31
33
35
37
3°
an
i ?
"3
uu
Lh
UA
bO
b?
r, u 1 1 0 TO G 1 1 1 b
C. u i i b TO FR1TCWARD C
P R I T r. H i R 0 C TO SEAVR
R i A y E H Tn r, j s H « C »K
r, K 4 M * TO STtA^GM
SUAUL!OAT TO S-F-C
S-F-r TO ^ f^R<<
N FOBK 70 r, j i 3 n
K«131 TO S t- 3* C"Al
"31 TO OFAD-AiN
Pt »o IP WU'ITH* ' s G
r.LLO TO SL4ur.MTj18
? L T " n u s TO c ' •< Y G >>' r
"6." 70 r,f< nuTULL
r.E" OUTULI TO 5 I-PK
r A' rn\ c jo r.Qi ii .s
<"• u i 3 ; . s TU L » * F r h K
L4Kt C TO MI vf '•TON
nu T K A LL TU r^Hi'ww o
OUTI «LI TO I-TL. cutr*
RIC riEiK TO -iluo r.
"lio in K 1 1 t P G r. OUT
OUT TO f, 01^3
r. » ; 3 3 TO Pl.t- r R ' t. K
p I ' L C f n C J 4 L^ ulv
S r h ' TO t>"l3b
r.ut 35 TO u-K'i Y C.RK
uJI'LY C TO ^51 TI.K
DL'^L C 7n KtLL'K"L"
"ILL 0 Tp [ILK* 1 4"t
OUT TO C0«[ t^-'t L*KF
Of- HEATH
?00
193
193
190
16?
17"
17fr
17?
10f
3"
30
26
21
h
a
19
17
17
1"
13
1'
7
^
S
?
167
1 t.U
Ibfr
Ib?
1 uf.
139
.0
.3
.1
.5
. 0
.9
, "
. a
,9
.?
.1
.5
.6
, u
.0
. 1
,R
. 1
. (1
.2
.?
.fl
,9
.0
. 3
. fl
.2
,3
.7
. 6
.11
RFACH FLO"
Lt^GTH RAIF
(MILES) CCFS)
6.7 117
.2 121
2.6 157
7.6 If?
u.O I5b
2.b 231
2.0 253
3.b 317
1.1 317
U.I 7
3.6 19
1.9 ?9
b.S 39
2.U 17
U.O 2b
1.3 66
. 1 6H
2. / 7u
1.3 7u
2.0 7u
3.U lib
.9 1 1 V
1.9 1 ? 7
2.7 127
2.3 133
3.b u^O
7.9 uSo
3.6 t5u
h. 1 U5B
l.d Ihi
0.0 "66
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.3
.0
. 0
.0
.0
.3
. o
. b
. 6
.0
. 0
.0
.0
.0
.0
.0
.0
.0
TANC-tT P.n. LFVFL = -.00
MEAN TEMPFNATl't'U (CEM)
(iir./L)
H.27
H. at
». 36
11.39
1) . UU
a. of
(i.bf
im"
H.55
11.9?
II . 7 1
(1.79
II .97
fl. 3J
U.frc
«.ac
ti. a F
(i . b"
H . 9°
a.h*-
6.61
'. . 7 c
£• . / c;
f-.l'c
{ .bi.
b .b7
a.7?
fr.7r
e.bf'
t, u1.
t .31
P.n.
199.3
193.1
192.8
162.9
10?. 5
l7K.fr.
17h.?
17?. 0
le«.1
33,8
2'.7
26.1
2".l
6.2
3.6
17. fl
17.1
16. «
1".3
11.2
1 n.9
7.7
6.7
1.7
. n
167. a
lei."
152.7
1^6.6
1 3°.u
131.1
(MG/L)
b.20
8. ull
e.3i
a.uu
fc.3"
6.b9
e.b«
e.uc
8.b5
£>. 1?
B .bi
e .71
b.ol
6.19
6.61
.2
1 ''O
? 1 ? 0 . 0
-------
S P 0 < A N t KlvER
NUM"EC OF RUN = 0
5FA SON OF Y 9 , = A UP
NO, I PEN T I f- T C A T I Ow R
OF OF
REACH PEACH
1
a
6
s.
10
I?
1"
15
lh
1 ft
c?0
27
cl3
75
f 6
^7
29
31
3?
35
37
3°
41
47
4j
4U
46
U fl
bO
b?
RU110 TO GUI lb
G Ml 1 5 TO PR1TC"ARD C
BHITCHAHO C TO 3EAVR
RtAyEH TO r.JlM" T.RK
T, KiH" TO STtAMRQAT
Mi-A^uOAT TO S-F-C
S-F.r TO N f f W K
N FCPK TO K.,130
n 4 1 3 1 TO s (.an CPAL
w 3 u TO u F A n •< 4 N
OEA0 TO ^U'JTtR'S G
r.ULB TO 5L4ur.MTto
SLTH^uS TO r.fiYC* r.
* b . 4 TO <~,t" 0 u T (• 4 L L
f; E M o ij T f. ft u TO 5 ^ P K
TAfyOK, c TO C- 41 31 .5
i". u 1 3 ; . 5 TU \_li\f C h K
L 4 H t t. TO 5 J 1 V f "• T 0 'J
0 u T f- A L I. TU r ^ H " r^w I"1
nuTt »Ll TO bib CDtcK
RIC nEtK To MLO r
^ 1 1 u TO K 1 1 L P r, r. OUT
HUT TU i; (11x3
n u i 3 3 T u Pl.t r u 1 1 K
PI' L C in r JAUKI ujv
5 F k x TO C> 4 i l b
f. U 1 3 S TO 4 J 1' I y CPK
iJji'LY C TO yHS' TkK
°yst C T" «tl_L«Kv'tY
"ILL 0 TC HL'C" I4"t
HUT TQ COSIt^'E L'RF
REGION NU
IvER MILE'
TO HF AP
Of HEATH
700
193
19(1
16?
17"
1 7ft
17?
1 e*
3"
30
26
£(1
6
U
19
17
17
1"
13
P
7
h
S
2
167
1 6"
1 bfc
1 b?
1 ufr
1 3°
.0
.1
!5
,9
.9
, U
. a
,9
.7
. I
.5
. 6
. u
. a
. i
.«
. i
. u
.2
.7
.A
. 9
.0
, 3
. fl
. 2
. 3
. 7
. 6
.11
"BE* ^ AuGUSI . (969
TAKRtT
MEAN T
RFACH FLOW D.O. RIVER
LtN(,TH R4TF HIM. AT "
(MILES) (CFS) (:ir./L) P.O
6.7
7.6
".0
2 .b
2.0
3.b
1 . 1
«.l
3.6
1 .9
b.b"
2. a
u.O
1 .3
, /
2. /
1 . 3
2.0
3.1
. 9
1 .9
2.7
2.3
3.6
7 .9
3.6
6. 1
/.2
tl.O
117
171
157
187
1°5
21 1
2^3
317
317
7
19
79
39
17
7b
66
hH
7u
7u
7u
1 lb
1 19
177
177
133
0 = 0
USo
L S U
USH
Uh2
466
.0
.0
.0
.0
.0
. 0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.3
. 6
. b
. 6
.0
. u
.0
.0
.0
. 0
.0
. 0
.0
H
II
II
H
I)
U
1!
U
II
II
II
11
II
(i
i)
II
;i
(1
H
U
b
t,
{,
f,
f
b
t
C
*•
e
t
.27
, 1.1
.36
.39
. on
. bd
.b"
. 09
.S5
. 9?
.71
.79
.97
.3?
.6C
.«c
.t)f
. »<•
. 99
.(,*.
.bl
.7C
. I'-
. L»?
.bl.
.57
.7?
.7?
. b'
. U7.
.31
199
192
te?
IB?
I7P
1 76
17?
IB«
33
2°
26
2"
6
3
17
17
1 6
1"
1 1
1"
7
6
4
167
1 b^
15?
1 ^6
139
131
P.O.
MILE
IN.
.1
.8
.9
.5
. t,
,?
.0
.1)
,S
.7
, I
.1
.2
.6
.ft
. 1
t ft
, 1
.2
,9
. 7
.7
.7
. o
, 4
.4
,7
. 6
. a
.4
LFVFL = -.00
AH'kE =1(<.35
DlSSuLvFp OXYGEN irnp
AT START 4T t«u (CENT)
b
B
6
e
b
6
t
b
8
e
b
6
b
e
6
g
b
B
P
9
6
b
0
H
e
b1
t
e
b
fr
t)
.20
.3"
.41
.3"
.b9
.bH
.-P
.b5
.1?
.bl
.74
.81
.19
.bl
.O'
.tl°
. ol
.9«
.0"
.b7
. /"
. 'If-
. P?
. /"
.b5
.bf>
.09
.7?
.bfr
.47
fl t
? .
fl f
K ,
6.
B.
P.
9.
f .
P .
<5.
» .
A t
fl .
fl f
<3 .
a ^
« .
c.
R .
R .
fc .
« .
R ,
fi p
q t
P .
f
f- .
44 19 /
44 19.0
49 19.0
39 70.0
70 1 / . b
62 : o . b
f- 0 • !: . 5
5b ;°.0
C6 19,0
00
99
94
?3
77
"b
P Q
« T
" 1
"3
"6
?1
0 1
f-f
t . 4
O . 4
n . 0
b.J
'•.0
7.0
- .3
/ .b
t) , ^
c.C
tf .0
o.7
/.('
o . r
'3 1 / . 3
* « ? ) • l'
6b 1 /.b
°U !o.7
72 '9.2
5b 1^.2
43 •> . b
31 ?u.O
-------
N
i
Nn.
OF
°tiCH
1
2
u
h
R
11
12
1"
15
lh
11
21
2'
23
as
a*>
fi
c.°
31
31
35
37
3"
"•>.
U3
a u
M»>
un
bO
S2
SPOKANE KTVEK R
HNXt" OF N"N = 0
EA^O'J f> Y". = A.ur.
TtENTIFJCATIOM
OF
Rt AL'«
Cull" TU l> u 1 ; 3
R u I l 5 TU PB:TCW»PI> c
PPI jr^ixn c TO UEAVH
P t. i v F R Tn r. tf a H w r ^ K
C.wAp1' TO :>TL«:'PL1»T
STFt"n"AT TU S-F-C
S-F-C TO N f 0-ix
u F u o n in i-, u 1 5 0
P u 1 3 n TO S ia* C n A 1
w 3 U TO U F A ^ V 4 M
TL»0 TH *•)>.;<:..< I S ,
r. LA n T n 5 1 4 >.; r. ,- T t "
SLT"nOS TO C A .vjN C
"b.1 Tl PL11 "un 4LL
r, 1 1 r>uTMi L TO s r UK
rt'ly^N C TH r"J 1 31 .5
r- j i 3 1 . 5 Tn L ' • ' r - K
1 i't L T" Ml. VF«T(,M
niTf- JLI Tu i'S '•'-" i
CuTfiLL TO oil, cctfK
^1^ CkFt^ TU MIL^ r
"•'ItO Tn «tLLP(,r. HUT
^LlT To L,'M13
c i • t c T^ COAL- °iv
f f k K TO t u 1 1 3
r-ul35 TU u.lf'y C°K
U JULv C Trj >,0b
. «<"
.«7
.°2
."1
.°1
1 .70
2. "7
1 .°1
1 .»/
2.ua
? • S H
2.10
1 .20
1.13
1 .23
i .2b
1 .11
1.11
.13
.32
.31
.11
.2b
.25
.20
.37
.17
2.22
1 .23
. b
. 2
. /
. a
. 9
. n
• b
. •>
. »
• b
.Jb
.23
• "?. "4
.24
. ! ?
. ! 7
.1 ;
. " ^
.^b
.PS
T4kPLT n.O. LFVFL = -.00
nFiK) Tt"PFK«]llkF =16.35
T^AVtL t-'EiN "t»S i'1'1" ;lcik-' SL^kc
T)"f. VFL^CTTV DFPTH U1 1 n T u (z;,^;." r-.^
(uirS) (FT/SEtJ [FTI (^j M4,v.;- u:.ra)
.-02
.112
.118
. J«0
.at i
.111
.l"i
.1"!
. Oliil
.2?i
. :iv
.070
. 179
,0°2
. 128
.O'lb
.02U
,UnO
. 0 « i
.t^b
.112
. U29
.07',
1 ("I ^
• 1 i' O
,kiO«
.175
• >»-b
.c'-l
.«'IQ
.bib
.b71
1.C2
1.C3
1 .15
1 .2?
1 .15
1.17
1.1"
l.b?
I.b2
1 . 1*-
l.b"
1.67
l.t«
1 .bO
1.91
1 .7h
1 ./B
1 .(•"
! .»u
1 .65
1 . (i 5
1 .KB
1 .3"
1 - b 6
l.S"
1.25
1 .25
.his
.h5
.h5
.BS
. o •) 1 1 e . L- 1
.Qli i 1 e. 7^
1.11 122. *.*.
1 . ' " J2v.br
l.lj lit,.?!
!.'->, 1 1 i . 7 "
l.hu 11c'.--r2
1 . 2 l 1 7 c . b u
1.22 171.7?
,2b jM.il
. a 9 a - . b?
.to 2 C . < <•
. ** ?..,-•-
.S3 lO.jH
.77 1 e . v ^
1 ,ce 1b.ec
1 . 1 f i j . - ."•
1 . ' 3 I •> . - 1
1 . 1 u i b . " -
1.1- 15.M*.
l.lrt ^b-.S''
t ."b li. /°
1 ,S0 52. 1
1 *• S u ', < S
1 . 79 U; . /a
V."2 !=•>."'
1 . H j 1 <= i . r ^
2.7i ! n t. . ni
>.Tc I"".!-
a. 73 I«V.CM
2 . 7 a i = v . 7 1
^ « l 7 ,, *.
. 0 "• .' 7 7 1"
. u " 0 J i !
,C''u;o5
.COL,'' /I
,,'rc5r =
. U " u 5 M =
. u 'M ' n 0
.C"! '(-2
. i in' /2
. '.' ••-,'11
. u r 3 ° .J u
.1 "c3 ,0
. i, r i i „ 7
.C"J717
. •"• 0 2 '• u 1
. L " <- ' ! =
. ll •' t " D u
. i. o i r r ;
. u ' f • l n
. C "' i ^ .'i -
. •„ r : r- b o
. v r (. '- 3 1
.1 " 'j '• - .1
. ; " c. - 1 •=
. f. c J - C (•
. V, V 1 J--
. i1 ^ 0 1 J >•
. b - ;• 1 j 7
. o n u i i 7
-------
s P o K 4 M E. H i v E n K n; i o N N u u e >•- H d «ur.usT«i«6
27
2"
31
33
35
(7
3"
u"
i.2
«3
u u
ab
ufl
bi
S?
BL4C»
r.uiifl TC r.ani
r. u 1 1 5 TO pommel) c
om7CH»wn r TO HFAVH
n E. » v E K T n r-HtH" c^x
ft k a H M T a :-iTtAMHQA7
M r 4 " l; n A T TO S - F - C
S-f-C T0 , MH<
>• FORK Tn r. u 13 o
r-uiio TO i (.!»«, cr>4(_
" 3 ^ TC l)F £ n" fi N
lUlD T" "r'.' T f .->' b ti
r.OLO TO 5Lij",.TLu
Sl_Tt-^u1 Ti) f.VjM r
"o.u Tn r'tM HUT-SLL
r.t'i nuT(-iLi TO 5 fa*
ri,NYn>^ c Tn f. .j i 3 ' . 5
r-uiii.1. TO L ' « ' CM"
It-t L Tn iHvr"TC,K'
Tu'MuL TU uSiii'^t. n
ru'(-«i.L T yio CULEK
nir, c^rLK T0 ..,7Lr C
"li n in »i\ tit.r, nyT
nuT TU t. '. i is
r.a.j3 To Tjt r^FtK
"pit L Tp rL, iL'« »lv
S FiK TU III! 1 3b
r.ui js TC * JjLy C0l<
£|J''LY C T;; s'lf1: Ch*
D[.it (. Tn '
KILL 0 TO •*L-L'< LiiKt
OUT TC C n 1 1. 1 "• I L » * E
ffll./1-)
.030
.012
.010
.007
. 0 1 ^
.0014
.Oii
.ecu
. 00
.030
,C2B
,0?3
,0?2
.0?2
.0?0
.019
.012
.012
.012
.013
.013
. OOt
. OOb
.o«s
.C"i
.01"
. 01"
ft/L) tuO/Ll
.012 . 1 u (•
.013 . 1 (.• t-
.008 .1^?
. 01 14 . 1 ( u
. 0 1 u . 1 i U
.003 .103
.oni .i"3
. 0 0 u .nO?
• OO-i .1(2
,U?b .Oc*«
.02s ,l( f.
,0?b .1C h
.022 ,nc«
.027 . 0 (. S
.02- .0£fc
.022 .10^
. 0 ? 1 •' 1 C f>
.020 .1C»>
.Of- . OCft
. 0 1 S . 0 C (-
.012 ."Oa
.0)2 .1C"
.013 . 0 ( "
. n 1 2 .f>i,o
. u M . n C "
,CO^ .01.3
. u 1 b . n o 3
, 0 0 u . 1 1' 3
.CO*. .003
.OOu . PO?
. J 0 i .00?
* * » »
|4KC-tT D.P.
M f A '. T t "• P F p
t ^ T R A 7 : P ^
47 £.NO
Cf-f-/L)
.CD*
.106
.00"
.003
.0^3
.003
.003
.00?
.00?
i ^ 0 f*
.00*-
.Ou^
. 1 0 ft
. 0 0 A
• PC»!
.PO*.
.no*
,0oh
.ftj*.
. POh
.10"
.00"
. "C"
.PO"
.10"
.103
.ro3
.103
.CO?
.00?
.102
LfvFL = -
/.M"k'F sljl
I
VQI-M rc'-x
4T STi"T
(Mf,/l )
.020
,C?1
,0?0
. 0 ?0
.1.21
.l.?l
.021
.u?l
.0?!
i 0 ? 0
.C?u
.C?0
. C'?0
.o2o
.020
.020
.o?w
,0?0
.c?o
.0?0
.c?o
.0?0
.0?D
.020
- . 0 2 0
.021
.C'l
.C21
.021
.O'l
.C21
.01
.iS
Ei>TH]lo1' pnu-h- LO^TLf TWtTjr-.
;T F'.n AiFTf^T JTi'u
l"i-/u' C'\-/.! (-'/D
.l-?l .IS" . " 5 J
.t?l .r=i . .-so
,0?1 »i^E ,fu°
.0?: .''-7 . fiur
.0?l . •• . " .i-l
.v.?l . r-7 .1."
. U ? 1 .1*7 . ^ 7
. u3 1 .1-" . "ac
. r ?i ."'•'' . "uK
r?i' ^L? ^3^
. 1, 2 ,'. . ' 3 " ."'3^
.02;: . 1 3 - .I;-K
. i ' r . i - " . " •> 7
. i, 2 0 • ' - 3 . r. u "J
. 0 "1 .13° .1^r'
.t,?l' .'3'' .'i5
.u?i. .r-.f .."3°
. L -' ,'• . '• > 1 . 1 J "
. u -1 1 . r J n . 1 s >"
.!;>!• .1^7 .P.<7
.H'P .'37 .'v7
• l/2ly .1^7 .-.-7
. r. ' u .'37 . n _i e
.u'l ."i" -1-5C
. i! ' i .1-3 ,p»li
,C2l . " - « .1^7
. c 2 J . r, - 7 . r u =
.c'l .'-- -osr
. 1 ' •! ,lii .'^•:>
. .' 2 j . ' i a . r s •-
-------
'KlVFH HFGIQM K'UMbFM
iU<".UST.
NUMBER 01- RUN c 0
SEASON P> Y". a AUS
NO. IDENTIFICATION
OF OF
REACH Dt'cw
1
2
U
h
A
11
12
1"
15
I'-
ll
21
25
23
25
2*
27
29
31
33
3S
i7
4T
no
tl?
U3
oo
uh
0»
bO
53
GUI to TO BO iii
r. u 1 1 5 TO p9iTcwA°D C
PHTTCMAkO f TO BEAvK
B t » v E « Tn GUAM" r K K
r, w A H M TO S T c i •*! n 0 A f
MEA"br>AT TO 5-h-C
S-F-C TO f. KHUK
N FORK TO r, a 1 3 n
("•U130 TO S F"K C^AL
M30 TO DFjnr SM
PfO Tn Wn'i-yc-^is (,
r.ULu TO Si'ur.MTt-u
SLTnCUS TU CAl;YiJ" C
v b . 4 TO r, £ H c u T I- a L L
r,t'< OUTF»Ll TO S F"K
CAMTON C TO r.in 31 ,s
f' 013), 5 Tn i '• < F (" ^ *
uxt c m s i L v F * T o >J
0 u T I- A L L Tu o S « i • - M n
nt'TclLl TU rfll. CPLI-K
BIG CKftK TO KILO C
"It 0 Tn U]33
r.ui J3 TO PUE rut- LK
PT-'E C TP ro4LN RIV
5 F,vK TO GUj3b
r.ul 35 TO oJuly CRK
OJIlLY C TO HOSE OH*
"U^E C Tn KJI. LA"
1 .tlu
, .000
.OnO
.ooc
.000
. 0 " 0
..OOu
.000
.000
.000
. o 1 0
. u 0 j
.t'Oo
.ill.
3.7^b
2.9 = 3
l.JV,
1 . J 'J H
I . J'O
1.1""
l.i.°i
. c" 0 u
. /"b
b.b3?
b.-|7]
B.3^;
1 .0" /
I./PC
1 . 7 *> 1
1 .0B7
1 . b"b
l.i'u
nFT-2 cnsc
AT PTAhT
("0/D
.000
."CO
.TOO
."0"
, oco
. 100
,"00
• o u o
."01
* " L 1
. "^ 0
* 0 ;.! 1
.00"
, r C' 1
.COO
.000
* " 0 "
. '• 'J n
. <• i/"
.00"
."C"
. " 0"
. ";:0
.lui
."0"
.ouo
."CO
."CO
.00"
. Ou"
."00
ItNTOAfjCN
AT L^u
(h'./L)
.001
.100
.000
.000
.000
."00
. 000
. "CO
."00
."01
. " t "
.000
. "CO
.00"
.100
.000
, "Q1
.00"
* 1 C "
.oco
."01
.101
. oco
.It1 1
."3"
.00"
.000
. 0 C "
."00
."0"
.00"
-------
SPQKANt RIVER
NUMiJEW OF HUN = 0
REGION NUUBFR 2 AUGUST, 1159
TARGE.T 0.0
SEASON o^ r". = AUG
NO.
OF
P.MCH
1
2
u
fc
f,
10
1?
1"
15
16
1 9
21
2?
21
25
2*
27
2"
31
33
35
37
3"
uo
U2
uj
«u
Uli
afl
50
b?
IDENTIFICATION
OF
"EACH
G u 11 0 TO G 1 1 1 5
G u 1 1 5 TO PSITCHiRO c
PRITCHAKO C TO BEAVH.
"tiyFR TH GUAM" CKK
G w J M H TU STL»''HG4T
S1F.AMBOAT TU S-f-C
S-^-C TO N KN«
N FORK TP GU13"
GU130 Tu S (-"K CDAL
"3d TO DF ADMiN
nt JD TH HH^TFH i s c
GOLD TO 5L'U'iTtP
Sl_TfnUS Tn LAMVUN C
"6.'L'
.COQ
.coc
.000
. 0"0
,C"0
.000
.000
.000
.0"0
.000
. 0 "0
, C o o
• u 0 c
,000
, U 0 (j
,000
.(/oc
2 . - ^2
2 ."^3
2.1-13
1 . 'j 5 7
1 . iOi.
1.1-17
1.K17
1 .3C3
. JOO
."10
.39(1
• .i^i
.3"C
. J»n
in\-i cnsC
AT ? 1 A rt T
t"L-/L)
.oon
.000
• oo "
,0i,0
. "00
. "00
.00"
.000
."or,
."00
. r. 0 "
. " 0 "
. "I'O
.CO"
."CO
.OQO
.OM"
. r u 0
. " u"
. "0"
. " 0 o
. o i; 0
,0j,1
.000
. "00
. OI..P
. ." u n
.OuO
. f I. 0
. "0 "
. "I,"
tv/T)ATlP\
..-CO
.000
.000
.000
. r c o
. "00
. o o o
.on"
, 0 ; 0
.^00
. " 0 0
."CO
.00"
.00"
.000
."00
."00
.00"
.000
."or
.000
.000
.000
.00"
. 0 L' 0
-------
VD
Nt
Sf
NO.
QF
"t'CH
1
2
u
fc
A
10
1?
1"
15
It,
19
21
22
23
25
2*
n
29
31
33
JS
37
3'
UP
«?
U3
uu
«h
u«
51
b?
SPOKA^t RIvEK
JM9E>< OF KUN = 0
^ASO*' o^ YR. = AUK
IDFMTIFICATION
OF
BEACH
GullO TO r. K.
CAVyOf^ c JO (-U131 S
("•"131.5 TO LAK^ rwK
LAKt c in SlLvFkTQ'J
HuTFJLL TO O^ni'hK' 0
fuTMLL TO BT^ CBLKK
BIK CWEtx TO MTLO C
"HU TO KLLinur. OUT
f. U T TO G"l33
r. u 1 J 3 TU plM: r«ELk
Pl"t C TO fu'L1-1 "IV
S FRK TO G"1J5
CU135 TO aJuLy C»K
ilJULY C To W.~iE THK
PCjSt C T" KILLAPMtY
KILL U Tn flLiC" L»JTsATION
AT S7ART AT t«.'0
(^G/L)
.OOP
.000
,000
. "1)0
.OOP
.000
."OP
.OOP
.00"
,001
. 1 1) f
.000
.10"
.OOP
, "0"
,• 1 1) 0
.000
. o i, r
,0(1"
.OOP
.11)0
."UP
.11)0
. OuO
.01)0
.000
.000
.nor
.000
."0"
.000
(MS/1 )
.00"
,000
.coo
.ceo
.PO"
.OOP
.000
.000
.f&n
."0"
.OOP
. "00
.000
.POO
.000
,"0"
.00"
."CO
.000
.10"
.000
.000
. "0"
."OP
.OOP
."OP
.000
.000
."OP
.000
.000
u-
A7 STA
(i-r./u
,000
.000
,C "0
,cno
.UPO
, 0 1 U
.000
.010
.010
,0"0
.010
,ono
, 0 "0
.linO
.000
.010
.0"0
("(> /L)
.COO
.occ
.010
,cno
, c"6
.000
.000
.C"0
.010
.0"0
. 0 1 1.
.cio
* U '" C
,000
, o"l,
.010
. 0"C
2 • -i?2
2 . ^ ^3
2 . ". i 3
1,'jS?
1 , b ib
i , *• i r
1 , "• 1 7
1.3^3
, J 10
.'•"'0
. 3°fi
• 3«i
. 3^0
. 3sn
JOs-1 C^
AT ? 1 £HT
t"0/L)
* i c "
."OP
,oon
. Oo"
,"0n
, 1 0 1
.'01
.10"
, i C "
,"00
. f' 0 "
."•-•"
. "o"
.CU"
,101
. r 0 "
. Pu"
1 f U 1
. 1 j"
."pi
. 1 0 1
. r k. 1
, 1 u .1
, 1 C "
. 1C. 1
.If"
. •" u 1
. 0 u "
,f I'O
."0"
• '" I/"
,rc"
"00
"C "
.0 0
, " c"
, ' C r
,10"
.100
.10"
.10"
, r 0 1
. " 0 1
.10"
."01
. Oi'"
. .1 C "
-------
*»»»*FIMAL SUM^AKY*****
SPOKANE HlVEK HF.GION NtjMbEK I AUC.UST t J •)&">
NUMnE1' OK RUN = 0 - T'xrcET 0.0. LF.VF.L = -.00
SEA-'U*) Of YP-, » AUT,
MO.
Hi-
HFAHwiTEM
1
£
3
-UO,
OP
STPLTCH
1
2
3
INITIAL
K LOW
(CTS)
117,0
1,0
13.0
FILIAL 1
K LOW
CCFS)
117,0
1 .0
13.0
MJG'1EKJT'TI01>
REuUIftEO
( Q F ^ 1
.0
.0
.0
-------
»T KEAtH tNO
RE4CH
5?
50
ufl
uf,
uu
ul
10
8
&
1
HZ
«0
39
37
35
33
29
27
26
DEP. VA«.
.83071+01
.8=608+01
.87233+01
.8h77u+01
.8UU02+0 1
.BSl.33+01
,6f 108+01
.88773+01
.8«86b+01
.92301*01
8.300
VO
ON
-------
On
KEACH tND C*G/L1
PE«CH
53
50
it*
10
8
6
1
nz
10
39
37
35
33
SO
27
26
2?
OEP. VA".
.83071*01
.e«Zfl2»or
.85618*01
.87233*01
.85^5*01
.8SUS6+01
.870U6+01
.83936+01
.901 ?9+01
.88773+01
.92301+01
8.300
-------
»T
PEACH
33
27
22
DEC. VA».
.911'10-00
.91812-00
."??t!OB-00
.7JO
52
50
1»
16
11
11-
15
11
10
p
h
1
12
Ufl
3'
37
.130'2*01
.13078+01
.12616+01
.12313+01
. 1 1 7-00
. 7<-Oq^-00
.73911^-00
.22f-l + 01
.2I7"b+ol
.2U219+01
. In717+01
0
0
0
0
0
0
ft
0
0»
t
*
*
0
0
0
n
*
«
ENfl (MS/L)
2.ceo
o
o
o
o
0
0
0
0
0
0
0
0
0
» 0
* 0
o
o
«o
0
o
o
o
-------
NH3 AT
H ENO ChG/L)
REICH
sn
US
ilh
aa
«3
15
1"
11
«
6
1
0?
UO
31
37
DFP. VA".
.JOJS9-0?
.3ii"i-oa
.i4?-ia-o2
.371473-0?
,000
o
0
0
0
o
0
0
0
0
o
o
0
0
n
o
MD
OO
-------
N02 AT
END (MR/LJ
.016
RtiCH
5?
50
IS
1"
1"
fl
(,
1
«?
«o
3"»
37
3S
33
2"
il
ah
22
DEP. VAR,
.20317-0?
.tl" 117-0?
. t>fl<477-u2
.000
0
0
0
0
o
0
n
o
o
o
n
o
n
o
o
o
o
o
o
o
o
o
-------
NOJ AT BE4CH END (MG/L)
RE'CH
b?
bn
03
1"
10
R
1
U5
un
IS
UFP. ViR.
b f-O - 0 1
i ?7u-fll
£ o u n 3 - o
.020
0
0
o
0
n
o
n «
o *
o
o
o
o <
t> <
o *
o »
o <
o *
o »
o *
0 »
.030
o
o
-------
N03 *T
END (MP/U
b?
10
R
h
1
u2
uo
3«
3?
<;=>
UFP.
. e! 1 V 1 u - 0 1
,2"M fr-01
.i^M-ul
n 1-01
.020
0
0
0
0
fl
0
fl «
n t
o
0
n
o <
o '
o >
o *
P •
o «
n <
0 *
0 »
0 «
0 «
.030
o
o
-------
AT RtiCM END t^G/1.)
DFP. V4P.
5?
50
IS
1"
10
s
*>
1
3"
3?
,U«6(JO-01
. J 7 i u r, - 0 1
.1B7PU-01
.030
o
o
0
0
0
0
0
0
0
0
0
0
n
0
0
0
0
0
0
0
0
0
.060
-------
COL AT RtiC"
Pt*CH
ua
U3
15
1"
10
*
IS
1
n?
un
33
d°
27
2*
2?
.ill!0»03
. 1 o a '> i » 0 3
IP3.000
0
0
0
n
0
0
0
n
*
*
n»
o
o
0
0
o
o
0
0
0
0
0
31?.000
0
0
0
o
o
o
o
o
o
o
o
o
o
* o
» 0
0
0
»0
(I
0
0
0
-------
COL AT Rt»CH FNP (M(Vl)
Rt*C» DEC. ViR. 1P3.000 31?.000
b2.23701*030 * 0
b<1 .2? 1 (V*03 0 * 0
ue,2^<474*u30 * 0
-------
AT REACH RNt> CMG/U
RE4C"
OFP.
.000
V?
bO
lift
U6
iiu
«3
15
I"
10
A
h
1
«2
uo
3'
37
35
j3
i|9
27
26
. 1^278*01
.1605/+01
.16866*01
•17604+01
»1779ti*Ul
,!Ptj6B*01
. 0"000
.00000
.00000
f Q 0 0 H 0
.oncno
.00000
.6Tbu7»ol
.61707+01
• 6S367*0 t
. 7f"3nu-cO
.6" 01 f)-00
, 1 P !y tlt>+ 0 1
. 1 1 1 (•• 0 + u 1
. 1 1«7B+01
. 1 1 bOi+0 1
0
0
0
0
0
0
(1
0
0
0
n
0
0
0
n
.00000
6.600
0
0
0
0
0
0
0
0
0
0
0
0
» 0
» 0
*0
0
0
0
(I
0
0
0
o
OJ
-------
CL- AT RtAtH END
o
JS
50
ifl
15
1"
10
S
6
I
3°
37
33
.31301-00
,OOuOO
.00000
.00000
. 0 0 0 n 0
.11173+01
. l^OSQ + Ol
,00000
.00000
.00000
0
0
0
0
0
0
0
0
0
0
0
0
0
u
o
o
o
» 0
•0
0
0
0
-------
SECTION IX
PROGRAM LISTING
Program Page
Typical JCL 106
Subroutine BLEND 10?
Subroutine CMIN 107
Subroutine DOEQU 107
Main Program 109
Subroutine DOSAT 116
Subroutine GETCON 116
Subroutine K2CAL 123
Subroutine NEWIN 123
Subroutine PLOT 127
Subroutine PLTSET 128
Subroutine REINI 130
Subroutine REPRE 131
Subroutine RMATC 131
Subroutine RUNON 133
Subroutine SCAN 137
Subroutine SETOPT 138
Subroutine TRIED 139
105
-------
TYPICAL JCL
A typical JCL deck for executing a DOSCI run is listed here. (JCL
variations may occur according to machine or installation).
// EXEC FORTGCL,REGION=160K
//FORT.SYSLIN DD UNIT=SYSDA,SPACE=(CYL,(5,5)),
//DSN=&LOADSET,DISP=(MOD,PASS),DCB=(RECFK=FB,LRECL=80,BLKSIZE=3120)
//FORT.SYSIN DD *
DOSCI FORTRAN SOURCE DECK INSERTED HERE
/*
//LKED.SYSLMOD DD DSN=CN1752.ZAS.DOSAG(DOSCI),
// UNIT=2314,VOL=SER=WYLB02,SPACE=(1024,(140,10,1),RLSE),DISP=(,KEEP)
// EXEC PGM=DOSCI,REGION=160K
//STEPLIB DD DSN=CN1752.ZAS.DOSAG,UNIT=2314,
//VOL=SER=WYLB02,DISP=SHR
//FT06F001 DD SYSOUT=A
//FT05F001 DD *
DOSCI DATA DECK INSERTED HERE
/*
//
106
-------
t, SUBROUTINE HLEND
2, COMMON/ULOCKS/JJ.KK.II
3. COMMnN/HLOCKVCONDZ(10»20)iCONOI(20i20),CONDE(20i20)
4, COMMON/MISC/IO,CALCK2,NrLOW,NCON
5. C
6, C STEP-1
7, C ADO TWO UPSTREAM FLOWS AT A
8, C JUNCTION TO GET DOWNSTREAM
9, C FLOW,
10, C
11, CONDJ(KKtNFLOX)=CONDE(JJ,NFLOW)iCONDEUl(NFLOW)
12, C
13. C STEP-2
i«, c CALCULATE DOWNSTREAM BOD AND
is, c 0,0, BY USING A WEIGHTED MASS
16. C BALANCE FORMULA,
17. C
18, DO 3 I=1,NCON
19, CONDI(KK,I)=((CONDECJJ»NFLOK}*CONDE(JJ,n)+(CONDE(II,NFLOW)*CONOE
20, *(II.l)))/cONDl(KK.NfLO>N)
21, 3 CONTINUE
22, C
23, C STEP-3
2«. C RETURN TO CALLER (R M A T C)
25. C
26, RETURN
27, END
26, SUBROUTINE CMIM
29, COMMON/BLOCK2/JNITC20) ,F( 10) iC(ll)
30. COMMON/BLOCK5/JJ(43,RMLOW,MAX(CLOW
31. C
32. C
33. C FIND SUBREACH WITH LOWEST
ju. C DISSOLVED OXYGEN LEVEL AND
35, C RIVER KILE AT WHICH IT OCCURS.
36. C
38, CLOw=c(l)
39, IF(MAX.EG.l) RFTuHN
40. DO 1 1=2,MAX
41, !F(CU) .GE.CLOW) Go TO 1
42. CLOKsCU)
43, RMLtln = F(I)
fl«, 1 CONTINUE.
45, C
IC«) ,RMLOW,MAX,Cl 0*,NTRIH(10) ,QUF',FIN'L(5) t IA,
55, *JCK,f- INL'-t3) ,NJ,K?,VEl.
56, REAL K?
57. COMMO'i/HLOCKl/CKMIN(100),DATA(100iOO),FINIS(lOO,50),RIDENT
50, *C100,'j) , KHOPTC 100)
59, COM^•,0^/t^LnCK2/J^lT (c"0) , F'Cl 0) iC( 1 1 )
bO, COMHON/KCHVAR/CCLK(faOO),DOK^(100)
107
-------
61, COMMON/TF_"iPFR/TtMPAV,1(>'KHA ( 100)
62,
63,
64, COMMON/COSiE'v.n/nOEfPnDF. fNH3e»N02t,>,n3fiP04EiALGEfCOLEiHMIE»HM2F.«
65, *HM3EiHne'i TUTSt
66, COMMON/OPT I "N /IF N, IK?, ICOLiICOMB
67. COfMON/OPT2/IMf AVY.ITDTf.1
68, RtAL NH.I. r,iOr>t .vO.it NH SEt NU2E»N03E
69, DIMENSION ARIM(l) f AKOUT(l)
70, EQUIVALENCE (ARINtDO) i (AKOUTtOOF)
71, DATA BLM
-------
l?2. 16 ARINCIL)
123, 20 CONTINUE
12«, CAUL CflM
125, IFdCK.LO.O) *RIU(NJ.100) IA,OATA(IA,2),TOTFLO,CLOh,RMLOWi -~
126, *XCOM ( 1 ) ,i',H. ,=Ai-"X , C.H' OA , H"K A , MPO, , HI- -jnn
179, *t AHM?i All'
-------
185, COMMON/TEMPtR/TE>-'t'AV.TEMHEA( 100) iSATRFAf 100)
1BU. COMMON/ OF'TION/JFNtIK;>,ICOLiICrj>-lHiINH3iIK'02iISl03»IPOf|ilAl.G«IFIKST
185. COMMON /OF'T2/IhEAVY« IKlTNt 1CHLOH
186. COMMON/OPT3/IP.INH,IN2«IN3
167, COMMON/r>IUC/IO.CALCK2iNFLUWiNCONiTOTFLO
168, REAL NOHEFR
189, REAL NH3«flO?|N03.MPO/4«MlN03
190, REAL K2Nni.HNH3i*L«NR,NH3KiN02K.IO
191, REAL NH3F. «N02E« N03E
192, DIMENSION AMOMTH(IZ)
193, RFAL K2
190, DATA AMONTH/3HOCTt3HN'OV,3HDEC.?HJAN,3HFEBt3HMAR,3HAPR,3HMAY»
1'iYPF RIvtR ,
2ijO, * iiOhllASIN t/)
?q], M?, (T1TLKI) i I = l»l«)
2/13, 1 FORMAT (20A«)
110
-------
215,
216.
2«7,
2UP.
2«9.
250,
251,
252.
253,
P5«,
255,
256,
257.
258,
259,
260.
26J,
262.
263,
264.
265,
266,
267.
268,
269,
270,
271.
272.
273.
270.
275.
276,
277,
278.
279,
280.
231,
262.
283,
28«,
285,
286,
287.
288.
289,
290,
291.
292.
293,
29
READ CNI,ij oo"iIOUM^
IF (OUf, l , Nh. .(..'•IDF ) GO TO 777
902 FORMAT OSXi'llH* * * FILE B • PHYSICAL DESCRIPTION * * *,//)
'IRITE (NJ.602)
602 FORMAT U'JXIIOHCAND NO, OF NO. OF no, t
40HOF NO, OF INSERT i MEAN ,/,
ISXtHOMTYPE HEADWATERS JUNCTIONS UtACt
«OHHE.S STHKTCHE3 FOR FINAL ELEV.,/,
1'JXi'IOH MAX OF 10 MAX HF 10 MAX fl,
ilOMK 100 MAX OF 20 3UM&ARY C^LY (FT)/)
Rf AD (NI . 2M DUMl ,l>UM2,t'INiT iNJUNC «'.REAiNTF.0)
WRITE (NJtOOM) OUMl,DUM^,NIMT«NJUNCiNRtA,MRIBiICK»KF.V
READ (NI , i ) OUdl iDUM2
"RITE (NJtflOl) OUM1.DUM2
IF (OUMl .Ht.ENOF) GO TO 777
WRITE (NJ»902*Tt(l-"0(I"w)»J-l#20)
805 FORMAT ( 1 5X , 2Aa , 2Xt 15 , 5X , 20(1Xt 12))
10 CONTINUE
READ (NItl) UUM1,DUM2
WRITE (NJtfiOU OlJMltDUM2
IF (OUMl.NE.tNDF) GO TO 777
IF(NJUNC.KQ.O) GO TO 97
WRITE (NJ.905)
905 FORMAT («OX,30H* * * FILE D » JUNCTIONS * * *i//)
WRITE (NJtbO'l)
60U FORMAT (15X, 40HCARD NO, NO, OF ,
* UOH NO, OF NO, OF •/•
* ISXtlOHTYPE OF UPSTREAM ,
* 40HUPSTREAM DOWNSTREAM ,/,
* 15X,«OH JUNCTION STRETCH ,
* «OH STRETCH STRETCH t/)
00 11 KcJ.NJU'JC
READ (M,33) DUMl ,DIIM2,I, (JUNC(I t J) i J=l»3)
33 FORMAT (2A«, 13X i I2,2X(3(UXtI2iUX))
WRITE (^J,806) DUMl,DUM2.I,(JUNCtItJ)tJ=lt3)
806 FORMAT (15X,2A«, 12XtIbi3(bXi15))
11 CONTINUE
97 CONTINUE
READ (NItl) DUM1,DUM2
KKITt (NJiHOl) DUM1,DUM2
IF (OUMl .NE.ENOF) GO TO 777
WRITE (NJ.906)
906 FORMAT MUXiilH* * » FILE E « HEADWATERS * *
WRITE (NJ»f>05)
60b FORMAT(I CARD STCH AUGM 00 BOD NH3-N
N02-SJ N03-N POU-P
111
-------
305,
306,
307,
300,
309,
310.
311.
•312,
313,
3H.
315.
316,
31V.
318,
319,
320,
321,
322,
323,
32«.
325,
326,
327.
32B,
329,
330,
331,
332,
333.
334.
335.
226,
337,
330,
339,
3«0,
301.
302,
313,
3«0,
305.
346.
307,
308.
309,
350,
351.
352,
353,
35«.
355,
356.
357.
358,
359,
360.
361,
362,
363.
36U.
365.
S TOT N
(MG/L) C'T./U
(I-.G/L) ("C/L)
* PHYTO COIUORMS HMI HM2 H
» KMIJI/ I TYPE NO OPT (Mf./L)
»G/L) (MG/L) (MPN/100) (MG/L) (KC./L)
*G/L) (MCi/L)'/)
DO 12 K=l,MI'MlT
Rf.AL'fN'I »??) DUM1 ,OUM2i I t I AUC, ( I) I (CONOZfl t J) i J=l i 16)
IN I T t K ) = I
HTI
CT./t) (> = ,15*1,7
IF(OATA(Iift),t3.0.) DATA(Ii6) = .«H
FORMAT <2««i8Xil2,2Xi6(2XiFB,0))
"RITE (NJ.S08) DUH1iDUM2iIt(DATA(I IJ)tJ=l16)
FORr*AT(rJX»2AO,7X»17,F8,l,F10.1,FlO,3,FlO,6,F10.3,F10.3)
CONTINUE
00 HI K=l t100.
BOOK(K)=DATA(K,3)
BODKS(K)=DATA(KrO)
READ (Mil) DUM1.DUM2
*RITE (NJ.BOl) DUM1.DUM2
IF (oU'M .NE.ENDF) GO TO 777
hRITE (NJ.9083
FORMAT(31Xi i* * * FILE F(2) - DATAdi?) THRU DATA (I i 1 0 ) , DATA (I i 21)
*THRU DATA(I;32) * * *'/)
WRITE (NJ.607)
FORMAT(I CARD REACH FLOW 00 BOO NH3-N M02-N N03"N POO"P
* PHYTO COLIFORMS HMI HM2 HH3 TOTN C^LOU H ^ 11 H M I <;
* HMI3I/ 1 TYPE NO (CFS) (MG/L) (MG/L) (MG/L) ClG/L) (MG/L) (H
*G/L) (MG/L) (MPN/100) (MG/L) (MG/L) (MG/L) (MG/L) (MG/L) (MG/L) CM
*G/L) (MG/L)'/)
DO 150 Ks1iNREA
REAO(NI,127) OUmiOUM2iI|STOR
FORMAT(2AU,3X,I9,F10,0)
INPUT ONLY RtAChCS ftlTH NONZERO FLOWS
IFCDUMl.EU.tMOF) GO TO 151
DATA(I,7)=STOR
RKADfNI.128) (DATA(I,J),J=e,10),(DATA(I,J).J=20,32)
FORMAT(8F]O.U/BF10.U)
KRITE(NJ,816) DuM2,I,(DATA(I,J),Jc7,H)1(DATA(I,J),J=?«,32)
112
-------
366,
367.
360.
369.
370.
371.
372.
373,
374,
375.
376.
377,
37fl,
379,
300,
301.
382,
3B3,
381.
305,
38b.
387,
38B.
389.
390.
391.
392,
393,
394,
395.
396,
397.
396.
399,
aoo.
401.
(102.
103,
404,
405.
406.
407.
408.
409,
410,
411.
412.
413,
414,
415.
416.
417.
41B,
419.
420,
421,
422.
423,
420.
425,
426,
010
150
151
909
60fl
170
BIO
160
153
161
910
609
C
C
240
811
230
911
FORMAT(lX,A/j,I5iF6,liF8.2.F7.2t5F-7,4iMO,lfflF7.2)
CONTINUE
READ (Mtl) t>UMl»DUM2
CONTINUE:
WRITS: («j»«oi) OUMI.OUHR
IF (IHJMl .t.'E.I'NDK) GO 10 777
WRITE (NJflO'J)
FORMAT C30X,«hH* * * FILK F(3) « OATACI.15) THRU DATA (1.20) i
* 5Ht * *i//)
WRITE (Nj,iiOB)
FOfU'AT ( lr:>X.40HCARD NO. VALUE COEF, 0»
* 41HN EXP, ON EXP, ON COEF, ON EXP, ONi/t
* 1SX»«OHTYPE OF FOR K2 V FOH Kt
» 41H2 V FOK K2 D FOH K2 Q FOR K? 0 FOR K2t/»
* ISX.^OH RtACH OPTION-1 OF'TION^i
* U1H2 OPTIQ'-HS OPTIUNf^ OPTION. 3 OPTIONnJ./)
DO 160 K=l.HRfcA
READ (MI. 170) DUM1.01IM2, I , (DATACliJ) » J=15»20)
IFCDUH1 .EO.FJNOF) GO TO 153
FORMAT (2A4(8X.I?i2Xi6(HX»FH.O))
FORK AT (15X,2A4,7X|I5,6(2X,K6,3))
CONTINUE
READ (NI.l) DUM1.DUM2
CONTINUE
D0161K=liNREA
IF(DATACK.l) .EO.O.) GO TO 161
IFCOATACK, 16). [-0.0.3 DATA(Kil6) = 3,3
IKLiATA(K.17).En.O.) DATA(K.17) = 1,0
IF(OATA(K,1B).EQ.O.) OATA{Kilfl)=l.i3
WRITE (NJ«UIO) LA01iLAB2iK.(DATA(K.J)iJ=15t20)
CONTINUE
WRITE C^JtSOl) DUHliDUM2
IF (Ouni .'
-------
027,
020.
129,
030,
OM,
t|32,
133.
lit,
035.
036.
'137.
138.
039,
000,
I'll.
002,
003,
000,
005,
006.
007,
000,
009,
050,
051.
052.
053,
050.
OSS,
056,
057,
05ts,
059.
060,
061.
062,
063.
060.
065.
066.
067,
068.
069.
070.
071.
072.
073,
070,
075,
076,
077.
078,
079,
060,
081.
082.
063.
080,
085,
086.
087,
610 FORMAT U5X.OOHCAPO
*i FRACTION KEMOVfO HY TREATMENT I/»
* 1'JXtOOHiYHE D.O, LEVFL (>'G/L)
*/)
ReAD(NI.S) (HIM ,DUM2» (OOL(I) » 1 = 1 iUj I CTIU AC( I) , 1=1 »0)
5 FORMAT (2A CONDZ(K,15) =
*PIHM2*CONDZ(K,10)
IF(CONOZ(K, 1 1) .EQ.O. .AND, CONDZ(K,16) . EQ . 0 .) CONDZ(K,16) =
»PIHM3*CO'JDZ(K, 1 1 )
250 CONTINUE
114
-------
480,
4fl9,
490,
491,
492.
494,
495,
496,
497,
498,
499,
500,
501,
502,
503.
504,
505,
506,
507,
508,
509,
510,
511.
512.
513.
514,
515,
516,
517.
5l8r
519,
520,
521,
522,
523.
524.
525.
526.
527,
528,
529,
530,
531,
532.
533.
534.
535.
536,
537.
538,
539,
540.
541.
542.
543.
544,
545,
546,
547.
548,
C
CHtCK TO INSURE THAT THE BOD
P 0 D « T = fl 0 1) C + 1 2 . / R 0 0 P C
fi 0 D M' K = H D () N 4 1 « , / IJ 0 D * T
B 0 0 P 'A R = } ? , / » 0 D W T
DO 256 K = 1 i S. f) E A
•N EXCEEDS THE ALGAt."N
AlGN=nATA(Ki 1 4) MiODNhK/ ( APH •* BODPUR )
779
256
C
C
C
C
C
C
20
22
23
24
3335
334
If- (BODN'X.GE.ALGN) GO TO 256
hRITe (NJ • 77V) BODNXi ALKMK
(•ORMATP *<:*FIOQ NITROGEN =i,
* " If-1 R E A C H I , I 5 , 1 * * * 1 )
CONTINUE
NRUNcQ
DO 13 I=l»4
IF(DOLCI) .EQ.O, .AND, I.GT
DOLEV-DOL(I)
DO Id J=l,4
TF = TRFAC(J)/MOO.
IF(J,GT.1 .AND, TF.EQ.O.)
00 15 K=l , 12
IF (TEMMO(K))2?il5i22
CONTINUE
VAR=0,
DO 23 LL=lfNREA
IF (TEMREACLL) .EO.O.) TEMREA
V A R = V A R + T E M H t A 1 1 1_ l_ )
TEMPAV=VAR/NREA""
NSEAS=K
DO 24 L=l( 10
CONDZ(L(NFLOW)=HWFLOW(L(K)
CONTINUE
IF(ICK,EQ.l .OR, ICOHB.EQ,
WRITE (NJ(2055)
kvRlTE (Nj,334)
F10,6(i AND ALGAL NITROGEN =>(F10,6i
STtP-2
SET BASIN CONDITIONS TO RE
CONSIDERED (TARGET r> , 0 . LEVELt
TREATMENT EFK.f AND TEMPERATURE),
,1) GO TO 13
GO TO 14
aUsTEHMOCK)
18) GO TO 3336
FORMAT (26X»49H* ****INTERMEDIATE SUMMARY*
WRITE (Nj.eoil) (TITIECL) iL=l»18)
8011
335
C
C
C
C
3336
C
C
C
C
FORMAT (19X, 18A4f //)
NRUN=NRUNtl
WRITECNji 335) NRUNiDOLEVi TF(
FORMAT (15X, 15HNUMBEK OF RUN
F5 1 2 1 / «
15X( 15HTREATMENT (C)
F5.2(/.
15X. 17HSEASON OF YR.
F5.2.//3
CALL RUNON
TFN(SFASON(KSEAS)tTEMPAV
=. I5.36X, 19HTARGET 0,0, LEVEL =»
=»F5.2i36X( 19HTREATMENT (N) =i
= |A3I36X( 19H"EAN TEMPERATURE =i
STtP«3
CALL RUNON
STtP-4
REPEAT STEP-2 THRU STEP-I
UNTIL ALL coMHiNATinNs OF
115
-------
509. '
550.
551.
552 ,
Sb3 ,
550,
555,
556,
RS7
•
559^
560.
561.
562,
563.
£,f*U
J O H .
C i C
J C «
581,
582.
583,
58U.
585.
586.
587.
588.
589,
590.
591,
592.
593.
591,
595.
596,
597.
59fl.
599.
600,
601.
602.
603.
604,
605,
606,
607,
608,
609,
C
C
C
C
C
C
C
C
C
1
C
C
C
C
C
C
C
C
C
C
C
C
C
C
UASIN CONDITIONS HAVE
ANALYZED,
15 CONTINUE
11 CONTINUE
13 CONTINUE
GO TO 7777
777 wHITf. (NJ,778)
* 22Xi3«H+ EXECUTION WAS T E t
33HR HINATED BECAUSE *i//»
22Xi !H+»31Xr3HO Ft31XilH*i//i
2?Xi34H* ERRORS IN ,
33" I N P U T DATA *f//i
GO TO 9999
7777 CONTINUE
9999 CONTINUE
STOP
END
SUBROUTINE DOSAT
COMMON/BLOCK5/JJ((3)»NRt:AiNJNITfNJUNC(ELEV
COMMON/TF.MPER/TEMpAViTEMREA(100),SATRr:A(lOO)
STEP»1
UFEN
CALCULATE THE 0,0. SATURATION
LEVEL FROM TEMPERATURE
BASIN ELEVATION.
I'U 1 I - 1 f "H f. rt
TEHP = TF HREACI)
CSAT=(1«.62»(,3898*TF.MPJ+(,006969»TEMP*+2J.(,00005897*TE
* (1.0-C0.00000697*ELEV))**5,167)
SATRtA(I)=CSAT
STEP-2
RETURN TO CALLER CR U
RETURN
END
SUBROUTINE GETCON(IO»DELT»lNDtDEPTH(VEL)
PHYTOPLAMKTON IS THE ONLY TYPE OF ALGAE MODELED BY THIS
10 IS THE AVERAGE LIGHT ITENSITY IN LANGLEYS/MlN DURING
STEP
DfcLT IS THF TIME STEP LFNGTH IN HOURS
IND IS THE INDEX OF THE STREAM REACH
DEPTH IS THE AVERAGE DEPTH OF THE REACH IN FEET
VEL IS THE AVERAGE VELOCITY OF THE RFACH IN FT/SEC
COMMON/MISC/HLAMK,CALCK2
COMMOI./CO'-.'DtG/DOtHOO»fjH3,N02i''l03iPOU»ALr. iCOLiHMl,HM2,MH3
AND
MP**3) }*(
,
N ON),
ROUTINE
THE TIhE
,HM,TOTN
COMMON/CO '; FsD/POEiBODE»NH3 F. •N02tiK.'03F,PO«EtALGF. »COLt»HMlEtHH2t»
* H M 3 f ,HMf»TOTN'E
RFAL fJH3.»JDf'fM03(NH3EtNC2EtM03F
COHf'ON/Cn^ST/TwKCnLiAHO^tAwMtCHMOC'THKNH^iVOLlT^iTHVnL^t
tpQQKiRn^^CiHOnof7iNOREpRir>HMAXiTHriHHXiCn.L'OAif-iMKA|MPO/JtMlN
*MNHJifL»APR,NR,ASH,ANDtATDibRKBODiBRRPf)ili8RRNH3»HtNOO
SOOC i
03ii2N03i
116
-------
610,
611,
612.
613.
614,
615,
616.
617,
618,
619,
620,
621,
622.
623,
624.
625,
626.
627,
628,
629,
630,
631,
632.
633,
634,
635.
636,
637.
638,
639.
640.
£" ' •
642.
643.
64U,
645,
646,
647.
648,
649,
650,
651.
652,
653.
654,
655.
656,
657,
658.
659,
660.
661,
662.
663.
664,
665.
666.
667,
668.
669,
670.
C
C
C
C
C
C
C
C
C
C
C
C
C
c
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
» , Hl'KA? (HMKA^i TH'.'QSK t Tuf'fVJK
CPXhON/KCHVAR/CQI. KC100)iHr;OK(100),BnpKS(IOO)iK'H3K(100),N02K(100).
* iPO«K( ion) ,f.n3i< cioo)
COMMON/ TF. »'Pe R/T F *PAVf TF.WREA( 100) i SATPtA (100)
REAL HH.tK.NOZKiHPOil, MllJOi, Mr>NtM|MNHj, I 0 , '.'H , I H AR , ML , K 2 A20
REAL NOiK
REAL NOPFFR
COM'.'ON/nPnnM/lFNiIK?iTCOL»lCUMHflNh3iIN02iINU3tIP04iIALGtniRST
*OHODUR i ONu 5f;0 t O^H^ i n.SIHiV t ONH i,",(i i O'.'H^AP i Dklt':U'R i [)NO?M3 i HNC2 t O'-'OiBOt
* D E L A t n A R C S t f ) A L 5 N K , 0 A 0 T H , D E L 0 i 0 M E t D » 0 G I V E
DIMENSION nf-LTS(l)
EQUIVALENCE (DEI.TS.DELCOD
DATA ULAN, ASTER/i I , I* 1 /
DtLCOL
DLBODO
DLBGDS
OBODAG
DBODAR
DHODAS
DROOAO
DBOOOR
DNH3BD
DNH3
DNH3V
P ^ H 3 A G
0 N H 3 A R
0 N H 3 R R
DNO?H3
DN02
ON03BD
n N 0 3 H 3
DN0302
DN03AG
DNOiAR
ONLY 30
FAKIT=1
DN03HR
PP04UD
DP04AG
DP04BR
DELA
DARES
DALSN'K
DADTH
DELO
ONEED
OGIVE
IS
IS
IS
IS
IS
IS
IS
IS
IS
IS
IS
i 0
IS
IS
IS
IS
IS
IS
IS
IS
IS
CHANGE
CHANGf.
CHANGE
CHANGE
CHANGE
CHANGE
CHANGE
CHANGE
CHANGE
CHANGE
CHANGE
C r ' M •'•< u F
CHA'-'GE
CHANGE
CHANGE
CHANGE
CHANGE
C H A N G E
CHANGF
CHANGE
CHANGE
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
IN
CONSECUTIVE
,
IS
IS
IS
IS
IS
IS
IS
IS
IS
IS
IS
CHANGE
CHANGE
CHANGE
CHANGE
CHANGE
CHANGE
CHANGE
CHANGE
CHANCE
CHANGE
CHANGE
IFUFIRST. NT, 1)
NOUT=6
JBUGiO
U1RST =
0
IN
IN
IN
IN
IN
IN
IN
IX'
IN
IN
IN
GO
COLIFORM
BOD
BOO
HOD
BOO
BOO
BOD
000
NH3-N
NH3-N
N M J - N
rx n } ~ ni
NH3-N
N H 3 1 N
N02"N
N02-N
N03-N
N03-N
N03-N
NO S-M
N03-N
COMMENTS ARE
N03-N
P04-P
P04-P
P04-P
PHYTOPHNKTPN
PHYTO PLANKTON
PHYTOPLANKTON
PHYTOPLANKTON
DO
DO
DO
TO 8
COMC.
CONC .
CONC.
CONC ,
CONC,
CONC,
CONC,
CONC,
cove.
CONC,
CONC.
C 0 N' C ,
CONc,
CONC,
CONC,
CONC.
CONC.
CONC.
CONC.
CONC.
CONC,
ALLOHFD
CONC.
CONC,
CONC,
CONC,
co»-c.
CONC.
CONC,
CD-JC,
CONC.
CONC.
DUE
DUE
DUE
DUE
DUE
DUE
DUE
DUE
DUE
DUF
HUE
out
DUE
DUE
DUE
DUE
cut
OUF
DUE
DUF.
DUE
DUE
DUF.
DUF
DUE
DUE
DUF
nuF.
DUE
DUE
DUE
TO
TO
TO
TO
TO
TO
TO
TO
TO
TU
TO
TU
TO
TO
TO
TO
TO
TO
TO
TO
TO
TO
TO
TO
TO
TO
TO
TU
TO
TO
TO
DECAY(GROWTH) (-OR+)
DECAY(-)
SETTLINGS)
ALGAL GROWTH (+)
ALGAL RESPIRATIONS
ALGAL SETTLINGS)
ALGAL DCATHC+)
BENTHAI. RELEASED)
BOO DECAYC+)
NH3 DFCAY(i)
NHi VOLITIZATIQN(-)
ALGAL GRO'^THC-) .
ALGAL RESPT.RATIONC +
HENTHAL RELEASEC+)
NH3 DFCAY(t)
N02 OFCAY(-)
BOD DFCAYU)
NH3 DECAYt*)
N02 DECiY(+)
ALGAL GROWTH:-)
ALGAL RESPIKATION(t
BFNTHAL RFLEASEC+)
BOP DF.CAY(*)
ALGAL GRO«TH(»)
GROwTH(+)
RFSHIRATIOM(^)
SINKING(-)
OKATH(K'ATIIRAL*TOXIC
ALL RhACTIONS(+OR-J
ALL OX-OF -^ANOS(-)
ALL OX-GF.NFRATION( +
POD»T=BOPC*12,/BODPC
117
-------
671,
672.
673,
67 1 ,
675,
676.
6/7.
678,
679,
600,
601.
682,
683.
60U,
6fl5,
686.
607.
608.
689,
690.
691,
692.
693,
694,
695,
696,
697.
698.
699.
700.
701,
702,
703,
70ft,
705,
706,
707,
708,
709.
710.
711.
712,
713.
71ft.
715.
716.
717.
718,
71"?,
720.
721,
722,
723.
72ft.
725,
726,
727,
728.
729.
730.
731.
8
9
C
C
C
10
c
c
c
c
c
c
15
20
C
C
C
MOf)'J*R = i)ODN+ 1/J./HOOUT
MOD'1 •*!.• = v./nnoM
FACl = ilOnuO + NUH|FR/(APR*HODPMO
CONTINUE
If- ( 1IIUG.R], 1) i-'RITFCNOUT.10?3) I MO f I 0 i Oh LT f OF f'TH . VEL i COL i BOD i NH3 i
* Nl 0 ? t 'J 0 i i P 0 4 t A L G i H U i H M i , H M 2 i H M i i T 0 T N
Of) 9 1=1.35
Df-LTSf. I) = 0,
0 N 0 3 0 S = 0 •
OPO.'K)S = 0.
TCOH = TFi'.Rf-.A( I,MO)-?0.
FAC2 = .i.28*Dl.l.T/(ntPTH*1000.)
IFCICOL.t'Q.O) GO TO 10
CALCULATE COLIFORM CONCENTRATION CHANGE DUE TO DEC AY ( GROWTH)
XTtf'P = COLK(IN01*THKCOL**TCOR
FACHM1=H'<1-CH^OC
FAChy2=HM?.CHMQ2C
FACHHJrHM^nCH^OSC
DELK = AKOQ*B(10
IFCFACHM1 .GT.O.) OELK = OELK-f AHM *FACHMl
IKFACHMP.GT.O.) DCLK = oELK + AHM?*FACHM2
IFCFACWM3.GT.O.) OEt.K = DELK + AHM3*FACHM3
IF(ABS(DKLK),r,T.AR.S(xTEMp)).l)fcLK = SH,N(xT£MP,DELK)
XK=XTEMP+DELK
DtLCOL=COL*(EXP(-XK*DELT)"l.)
IF(IHuG.EQ.l) xKITF(NOuT,1008) XTEMP,FACHhl,FACHH2tFACHM3iOELK»
*XK,De-LCOL
CONTINUE
IFCICOMB.CT.U) GO TO 20
CALCULATE BnD CHANGE DUE TO DECAY AND SETTLING
XTFMP=BODK(IND)*THKCOL**TCOR
DLBODD=BOn*CEXP(-XTE^P*DELT)-l.)
DHOTnT=BnD*(KxPft-XTEMP,BODKS(IND))*OELT).l.)
DLBODS=DBOTOT-OLBODD
jF(jBUf',F0.1) *-'RITE(NOUT(1009) XTFMP i DLHOOD . PlLBoDS
ONEtDsONEED+DLBODO
CALCULATE CHANGES IN POft»NH3» AND N03 DUE TO BOD DECAY
BODMTL=A8SCDLBODD)/BOOOQ
BOD"C=BODMTL*NOREFR
IFdPOft.En.l) DPOaBD = BODMC*BODPi^R
IF (INM3,EO,0) GO TO 15
DNH3BO=BODHC*hODNKR
GO TO 20
CONTINUE
IF(IN03.CO,0) GO TO 20
DN03BO=BODMC*BODN*R
ONEtO = Of.'Et-D»DN03BD*1.33
CONTlNt't
IKICOWB.LE.il .AND, IBUG.t.Q.l) URITE (K'OUT f 1 0 1 0 ) DNHJHQ , DNO JBD i
*DP01BnfO'JKEO
NITRIFICATION DOES NOT OCCUR IF oo is LESS THAN 2 MG/L
IFCOO.LT.2.) GO TO UQ
lF(lNH3.Ef;.0) GO TO 30
NH3 DECAYS TO EITHER N02 OR N03
118
-------
732,
735. II- (INH.E-3. 1) ONMisNHl* (FXP(-Tr.MPl)-l ,)
' , + TtHPl*NH3)
73V.
738,
7"o!
7'U.
712,
7«0,
715,
716,
7U7,
718,
7*(EXP(-TF"MPn-l,)
lF(IN2.f.O.?) DN02=(-«TtMPl*N02*t2)/(l. + Tf.MPl*N02)
DN0302 = -n>''02
ONEE020^EEO+1,11*"^02
IFCIBUG.EQ.I) WRITEtNOUTt 1012) XT£HP , ON02(DN0302 i ONEEO
CONTINUE
1RIN03,EG,0) GO TO 33
DECAY (SETTLING) OF N03
XTFHP = NQ3K(lNO)*THt>io3K*»-TCOR*DELT
IKI^i.EO.l) OMO»OS = N03*(FXrC-XTFHP)il,)
IF (IN3,tQ,2) DNQ30S=(-'XTEMP*N03**2)/( j , + XTE^P*N03)
CONT INUE
IF(1PO
-------
793. IF
79/1, FH^=AMIM (F MS'] «F'l-""VtF'HM)
795. F
79(,. F
797, Fl.'t<3iNH3/(IU M J+NH3)
790, J F n F N. K>, 2) F" *> = A M ;
799, I f- (I f N . E 0 . n , o R ,
COO. IF ( IUUG.FO. 1 ) WliITt (NOUT , 1 010) FHM, f. p , KN03 , FMH3 , FN
flfll , XK^FXTM IMm ,00'457 + ALG
802, I!'AI<=(IO*Cl.-FXF'(~XK*DEP7H))/(XK*DEPTH))#(24,/9.)
803. DU=ALnGC100,)/XK
80fl, IF(PU.l.T.OEPTH) IBARs. 215*10*20, /9,
805, F L =10 A H/(M|
P06, FLlH-AMJNii
807, C
808, C
809, IF(IHUG.KO. 1) WRITtCNOUTflOlb) XK , I B AF)» DU r FL • GRLIh • DEL A
010. 50 CONTINUE
611, C
812, C CALCULATE CHANGES IN BOD(POOi^03i AND MH3 DUE TO ALGAL GROWTH
6S3, C
810. IF(A»S(DF.LA) ,LT. .00001) DFLAsfl.
816, IF(ICO^B.GT, 1 1) DaODAG=0,
Sis! IFCIHUG.tQ. 1) wRlTECNOUTi 1016) OBOOAG»DPOOAG
819, jFCABSCDPrjOAG) .LE.POO) Go TO 55
820. OELA=,9+DELA
821, GO TO 50
822. 55 CONTINUE
823, C IT), THE NITROGtN DEMAND, IS POSITIVE
820, n »o - ~ 0 P 0 a i r. * u n n ..i v R / n n p p u o
825, IF-CDN.GT, ,9»No3) GO TO 60
827. . . .
828, GO TO 70
629, 60 CONTINUE
830, IFCI.-JH3.MF..O .AND. DN.LE, ,9*N03+.B*NH3) GO TO 65
831. DELA=.9*OELA
832. GO TO 50
833, 65 CONTINUE
$30, DN04AG3",9*N03
835, DNjH3AG = «(D^ + DN03AG)
836, 70 CONTINUE
637. OGIvf=OGIvE+DELA*BODOO/(APR*BODPwR)
838. lF"(IBU5.EO,n '"« I TE t NODT i 1 0 1 7 ) DN03AG »D*H3AG iOGI VE
839. C
840, C CALCULATE ALGAL Rt'SPIRATION QUANTITIES
801, C
802. ARR=NR*TEHRFA(IND)
603, DARES=-ALG*ARR*OELT
805. IF(ICOMi3.r,T.l 1) DHOOAPsO,
gif,, DPO'4AR = »DARf-S/APR
807.
808. IF(INriJ.EQ.O) GO TO 75
809.
650, GO TO 80
851. 75
652.
853.
120
-------
85a,- no CONTINUE
855, IM IllUG.Ff). 1 ) hRITECNOUTi I0l8) DARES i DBODAR F OPO«A« F DNHTAH F DN03AK F
656. *ONEtU
857. C
058. C CALCULATE AI.GAL SINKING AND DEATH
859. C
860, OALSNK = -ASR*ALG*D(.l.T/OfPTH
661, IF ( ICOMB.Lt.. I I ) OBODAS=OALSNK*F AC 1
862,
663,
86ft,
865,
866, FACX=0.
667, IFfKACMl .fiT.O.) FACX=ATn*FACMl
860. IFCF-'AC'ia.GT.O.) r ACX = F AC X + A 1 [)2*F ACM2
669, lK(FAf.M3.RT,0,) FACX = F ACXt ATOi*F ACM3
870,
871,
87a, IF(IBuG.fcQ.l) wRITECNOUTF 1019) 0 ALSNK , OBORAS , DAIND , D ALTOX F OUOOAD
873, 90 CONTINUE
87«. C
875. C CALCULATE BENTHAL RELEASE TERMS
876, C
877.
878, IFdPO/l.EO.O) GO TO 91
879, DPOUBR=BRRPOU*FAC2
B80, 9J IF( IMH3.EO.O) GO TO 92
B81, ONH3BR=8RRNH3f FAC2
682, GO TO 95
883, 92 IFCIM03.EO.O) GO TO 95
880, DN03BR=BKRNH3*FAC2
866, 95 CONTINUE
887. IFCIBUG.EQ, 1) wRITt ( NOuT , 1 020) DBOD0R , DPO'tQR, ON03BRt DNH5HRF ONEED
868, C
669, C CALCULATE OXYGEN REAERATION TERM AND BENTHAL DEMAND TERM
890, C
891, XK = DOK2(INf))
892. IF(IK2,tQ.O) GO TO 100
893. K2A20=Ul.M*V(-L**.969*OtPTH**(»l,673))/2«»D*Of AC
911. ONH3Br = pNHjr\o*cf- AC
9J2, DN03liD = ')'JD3HD*PF AC
121
-------
915,'
91h.
917.
910.
919,
920,
921,
922.
923,
92" ,
925,
926,
927.
928.
929,
930,
931,
932,
933,
934,
935.
936.
937.
930.
939.
940,
941 ,
942,
943,
945!
/I d <
/ *4 O f
947,
948.
949,
950,
951,
952.
953.
95U,
955.
956,
957.
956.
959,
960.
961,
962,
963.
964,
965.
966.
967.
968,
969,
970,
971.
972.
973.
974.
975.
110
C
C
c
1001
c
c
c
c
c
c
c
c
1008
009
010
on
012
013
1014
1015
1016
1017
1016
1019
1020
1021
1022
1023
0 f J 0 3 H 3 = D '•' n 5 M 3 * n r A C
liMC'^ = [)f-'p? + OF AC
o '•' o 3 n ? = o '< n 3 n '? t o f A c
D A R K 8 = 0 A H F. 3 * 0 f- .'• C
0 F' 0 ;i A R .-. 0 P f. /It, H * 0 F' A C
D v 0 3 A <> = D V P 3 » R * P F ' A C
D P () /J [! 9 = D P n u f i p, * p F A C
DM03[*H = ON03'JR*OFAC
D £ L 0 = - D 0
CONTINUE
HEAVY METALS
DfLHMi=HMl*(EXP("HMiK(IK'D)tOELT)M,)
[)FLKM^3HM2*(txr>(«hM2K(IfD)*DELT)»l,)
DEL H'i3 = HM3+(ExP(-'l||JiK(IM})*DELT)-l.)
IF (1BUG.EO, 1) WHITE CNCHJT . 1001) OElHMi,DELHH2,PELHM3
FORMAT C I DtLHMl iOF-l,HM2iDFtHM3= 1 , 3E16.8)
UPDATE CONSTITUENT CONCENTRATIONS
COLf =COL + OELCOt.
BOOE-BOO*OLBU010 + OLBODS'H)ROnAG*-DBODAR + pBCOAS-fpROOAD*DBODBR
NH3F.='>H3-t-pNH3Wp + 0'''M3 + DNH3V+'DNH3AG + DNH3AR + ONH38R
Nf)^t=N02 + llN02H3-fDM02
NP3F. = K03 + DNO..'tto= i tt 16. a)
FORMAT; OALSSK.nHODAS.PAL^DrDALTOX.nnnrAO^I , 5 E I 6 , R )
FQRf'iT( pRPDHR«OPO'JHS«D'J03f;HiDVH\bRiOM:rO=l . SE1 6,8)
F"CRh'AT; XK,OOn,POHEN»OGIVl'iO'»t!-. n,OFLO=l ,6Elh,8)
FORMAT; FACTOR ='.ei6,R>> ON RfACHiiih)
FpRrAT;//i INITIAL CONDITIONS FOR REACMI.IS/I io»DELT»orPTH,vELOd
*TY=I «4F1P.«/
122
-------
976.
977.
978,
979,
980.
9fll,
982,
983,
905!
906,
987.
9RO,
909,
990,
991 ,
992,
993,
990,
995,
996,
997.
998,
999.
1000.
1001.
1002.
1003.
1000.
1005,
1006.
1007.
1008.
1009.
1010.
1011,
1012.
1013.
loia.
1015,
1016,
IOIB!
1019,
1020,
1022,
1023,
1020,
1025.
1026.
1027.
1028.
1029.
1030.
1031.
1032,
1033,
1030.
1035,
1036.
H
li
C
C
c
c
c
c
3
C
C
c
c
* I COL i BOO, MM 3 i 'J02.HO 3= I ir
* i PO/J , ALG.O'.IS i, v.i'i.n/
1020 FOR'IATt/l UNAL. CONDITIONS- FUR REACHI,JS/l COL f 000 » NH3 , M02 t ^03= I I 5
+ F ife.«/ i POU, Air- inn= i. u u-.n/
* I HM 1 i MM^ , »<", 3 , 1 PT^: I i Ut 1 iS . fl//)
FORHATC I ON'03r>SiOPOOOS=l i r? E 1 6 , 8 )
SUHROUTI'.'K K?CAL(H)
*K?OPT(1 Oft)
CO"MON/ULOCK5/JJ(21)»IA,ICK(2),OELQiM(2)iK2.VFL
REAL K2
STFP-1
CAI.CU1 ATF THE
COtFFICIENT FROM VARIOUS USER
OPTIONS.
IOPT=K2UPT(IA)
IF CIOPT.EO.'.;) RETURN
GO TO ( 1 »2»
-------
1037,
1038,
1039.
1 0 'I 0 ,
1 0 'J 1 .
1 0 '1 2 ,
1043,
1 0 '1 U ,
1 0 ') 5 ,
10«6,
10U7.
1 0 2Ail,li4iF10,UFFll.(|Ff-ll,«FF11.0FMl,lFFll,lFFll«2F
*F9.'JFF-8.y,F8.a,F8,2)
300 CONTINUE
tsRITE(ivjF2) OUMltOUMg
2 FORMATC 1X.2AO)
GO TO 116
113 CONTINUE
R E A D C N I , i ) nUM],DUM2
1 FORHATC2AO)
IFCPUM1 .EfJ.ENDF) GO TO 117
GO TO 777
1)6 CONTINUE
WRITE (NJ, 121 1)
1211 FORMAT ( 1H1 )
«RITE(NJ,121)
121 FORMAT (31X. i * * * FILE K - MISCELLANEOUS VARIABLES * * *'i//)
124
-------
1096.
1099,
1100.
1101,
1102,
1103,
1 1 0 4 ,
1105,
1106.
1107,
1108,
1109,
1110,
1111.
1112.
1113,
•111').
1115,
1116.
1117.
1118,
1119,
1120,
1121 .
1122,
1123,
tl2'j|
1126,
1127,
1128,
!! 29,
1130,
1131,
1132.
1133.
1134,
1135.
1136.
1137.
1138.
1139,
1140.
1142,
1143.
1144,
1146J
1147.
1148.
1149,
1150,
1151,
1152,
1153,
1.154,
1155,
1156.
1157,
1158,
122
123
10
400
405
410
415
i 123) OMM2itHI"2 i 10
MEAO(M.1?2) niJ'M if)iiM2, IFN ,ICOL,
READCui t 1 22) OUMi t!)U'V. IMf-.AVY i ITOTN, ICHLOR
REAO{fxI . 122) Dili* 1 iOUM?« IM1, IN
-------
1159. •
1160,
1161,
162,
163,
16'!.
165,
166,
167,
160,
169,
170,
171,
172,
173,
1 7'f ,
•1175.
1176.
1177,
1178,
1179,
1180,
1181,
1182,
1183.
1184.
1165.
1186,
1187.
1188.
1189,
l!9n,
1191,
1192,
1193.
1194.
1195.
1196.
1197,
1198.
1199,
1200,
1201.
1202,
1203,
1204,
1205,
1206,
1207,
1208.
1209.
1210,
1211.
1212.
1213,
1214,
1215,
1216,
1217.
1218,
1219,
420
1002
1001
1000
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
10?0
1021
1022
1023
1024
1025
1026
1027
w R n c ( N j 1042) ML
w H I T e C '•' J 1033) APR
WRITFC-iJ 1034) f.H
V« R IT e ( '*•' J 1035) A 5 R
WKITf (NJ 1056) AMD
IFC I"KAVY.r,T .0) VRTTE (NJ, 1037) ATD
IFC IHEAVY.GT. l) «'
-------
1220,
1221,
1222.
1223.
122<4 ,
1225.
1226.
1227,
1226,
1229,
1230.
1231.
1232.
1233.
1 231 ,
1235.
1236.
1237.
1230,
1239.
1210.
1211,
1212,
1213.
1214.
1215,
1216.
1217,
1218,
1219.
1250.
1251,
1252,
1253.
1251,
1255,
1256,
1257,
1258.
1259.
1260.
1261,
1262,
1263,
1261,
1265.
1266,
1267,
1268.
1269,
1270,
1271.
1272.
1273.
1271.
1275,
1276.
1277.
1278.
1279.
1260.
1028
1029
lf>30
1031
1032
1033
103"
1035
103(.
1037
1038
1039
1010
1011
1012
1013
1011
1015
10«6
10U7
1018
1019
1050
777
778
* A 1 1 0 N El,3|XrFlO,5)
FORMATS l KTC-M4FI.IS-MF. NTiJM CONSTANT (FT. P/L) FOR PHOSPHORUS LIHITAT
*ION OF PHY TflPL ANKTON GROWTH = 1 i 9XiF10.5)
F'ORlAT(' f'lCHAFLIIJ-'iENToM Cc "iS T A ';T d'C, >. /U FOR PHOSPHnRUS LI"ITAT
*ION OF PHYTOPLANinO1; r,i)0*-TH =l, 9 X i F 1 0 . 5 )
FORKATU s' ICHAEL I ?-•". KNTOM C INSTANT (MG ^'(-N/L) FOR NITROr.EN LIMIT
+ ATIO'I OF PHYTf'PL ANKTOM GRfUIH s \ , 7 X . F 1 0 , 5)
FORl'AT(l !•: ICHAfLIS",''.OJT[lN CON STA'^T ("-T. KHJ-N/L) FOR NITROGEN LIMIT
*ATION OF PHYTOPL.'.'-.'KTON GROWTH =li7X,F10.5)
FORnATCI LIGHT INTFUSITY CALCULATION FACTOR ( L AhGLFY 3 /M I N) =lf«UXi
+F10.5)
F 0 R ' 1 A T ( 1 P H Y T 0 P L A N K T 0 N TO PHOSPHORUS RATIO = I t 6 0 X i F 1 0 . 5 )
FURM£T(I PHYT UPLAND TON RtSP I *< A T 10') FACTOR r 1 , h 1 X , F 1 0 , ri )
FORMAT(I PHYTOPLA'.'K TO'.' SINKING RATL (FT/HW) = ' f 59X i F 1 0 , 5)
FORI'AH' PHYTOPl.A'JKTUN NATURAL DEATH RATE C /HR) = • » 51X , F 1 0 . 5 )
FOR^'AT(I PHYTOPLA'iKTON TOxIC OF.ATH COEFFICIENT FOR HEAVY f'FTAL 1 =
* " i38X,FlO.'j)
FORHAT(I PHYTOPI.ANKTON TOXIC DEATH COEFFICIENT FOR HEAVY ^ETAL 2 ^
*' .38X.F10.5)
FORMAT(l PHYTOPLANKTON TOXIC DEATH COEFFICIENT FOR HEAVY METAL 3 =
*' i38X,F10.S)
FORKATC1 HOD 6ENTHAL RELEASE RATE (MG/SOijiRE MFTFR-HR) =i»18XiF10.
*5)
FORHATM PHOSPHORUS DENTHAL RELEASE RATE (MG/SQUARE METER-HR) ='»4
*1X»(- 10,5)
FORMAT(l NITROGEN BENTHAL RELEASE RATE (MG/SOUARE METER-HR) =l(13X
*iF10.5)
FORr'.AT(l HENTHAL OXYGEN DEMAND (MG/SUijARE WETF.R-HR) = 1 , 5 1 X t F 1 0 ,5)
FORHATCI FRACTION OF HEAVY METAL 1 IN ION FORM = « i5bX i F10 .5)
FORf'.ATC' FRACTIO'J OF HEAVY HFTAL 2 IN ION FOR^ = l » 56X , F 1 0 ,5)
FORMAT(I FRACTION OF HEAVY METAL 3 IN ION FORM = 1 , 56X i F 1 0 .5)
FORMATM TEMPERATURE COKRtCllUN C 0 N S T A N T FOR BOD REACTION COEFFICI
*ENT =l .33XtF10,5)
FORHATC TF:UPFRATURE CORRECTION CONSTANT FOR Noa DECAY COEFFICIENT
* =l t36X.Fl 0,5)
FORMATC' TEMPERATURE CORRECTION CONSTANT FOR N03 DECAY COEFFICIENT
* =l ,36X(F10.5)
FCRMATM TEMPERATURE CORRECTION CONSTANT FOR POI DECAY COEFFICIENT
* =1 ,36X,F10,5) '
RtAD(Nl,123) DUHliDUM2
IF (DU^l .EO.ENDF) l"R I TE ( N J ( 23 DUMi,DUMa
IF (D'J^l ,\'E .ENDF) GO TO 777
RETURN
WRITE CNJ.778)
FORMAT ( 1 DATA CARD ERRORI )
STOP
ESD
SUBROUTINE PLOT(N,A,TiLAB)
DIMENSION FINCR(12).A(nfICHAR(3),T(l)»IPLOT(100)
INTEGER T
DIMENSION LABC6)
DATA FlN'CR/,01(,05f.lt.25t,5il.i2.f3,,1,i5,<10,ilOO,/
DATA ICHAR/i l , i* l i i 0 l /
ILAST=100
FMIN=1,E10
FKAX=-1,E10
DO 1 1 = 1, N
IF (A( J ) ,LT .FMIN) FMIN=A(I)
IF(A{ I ) .GT.FMAX) FMAX = A(I)
1 CONTINUE
DX=(FMAX-FMIN)/«o,
127
-------
1301,
12R2,
120b,
1206,
1307,
1288,
1289,
1290,
1291,
1292,
129i.
1295.
1396.
1297,
1298,
1399.
1300.
1301,
1302,
1303,
130H,
1305,
1306,
1307.
1308.
1309,
1310,
1311.
1312.
1313.
1311.
1315.
1316.
1317.
1318.
1319.
1320.
1321.
1322.
1323,
1324.
1325.
1326.
1327.
1328.
1329,
1330.
1331.
1332.
1333.
1331,
1335.
1336,
1337.
1338,
1339.
1300.
1311.
on 10 J=lilP
K = J
If (Fir-CR(J).r,T.DX) r;o TO 11
10 COMTlNUh
11 IKK. to. I)
I M A X e F K A X / 0 X -
600
601
20
603
30
S*AX=IMAX*DX
* R I T f- ( h i (, 0 0 ) LAB
FORMAT ( I'M iHOX»6A«)
wRITE(ht601) SM1N»SMAX
KOKMATC/il RfACH
00 30 J=lffl
00 20 IP-fdLAST
IPLOT(IP)=ICHAR(l)
IPLOT(13=ICH*R{3>
IPLOT(It*ST)=ICHA«(3)
IP=(A(J)»SMIN)*ILAST/(SMAX»SMIN)
IF(IP.LK.O) IP=1
IFOP.GT.1LASTJ IP=ILAST
IPLOT(IP3=ICMAR(2)
wRIU(6,602) T(J) ,A(J) , IPI.OT
FORMATUH . 1 1 0 , E 1 2 , 5 i 7X i I 00 A 1 3
CONTINUE
RETURN
DEP, VAR , I , F 1 0 , 3 » 85X i F 1 0 . 3/3
SUBROUTINE PLTSET (0 , IF-'UAG)
DATA tNDF/'ENDF'/
DIMFNSION 0(100)13
fJl^ElvSIOM INR(100)iINt)S<«3)iXPLT(t003fYPLTUOO)
INTEGf-tt XPUT
DIMENSION TITLE(6f«3)iTl(60)»T2C60)iT3(60)iTU(<>0)»T5(18)
DIMENSION fl(2!583
EOUIVAUfNCF (TIiTITLE)
FnuiVAuEfyCE (TI»Tl)t(TI(61)iT2)i(Tl(l?niT3)i(TI(181),Tfl3iCTI(2fll
*3iT5)
DATA Tl/
* REACI
* RIV i
* REACI
* FLO* I
* MINI"
* MINI '
* 00 A I
* K 2 F' 1
* TRAVi
* 1 MEAN i
DATA T
*IBOD I
i i
1 MEAM
I DO A I
teon i
1 NH J 1
1 M K 3 1
1^02 '
IK.Q2 1
IN03 I
1 H NU' »
1 MJLt ' »
1 H L E ' .
1 AT
1 MU"4
MUM
T RE
OR H
EL T
VEL
/
AT R
i
t
f
i
i
i
»
t
t
DEPi ,
T RU ,
AT Ri •
AT Hi ,
AT HI ,
AT RI i
AT R' i
I AT R i i
M B t. R I f I 1
TO HI . (EACH"
fvGTHl . 1
END 1 , IOF R
DO I I i IN RE
DO 11 ' i I I L E
ACH i . I EMO
EACH" « I
IME i . i
OCITl t i Y
EACHI.I END
"ft
TH 1,1
ACH 1 t 1 STAR
FACHI ,i STA
E t C K i , 1 STA
E « C H I • ' E fj 0
F i r. H 1 , I STA
EACH 1,1 ENH
EACWi,! STA
t
HtADl
1
EACHi
ACH
(MG/
(HG
T CM
RT(*
RT (K
(HT,
HT(M
(PC
RT(r
1
1
1
i-CFS
1
1
'L3
1
t
i
I/L3
i
i
ir./L)
IG/L)
'G/L3
i/L)
'G/L3
'/L>
'G/L3 /
DATA TV
128
-------
1 3 « 2
13l ,
AT R i i EACH i ,
AT R I i f A C H I »
AT R ' . f ACHI ,
AT R ' . E A C H 1 i
AT Rl , EACHI ,
AT R 1 ( EACH",
AT R 1 t f ACHI ,
t
AT RlflMCHIt
AT R i . IEACHI ,
AT R" i IEACHI ,
IN = 5
NJ = 6
f M 0 " t
STA 1 i
K M) 1 t
y T A i i
f.Nfil t
STA 1 i
t"N[i " ,
S T A " t
t NO 1 ,
STA 1 .
ENOI t
STA 1 ,
t.NOI ,
STA" ,
L'NIH ,
STA " .
END' t
STA I ,
EIH'l f
STA' »
END" i
STA I ,
END" .
' ( M G
I R T ( M
" ("r.
' R T ( :'
1 ( 11 G
' R T ( M
" (IT,
IRT(M
" (T,
IRT(M
' (MG
"RT(M
' (MG
'RT(M
" (Mfi
' RT (M
' ("G
I R T ( M
" (MG
IRT(M
1 (MG
1 R T ( M
' (MG
• '/I)
• IG/D
• l/L)
i I C. / L )
. '/D
• "'•/!.)
. '/D
i i (; / L )
. i/L)
f iG/t.)
i '/D
> 'G/L)
i i/D
• ' G/L)
• i/L)
I'G/L)
• '/D
t 'G/D
• l/L)
• 'G/D
t '/D
,'G/I.)
i l/L)
,
i
,
,
t
r
,
i
f
i
f
i
i
i
f
,
/
t
t
/
IFdFLAG.EO.O) GO TO 20
READCINilOl) DUMliDUM2t
IOUT=0
IFtDUMi.EQ.KNOF) IOUT=1
NR t MIND
1 VM , TFfnllMljCfj.f-Nnt:) t} £ T l_l W Kl
1371
1375
1376
1377
1378
1379
1380
1381
1382
1383
1384
1385
1366
1387
1380
13B9
1390
1391
1392
1393
13941
1395
1396
1397
1398
1399
14100
14J01
1002
•
,
•
.
.
*
,
,
,
•
,
,
,
m
,
,
»
,
,
,
,
,
»
,
,
.
•
.
•
60
410
20
S
10
100
101
102
103
1 041
READ(IN,100) (INR(I),I=
R t" A 0 (I N , 1 0 0 ) ( I N D :> ( 1 ) i I
Rf.ADdN, 101 ) DUM1,DUM2
IF(DUM1 .NE.tNDF) GO TO
rtRITE(NJ, 10JI)
"JRITECNJ, 105) (INR(I)i
WRlTE(Njt 106)
DO 60 1 = 1, NINO
IM02=INDSd)
^RITE(NJ, 107) (TITLE(J
WRITECNJ, 1 02) DUM1,OUM2
RETURN
*RITE(NJ( 103)
STOP
CONTINUE
IFdOUT.KO.l) Rf.TURN
DO 10 1=1 i NINO
IND2=INDS ( I )
DO 5 J=1»NR
IN01=IMR( J)
XPLT(J)=IN01
YPLT( J)=0(IN01 i IND2)
CALL PLOT(NR,YPLTiXPLT(
RETURN
FnR"Ay ( 1 OX » 1 41 15 J
FORMAT C2AU.HX, 215)
FORMiT(lXi2A
-------
1403,
1404,
1405,
1406.
1407,
1408.
140V.
1410.
1411.
1412.
1113.
1415,
10'j
106
107
FCfO'AT(l OUANTITItS URL HF. PLOTTFD FOR THf FOLLOr I kIC,
* ( 1 0 X • 1 4 I S ))
FORMAT (//I TH£ FOLLOWING UUANTlTItS WILL ttf. PLOTTFOi/)
H" 0 M r A T C 1 0 X ( 6 A q )
E N ()
SlMUHiuTINC KFIH!
CO'iMON/lH."CK?/JMT (r>9l) r IN'ITC 10) i IAUG(?0)
/ilLOCKi/Tfc ("-'H 1 25 i I DM CH (?n . 10)
f 4 ) .'.'IMIT •"••JUNf. (/i) ."(Qut^F IML(H) i JA
COMMON /^ISC/IO,C*LCK2»NFLOW|NCC1N
1116,
1417.
1110,
1119.
1420,
1421,
1422.
1123.
1121.
'1425.
1126,
1127.
142fl.
1129,
1130,
1131.
1132.
1433,
14}4,
1435.
1436.
1437.
1438,
1439,
1140.
1441,
1442,
1443,
1414,
1415,
1446,
1117.
1448.
1119,
1450,
1451.
1452.
1453.
1454,
1455.
1156.
1157.
1156.
1159.
1460.
1461 .
1162.
1463.
C
C
C
C
C
00 i I=I.NINIT
J A = I N I T C I )
3
2
C
C
C
U
10
C
C
C
c
I F ( I A u G (
IDHCH (JA
IDMCH (JA
GO TO H
IDMCH (JA
CONT I NUE
QUP=CONDI
00 10 K=l
XCON£(K)=
JA) ) 2, ?. , 3
,D=JA
• 2) = 0
f l)nO
(JAiNFLOW)
(NCOS)
COMO! ( J A »K)
CALL TR1BD
C
C
C
C
C
C
C
C
c
20
C
C
C
c
c
1
c
c
c
c
CALL SCAN
CONDE( JA,
00 20 K=l
CONPECJAt
I F ( M ) 1 , 1 ,
CONTINUE
.
'
NFLOW)=OUP
• MCON
K)=XCONt(K)
5
STEP-1
SFT '•'ACHINF. IOF^TIFIFR s FOR
COUNTI'Jf. UPSTREAM AUGMENTATION,
STEP«1
CALL T R I B D
CALL SCAN
SET FINAL CONDITIONS TO
CALCULATED CONDITIONS
SUP-6
REPEAT STEP-! THRU STFP-5 FOR
EACH HEAOKATF.R STHETCH,
STEP.7
RETURN TO CALLFR (R U N 0 N),
130
-------
RETURN
1 1 (' 7 ,
111.8,
1*169,
1«70.
1171,
117?,
1 '1 7 5 ,
1 17 1 ,
1175.
1176,
H77,
1176,
H79.
1130.
1181,
1182,
1183,
1181.
1185,
1186.
1«87.
1188,
1189.
1190,
1191,
1192.
1193.
1191,
1195 ,
1196,
1197.
1198,
1199,
1500,
1501.
1502.
1503.
1501.
1505,
1506,
1507.
1506,
1509.
1510.
1511,
1512.
1513,
1510.
1515.
1516.
1517,
1518,
1519,
1520.
1521.
1522,
1523.
1521,
C
C
C
C
C
5
15
10
C
C
C
C
C
1
333
331
2
3
16
OATAU4.7) IS INfRf MfNTAl. FLOW IF REACH IfiMGTH IS NON7F--RO.
IT IS A POIt'T SOUPCF:((1R ,1 1 T^OH A r. At. ) IF THK «F4CH LENGTH IS ZF.RO,
IF THK Tf'FATMHJT FACTOR IS TO HE APPLIED. THE 00 CONCFN THAT ION
IS INPUT AS A M-r.AFIVE MUKflf.K,
COMMON/ WORK] /xCDNil (20)
CO''>'ON/m(iRK2/xcONt (20)
C 0 Ll. '«OK' / lit OCM/C.w "!'•'( 100) . OATA( 100. 10)
CO'tMON/h'LOCK'j/JJm) • TF .Tf->Pm iQUP.F INI, (3) ) IA
COMMON/MISC/IOiCAl.CK2»sIFLOK»NCOfJiTOTFLO
COVKOU/OPTION/JFNC3) t I COMB
se.T UP INITIAL CONDITIONS IN REACH IA
X^AX=l ,
IF (DATA(IA, ) ) .Nt .0.) XMAX=10,
IF ( !C(V'H,f- a, in) XMAX = I
FACFLO=OATA(IA»7)/XMAX
TOTFLOcOuf'+FACFLO
FAC=1,
IF('.UTA(IA» 1) .Nfc.O.) GO TO 5
IF (OATA(IA.8) .LT.O.) FAC=1.-TF
IF(OATA(IA,8) .LT.O.) DATA(IA,0)=-DATA(IA»8)
CONTINUE
DO 10 I~l rNCON
IF(DATA(IA,7).GE,0.) GO TO 15
VAR=XCONE(I)
GO TO 10
IF(I.LT,8) VAR=DATA(IA,7+I)*FAC
IFd.GE.8) VAR = r>ATA ( I A, lh+I)*FAC
XCOM(I)=(XCONE(I)*OUPtFACFLO*VAR)/TOTFLO
RETURN
ENO
SUBROUTINE RMATC
COHHON/BLOCK2/JNIT(20)iF(21)tIONE(20),G(170).JUNC(20.3).INIT(10).
*IAUG(20)
COMMON/PL DC K3/CONDZ(10,20)»CONDIC20,?0)iCONDE(20i20)
COMMON/ BLOC Ka/TEMMO( 12) ,IDMCH(20i 10)
COMMON /BLOCK5/JJ,KK, II»D(1) iNTRIB,NRFAiNINIT,NJUNCiELEV(5) iM,OUPi
*FIML(2) t JA
COfMOK/«ORKl/XCONl(20)/WORK2/XCONE(20)
COMHON/MISC/IO.CALCKa.NFLOW.NCON.TOTFLO
STtP-1
FIND THE HEADWATFR STRETCHES
THAT tNTER A JUNCTION.
DO 1 1=1.20
IONE(I)=0
IF(NJUNC)333.333.331
IONE(1)=1
DO 2 1=1 .20
JNIT(I)=0
I1'JIT = NINIT
DO 3 1 = 1 . UNIT
J^'IT(I) = I'JIT(I)
J=INIT(I)
IONE ( J) a 1
DO H I = 1»NJIINC
DO S J=l .UNIT
IF(JUNC(Itl)->JNIT(J))5,6i5
131
-------
1525. •
1526,
1527.
1528,
1529.
1530.
1531.
1532.
1533,
1530,
1535,
1536,
1537,
1530.
1539.
1500,
•1501,
1502,
1503,
1500,
1505,
1506,
1507,
1508,
1509,
1550,
1551,
1552.
1553,
1550,
1555,
1 S 5 h .
1557,
1558,
1559,
1560,
1561.
1562,
1563.
1560,
1565.
1566.
1567.
1568.
1569,
1570.
1571.
1572,
1573,
1570,
1575,
1576,
1577,
1576.
•1579,
1580.
158J.
1582,
1583,
1580,
1585.
5
6
7
6
C
C
C
C
20
C
C
C
C
C
C
IS
10
9
12
11
C
C
C
C
C
C
C
C
C
C
C
C
C
30
C
C
C
C
C
C
CONTINUE
no TO o
DO 7 j=i i UNIT
IKCJUNC( I »2)-JNIT
-------
1507,
1508,
159oi
I F (M M 0 i 1 " i 1 5
10 IONc(KK)=l
JMT( I1NIT)=KK
o coin INUF
DO 17 1=1 iNTRIH
1592.
1S93,
1590,
1595.
1596,
1597,
1598,
1S99.
1600,
1601,
1602,
1603,
1600,
1605,
1606,
1607.
1608.
1609.
1610.
1611,
1612.
1613,
1610.
1615.
1616,
1617.
1618.
1619,
1620,
1621.
1622.
1623.
1620,
1625,
1626.
1627.
1628,
1629,
1630.
1631,
1632.
1633.
1630,
1635,
1636,
1637.
1638.
1639,
1600,
1601.
1 A U P
1 O u C. .
1603,
1600,
1605,
1606,
17
C
C
C
C
15
IFdONEd)) 16, 16t 17
CONTINUE
RETURN
FNO
SUBROUTINE RUNON
COMMON/BLOCKl/fRMINC 0100) i F I N IS (
COMMON/PLOCKrVJMiTt HI) . TITLF (20
COMhON/»I.OCK3/cONOZ ( lOi 20) fCUf-OI
C 0 M M 0 N / (i L 0 C K « / T E M M 0 ( 3 it 2 ) i a t. A 5 0 U (
COMMOK/RLUC*'5/JJ(fl) , NRfc-Af >.'IMT |N
STEP-0
RETURM TO CALLER (R U N 0 N),
100,50) tRIDENTf 100(5)
).IP«0(20.?0)iJUNC(20t'5)iIMIT(l(
t20(2P) iCO.MOE(20,20)
12)
))
JUNCC2) iDOLEViTF»TF.Mp,cSATiM,QUP(5
*) i ICK,FIrtLN(3) ,hJ»K2(2) »NSF.AS»:iRUN.TFN
C
C
C
C
C
1
C
C
C
C
C
3
C
C
C
C
C
3335
388
333
COMMON /TMPER/TtMpAV,Tt.MRtA( 100)
COMMQN/Ml3C/IO(2)iNFLPW
COMMON/OPTION /IF fJ,IK2,ICOL,ICOMH
COMMON/OPT2/IHtAVY,lTOTN,ICHLOR
DATA ASTt.R/l*!/
00 1 IE 1,10
00 1 J=l»20
CONDId . J)=CONDZ d»J)
CALL OOSAT
McO
IFdCK.EQ.l ,OR. ICOMD.fQ.18)
,SAFR£A( 100)
(IMH3,IN02fIM03,IPOO,IALG
STEP-i-1
SET ALL HEARWATE.R CONDITIONS
tQUAL TO ZERO,
STEP-2
CALL OOSAT
GO TO 3333
STEP-0
URITE HF.ADINT, FOR INTERMEDIATE
CRITIC (NJ.388)
FORMAT (1HO)
WRITE (NJ.333)
FORMAT ( 7X|OOH NO. RIVER MILE
* IOLVED OXYHEN
* 7X.08H OF TO HEAD
* "TART AT END AT
* 7X,08HREACH OF hFACH
REACH SUMMARY,
FLOK 0,0. RIVER MILE DISS
BODi ,/,
»
RATE MIU, ATMIN, AT3?
START AT ENpli/t
(CFS) (MG/L) 0,0. (MC
»
* i/L) (MG/L) (MG/L) (Mf;/i.)'»/i
C
C
C
STf.P-5
CALL R E I N I
133
-------
1607.
16/19,
1650.
1651 .
16S2.
1653,
16'>u,
1655,
1656,
1657,
1656,
1660,
1661,
1662.
1663.
166U,
1665,
1666,
1667.
1668,
1669,
1670.
1671,
1672,
1673.
1674.
1675.
1676.
1677,
1678.
1679.
1680,
1681.
1682.
1683.
168«.
1685,
1686,
1687,
1688.
1689,
1690,
1691.
1692,
1693.
1694,
1695.
1696,
1697.
1696,
1699,
1700.
1701.
1702.
1703.
1700,
1705.
1706.
1707.
3333 CALL PF.IM
IF ( ") U t « . 3
CALL RHATC
IF(M)5i5i 3
STEP-6
CALL R M A T C
STEP-7
V-RITE HEADING FOR FINAL
SUMMARY ,
5 WRITE (NJ.2055)
2055 FCJHMAT (1H1)
HHITK (NJ,2U)
21 FORMAT (33Xt<45H* ****F1NAL
WRITE (MJ.336) CTITLE(I) . 1 = 1 » 18)
336 FORMAT (19X.UA/1,//)
SUMMARY**** *»//)
335 FORMAT ( 15X , 15MNIJMUER Of- RUN =, 15 • 36X , 1 9HT ARGET D.O. LEVEL
15X.17MSEASON OF
F5,2f//)
.A3,36Xi 19HMFAM TEMPERATURE a,
25 FORMAT ClOXiU'jH NO, IDENTIFICATION RIVER MILE RtACH i
*50H FLOW D.O, RIVER MiLf-- DISSOLVED OXYRFN Tt-np,/,
* iox,^5H OF OF TO HE-AD LENGTH t
* I RATfc MIN, AT HIN. AT START AT f-.M? (CENT) I/
* 10X(U5HREACH HEACH OF REACH CMILESJt
(MC/L> s/r
30
SUP-B
WRITE FINAL SUMMARY,
00 88 I=1,NHEA
IF(FI'JIS(I.3).EO,0.) GO TO 88
IR = F INISt I»l)
WRITE (NJ,28) IR, (RIDfc'NT(I.J) iJ=1.5)» (FINIS(IiJ)iJ = 2i6)i
* FlNISCl»10).FINISCIi7)iTtHKtA(I)
68 CONTINUE
28 FORMAT tlOX.I3i3X.5Ati2X,F7.1«i
-------
1708,
1709,
1710.
1711.
171?.
1713.
171 «,
1715,
171&,
1717.
1710,
1719,
1720,
1721,
172?,
1723,
J«aHNIS(I,l)
HDK = HM.S(I»«)/{FINIS(IflO)*Fl»j IS (1,13))
h"ITK (nJ,V) I*. (KIIH.NT(I»J) t J=l .
» (FIMIS(T.J)i.J = OtlO).
*FlNI5(I,13)tV,inK,x.Sl_OP
MlS(I • 13)
iF INISU, IS)
IF (CHK.GT , 3?,) wRlTE(Nj,X,I5»5XibAU,lX.F7.2.F10.2.FlO,2,F10.3.2l-9,2iF7.2.F12,6)
IFf INHi.erj.O .AMD. IN02.EO.O .AND. INOi.CO.O .AMD.
* IPOO,HJ.O) GO TO 2051
KRITE(NJ,20r>5)
1725.
1726,
1727,
1728,
1729.
1730,
1731,
1732.
1733.
1730,
1735.
1736,
1737,
1738,
1739.
17«0.
1712)
1703,
17«s|
170&!
1707,
1708,
17«9,
1750,
1751.
1752.
1753.
1750.
1755,
1756.
1757,
1758,
1759,
1760,
1761,
1762.
1763.
1760.
1765.
1766.
1767.
1768.
210
205
205
206
220
220
207
230
206
Mj.JUfe) (TITLE(I)»I=l»ie)
uHITE(NJ,335) NRijN . DOLE v i SKA SON ( NSEAS )» TEMPAv
00 210 1=1,HREA
IF(FIN'IS(Ii3) .ECJ.O.) GO TO 210
IR=FINISCI i 1)
fcRIlE(>JJ,20b) IRi(PIDENT(l»J)«J=li5),(FINlS(l'tJ)»J=16,23)
COHTINUE
FOHMAT(IJ3»5X,5A«,1X|) NHUN.DOI.FVtSeASONCNSEAS) iTEHPAV
WRlTt (NJ.201)
00 220 I=1,NREA
IF(FINlS(Ii3) .F.0.0.) GO TO 220
FnRMATCI13i5Xi5A«,lXtF7,3,5XiF7,3,F9.1i3XfF9.1.2Xf2(F7.3i5X»F7.3i
*2X))
CONTINUE
CONTINUE
IKIHEAVY.EO.O .AND, ITOTN.EQ.O .AND. ICHLOK.EQ.O) RO TO 2301
(-RITE(NJ.2055)
WRITE(NJ,?0)
WRITEtMi 336) (TITLE(I)t I = lt IB)
KRITE(NJ,335) NKuNiDOLEV.SEASUN(NSEAS),TEMPAV
tiR!Tt(NJ»202)
IAS = 0
DO 230 I = 1»MRF-A
IF(FINlS(Ii3) ,F.Q.O.) GO TO 230
IRcFINiSd. 1)
wHITL(KJ,?n7) IHi(RIDENT(I.J)tJ=l,5),FINIS(I,32),FINIS(I,33)i
*FINlS(Ii3«)iFIMIS(I,ou),(FIMS(I,J),J=35»39)
FQWMAT(I13»bX.5AO,lXrF7.3t5x«F7.3t2XiF7,3«2XtAl,2XiF7.3»2Xt
*F7,3tbX.F7.3i2XiF7,4,5XiF7.3)
IF(FIMS(IiU<4).tO.ASTER) IAS=1
CONTINUE
IFflAS.EQ.l) WRITE(NJ,208)
FORMAT(/I * INDICATES SUM TOTAL OF TMF MTROGFN n CONSTITUFSTS BE
*ING HODFLEO f-XCfctDS THE TOTAL NITROGEN HEING MODELED AS 4 CONSEHVA
*TIVE,D
135
-------
1770.
1771,
1772,
1773,
177«.
1775,
1776.
1777,
1770,
1770.
17RO,
1781,
1702,
1703,
17CU.
V/05,
1786,
1707,
1788,
1789,
1790,
1791,
1792.
1793,
1794.
1795,
1796,
1797,
1798,
1799.
« O A A
ieoi!
1802,
1803.
1800,
1805.
1806,
1807.
1808.
1809,
1810,
1811,
1812.
less,
2301 CONTINUE
ii"( IHF.AVY.I.T.?) GO TO 2/401
1815,
1816.
1617.
1818.
161V,
1820,
1B21.
1822,
1823,
1824,
1825,
Ift26,
1827.
1826,
1829.
uwiTC(NJ,33M (TITLE ( I)t 1 = 1 i10)
MMTfc CNJ .33b) NHUNiOOLKV tSEASON(MSEAS)|TF MPAV
GO TO 240
210
2001
, (FINISH ?
RKHHRED »/?
(CFSJ ?/?
200
201
DO 2'10 I=1?NR£A
IF(F IMS(I i3) .CQ.O.)
IR = F INISd? 1)
n«ITF(NJ?205) IR. (RIOENTU t J) i J=li
CONTINUE.
CONTINUE
WRITE (NJ.2055)
KRITE (NJ,2/I)
V.RI1E (NJ.3U) (TITLE(I) i 1 = 1 r 10)
WRITE(NJi33S) NRUN.nOLEV i SEASON(NSEAS)iTEMPAV
WRITE (NJ,55)
55 FORMAT (3ix(unH NO. NO, INITIAL. FINAL AUGMENTATION)/?
* 3]X(08H OF OF FLOW FLOW
* 31Xf«HHHfADwATER STRETCH (CFS) (CFS)
* 31X?«6H , , ., ?/)
DO 100 I=1?NINIT
K = INIT(I)
OAUG=CONDI(K?NFLOH)-CONDZtKiNFLOW)
'*'RITE(NJ?56) J i IMT (I ) ?CQl>D/ (KiMFLOK) iCONOI (K?NFLO«) ?OAUG
56 FORtiAT (3lX,I5.5X,l5i«X,F7.1ilX,F7.1,3X,F7,l)
100 CONTINUE
STf-P-9
RETURN TO CALLER (o o s A ro,
PLOT DESIRED QUANTITIES FROM FINIS ARRAY,
CALL PLTSfT(FINISiO)
RETURN
FOR|JAT(1IX?"NO. IDENTIFICATION NH3-N CONCENTRATION N02»
*N CONCENTRATION N03«N CONCENTRATION P04-P CONCENTRATION"/
nx?"OF OF AT START AT END AT START
AT END AT START AT END AT START AT ENOI/
lOXi"PEACH REACH (MG/L) (MG/L) (MG/L)
(MG/L) (MG/L) (MG/L) (MG/L) (MG/L)'/
lOXi"
, '/)
IDENTIFICATION ALGAE CONCENTRATION COH
1 CONCENTRATION MET-? CONCENTRATION"/
OF AT START AT END AT START
AT END AT START AT END"/
REACH (MG/L) (MG/L) (MPN/100)
(HG/L) (MG/L) (MG/L)1/
202
FORMAT( 1 Ix? "NO,
*F CONCENTRATION M
* 1IXi I OK
* AT END AT START
*iOXi "REACH
* (MPN/100) (MG/L)
*
-------
1630,
1831 ,
1832.
1H35.
1836,
1037.
Ifl38,
1839.
1 0 '1 0 ,
* AT
*10X i
* (
*10Xi
tf
END I/
' R t. A C H
MG/L) ' /
SUBROUTINE; SCAN
C 0 M l1 0 N / H L U C K 2 / J N I T (
C 0 M >U) S / 11 L (! C K 3 / C 0 N 0 /
OF AT START
REACH (MG/L)
1(100)
'( in, 20) •CO^iD((20,20)
AT tNO AT START
(MG/L) (MG/L)
,IAUG(20)
i I 0 M C H ( 2 0 i 10)
COMMON/ HLOCK5/JJ(12)il.'OLKV,TK (3), M,QUP,HNU(2)iJA,IA
CHECK FOR AVAILARLE
1843.
18«« ,
1 8 1 'j ,
10«6.
1817.
IB'tB.
18/19.
1850,
1851,
1852.
1853,
1351,
Ifl55,
18S6,
1857,
1858.
1859,
I860,.
1661 ,
1862.
1863.
1060,
1865,
1866,
1867,
1868,
1869,
1870,
1871.
1872,
1873.
1870.
1875,
1876,
1877,
1878.
1879.
1880,
1861,
1882,
1883,
1 88(1 ,
1885,
1686,
1887,
1888,
1889,
1890,
C
C
c
c
c
2
C
c
c
c
c
3
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
7
6
e
c
c
c
c
c
STEP'
CHECI
A U G M 1
1 = 0
IFUAUG(JA)) 1, 1,2
1=1 + 1
STEP<
CHtCI
STRE"
IF(lORD(JAiI)Hilt3
IA=IORO(JA,I)
STEP'
CAl.Cl
HEACI
Z=DOLEV -CRMIN(IA)
IF(Z) 2,2,5
STfP
CALC>
OEFI
5 QADD = QUP*(Z/DOLEV * 0 ,25* (Z/DOLEV) **2)
STEP
CHEC
SOUR
AVAI
L=l
LL°0
Ms 1
IF(IDMCH (JA»L))1»8(6
LL=LL+i
L = L+ 1
GO TO 7
•13 LI.
STFP
DIVI
AMON
AUGM
CHtCK FOR Rt'ACHF.S IN A
CALCULATE D.O. DEFICIT FOR
CALCULATE AMOUNT OF DILUTION
NEEDED TO SATISFY THE
CHECK TO SEC HOW MAMY UPSTREAM
SOURCES OF AUGMENTATION ARE
DIVIDE DILUTION WATFR F.OUALLY
AMONG ALL SOURCES OF
137
-------
109J.
1092,
1893.
1094.
1895.
1096,
1097,
1098,
1899.
1900.
1901,
1902.
1903,
1904,
1905.
1906.
1907.
1900,
1909,
1910,
1911,
1912,
1913,
1914,
1915,
1916,
1917,
1910,
1919,
1920,
1921,
1922,
1923.
1924.
1925.
1926.
1927.
1920,
1929.
1930.
1931.
1932.
1933.
1934,
1935,
1936,
1937.
1938,
1939,
1942^
19«3,
1945J
1946,
194ft!
1949,
1950,
1951,
C
10
C
C
C
C
C
9
C
C
C
C
C
1
800
L=l
IFUDMCH (JA.L))1.1.9
STFH-7
Aon DILUTION WATER TO
AUGMENTATION SOURCES,
tJA.L)
CONDI (IB i I TO f V f ICHI.OR
COMMON /OPT3/IP.lNH,lN2rlN*
DATA IFlRt I2NOiN.H3f N02»N03tLPOO|N,LP/l 1ST 'tl2ND I » ' NH3- ' • ' N02- I •
* IN03-I i IP01» I i IN ' » 'P '/
DIMENSION LAB(IO)
DATA LAD/l l(l I » I MODE I t I LED '»IBY 't' '•'
*ltR Rl i 1EACT' i I ION I/
JKIR5T=1
NOUT=6
!NH3=0
IN02=0
I P 0 « = 0
!ALG=0
IF(ICOMB,EO,18) GO TO 800
IMHJ=1
IFdCOMB.GE.7 .AND, ICOM8.LE.il .OR. ICOMB.GE.16) INH3=0
IFfICOMB.tO.19 .OR. ICOMB.EO.aO) IN«3=1
IN02=0
IFf ICOMR.LE.2 .OR. ICOM0,EQ,« .OR, ICOM0.E0.12 .OR.
ICOMB.ET, 13) IN02=1
IFUCOMB.EO. 19) IN02=1
.OR.
IFfICOMB.EO.10 .OR. ICOMB.EO.il) IN03=0
JFflCOMB.GT.Pl) IN03=0
IPO
-------
1952, .
1953,
1951,
1955.
1956,
1957.
1950,
1959,
1960,
1961,
1962,
1963,
1961,
1965,
1966,
1967,
.I960,
1969.
1970.
1971,
197R,
1973.
1971,
1975,
1976,
1977,
1978,
1979,
1980,
1981,
1982,
1933,
198U,
1965,
1966.
1987.
1988.
1989.
1990,
1991.
1992,
1993,
1991 ,
1995,
1996.
1997,
1998,
1999,
2000,
2001,
2002,
2003,
200U,
2005,
2006,
2007,
2008,
2009,
2010,
2011,
2012.
10
20
30
10
999
1000
1001
1002
1003
1007
loos
1009
1010
1111
1112
LA!.\(6) = IFIR
IFUNH.f.0.2) LAtU6) = I2NO
wRITECNOUT) 1003) LAB
CONTINUE
IF-(IM02,EQ,0) CO TO 20
LAH( 1 )=*02
L A B t 2 ) 3 N
LA5(M = IFIR
IF( IN2,KT,2) LAB(6) = I2NO
v.RITf-(NOUT»1003) LAB
CONTINUE
If- (IN03.EO.O) CO TO 30
LAB(1)=N03
L AB ( 2) sN
LAB(6)=IFIR
I F C I M 3 i F. 0 . 2 ) L A H ( 6 ) = I ?. N 0
k'UTECNOUTi 1003) LAB
CONTINUE
IFdHO^.EQ.O) GO TO 10
LAB ( 1 3 =LP04
LAB(2)=LP
LAB C 6) = IF IR
IF C I P. E 0.2) I.AQ(6) = I2ND
V'RITE(NOUT . 1003) LAB
CONTINUE
IF(IALG.EQ.l) WRITE(NOUT,1007)
IF(IHf-AVY.tO.l) wRITF (MOOT. 1008)
IF(IHf.AVY.GT.l) WRITE (NO UT. 1111)
1FCITOTN.NE.O) WRITE (NOUT, 1009)
IFCICHLOR.NE.O) WRITE(NOUT.IOIO)
•«RITE(NGUT« 1112)
RETURN
IHEAVY
IHEAVY
FORMATC" THE FOLLOWING CONSTITUENTS ARE BEING MODELED 1)
FORMATC' T^E F-OLLOwING CONSTITUENTS ARF BEING MODELED DISSOLVED
*OXYGENl )
FOR"AT(1RX» 'CoLIFORMS')
KORMAT(U8X»'DODI)
FORMAT(W8X. 10A«)
FORMAT(48X, IPHYTOPLANKTONI)
FORMAT(U8XiIl,l HfAVY METAL (AND
FORMAT («3X? ' TOTAL NITROGEN')
FORMAT (U8X« i CHLORIDES I )
FORHATty8x.il.1 HEAVY METALS (AND
FORMATC//)
E NO
SUBROUTINE TRIBD
COMMON/bLOCKl/CwHIN(100).DATA(100
•
ITS ASSOCIATED ION) 1)
THEI-R ASSOCIATED IONS) ' )
."0)
COM"ON/BLOCK2/jNlT(6n»G(bO)«TlTLE(20)iIORD(20»20)
COMHON/BLOCK3/CONOZ(10.20).COi< 01(20,20)
COMMON /BLOCK5/JJ(6)iCLOW,NTRIbUO)(QUP,!-IHL(2).JA.IA
C
C
c
c
c
c
2
C
c
1 = 1
IFCIORD(JAiI))l»lt3
STEP-1
CHECK TO SEE IF THCRF. ARE ANY
DOWNSTREAM REACHES TO BE
CONSIDERED IN THE STRETCH,
STEPr2
139
-------
2013.
2011,
2015,
2016,
2017,
2018,
2019.
2020,
2021,
2022,
2023,
2021,
2025.
2026,
2027.
2028,
2029,
2030,
2031,
2032,
2033,
203U,
2035,
2036,
2037,
2038,
2039.
2010,
2011.
2012.
2013,
2011,
2015,
2016.
2017.
2018,
2019,
2050.
2051,
2052.
C
C
C
3
C
C
C
C
C
C
C-
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
1
IA=IOKD(JA.I)
CALL REPKF.
CALL DOEQU
OUP=QUP+OATA(IA,7)
IF(OUP) 10ilOi20
10 CRMlN(IA)=CONDI(JA, 1)
GO TO 30
20 CRMIM(IA)=CLOW
30 CONTINUE
1-1 + 1
GO TO 2
RETURN
END
CONSIDER MFXT DOWNSTREAM
REACH,
STH'-3
CALL R E P R E
STEP-1
CALL D 0 E n U
STfP"S
STORF MINIMUM DISSOLVED OXYGFN
OF REACH,
STEP-6
REPfAT STKP-1 THRU STEP-1)
UNTIL ALL REACHES HAVE HEEN
CONSIDERED,
STF.P-7
RETURN TO CALLER CR E I N I)
OR (R H A T C) ,
140
-------
SECTION X
REFERENCES
1. Churchill, M.A., Elmore, H.L., and Buckingham, R.A., "The Prediction •
of Stream Reaeration Rates," Journal, Sanitary Engineering Division,
A.S.C.E., vol. 7, 1962.
2. Eckenfelder, W.W. and O'Connor, D.J., Biological Waste Treatment,
Pergamon Press, London, 1962.
3. Langbien, W.B. and Durum, W.H., "The Aeration Capacity of Streams,"
U.S.G.S. Circular No. 542, U.S. Department of the Interior, Washing-
ton, D.C., 1967.
4. O'Connor, D.J. and Dobbins, W.E., "Mechanism of Reaeration in Natural
Streams," Transactions, A.S.C.E., vol. 123, New York, 1958.
5. O'Connor, D.J., "An Analysis of the Dissolved Oxygen Variation in a
Flowing Stream," paper presented at Symposium on Advances in Biologi-
cal Waste Treatment, University of Texas at Austin, Austin, Texas,
1966.
6. Owens, M., Edwards, R.W., and Gibbs, J.W., "Some Reaeration Studies
in Streams," Inter. Journal of Air and Water Pollution Res., vol. 8,
1964.
7. Streeter, H.W. and Phelps, E.B., "A Study of the Pollution and
Natural Purification of the Ohio River," Bulletin 146, U.S. Public
Health Services, 1925, reprinted 1958.
8. "DOSAG-I Simulation of Water Quality in Streams and Canals. Program
Documentation and Users Manual," Texas Water Development Board,
Austin, Texas, 1970.
141
-------
SECTION XI
ABBREVIATIONS
ASCE
EPA
SCI
USGS
BOD
Cent0
cf s
cms
deg
DO
°F
ft
g
HM
HM1
HM2
HM3
hr
JCL
m
mb
mg
min
mL
MPN
N
NH3 -N
NO -N
PO.-P
4
sec
American Society of Civil Engineers
Environmental Protection Agency
Systems Control, Incorporated
United States Geological Survey
biochemical oxygen demand (5-day)
chloride
centigrade degrees
cubic feet per second
cubic meters per second
degrees
dissolved oxygen
Farenheit degrees
feet
acceleration due to gravity
heavy metal
heavy metal one
heavy metal two
heavy metal three
hour
job control language
meters
millibars
milligrams
minutes
milliliter
most probable number
nitrogen
ammonia nitrogen
nitrite nitrogen
nitrate nitrogen
phosphate phosporus
seconds
143
-------
SELECTED WATER
RESOURCES ABSTRACTS
INPUT TRANSACTION FORM
w
Spokane River Basin Model Project
October, 1974
T '''it.-. - ,,•„" 's.
i' ,. 5
M. ' of
/-'-,.'»
,? r'-'c
Send To :
WATER RESOURCES SCIENTIFIC INFORMATION
US DEPARTMENT OF THE INTERIOR
WASHINGTON. D C. 2O24O
CENTER
, E. John Finnemore Systems Control, Inc.
-------
|