SPOKANE RIVER BASIN MODEL PROJECT
Volume III - Verification Report
by
E. John Finnemore, Ph.D.
John L. Shepherd
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 EPA and approved for publication.
Approval does not signify that the contents necessarily reflect the
views and policies of the Environmental Protection Agency, nor does
mention of trade names or commercial products constitute endorsement
or recommendation for use.
11
-------
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
Background 1
Objectives 1
Scope of Present Report 1
II. PROCEDURE 5
Approach 5
Model Modifications 5
III. RESULTS 7
DOSAG Verification 7
SWMM (RECEIV) Verification 19
DRM Verification 39
IV. DELIVERABLES 51
Appendix
A. MODIFICATIONS TO DOSAG 53
B. MODIFICATIONS TO SWMM (RECEIV) 57
C. MODIFICATIONS TO DRM 63
-------
FIGURES
NO. PAGE
1. Spokane River Basin 2
2. Layout of Spokane River System 3
3. DOSAG Verification for DO on River Region 2,
August 1969 11
4. DOSAG Verification for BOD on River Region 2,
August 1969 12
5. DOSAG Verification for DO on River Region 2,
September 1969 15
6. DOSAG Verification for BOD on River Region 2,
September 1969 16
7. RECEIV Verification for DO on River Region 2,
August 1969 26
8. RECEIV Verification for BOD on River Region 2,
August 1969 27
9. RECEIV Verification for Zinc on River Region 2,
August 1969 28
10. RECEIV Verification for DO on River Region 2,
September 1969 30
11. RECEIV Verification for BOD on River Region 2,
September 1969 31
12. RECEIV Verification for Zinc on River Region 2,
September 1969 32
13. DRM Verification for Temperature in Long Lake,
1971 (SURFACE) A2
14. DRM Verification for Temperature in Long Lake,
1971 (1 METER DEPTH) 43
15. DRM Verification for Temperature in Long Lake,
1971 (6 METER DEPTH) 4A
VI
-------
FIGURES (continued)
NO. PAGE
16. DKM Verification for Temperature in Long Lake,
1971 (14 METERS DEPTH) 45
17. DKM Verification for Temperature in Long Lake,
1971 (21 METERS DEPTH) 46
18. DRM Verification for Temperature in Long Lake,
1971 (24 METERS DEPTH) 47
19. DRM Verification for Temperature in Long Lake,
1971 (30 METERS DEPTH) 48
20. DRM Verification for Temperature in Long Lake,
1971 (OUTFLOW) 49
VII
-------
TABLES
NO. PAGE
1. DOSAG Verification, River Region 1 for
July-August 1971 8
2. DOSAG Verification, River Region 1 for
August-September 9
3. DOSAG Verification, River Region 2 for
Augus t 19 69 10
4. DOSAG Verification, River Region 2 for
September 1969 14
5. DOSAG Verification, River Region 3 for
August 1969 17
6. DOSAG Verification, River Region 3 for
September 1969 18
7. DOSAG Verification, River Region 4 for
July-August 1968 20
8. DOSAG Verification, River Region 4 for
August-September 1968 21
9. DOSAG Verification, River Region 5 for
August 1971 22
10. DOSAG Verification, River Region 5 for
September 1971 22
11. RECEIV Verification, River Region 1 for
August-September 1971 23
12. RECEIV Verification, River Region 2 for
August 1969 25
13. RECEIV Verification, River Region 2 for
September 1969 29
viii
-------
TABLES (continued)
NO. PAGE
14- RECEIV Verification, River Region 3 for
August 1969 33
15. RECEIV Verification, River Region 3 for
September 1969 34
16. RECEIV Verification, River Region 4 for
July-August 1968 36
17. RECEIV Verification, River Region 4 for
August-September 1968 37
18. RECEIV Verification, River Region 5 for
September 1971 38
19. DRM Verification, Long Lake 1971 40
20. DRM Verification, Coeur d'Alene Lake 1971 50
IX
-------
SECTION I
INTRODUCTION
BACKGROUND
The Spokane River Basin Model Project, initiated on August 31, 1972,
essentially consists of the modification and application of mathe-
matical models for water quality to rivers and lakes in the Spokane
River Basin in Washington and Idaho (see Figure 1).
The lakes- to be modeled naturally divide the river system into five
sub-areas. For the purposes of this project, these have been named River
Regions, and numbered on Figure 2.
Three mathematical models are included. Both DOSAG (Texas Water Develop-
ment Board) and the Receiving Water Module ("RECEIV") of the Storm Water
Management Model (EPA) are to be applied to each of the five above mentioned
river regions. The Deep Reservoir Model (EPA) is to be applied to Coeur
d'Alene Lake, Long Lake, and if possible the Spokane River Arm of F. D.
Roosevelt Lake.
The collection and assessment of the various data required by these models
was completed with the Project Officer's acceptance of the Phase I Data
Report. The Data Report also includes the selection and recommendation
of simulation periods to be used for each lake or river region.
OBJECTIVES
The objectives of this verification report are to:
• report on the capabilities of the previously existing models
• itemize the results of the verification simulations
• record modifications to the models (i) necessary before their
execution in Phase II, and (ii) desirable before their use in
Phase III.
SCOPE OF PRESENT REPORT
This report documents the procedures used in, and results obtained from,
verification executions of the three specified mathematical models simulat-
ing the prescribed rivers and lakes in the Spokane River Basin. It consists
of two major sections.
Section II including the Appendices, describes the approach used and modifi-
cations to the models which were needed or suggested.
Section III itemizes and discusses the results of the verification simula-
tions .
Since this report was a deliverable required during the earlier stages
of the project, it was written in May 1973.
-------
WASHINGTON j IDAHO
I'PSTP.W! L.'MT C.-
oi^:,£ RIVER :A;.'N
A?PnCX!,";>.7E \
STUDY EOIKARY.^ \
Scale
20 30 40 Xi'les
Figure 1. Spokane River Basin (portions to be modeled in bold).
-------
u>
FIGURE 2. UYOUT Of SPOKANE RIVER SYSTEM
-------
SECTION II
PROCEDURE
APPROACH
The Phase II version of DOSAG was verified by simulating BOD and DO
in the five river regions. The Phase II version of RECEIV (Receiving
Water Module of SWMM) was verified by simulating DO, BOD, and zinc
in the five river regions.
Zinc was modeled as conservative and coliforms were not modeled at
all, even though some observed coliform data .were available. This was
due to the shortage of coliform data and the fact that programming
changes would have been needed in RECEIV to provide a different conver-
sion factor and input and output format for coliforms on account of the
different units required to describe them. No TDS data were available;
and hence, TDS were not modeled.
For all runs the BOD decay coefficient was set at 0.2 DAY . For all
RECEIV runs, the reaeration coefficient was set at 1.0 DAY
The Phase II version of DRM was verified by simulating internal
isotherms and discharge temperatures in Long Lake and Coeur d'Alene
Lake.
Tables and brief explanations of the verification run results are included
in Section III.
MODEL MODIFICATIONS
Article I of the contract specifies that program modifications made
under Phase II shall be limited to those necessary to the Spokane River
Basin, such as changes in DIMENSIONS, etc. However, SCl's experience
with these models suggested that (a) they still contained basic errors,
(b) some of the options they offered had not been tested, and (c) some
of the types of applications made under this contract were new and un-
tested ones for these models. Furthermore, large sections of the DRM
code were non-general, and applicable to only one specific reservoir.
All the significant modifications (other than those obviously needed,
such as changes in DIMENSION, etc.) necessary to overcome the above
mentioned difficulties, and thus enable execution of the 'as-delivered'
model versions, are collected together under Class (i) modifications
for each of the three Appendices.
Class (ii) modifications for each of the Appendices contains suggested
further modifications desirable before the use of the models in Phase
III.
-------
SECTION III
RESULTS
DOSAG VERIFICATION
Tables, figures and brief discussions of the verification simulation
results follow under separate sections for each River Region as defined
in Figure 2.
In the tables, the units for concentrations are all given in units of
mg/L; RM represents River Mile.
River Region 1: St. Joe - St. Maries (DOSAG)
The only observed data available for this region during July 16 - August
15, 1971, consist of coliform data and temperature data. A run was made,
however, with the DO concentration in the headwaters set at 100% of satur-
ation and the input temperature for the region set at 20.1° (see Table
1). At an elevation of 2150 feet and a temperature of 20.1°, the saturation
level for DO is approximately 8.45 mg/L. Since the value of the reaeration
coefficient, which is calculated in DOSAG as a function of velocity and
depth, is such that DO levels remain near saturation, and since there
is no BOD to decrease the DO level, the DO concentration remains at or
near saturation throughout the region.
For the August 16 - September 16, 1971 period, DO data and temperature
data are available. The input temperature for this period is 20.5° and
the resultant saturation level of DO is approximately 8.35 mg/L (see Table
2) . The fact that the observed values of DO are higher indicates that
the regionwide input temperature is probably too high.
River Region 2: Coeur d'Alene (DOSAG)
Extensive observed DO, BOD and temperature data are available for both
August and September, 1969 for this region. In general, the observed
values of DO are higher throughout the region in September than in August,
while the observed values of BOD are lower in the South Fork and higher
in the main stem in September. Small sewers enter the river system at
three locations, and there is extensive leaching and seeping from various
mining operations along the South Fork of the Coeur d'Alene River.
For the August simulation, DO concentrations were set at approximately
8 mg/L throughout the system and BOD concentrations were set at approxi-
mately 1 mg/L. BOD and DO concentrations in the sewers were set at 200
mg/L and 1.0 mg/L respectively. The temperature for August was input
at 19 degrees. As can be seen from Table 3 and Figures 3 and 4, the DO
stays near saturation at about 8.5 mg/L while the BOD stays near 1 mg/L
except in the vicinity of the outfalls.
-------
TABLE 1.
DOSAG VERIFICATION (PHASE II)
RIVER REGION 1 (ST. JOE - ST. MARIES)
JULY 16 - AUGUST 15, 1971
w
o
1-)
H
CO
H
CO
DO
*
RM OBS MOD
42.9 8.44
33.5 8.44
31.1 8.44
24.5 8.44
15.7 8.44
15.4 8.44
15.0 8.44
10.0 8.44
0.7 8.44
27.8 8.44
14.8 8.44
12.0 8.44
10.0 8.44
3.9 8.44
River Mile
-------
TABLE 2
DOSAG VERIFICATION (PHASE II)
RIVER REGION 1 (ST. JOE - ST. MARIES)
AUGUST 16 - SEPTEMBER 16, 1971
w
o
>->
H
CO
i
H
CO
RM
42.9
33.5
31.1
24.5
15.7
15.4
15.0
10.0
0.7
27.8
14.8
12.0
10.0
3.9
DO
OBS
9.5
9.2
9.7
9.0
9.6
9.4
9.4
10.0
9.5
MOD
8.37
8.37
8.37
8.37
8.37
8.37
8.37
8.37
8.37
8.37
8.37
8.37
8.37
8.37
-------
TABLE 3
DOSAG VERIFICATION (PHASE II)
RIVER REGION 2 (COEUR D'ALENE)
AUGUST 1969
erf
w
w
w
a
w
o
u
g
o
CO
RM
138
145
148
154
160
166
168
170
177
182
188
194
1
4
6
7
10
12
15
17
18
22
DO
OBS
7.5
7.6
7.6
8.3
8.3
8.5
8.6
7.6
8.1
8.3
8.0
8.2
8.2
7.9
7.9
8.2
8.0
8.2
8.7
7.8
7.9
8.0
MOD
8.40
8.42
8.46
8.51
8.54
8.53
8.54
8.53
8.52
8.56
8.56
8.50
8.49
8.49
8.47
8.48
8.48
8.53
8.56
8.56
8.56
8.58
BOD
OBS
1.3
1.0
0.85
1.3
1.0
1.0
0.9
1.0
1.0
1.0
0.8
1.7
3.1
0.9
1.4
3.7
1.7
2.9
1.7
1.75
1.3
MOD
1.30
1.29
1.26
1.20
1.10
1.18
0.786
0.792
0.762
0.741
0.762
0.741
2.27
2.34
2.39
1.85
1.89
2.44
0.883
0.890
0.894
0.915
10
-------
1U
9
— i
C3
o 8
1 —
CcL
I—
UJ -.
o 7
o
C_J>
o
Q
6
t;
1 i i i i i i i i i i i i
0 0 0 o 0 s^e 0 Q 0 0
V V V
v v
v v v v
-
-£
o
« Modeled g
VObserved oo
i , T i , i i i i
iii iii
V
0 000 0° 00° °
V V V
^ V vv V
-
.
l/l C_3
n3 ro o
4-> -4-> c
^ Z3 ^3
0 O 0
, 1 , tt, T , ,
O o Ul
MAIN STEM
SOUTH FORK
RIVER MILE
FIGURE 3. DOSAG VERIFICATION FOR DO ON RIVER REGION 2, AUGUST 1969
-------
i
A
i
C3
n
i — i
| —
UJ 9
O *•
O
O
O
ca i
n
i I 1 1 1 I 1 I i 1 i i
-
O
u_
_c:
0
0 0
V $ ^ _ V V V
0 Modeled
V Observed
i , , T
MAIN STEM
„ o
v
4-> E
rs ro
oo
V
V
. i ,
v
00° 0
o ui
SOUTH FORK
RIVER MILE
FIGURE 4. DOSAG VERIFICATION FOR BOD ON RIVER REGION 2, AUGUST 1969
-------
For the September simulation, DO concentrations were set at approximately
9.5 mg/L throughout the system and BOD concentrations were set at approxi-
mately 1.5 mg/L. The BOD concentration in the sewers was set at 75 mg/L.
The temperature for September was input at 17 degrees. The resulting
values match the observed values very well (see Table 4 and Figures 5
and 6) with the exception of DO in the upper reaches of the South Fork.
This descrepancy is probably due to the fact that the temperature of these
upper reaches was less than 17 degrees; and, thus, the DO concentration
level there was greater than 9 mg/L.
River Region 3: Upper Spokane (DOSAG)
Sparse DO and temperature data are available for both August and September
1969 for this region. Small diversions from the river occur at two
locations, and there is considerable groundwater flow into the system
throughout the region.
For August the region temperature was set at 21.1 degrees and the DO level
in the main inflow (below Post Falls Dam) was set at 85% of saturation.
DO in the groundwater inflow was set at 7 mg/L. The resulting DO concen-
tration (see Table 5) agrees with the observed values in the upper part
of the Spokane River, but is too high in the downstream part. This is
perhaps caused by the large groundwater inflow which may have a DO concen-
tration less than the 7 mg/L which was used. Reactions which consume
DO, such as BOD decay, are also taking place, although no observed data
is available. Also, the value of the reaeration coefficient, which is
calculated as a function of depth and velocity, may not be suitable for
this area of the region.
For September the region temperature was set at 17 degrees and the DO
level in the main inflow was set at 90% of saturation. The DO level in
the groundwater was left at 7 mg/L. Results similar to the August results
were obtained (see Table 6).
River Region 4: Little Spokane (DOSAG)
Observed DO, BOD, and temperature data are available for July 11 - August
10, 1968 for this region. Observed DO and temperature data are available
for August 11 -September 10, 1968. There is considerable groundwater
flow into the region.
For the July 11 -August 10 run, the input temperature was 14.6 degrees.
The DO level in the headwater inflow was set at 94% of saturation. The
DO level of Deadman Creek (which has a flow approximately equal to that
of the Little Spokane River above their junction) was set to 9 mg/L.
The DO concentration of the groundwater (which triples the Little Spokane
13
-------
TABLE 4
DOSAG VERIFICATION (PHASE II)
RIVER REGION 2 (COEUR D'ALENE)
SEPTEMBER 1969
Pi
w
1 1
Pi
w
w
Q
Pi
W
O
g
O
C/3
RM
138
145
148
154
160
166
168
170
177
182
188
194
1
4
6
7
10
12
15
17
18
22
DO
OBS
9.0
9.5
9.3
8.9
9.2
8.8
8.8
9.3
8.9
9.1
8.6
8.7
9.3
10.2
7.9
10.3
10.2
10.0
9.8
9.5
MOD
8.78
8.80
8.83
8.88
8.91
8.91
8.94
8.94
8.94
8.93
8.94
8.89
8.92
8.94
8.97
8.99
8.99
8.94
8.96
8.96
8.97
8.95
BOD
OBS
1.6
1.6
1.4
1.2
1.2
0.6
0.5
1.7
1.0
1.0
1.1
1.6
1.3
1.8
1.9
1.7
1.5
1.7
1.4
1.3
MOD
1.38
1.35
1.29
1.21
1.07
1.15
0.914
0.921
0.762
0.741
0.759
0.743
1.88
1.91
1.95
1.71
1.73
1.92
1.25
1.25
1.26
1.38
14
-------
1 1
10
1
E 9
t— (
<:
i—
UJ O
o
O
o
Q 7
6
i i i i i i i i i i i i
V
*
^ 0 0 V © 0 00^0
0 V V
-
o
u_
.c:
©Modeled ^
o
V Observed ^
i i i i i T i i i i i I I
1*3 -f -Pen LnOlO^-vJ «^J CDCOuD ijO C
Cno cno cnouio ui ouica en c
MAIN STEM
3
3
1 1 1 I ' '
w V
V
V V
0 0 00 o 0 000 V
V
V
S-
00 O
^ •— c
ro (13 O ~
000
, * , iV, T, , ,
O LH C3 Ul CD O
SOUTH FORK
RIVER MILE
FIGURE 5. DOSAG VERIFICATION FOR DO ON RIVER REGION 2, SEPTEMBER 1969
-------
D
_j 4
(J5
^:
ri
0 3
i — § -1
i—
=a;
a;
h-
z:
uu
k? 2
0
o
o
0
m 1
0
1 1 1 1 1 1 I 1 1 1 II
-
o
u_
-C
4->
V V J *
6 S « Y ¥
© V V
0 Modeled Q 0 0 0
V Observed
i i i i i i T i i i i i i
(
^ ^ ** t~* ^ *-• N)
gggS gg§C
MAIN STFM
1 1 1 r^ r i
t_5
<— i — C
ro ro O ~
M- <4- >,
^ 3 fO
O O CJ
, 0 © V ©
v ¥ «v vv
„ V 001 ^
V
i t i tti Ti i i
> - S K g s s ,
SOUTH FORK
RIVER MILE
FIGURE 6. DOSAG VERIFICATION FOR BOD ON RIVER REGION 2, SEPTEMBER 1969
-------
TABLE 5
DOSAG VERIFICATION (PHASE II)
RIVER REGION 3 (UPPER SPOKANE)
AUGUST 1969
Pi
w
>
1 — 1
&
w
• o
p-l
en
pi
O
O
^
DO
RM OBS MOD
110.7 7.18
106.6 7.68
102.1 8.00
101.8 8.2 8.03
98.7 8.1 8.13
96.4 8.19
93.9 8.0 8.22
88.7 8.11
84.8 8.04
80.2 8.04
77.9 8.03
76.2 8.02
74.2 8.01
72.9 8.02
72.4 7.93
64.2 7.97
58.1 8.02
56.7 7.4 8.04
39.0 7.17
32.9 8.33
20.2 8.36
14.5 8.36
0.8 8.36
17
-------
TABLE 6
DOSAG VERIFICATION (PHASE II)
RIVER REGION 3 (UPPER SPOKANE)
SEPTEMBER 1969
pi
w
1— 1
w
o
P-4
CO
g
|
0
1
DO
RM OBS MOD
110.7 8.22
106.6 8.49
102.1 8.71
101.8 8.8 8.73
98.7 8.8 8.81
96.4 8.87
93.9 8.90
88.7 8.86
84.8 8.78
80.2 8.77
77.9 8.76
76.2 8.73
74.2 8.73
72.9 8.72
72.4 8.62
64.2 8.64
58.1 8.70
56.7 7.7 8.71
39.0 8.43
32.9 9.07
20.2 9.09
14.5 9.09
0.8 9.09
18
-------
flow between RM11 and RM4) was set at 7 mg/L. The headwater BOD level
was set at 0.8 mg/L. All other BOD concentrations were set at 0.6 mg/L.
As can be seen from Table 7, the DO concentration is driven toward satur-
ation by reaeration and diluted toward 7 mg/L by the groundwater flow.
The BOD values agree fairly well with the observed values.
For the August 11 - September 10 run, the input temperature was 14.1 de-
grees, all BOD concentrations were zero, and the DO level of Deadman Creek
was set at 9.5 mg/L. Groundwater DO remained at 7 mg/L. A DO profile
similar to the July 11 - August 10 run was obtained (see Table 8).
River Region 5; Lower Spokane (DOSAG)
Only one observed value of DO is available for this region during each
of the months of August and September, 1971. No BOD data are available.
For August, the input temperature was set at 20.1 degrees and the DO level
of the inflow was set at 54% of saturation. Simulation results are pre-
sented in Table 9.
For September, the input temperature was set at 17.1 degrees and the DO
level of the inflow was set at 54% of saturation. Simulation results
are presented in Table 10.
SWMM (RECEIV) VERIFICATION
Tables, figures and brief discussions of the verification simulation results
follow under separate sections for each River Region as defined in Figure
2.
River Region 1: St. Joe - St. Maries (RECEIV)
The only observed data available for this region during July 16 - August
15, 1971 consist of coliform data and temperature data. Since RECEIV
simulates neither of these quantities (see Section II), the region was not
modeled for this period.
For the August 16 - September 16, 1971 period, DO data and temperature
data are available, and for this period the DO-BOD constituent was simu-
lated, with all BOD levels set to zero. The reaeration coefficient K_
was set to keep DO levels near saturation and DO saturation was set at
10 mg/L. The initial concentrations of the various junctions were set
at values close to the observed data, and DO concentrations of junction
inflows were set accordingly. Simulation results are presented in Table
11.
19
-------
TABLE 7
DOSAG VERIFICATION (PHASE II)
RIVER REGION 4 (LITTLE SPOKANE)
JULY 11 - AUGUST 10, 1968
RM
37.6
34.6
32.9
31.0
21.3
13.5
13.1
11.4
10.8
7.9
3.9
0.1
DO
OBS
9.7
10.4
9.9
9.4
8.8
8.6
8.7
9.8
MOD
9.14
9.49
9.51
9.51
9.51
9.51
9.23
9.40
8.30
8.53
9.24
9.42
BOD
OBS
0.6
0.7
0.9
0.6
0.6
0.6
0.7
0.4
MOD
0.80
0.79
0.78
0.77
0.73
0.70
0.64
0.63
0.61
0.61
0.60
0.59
20
-------
TABLE 8
DOSAG VERIFICATION (PHASE II)
RIVER REGION 4 (LITTLE SPOKANE)
AUGUST 11 - SEPTEMBER 10, 1968
RM
37.6
34.6
32.9
31.0
21.3
13.5
13.1
11.4
10.8
7.9
3.9
0.1
DO
OBS
9.8
10.6
9.4
9.6
9.8
8.2
8.6
8.6
MOD
9
9
9
9
9
9
9
9
8
8
9
9
.14
.54
.59
.62
.62
.62
.55
.58
.41
.49
.16
.40
21
-------
TABLE 9
DOSAG VERIFICATION (PHASE II)
RIVER REGION 5 (LOWER SPOKANE)
AUGUST 1971
DO
RM
OBS
MOD
33.9
31.8
28.2
24.0
4.6
4.61
4.79
4.98
5.36
TABLE 10
DOSAG VERIFICATION (PHASE II)
RIVER REGION 5 (LOWER SPOKANE)
SEPTEMBER 1971
DO
RM
OBS
MOD
33.9
31.8
28.2
24.0
4.8
4.99
5.16
5.32
5.64
22
-------
TABLE 11
RECEIV VERIFICATION (PHASE II)
RIVER REGION 1 (ST. JOE - ST. MARIES)
AUGUST 16 - SEPTEMBER 16, 1971
w
o
H
3
i
H
CO
RM
42.9
33.5
31.1
24.5
15.7
15.4
15.0
10.0
0.7
27.8
14.8
12.0
10.0
3.9
DO
OBS
9.5
9.2
9.7
9.0
9.6
9.4
9.4
10.0
9.5
MOD
9.58
9.70
9.70
9.77
9.84
9.84
9.84
9.86
9.87
10.0
10.0
9.95
9.95
9.94
23
-------
River Region 2; Coeur d'Alene (RECEIV)
Extensive observed data for DO, BOD, zinc and temperature are available
for this region for both August and September, 1969.
Initial conditions similar to those used by DOSAG (see DOSAG VERIFICA-
TION) with regard to DO and BOD were input and the saturation value of
DO was set at 8.3 mg/L for the August simulation. In addition, zinc levels
of approximately 2 mg/L were input for Canyon Creek and the zinc concentra-
tion of the inflow in Junction 11 (representing the discharge of a zinc
processing plant) was set at 125 mg/L. Good agreement with the observed
data was obtained (see Table 12 and Figures 7 through 9).
For September, the DO saturation level was set at 9.3 mg/L, the zinc plant
inflow concentration for zinc was set at 250 mg/L, and the sewage BOD levels
were reduced (as with DOSAG). Good agreement with the observed data was
again obtained (see Table 13 and Figures 10 through 12).
River Region 3: Upper Spokane (RECEIV)
Sparce observed DO, zinc, and temperature data are available for August
and September, 1969, for this region. Small diversions from the river
occur at two locations and there is considerable groundwater flow into the
system throughout the region.
For August, the DO saturation level was set at 8.2 mg/L and the DO level
of the groundwater inflow was set at 7.0 mg/L. The zinc level was set at
0.15 mg/L in the main inflow (below Post Falls Dam). The results agree.
with the observed data reasonably well (see Table 14), except for DO in
the lower reaches of the Spokane River where the simulated value is higher
than the observed data. This may be due to unmodeled DO consumption, such
as BOD decay, or the input reaeration coefficient may not be realistic for
this area of the region.
For September, the DO saturation level was set at 9.0 mg/L and the zinc
level in the main inflow was set at 0.19 mg/L. The DO concentration of the
groundwater was left at 7.0 mg/L. Results similar to the August results
were obtained (see Table 15).
River Region 4: Little Spokane (RECEIV)
Observed DO, BOD, and temperature data are available for July 11 - August
10, 1968, for this region. Observed DO and temperature data are available
for August 11 - September 10, 1968. There is considerable groundwater
flow into the region.
For the July 11 - August 10 run, the DO saturation level was set at 9.8
mg/L. The DO level in the headwater inflow and in Deadman Creek was set
at 9.8 mg/L. The DO level in the groundwater inflow was set at 7.5 mg/L.
24
-------
TABLE 12
RECEIV VERIFICATION (PHASE II)
RIVER REGION 2 (COEUR D'ALENE)
AUGUST 1969
<&
w
M
W
W
n
COEUR
^
0
CD
RM
138
145
148
154
160
166
168
170
177
182
188
194
1
4
6
7
10
12
15
17
18
22
DO
OBS
7.5
7.6
7.6
8.3
8.3
8.5
8.6
7.6
8.1
8.3
8.0
8.2
8.2
7.9
7.9
8.2
8.0
8.2
8.7
7.8
7.9
8.0
MOD
8.10
8.11
8.11
8.13
8.16
8.23
8.26
8.28
8.30
8.30
8.30
8.30
8.26
8.30
8.30
8.30
8.27
8.30
8.26
8.26
8.30
8.30
BOD
OBS
1.3
1.0
0.85
1.3
1.0
1.0
0.9
1.0
1.0
1.0
0.8
1.7
3.1
0.9
1.4
3.7
1.7
2.9
1.7
1.75
1.3
MOD
1.15
1.30
1.30
1.45
1.39
1.16
1.21
0.81
0.83
0.85
0.91
0.93
1.21
2.46
2.54
2.54
2.69
2.46
3.04
3.04
1.41
1.43
ZINC
OBS
2.6
1.8
1.65
4.7
4.7
1.9
0
0
0
0
0
0
5.1
6.65
0.5
0.3
0.3
0.6
0.65
0.8
1.0
0
MOD
2.03
1.96
1.96
1.91
1.90
1.89
0
0
0
0
0
0
1.88
6.64
0.59
0.59
0.66
0.66
1.01
1.01
1.05
0
25
-------
10
< 9
CD
_,T
1 8
o:
i—
LU
2: 7
o
O
0
o
6
S
V
V V 000 ©"VQ 0 ®-
- 0 00 ^
v v v v
-
-^
o
u_
©Modeled "3
o
v Observed °°
i il ii i V i i i ] ii
MAIN STEM
I I I I
RIVER MILE
0 ©£, 0 ^ 0 00 0
V ™ V
V
4-
4->
3
O
re
i+-
-i->
3
O
c
03
>— l-l N> tJ LU UJ
tn o ui
SOUTH FORK
FIGURE 7. RECEIV VERIFICATION FOR DO ON RIVER REGION 2, AUGUST 1969
-------
0
4
— i
CJ3
SI
o 3
i — i
h-
i—
uj 2
0
o
0 -,
CO 1
n
i i i i i i i i i i i i
-
s_
0
U-
g
v v ' ¥ 0 0 0 f
0 Modeled
v Observed
ii 1 i i i V l ill ii
r- V " °
ro ro O
^- H- >,
W -M -t->0 C
V ^ ZJ ro
o oV t-3
0
® 00
e 0
-
V V rf&
TO 00 ©
V
i T i Iti ^i i i
o tn o
O l/l
MAIN STEM
RIVER MILE
SOUTH FORK
FIGURE 8. RECEIV VERIFICATION FOR BOD ON RIVER REGION 2, AUGUST 1969
-------
ro
CO
V V
-F
00 0
V
w
»
® Modeled
^7 Observed
i.
o
o
°°
MATM CTCM
MAIN STEM
*•
V
CO
O-i—
04-
r-^lO
c
0
>,
ro
^°©f
o w
V w
RIVER MILE
FIGURE 9. RECEIV VERIFICATION FOR ZINC ON RIVER REGION 2, AUGUST 1969
-------
TABLE 13
RECEIV VERIFICATION (PHASE II)
RIVER REGION 2 (COEUR D'ALENE)
SEPTEMBER 1969
RIVER
ALENE
«
COEUR
%
o
fn
C/2
RM
138
145
148
154
160
166
168
170
177
182
188
194
1
4
6
7
10
12
15
17
18
22
DO
OBS
9.0
9.5
9.3
8.9
9.2
8.8
8.8
9.3
8.9
9.1
8.6
8.7
9.3
10.2
7.9
10.3
10.2
10.0
9.8
9.5
MOD
9.01
9.04
9.04
9.09
9.12
9.17
9.19
9.19
9.20
9.20
9.23
9.24
9.19
9.26
9.30
9.30
9.28
9.30
9.24
9.24
9.28
9.30
BOD
OBS
1.6
1.6
1.4
1.2
1.2
0.6
0.5
1.7
1.0
1.0
1.1
1.6
1.3
1.8
1.9
1.7
1.5
1.7
1.4
1.3
MOD
1.57
1.71
1.71
1.36
1.31
1.07
0.99
0.79
0.82
0.84
0.90
0.92
0.99
1.65
1.69
1.69
1.73
1.63
1.75
1.75
1.41
1.42
ZINC
OBS
2.6
4.1
3.1
3.0
6.0
4.7
0
0
0
0
0
15.0
14.0
2.0
2.0
1.2
1.5
1.8
2.8
1.5
0.1
MOD
4.06
4.00
4.00
3.97
3.97
3.94
3.94
0
0
0
0
0
3.94
14.50
2.07
2.07
2.08
2.08
2.12
2.12
2.13
0
29
-------
CD
o 9
t— 1
1—
o;
i—
O O
2:
o
C_J
0
7
fi
i i i i i i i i i i i i
V
^ w
© n ^ 0 © G) 0 _G^
V V ^
-
_^
o
-C
0 Modeled o
v Observed °°
i i i i i i
V V
v® 00@ 0^ 0
0 '
V ^
V
i-
c_>
T3 f^3 O
M-j M- >,
O O O
i t , Mi Ti , i
f° f° u> o»
MAIN STEM
SOUTH FORK
RIVER MILE
FIGURE 10. RECEIV VERIFICATION FOR DO ON RIVER REGION 2, SEPTEMBER 1969
-------
b
_j 4
0
ft
° 1
i — i 0
1—
OL
I
1
UJ
C ") ••)
O
O
O
CO
n
i i ' i i i i i i i it
^
o
u_
5
^ $ 1 *
^Si
^ V
® V V
0 © @
® Modeled v
V
V Observed
i i i i i i
s_
(/) CJ
03 tO O
-4-J 4-> C
3 ^ m
o o c_>
_ V
\j7 fa
V ® v
V
i t , H, T,
MAIN STEM
RIVER MILE
SOUTH FORK
FIGURE 11. RECEIV VERIFICATION FOR BOD ON RIVER REGION 2, SEPTEMBER 1969
-------
U)
K)
15
14
13
12
o 9
i—
•* 8
a:
o 6
o
o 5
^
i — i
M 4
3
2
1
0
i i i i i i i i i i i i
-
V
-
V
V
V v -*
i_
V o
u_
® Modeled ^
o
^ Observed °°
1 1 1 1 1 iT^B 1 » 1 rp 1 «p 1 -J
^^
5 ®' ' ' ' ' '
V
-
-
-
D-
i.
U 0
c
hJ O
to >>
i — i — fO ~
<+- ^t-
o o ^
^5^> 00 rs ®s _
V
vv v
i
T 1 i Hi Ti S7 i i
S
tnocnocn
MAIN STEM
<;nilTH
ouum
RIVER MILE
FIGURE 12. RECEIV VERIFICATION FOR ZINC ON RIVER REGION 2, SEPTEMBER 1969
-------
TABLE 14
RECEIV VERIFICATION (PHASE II)
RIVER REGION 3 (UPPER SPOKANE)
AUGUST 1969
w
s
w
o
P-I
c/3
pi
u
o
*
DO
RM OBS MOD
110.7 8.2
106.6 8.2
102.1 8.2
101.8 8.2 8.2
98.7 8.1 8.2
96.4 8.2
93.9 8.0 8.2
88.7 7.98
84.8 7.99
80.2 8.03
77.9 8.03
76.2 8.03
74.2 7.96
72.9 7.98
72.4 7.98
64.2 7.84
58.1 7.96
56.7 7.4 7.96
39.0 7.33
32.9 7.58
20.2 7.69
14.5 7-78
0.8 7.98
ZINC
OBS MOD
0.150
0.150
0.150
0.150 0.150
0.145 0.150
0.150
0.200 0.150
0.115
0.111
0.104
0.104
0.091
0.091
0.091
0.091
0.084
0.083
0.080 0.083
0
0
0
0
0
33
-------
TABLE 15
RECEIV VERIFICATION (PHASE II)
RIVER REGION 3 (UPPER SPOKANE)
SEPTEMBER 1969
ft
W
>
&
W
0
PM
en
HANGMAN CR.
DO
RM OBS MOD
110.7 8.87
106.6 8.87
102.1 8.93
101.8 8.8 8.93
98.7 8.8 8.94
96.4 8.95
93.9 8.95
88.7 8.71
84.8 8.73
80.2 8.80
77.9 8.80
76.2 8.80
74.2 8.77
72.9 8.79
72.4 8.79
64.2 8.63
58.1 8.63
56.7 7.7 8.74
39.0 7.53
32.9 7.96
20.2 8.14
14.5 8.30
0.8 8.42
ZINC
OBS MOD
0.190
0.190
0.190
0.190 0.190
0.190 0.190
0.190
0.190
0.159
0.158
0.154
0.154
0.146
0.146
0.145
0.145
0.137
0.137
0.140 0.137
0
0
0
0
0
34
-------
Input BOD concentrations were approximately 0.7 mg/L except for the ground-
water BOD which was input as 0.55 mg/L. As can be seen from Table 16,
the DO level stays near saturation until the groundwater flow enters
the river (between RMll and RM4) and reduces it to approximately 8.2
mg/L. A better match to the observed data could probably be obtained
by increasing the DO level in the groundwater or by increasing the reaer-
ation coefficient.
For the August 11 - September 10 run, the DO saturation level was set
at 9.5 mg/L. All BOD levels were set to zero. Other concentrations
remained the same. A DO profile similar to the July 11 - August 10 run
was obtained (see Table 17).
River Region 5: Lower Spokane (RECEIV)
For each of the months of August and September, 1971, one observed value
of DO and one observed value of zinc are all that is available.
For September, the DO saturation level was set at 7.5 mg/L and the DO
concentration of the inflow was set at 4.6 mg/L. All zinc concentrations
were set at .04 mg/L. Simulation results are reported in Table 18.
Because of the great similarity of the August and September data, an August
run was not made.
35
-------
TABLE 16
RECEIV VERIFICATION (PHASE II)
RIVER REGION 4 (LITTLE SPOKANE)
JULY 11 - AUGUST 10, 1968
RM
37.6
34.6
32.9
31.0
21.3
13.5
13.1
11.4
10.8
7.9
3.9
0.1
DO
OBS
9.7
10.4
9.9
9.4
8.8
8.6
8.7
9.8
MOD
9.8
9.8
9.8
9.8
9.7
9.8
9.8
8.3
8.3
8.1
8.2
8.3
BOD
OBS
0.6
0.7
0.9
0.6
0.6
0.6
0.7
0.4
MOD
0.69
0.67
0.67
0.65
0.60
0.73
0.73
0.60
0.60
0.58
0.58
0.56
36
-------
TABLE 17
RECEIV VERIFICATION (PHASE II)
RIVER REGION 4 (LITTLE SPOKANE)
AUGUST 11 - SEPTEMBER 10, 1968
RM
37.6
34.6
32.9
31.0
21.3
13.5
13.1
11.4
10.8
7.9
3.9
0.1
DO
OBS
9.
10.
9.
9.
9.
8.
8.
8.
8
6
4
6
8
2
6
6
MOD
9.5
9.5
9.5'
9.5
9.5
9.4
9.4
8.3
8.3
8.1
8.2
8.3
37
-------
TABLE 18
RECEIV VERIFICATION (PHASE II)
RIVER REGION 5 (LOWER SPOKANE)
SEPTEMBER 1971
DO
ZINC
RM
OBS
MOD
OBS
MOD
33.9
31.8
28.2
24.0
4.8
4.78
5.09
5.53
5.84
0.04
0.04
0.04
0.04
0.04
38
-------
DRM VERIFICATION
This section presents brief discussions with tables and figures, where
appropriate, of the verification simulation results obtained from appli-
cations of the Deep Reservoir Model to both Long Lake and Coeur d'Alene
Lake (see Figures 1 and 2).
Long Lake (DRM)
Considerable temperature data are available for Long Lake for June - Novem-
ber, 1971. The lake inflow rates and inflow temperatures were estimated
using measurements from USGS gaging stations 4225, 4260 and 4310. Outflow
rates and outflow temperatures were determined from USGS gaging station
4330. Meteorologic data from the City of Spokane were used. The lake
was modeled as an inverted trapezoid 22 miles long, 32 meters deep, and
with the bottom width equal to .01 times the top width. This results in
a full volume of 3.19 x 10 cubic meters and a maximum surface area of
1.98 x 10' square meters Based on available data, all of these values
appear reasonable. The lake was assumed to be completely mixed at a
temperature of 12 degrees on June 1.
The simulation indicated that Long Lake "turned over" during the first
week of September, 1971. Simulation results are compared with observed
data in Table 19 and Figures 13 through 20; the comparison in general
is very good.
Coeur d'Alene Lake (DRM)
Surface temperature data are available on only four days during June
- November, 1971 for Coeur d'Alene Lake. No outflow or inflow tempera-
tures are available. The lake inflow rates were estimated using measure-
ments from USGS gaging station 4135, 4145, and 4149. Outflow was
determined from USGS gaging station 4190. Inflow temperatures were
estimated from surface temperature data. Meteorologic data from the
city of Spokane were used. The lake was modeled as being 22 miles long
and 25 meters deep, with a full volume of 1.5 x 10" cubic meters and
a maximum surface area of 1.77 x 108 square meters. The full volume
value was roughly estimated from available data for the top ten feet
only. The lake was assumed to be completely mixed at a temperature of
14 degrees on June 1.
The simulation indicates that Coeur d'Alene Lake "turned over" about
September 15th. Simulation results are compared with the observed data
in Table 20.
39
-------
TABLE 19
DRM VERIFICATION (PHASE II)
LONG LAKE, 1971
Date
June 1
June 13
June 22
July 8
July 11
July 25
July 27
August 1
August 10
August 15
August 18
August 23
Depth (m)
0
0
0
0
0
0
1
12
15
21
0
1
6
12
14
24
0
1
6
12
24
1
6
12
14
21
Lake Temp.
OBS
12.0
13.7
17.1
14.8
15.8
20.2
23.6
17.0
16.7
15.5
21.0
25.8
20.9
19.8
19.0
16.0
19.5
23.4
21.0
19.5
16.0
21.4
20.4
19.0
18.5
15.5
(°C)
MOD
12.2
14.1
18.9
16.3
16.1
21.2
21.7
17.2
16.5
16.1
22.6
22.6
21.6
19.1
18.3
16.7
21.0
20.0
20.0
19.3
17.1
19.2
19.2
19.2
18.8
17.6
Outflow Temp. (°C)
OBS MOD
12.9 12.0
14.1 13.5
16.0
16.0
16.4 16.0
19.4 18.8
19.2
20.2 19.9
20.5
20.0 20.3
19.9
19.2
(continued)
40
-------
TABLE 19 (CONTINUED)
Date Depth (m)
August 24
September 1
September 6
September 15
September 19
September 21
September 22
October 5
October 19
October 27
November 3
0
1
6
12
14
15
21
30
0
1
6
12
15
21
30
0
1
12
14
15
6
0
0
0
0
Lake Temp.
OBS
20.6
18.9
18.8
19.3
18.0
19.5
17.2
16.5
16.8
17.6
17.3
16.8
17.0
16.5
16.0
14.4
16.7
13.5
17.0
15.4
14.5
14.1
10.0
8.6
7.5
(°C)
MOD
19.3
18.9
18.9
18.9
18.8
18.7
18.1
17.9
17.9
16.8
16.8
16.8
16.8
16.8
16.8
16.2
15.7
15.8
15.8
15.8
15.8
14.1
11.3
9.2
7.4
Outflow Temp. (°C)
OBS MOD
19.2
19.1
17.5 18.0
16.9
16.7 16.2
15.8
15.8
14.0
11.5
9.4
7.5
41
-------
30
28
26
24
. 22
o>
-a
. 18
16
a:
14
12
10
V
0
0 Modeled
v Observed
V
V
0
V
0
v
©
0
V
JUNE
JULY
AUGUST
SEPTEMBER OCTOBER
NOVEMBER
FIGURE 13. DRM VERIFICATION FOR TEMPERATURE IN LONG LAKE, 1971 (SURFACE)
-------
28
.26
o
^24
QJ
"22
UJ
P20
18
16
14
V
v
0 Modeled
v Observed
JUNE
JULY
AUGUST
SEPTEMBER
OCTOBER
NOVEMBER
FIGURE 14. DRM VERIFICATION FOR TEMPERATURE IN LONG LAKE, 1971 (1 METER DEPTH)
-------
301-
28
26
24
22
20
16
14
UJ , 0
D- 12
10
8
6
4
2
©
V
V
©
® Modeled
v Observed
JUNE
JULY
AUGUST
SEPTEMBER
OCTOBER
NOVEMBER
FIGURE 15. DRM VERIFICATION FOR TEMPERATURE IN LONG LAKE, 1971 (6 METER DEPTH)
-------
en
O)
UJ
Q-
L"t
22
20
18
16
14
19
-
. v
V V 0 ®
V V
0 0
V
0
© Modeled
v Observed
i I 1 1 1
JUNE
JULY
AUGUST
SEPTEMBER
OCTOBER
NOVEMBER
FIGURE 16. DRM VERIFICATION FOR TEMPERATURE IN LONG LAKE, 1971 (14 METERS DEPTH)
-------
o
c^20
OJ
-o
, 18
UJ
1 16
1 14
1 IP
-
V
V a ^
^ ®
® Modeled
^ Observed
JUNE JULY AUGUST SEPTEMBER OCTOBER NOVEMBER
FIGURE 17. DRM VERIFICATION FOR TEMPERATURE IN LONG LAKE, 1971 (21 METERS DEPTH)
-------
26
24
22
£20
2 16
LU
Q-
F, 14
12
10
V
0
® Modeled
V Observed
JUNE
JULY AUGUST SEPTEMBER OCTOBER NOVEMBER
FIGURE 18. DRM VERIFICATION FOR TEMPERATURE IN LONG LAKE, 1971 (24 METERS DEPTH)
-------
00
26
24
"22
|20
£ 18
§16
o;
UJ
10
Modeled
Observed
JUNE
JULY
v
9
AUGUST
SEPTEMBER OCTOBER
NOVEMBER
FIGURE 19. DRM VERIFICATION FOR TEMPERATURE IN LONG LAKE, 1971 (30 METERS DEPTH)
-------
£0
24
o 22
^20
-a
uj" 18
a:
I16
UJ
P- i n
UJ
H- (
12'
1 n
_
-
A ^
— V7 ©
® V
V
V
0
ra ®
V
e
V
® Modeled
V Observed
i i i i i
JUNE
JULY
AUGUST
SEPTEMBER
OCTOBER
NOVEMBER
FIGURE 20. DRM VERIFICATION FOR TEMPERATURE IN LONG LAKE, 1971 (OUTFLOW)
-------
TABLE 20
DRM VERIFICATION (PHASE II)
COEUR D'ALENE LAKE, 1971
Date
June 16
July 14
August 20
August 27
Depth (m)
0
0
0
0
Lake Temp .
DBS
16.0
19.0
24.0
21.8
C°c)
MOD
15.6
20.1
19.7
20.6
Outflow Temp. (°C)
OBS MOD
15.1
17.8
20.1
19.9
50
-------
SECTION IV
DELIVERABLES
The deliverables required from Phase II of this contract are a
Verification Report (this document) and source and data decks for
all the verified models.
The following three source decks (in level G FORTRAN IV language,
executable on the IBM 370/155 system), and 20 data decks corresponding
to the verification runs described in Section III are being separately
supplied to the Project Officer.
Source Deck
Data Deck (for the indicated Source Deck)
Region Period
DOSAG
RECEIV
DRM
RR1
RR1
RR2
RR2
RR3
RR3
RR4
RR4
RR5
RR5
RR1
RR2
RR2
RR3
RR3
RR4
RR4
RR5
Long Lake
Coeur d'Alene
Lake
July - August 1971
August - September 1971
August 1969
September 1969
August 1969
September 1969
July - August 1968
August - September 1968
August 1971
September 1971
August - September 1971
August 1969
September 1969
August 1969
September 1969
July - August 1968
August - September 1968
September 1971
June - November 1971
June - November 1971
51
-------
APPENDIX A
MODIFICATIONS TO DOSAG
53
-------
APPENDIX A
MODIFICATIONS TO DOSAG
Class (ii): Desirable before use in Phase III
Modification Dl.
Purpose:
To enable DOSAG to handle a region with no junctions (branching
points) directly, i.e., without the need for dummy tributaries, etc.
Description:
When there are no junctions included, the section of the code which
reads in the File D data is skipped.
Modification D2.
Purpose:
To facilitate use of program.
Description:
Default values were added for K and for the exponents and coefficients
of flow used to calculate depth and velocity.
Modification D3.
Purpose:
To facilitate use of program.
Description:
Only reaches with non-zero lateral inflow need to be input to file F-2.
Modification D4.
Purpose:
To facilitate interpretation of output results.
Description:
A plotting capability (subroutines PLTSET, PLOT) was added to plot
constituent concentration versus reach.
55
-------
APPENDIX B
MODIFICATIONS TO SWMM (KECEIV)
57
-------
APPENDIX B
MODIFICATIONS TO SWMM (RECEIV)
Class (i): Necessary to execute Phase II model
^Mod if icat ion JR1.
Purpose:
To enable RECEIV to simulate shallow, steep, non-tidal river reaches
(particularly in River Region 2).
Description:
(a) For the non-tidal case, the check for a dry junction (negative
depth) at the integration half step was removed.
(b) For the non-tidal case, the junction volumes written on disk
for use by the quality portion of the program were made
consistent with the channel volumes, i.e., if the channel
volumes are known, the junction volumes are also known
(neglecting overlap, which will be small for river applications)
Modification R2.
Purpose :
To correct quality errors in SWMM subroutine LOOPQL.
Description:
The error found by SCI in subroutine LOOPQL is the same error documented
by Bob Shubinski of WRE in his letter of November 10, 1972 to Mr. Torno of
EPA, and subsequently supplied to SCI, on February 27, 1973. SCI recom-
mends changing only the statement labeled 288 in LOOPQL from
288 C(J,KC)=C(J,KC)+(MADD(J,KC)/VOL(J))*DELTQ
288 C(J,KC) = (1.-QIN(J)*DELTQ/VOL(J))*C(J,KC)+(MADD(J,KC)/VOL(J))*DELTQ
This results, in effect, in the same new statement 288 as that suggested
by WRE in their letter.
The remainder of the changes proposed by WRE prevent the above-mentioned
correction from being made at the downstream boundary junction (where
J = JGW). Since in the non-tidal case the program sets QIN(JGW) equal
to zero, no matter what value is input, the correction is not required
at the boundary junction, and the additional WRE changes are satisfactory
They appear not to be necessary, however, since statement 288 as recom-
mended above, produces correct results whether QIN(J) is zero or not.
59
-------
Modification R3.
Purpose:
To prevent possible storage of unwanted plot data in an invalid core
address.
Description:
The array JPLT in common is not initialized. If the input value of
NPLT is zero, i.e., if no plots of junction surface elevations are
required, an IBM "error 240" may occur, depending on what random
value is in JPLT(l) . The following fix is recommended.
Insert the following after card SWFL203 in subroutine SWFLOW.
IF(NPLT.EQ.O) GO TO 350
Insert the following after card SWFL522 in subroutine SWFLOW.
IF(NPLT.EQ.O) GO TO 1040.
Class (ii): Desirable before use in Phase III
Modification R5.
Purpose:
To facilitate use of the program.
Description:
(a) Channel length is input in miles, not feet.
(b) Constituent concentrations of inflows are input in mg/L, not
Ibs/day.
(c) Junction surface area and depth are calculated internally
instead of being input (in steady-state case only).
(d) Input channel depth is actual depth, not distance from datum
plane to channel bottom.
(e) Convergence checks were added to subroutines SWFLOW and SWQUAL
to terminate integration if convergence is reached before the
input stop time is reached (in steady-state case only).
(f) Average daily values of channel flow and velocity are calcu-
lated and printed for each channel.
60
-------
(g) The capability was added to write the hydrodynamic data
onto tape for the single day (initially unknown) on which
hydrodynamic convergence occurred, as determined by item
(e) above. (Note: previously, if only one day's data were
required to be written on tape, as for the case of steady-
state in which the quality model SWQUAL would use the same
day's data repetitively, this data had to be for the
initially prescribed stop date.)
Modification R6.
Purpose:
To provide better hydrodynamic convergence, reduce input preparation
and run time.
Description:
As an optional alternative to prescribing Manning's n values for each
river reach, compute these from the channel configurations and initially
prescribed flow conditions. (This results in consistency in the
energy losses amongst the various reaches, whence more rapid convergence
in steeper streams.)
61
-------
APPENDIX C
MODIFICATIONS TO DRM
63
-------
APPENDIX C
MODIFICATIONS TO DRM
Class (i): Necessary to execute Phase II model
Modification Ll.
Purpose:
To enable the as-delivered model version to be applied to other cases
than Dworshak Reservoir/Libby Dam.
Description:
Generalize the code and documentation, supplement poor documentation,
including the following:
(a) Numerous hard-wired values were changed to be input
variables. Among the values generalized were first and
last day of run, fitting coefficients, tape (disk) units
used for interfacing, lake surface elevation, lake bottom
elevation, dam outlet elevations, and the number of
outlets.
(b) An option was added to allow an outlet to always be at the
lake surface, i.e., to simulate a spillway.
Modification L2.
Purpose:
To prevent referencing data from an invalid core address.
Description:
The variable NUME has not been initialized when subroutine SUBC is
called. This may also cause an IBM "error 240". The following fix
is recommended.
Change the 5th line after 'statement lable 105 in subroutine SUBC from
IF(T(NUME)-T(1).GT. .1) IMIX = NUME
to
IF(T(MAXE)-T(1).GT .1) IMIX = MAXE
65
-------
Class (11); Desirable before use in Phase III
Modification L3.
Purpose:
To facilitate the use of the program.
Description:
(a) The three separate subprograms BAL, MIPP, and TSIP were
integrated into one program.
(b) The format of the input cards defining inflow, outflow, and
inflow temperatures was changed so that the quantities are
input grouped together by day number, instead of being input
separately in blocks to match the months being simulated.
(c) The initial temperature is printed for each level of the lake
in the output.
(d) Automatic interpolation is done if the temperature history
of the lake inflow is incomplete.
(e) An option to apply a factor and a bias (representing ground-
water flow) to the inflow was added.
66
-------
SELECTED WATER
RESOURCES ABSTRACTS
INPUT TRANSACTION FORM
Spokane River Basin Model Project
October, 1974
i; /•>.- 'nrriV'- Or^r.r. :;; Environmental Protection Agency
Set of six volumes: Volume I - Final Report, Volume II - Data
Report, Volume III - Verification Report, Volume IV - User's Manual for Steady-
state Stream Model, Volume V - User's Manual for Dynamic Stream Model, Volume VI -
User's Manual for Stratified Reservoir Model.
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 quality 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 con-
ducted 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.
Send To:
WATER RESOURCES SCIENTIFIC INFORMATION CFNTFD
US DEPARTMENT OF THE INTERIOR
WASHINGTON.D C 2O24O
E. John Finnemore
Systems Control, Inc.
-------
|