-------
Table D-l (Cont'd)
Sample No.* pH Cond Na Ca Mg K HC03 C03 S04 N03 Cl
SS-12/2045-820
T-13/0970-821
S-13/1017-821
S-13/1017-821
SS-13/1050-821
S-13/1055-821
S- 13/1117-821
SS-13/1145-821
SS-13/1215-821
SS-13/1304-821
SS-13/1355-821
SS-13/1508-821
7.7
8.0
7.7
7.6
7.5
7.9
7.9
7.4
7.8
7.5
7.8
7.7
3600
300
400
400
2600
300
300
2800
3000
3200
3200
3200
4.8
0.2
0.2
0.2
2.7
0.2
0.2
2.8
3.0
3.2
3.4
3.4
20.6
1.5
1.8
2.6
17.0
1.8
1.6
17.4
18.2
18.7
19.4
19.7
23.6
1.0
1.1
1.4
16.0
1.1
1.0
.18.1
18.3
19.8
20.1
20.0
0.7
0.2
0.2
0.2
0.5
0.2
0.2
0.6
0.6
0.6
0.6
0.6
5.0
2.3
2.4
2.3
4.8
2.4
2.1
4.7
4.5
4.9
4.9
5.0
0
0
0
0
0
0
0
0
0
0
0
0
44.3
0.8
1.0
2.5
31.2
1.2
1.0
35.2
36.4
39.1
39.1
37.8
0.07
0.01
0.01
0.01
0.07
0.01
0.02
0.14
0.16
0.19
0.21
0.21
0.1
<0.1
0.1
0.1
0.1
0.1
<0.1
<0.1
0.1
<0.1
<0.1
<0.1
Sample No. indicates type of sample (S=surface, SS=subsurface, T=water supply in
tank), run number, time, and date (month and day).
242
-------
Appendix E
Discharge Data
Table E-l
Date
Mar 28
29
30
31
Apr 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
May 1
2
3
4
5
6
7
8
9
10
-C2
mj/sec
_
-
-
0.29
0.29
0.30
0.31
0.33
0.34
0.35
0.37
0.37
0.37
0.34
0.31
0.34
0.37
0.37
0.37
0.37
0.36
0.36
0.36
0.38
. 0.44
0.53
0.62
0.65
0.80
0.82
0.74
0.60
0.52
0.49
0.45
0.49
0.60
0.69
0.63
0.55
0.52
0,54
0.65
0.84
Mean Daily
C3
m3/secx!03
2.5
2.5
2.5
2.5
2.5
2.8
3.1
3.4
4.0
4.5
4.8
4.2
3.4
3.4
4.2
4.8
5.7
6.2
7.6
10.2
12.5
11.6
9.9
9.4 -
12.5
14.2
18.7
22.4
26.3
32.9
39.1
41.6
42.5
42.2
41.3
39.4
37.4
37.4
38.5
39.1
40.8
43.0
38.8
42.5
Discharges
-C6*
nr/sec m
_
-
-
0.52*
0.52
0.52
0.54 .
0.55
0.57
0.58
0,60*
0.54
0.49*
0.51
0.52*
0.58
0.64*
0.67
0.68*
0.72
0.74*
0.67
0.60*
0.69
0.84*
0.94
0.99*
1.16
1.45*
1.38
1.25*
1.13
1.06*
0.95
0.88
0.85*
0.99
1.11*
1.05
0.95*
0.90
. 0.90*
1.03
1.27* .
, 1975
C9
3/secxl03
1.1
1.1
0.9
0.9
0.9
0.9
0.6
0.6
0.6
0.6
0.9
0.9
0.9
0.9
1.1
1.7
1.7
1.1
2.0
2.3
4.1
0.9
0.9
2,0
7.1
12.2
19.3
36.8
62.9
83.8
81.3
62.0
43.3
31.2
26.3
25.2
29.5
39.1
44,5
36.5
25.8
22.4
20.7
„ 24.6
3 CIO .
nr/secxlOJ
3.7
3.4
3.4
3.1
3.1
2.8
2.8
3.1
3.4
4.0
4.8
5.1
5.4
5.4
5.9
6.8
7.6
9.1
11.0
17.0
29.2
32.9
27.2
23.5
22.9
28.6
39.1
53.0
58.3
68.8
74.8
60.9
44.5
32.9
25.8
21.8
21.2
30.9
45.0
39.4
30.9
25.8
23.8
24.9
243
-------
Table E-l (Cont'd)
Date
May 11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
Jun 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
-C2
nr/sec
1.12
1.20
1.18
1.70
2.31
2.64
3.12
2.93
2.93
3.71
3.63
3.43
3.12
2.92
2.69
2.46
2.18
1.93
1.69
1.49
1.66
2.12
2.55
2.83
3.26
3.64
4.06
7.13
9.92
8.62
4.46
3.27
3.47
5.34
9.56
11.67
11.52
8.33
7.47
7.62
6.51
5.24
4.67
4.39
4.96
5.28
4.40
4.46
o C3
nr/secxlO0
56.9
70.8
51.0
42.8
42.5
42.5
45.0
47.9
52.4
60.9
69.7
77.6
49.3
45.6
41.9
39.1
36.0
33.4
30.9
23.5
13.6
20.4
20.4
20.1
19.0
18.1
18.4
18.4
18.4
18.4
18.1
18.1
18.7
18.7
18.7
18.7
19.0
19.3
18.4
18.1
18.4
18.7
19.0
19.5
18.7
1-7.0
15.0
12.5
_C6*
nr/sec
1.52
1.57*
1.53
2.01*
2.66
3.00
3.51
3.35
3.38
4.19
4.13
3.95*
3.64
3.44
3.21
2.99
2.71
2.46
2.22*
2.02
2.19
2.66
3.08
3.37
3.80
4.18*
4.77
8.05
10.98
9.87
5.85
4.76*
4.89
7.75
10.72
12.66
12.37
9.06
8.09*
8.27
7.18
5.92
5.37
5.09
5.72
6.10
5.31
5.41*
3 C9 1
m°/secxlOJ
32.0
41.1
40.5
38.8
36.8
36.5
34.6
34.6
36.8
34.6
30.3
26.1
24.9
24.6
24.1
23.8
23.2
22.7
22.1
20.1
19.0
17.8
16.7
15.9
15.6
15.0
14.7
15.0
15.6
14.4
13.9
12.7
12.5
11.9
11.6
11.0
11.3
10.8
10.8
10.5
10.2
9.9
9.6
9.3
9.3
8.8
8.5
8.5
3 C1° 3
nr/secxlO
32.9
39.1
41.6
44.2
57.5
62.0
60.3
56.6
47.3
44.5
42.5
41.6
38.5
33.1
30.3
28.9
28.3
28.3
28.3
26.3
25.2
24.4
24.9
24.6
22.1
19.5
19.3
18.4
20.4
20.4
19.5
17.6
16.1
16.4
17.6
19.0
21.5
26.3
32.9
28.6
22.4
16.7
11.9
9.1
8.2
7.9
7.6
7.4
244
-------
Table E-l (Cont'd)
0C2 . C3 -
Date m3/sec nr/secxlO"3 m
Jun 28
29
30
Jul 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
Aug 1
2
3
4
5
6
7
8
9
10
11
12
13
14
5.45
5.24
4.96
4.81
4.46
3.98
5.32
4.67
3.96
3.40
2.93
3.03
2.71
1.97
1.68
1.48
1.60
1.47
1.55
1.08
0.83
0.75
0.80
0.69
0.56
0.51
0.47
0.46
0.47
0.49
0.51
0.67
0.69
0.61
0.59
0.47
0.47
0.44
0.35
0.37
0.34
0.29
0.30
0.43
0.47
0.54
0.59
0.55
9.9
7.1
5.1
5.9
7.6
6.5
5.1
5.7
5.4
5.1
5.4
5.7
5.4
4.8
4.5
4.5
5.7
5.7
6.2
5.9
5.7
5.1
4.8
4.5
4.2
4.0
4.0
3.7
2.8 _
2.8
2.3
2.8
3.4
3.1
2.8
2.5
2.3
2.3
2.3
2.3
2.0
2.0
2.0
1.7
1.7
1.7
2.5
2.0
C6* 3
6.30
5.95
5.56*
5.31
4.81
4.23*
5.61
4.97
4.28
3.73*
3.28
3.40
3.10
2.37*
2.19
2.05
2.23*
2.08
2.10
1.61*
1.33
1.23
1.25
1.09*
0.95
0.89
0.82
0.80*
0.82
0.85
0.86
1.04
1.07*
0.94
0.85
0.72
0.73
0.71*
0.61
0.63
0.59*
0.53
0.52
0.66*
0.71
0.79
0.88
0.86*
C9 , , CIO -
secxlOJ nrvsecxlO-3
8.2
7.9
7.9
7.9
7.4
7.4
7.1
7.1
7.1
7.1
7.4
7.4
7.1
6.8
6.5
6.5
7.1
7.1
7.1
7.1
7.1
7.1
7.4
7.4
7.4
7-4 „
7.1
6.8
6.8
6.5
6.8
6.8
7.1
7.1
7.1
6.8
7.1 .
6.5
6.5
6.8
6.5
6.5
6.5
6.2
5.9
5.7
5.4
5.4
7.4
7.1
7.1
5.7
40
.5
4.5
4.0
4.0
4.8
5.4
4.8
4.5
4.0
3.7
3.1
2.8
2.3
2.3
2.3
2.3
2.0
2.0
2.0
2.0
2.0
2.0
2.0
2.0
2.3
2.0
2.0
2.3
2.0
2.0
1.7
1.7
1.4
1.1
1.1
1.1
1.1
1.1
1.1
1.1
1.0
1.0
1.0
1.0
245
-------
Table E-l (Cont'd)
-C2
Date nr/sec
Aug 15 0.52
16 0.50
17 0.48
18 0.53
19 0.53
20 0.52
21 0.55
22 0.54
23 0.44
24 0.45
25 0.43
26 0.42
27 0.43
28 0.39
29 0.35
30 0.33
31 0.32
Sep 1 0.30
2 0.31
3 0.30
4 0.29
5 0.29
6 0.30
7 0.30
8 0.31
9 0.31
10 0.31
11 0.32
12 0.32
13 0.32
14 0.33
15 0.33
16 0.34
17 0.35
18 0.35
19 0.36
20 0.36
21 0.37
22 0.37
23 0.37
24 0.36 .
25 0.36
26 0.37
27 0.37
28 0.36
29 0.36
30 0.31
Oct 1 0.31
3 C3 3
nr/secxlO0
2,0
1.7
1.7
1.7
1.7
1.7
4.0
5.1
2.3
1.7
1.7
1.7
1.4
1.4
1.4
1.4
1.4
1.4
1.4
1.4
1.1
1.1
1.1
1.1
1.1
1.1
1.4
1.7
1.1 ~
1.4
1.4
1.1
1.1
1.1
1.1
1.1
1.1
0.9
1.0
1.0
1.0
1.0
1.0
1.0
1.0
1.0
1.0
1.0
C6*
nr/sec
0.82
0.78
0.77*
0.80
0.78
0.76
0.77*
0.76
0.67
0.69*
0.67
0.66
0.68
0.64*
0.59
0.58
0.57.
0.55
0.55
0.54
0.53
0.54*
0.54
0.54
0.55
0.55
0.56
0.56
0.56
0.57*
0.57
0.58
0.59
0.60
0.61
0.62
0.63
0.64*
0.65
0.64
0.63
0.64
0.64
0.65
0.63*
0.64
0.59
0.59
T °9 T
nr/secxlOv5
5.4
5.4
5.4
5.4
5.4
5.4
5.4
4.8
4.2
4.0
3.4
' 2.8
2.5
2.3
2.3
2.3
2.3
2.3
2.0
2.0
2.0
2.0
2.0
2.0
2.0
270
2.0
2.0
2.0
2.0
2.0
2.0
2.0
1.7
1.7
1.7
1.7
1.7
1.7
1.7
1.7
2.0
2.0
2.3
2.3
2.3
2.3
2.3
. CIO
m /secxlO
1 0
4. * \J
1 0
X » w
i n
x • w
1 0 '
JL « W
1 0
* • w
1 0
X • W
1 0
X « W
1 0
*• • \J
1 0
X * \J
1 0
X • W
1 0
X • \J
1 1
X • X
1 1
X • X
1 1
X • X
1 1
+ * X
1 1
X • X
1 1
X • X
1 0
x • w
1 0
X » \J
1 0
X • \J
0 9
w • J
0 9
w • «/
0 9
W • *J
0 9
w • •/
0 9
w • ^
0 9
w * ^
0 9
w * «>
0 7
W • /
0 7
v • t
0 7
W • /
0 7
W • /
0 7
w • f
0 7
w • /
0 7
\J • f
0 7
W t f
0 7
w • f
0 6
w • v
0 6
w • v
0 6
W t W
0 6
w • \j
0.6
w * w
06
w • w
0 6
w • v
0 6
W • U
0 6
w • w
0 6
V • w
0 fi
\I * \J
0.6
246
-------
Table E-l (Cont'd)
C2
Date m3/sec rrr
Oct 2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
Nov 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
0.31
0.31
0.25
0.26
0.28
0.29
0.28
0.30
0.31
0.31
0.31
0.31
0.32
0.32
0.30
0.29
0.28
0.28
0.26
0.25
0.25
0.24
0.24
0.46
0.46
0.45
0.44
0.44
0.44
0.44
0.43
0.42
0.42
0.42
0.42
0.42
0.42
0.43
0.40
0.43
0.42
0.44
0.45
0.47
0,48
0.63
0.52
0.77
03 3
/secxlO-3
1.0
1.0
1.0
1.0
,1.0
1.0
1.0
1.1
1.1
1.1
1.1
1.1
1.1
1.1
1.1
1.1
1.1
1.1
1.4
1.7
2.0
2.0
2.3
2.0
2.0
1.8
1.7
1.7
1.7
1.7
1.7
1.7
1.7
1.7
1.7
1.7
1.7
3.4
3.1
2.5
2.0
2.0
1.7
1.7
1.1
1.1
1.1
1.1
C6* . C9 CIO
m3/sec m3/secxlOJ nvVsecxlO-3
0.59
0.59*
0.52
0.52
0.53
0.54
0.52
0.52
0.54*
0.54
0.55
0.57 .
0.58
0.59
0.58
0.58*
0.58
0.58
0.56
0.55
0.54
0.54
0.53*
0.74
0.75
0.74
0.73
0.73
- 0.74
0.74
0.71
0.71
0.70
0.70
0.69
0.69
0.69*
0.68
0.65
0.67
0.65
0.66
0.65
0.65
0.66*
0.81
0.71
0.99
2.3
2.3
2.3
2.3
2.3
2.3
2.3
2.3
2.3
2.3
2.3
2.3
2.3
2.3
2.3
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
3.1
3.1
3.1
3.1
3.7
3.7
3.7
3.1
3.1
3.1
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
0.6
0.6
0.6
0.6
0.7
0.7
0.6
0.6
0.6
0.7
0.7
0.6
0.6
0.8
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.6
0.3
0.3
247
-------
Table E-l (Cont'd)
Date
Nov
19
20
21
22
23
24
25
26
27
28
29
30
,C2 C3
m-Vsec m°/secxlOJ
0.53
0.53
0.53
0.53
0.53
0.53
0.50
0.50
0.50
0.50
0.50
0.48
0.9
0.9
0.9
0.9
0.9
0.9
0.6
0.6
0.6
0.6
0.6
0.6
,C6* C9 CIO .
nr/sec nrVsecxlO-5 nr/secxlOJ
0.76
0.75
0.75
0.75
0.75
0.75
0.72
0.72
0.72
0.72
0.72
0.70
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.5
0.3
0.3
0.3
0.3
0.3
0.3
0.0
0.0
0.0
0.0
0.0
0.0
Indicates actual measurement taken at Station C6.
248
-------
Table E-2
Mean Daily Discharges, 1976
-C2
Date m3/sec m
Mar 22
23
24
25
26
27
28
29
30
31
Apr 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
May 1
2
3
4
5
6
7
0.68
0.47
0.46
0.44
0.39
0.55
0.51
0.51
0.49
0.49
0.49
0.51
0.55
0.52
0.54
0.54
0.48
0.50
0.60
0,65
0.64
0.84
0.83
0.76
0.75
0.72
0.71
0.68
0.65
0.63
0.64
0.70
0.71
0.61
0.59
0.61
0.60
0.66
0.76
0.79
0.71
0.79
0.86
0.92
1.04
1.26
1.22
C3 C6* C9 CIO
J/secxlOJ nr/sec mJ/secxlOJ . nr/secxlO°
_
1.7
1.9
2.7
2.7
2.8
3.1
2.9
2.4
2.7
6.3
5.7
6.4
7.2
9.5
10.7
9.4
8.7
9.4
9.6
10.4
12.8
16.9
18.4
13.5
17.8
25.3
30.6
43.0
39.9
29.6
27.0
24.5
21.2
19.1
18.5
16.3
17.0
17.6
17.8
14.4
10.1
9.0
8.0
11.8
11.4
9.5
1.12
0.95
0.97
0.98*
1.18
1.54*
1.42
1.27
1.14*
1.20
1.23*
1.30
1.37*
1.31
1.32
1.28*
1.33
1.44
1.59*
1.61
1.67*
1.72
1.64*
1.52
1.45*
1.48
1.49*
1.56
1.59
1.64*
1.44
1.33*
1.36
1.28
1.62*
1.52
1.28*
1.31
1.39
1.38*
1.39
1.55*
1.63
1.69*
1.78
1.98*
1.93
_
0.1
0.1
0.1
0.1
0.1
0.1
0.1
0.1
0.2
0.2
0.6
0.8
1.1
1.4
1.7
5.7
8.5
10.9
19.8
25.7
25.4
30.7
28.3
2"7.3
26.5
25.3
21.8
17.8
15.4
14.9
14.9
15.6
15.0
14.9
15.9 *
16.3
15.9
16.7
17.4
17.4
16.7
16.3
16.2
15.9
16.7
17.4
_
4.1
4.4
5.7
6.1
7.2
8.2
8.9
8.5
8.6
.9.8
10.5
11.0
11.6
12.2
14.4
17.0
19.3
21.0
20.7
20.1
21.8
23.6
22.4
19.8
17.5
16.6
15.1
13.7
11.0
9.9
9.3
8.6
9.7
12.8
10.2
8.6
8.2
7.9
7.6
7.6
7.4
7.1
6.8
6.8
6.6
6.8
249
-------
Table E-2 (Cont'd)
Date
May 8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
Jun 1
2
3
4
3 C2
nr/sec
1.56
1.38
1.46
1.67
1.56
1.43
1.95
1.96
1.60
1.91
2.17
2.36
2.87
2.18
2.07
2.31
3.45
3.06
2.46
2.87
3.60
4.36
4.92
3.18
3.56
5.04
5.19
7.25
3 " 3
nr/secxlO0
12.7
12.1
11.1
10.7
9.7
9.9
8.9
6.1
8.5
10.8
10.1
9.4
10.2
10.0
10.7
10.7
10.3
10.1
9.7
8.9
7.8
7.1
6.8
6.6
6.3
6.1
5.6
5.6
_C6*
m /sec
2.24*
2.12
2.23
2.49*
2.38
2.21*
2.75
2.78
2.45*
2.59
2.74*
2.97
2.51
2.85*
2.72
2.93
4.06*
3.91
3.51*
3.77
4.34*
4.96
5.38
3.51*
3.97
5.42*
5.64
7.79*
3 C9 3
nr/secxlOJ
17.6
18.9
15.6
16.1
15.3
15.1
13.3
12.6
12.0
11.3
10.8
• 16.3
22.3
21.0
21.6
22.1
20.5
18.6
11.7
11.1
10.6
9.2
7.1
7.1
7.1
7.1
6.6
6.6
3 C1° 3
nT/secxlO-3
7.1
8.4
10.1
9.7
9.0
9.0
8.5
7.7
7.0
6.6
6.6
8.5
8.4
9.2
8.9
9.2
10.9
8.3
8.1
7.8
7.6
7.2
6.9
6.3
6.1
6.1
6.5
6.6
Indicates actual measurement taken at Station C6.
25fl
-------
Table E-3
Discharge at the Time Monthly Water Quality Samples Were Taken
Month nr/sec
C3 C6 C9 I OO
nr/secxlQ0 nr/sec nr/secxlO0 nr/secxlO
4/75
5/75
6/75
7/75
8/75
9/75
10/75
11/75
12/17/75
0.31
0.57
3.88
4.14
0.78
0.32
0.29
0.37
0.32
2.4
38.2
18.7
6.5
2.6
0.9
1.1
1.7
0.9
0.55
0.82
-
5.55
0.96
-
0.57
0.60
-
3.3
25.2
14.2
8.6
7.7
2.3
1.7
2.6
-
3.3
29.9
21.0
5.9
2.3
4.7
0.6
0.4
-
1/76
2/76
3/ 2/76
3/31/76
5/ 5/76
6/ 2/76
0.28
1.07
1.01
2.0
12.2
5.9
1.12
1.59
5.24
0.0
15.6
6.5
1.7
7.7
7.7
5.9
251
-------
Appendix F
Saturated Paste Data - Soil and Spoil Analyses
Table F-l
1976 Samples
Sample No.
SS100
0
30
60
SS101
0
45
120
SS102
0
30
60
SS103
0
30
••
SS104
' 0
30
Cond.
(umhos/cm)
2700
2900
3500
500
700
200
1200
2600
2600
300
3100
_
.
500
300
Sample No.
SS105
0
30
60
SS106
0
30
60
SS107
0
30
120
SS108
0
30
60
120
SS109
0
45 -
Cond.
(ymhos/cm)
300
600
300
500
3000
4000
200
300
500
300
200
400
1400
1800
2700
252
-------
Appendix G
Leaching Data - Composite Spoil Sample
Table G-l
Sample Size Distribution
Size Weight (gm)
Greater than 8 MM
Greater than 5 Mesh
Greater than 9 Mesh
Greater than 16 Mesh
Greater than 32 Mesh
Greater than 35 Mesh
Greater than 60 Mesh
Greater than 80 Mesh
Greater than 100 Mesh
Greater than 120 Mesh
Greater than 200 Mesh
Less than 200 Mesh
248.3
247.4
256.0
236.8
201.4
23.2
74.6
35.0
14.1
7.9
22.3
19.4
253
-------
Table G-2
Results of Leaching Test
Volume
Leachate
(At) pH
100
200
300
400
500
600
700
800
900
1000
1100
1200
1300
1400
1500
1600
1700
1800
1900
2000
2100
2200
2300
2400
2500
2600
2700
2800
2900
3000
3100
3200
3300
3400
3500
3600
3700
3800
3900
4000
4100
4200
4260
7.75
7.70
7.64
7.61
7.63
7.68
7.50
7.92
7.90
7.82
7.82
7.87
7.92
7.95
8.05
8.02
8.02
8.09
8.11
8.11
8.01
7.99
8.10
8.10
8.11
7.98
8.08
8.10
8.10
8.09
8.10
8.15
8.10
8.12
8.10
8.11
8.10
8.11
8.10
8.12
8.13
8.14
8.15
Spec. CaC03
Cond. Total Alkalin
025'C K Na Mg Ca Zn C1 Hard. Ity
(umho/cm) (mg/0 (nig/t) (mg/0 (mq/fc) (mg/i) (roq/J.) (mq/0 (mq/O
First Run (Initial Conditions)
1470.
1394
1314
1247
1220
1200
1000 J
941
867
787
728
677
630
570 '
500
480
450
437
413
413
406
372
336
318
300-
300 ,
281
270
270
262
250
22.0 20.8 80.3 231 0.058 0.12 995
27.0 9.9 49.2 158 0.21 12.5 660
18.0 5.4 31.0 99 0.020 3.0 425
0.96 321
9.0 4.1 19.8 81 0.018 0.48 231
""
6.7 3.0- 14.7 61 0.017 0.27 168
248 J . '
248
240
313
296
277
272.
257,
255
257
255
t
5.1 2.6 13.0 50 0.012 0.29 166
4.8 2.6 11.0 48 0.012 0.27 165
306
94.7
89.3
88.6
76.6
72.6
69.1
67.5
64,4
254
-------
Table G-2 (Cont'd)
Volume
Leachate
(•*)
J!tL
Spec.
Cond.
025'C K Na Kg Ca In
(pfnho/cm) (mg/t) (mg/t) (mg/t) (mg/t| (rcg/t)
CaCOi
Total Alkalin-
Cl Hard. 1ty
(mg/t) (mg/t) (mg/Q
4360
4460
4560
4660
4760
4860
4960
5060
5160
5260
5360
5460
5560
5660
5760
5860
5960
6060
6160
6260
6360
6460
6560
6660
6760
6860
6960
7060
7160
7260
7360
7460
7560
7660
7760
7860
7.80
7.88
7.90
7.93
7.95
8.02
8.11
8.06
8.98
8.02
8.12
8.11
8.18
8.22
8.31
8.27
8.24
8.31
8.43
8.47
8.42
8.33
8.40
8.34
8.52
8.55
8.44
8.51
8.49
8.50
8.45
8.42
8.58
8.55
8.60
8.49
450
290
257
206
168
146
140
122
118
107
92.4
85.3
85.3
83.1
83.0
83.1
82.5
72.'0
67.5
61.9
61.9
61.5
59.5
59.5
59.5
62. O
Second Run (After 60 hr aeration)
3.4 1.06 7.2 20.4 0.0105
2.1 0.620 4.1 12:0 0.009
1.7 0.459 3.0 8.0 0.007
2.0 0.380 2.3 6.5 0.005
1.5 0.330 1.9 6.2 0.004
1.6 0.288 1.6 5.4 0.0025
U2 0.280 1.5 5.8 -
33.9
0.34 85.6 32.2
0148 61.8 30.1
0.27 47.9 29.3
°'26 41'3
26'°
0-23 36.1 25.3
0.26 36.9 26.9
255
-------
Table G-2 (Cont'd)
Volume
Leachate
(IM)
7985
8010
8035
8060
8085
8110
8135
8160
8185
8210
8235
8260
8285
8310
8360
8860
9260
9360
9460
9560
9660
9760
9860
9960
10060
10160
10260
10360
10460
10560
10660
10760
10860
10960
11060
11160
11260
11360
11460
11560
11660
11760
11860
11960
12010
pH
7.80
7.90
8.03
8.30
8.37
8.31
8.36
8.33
8.40
8.44
8.19
8.24
8.40
8.41
8.42
8.00
8.10
7,18
7.39
7.80
7.93
8.21
8.28
8.33
8.30
8.30
8.18
8.20
8.19
8.22
8.25
8.05
8.01
8.01
8.10
8.15
8.24
8.10
8.14
8.19
8.13
8.25
8.12
8.14
8.52
Spec.
Cond.
025'C K
(ymho/cm) (mg/i)
Third
1293
990
943
901
903
889
862
858
830
847
820
8ZO
801
797
6.9
11.2
12.8
11.0
796 12.0
676 8.4
526 4.0
Fourth
518
441
346
271
6.4
248
233.
198
175
164
4.1
161
147,
139
132
127
3.5
115
110
107
102
93
2.8
89
90
89
87
85
2.3
81
87
87 2.1
88.5
Na
(mg/i)
Mg
Run (After 120
10.3
8.2
9.6
6.7
7.1
4.9
2.4
62.9
39.2
38.0
34.6
35.0
27,0
12.7
Ca
) (iwj/i)
Zn
hr aeration)
211.0
130
.
124
112
106
80
43
0.049
0.60
0.022
0.016
0,011
0.018
0.025
CaOh
Total AlkaHn
Cl Hard. Uy
(mq/O (mg/i) (mq/4)
0.72
0.21
24.2
9.1
'
Run tAfter drying & crushinq)
4.4
0.860
._L._
0.83
0,42
0.30
0.31
22.8
9.7
6.5
4.7
3.9
3.8
64,0
28.9
20.6
15.7
13.7
14.0
0.021
^
0.016
0.011
0.011
0.010
0.020
-
227 135
.
103 131
.
71 109
-
ce gfi
9W ww
-
48 76
^O * V
47.8 76
256
-------
APPENDIX II
SEDIMENT TRANSPORT MODEL
FORTRAN label
AOF
AGB.AWP,
BEX.BEW
AIM
BIM
CW
OR
DT
DMW
GBOUT
GMAX
GWOUT
ITCOM
ITMAX
NDX
NSTOM
PB
PW
LIST OF FORTRAN LABELS
Definition
Detachment coefficient of runoff
Parameter describing sediment transporting
capacity
Coefficient in raindrop soil detachment equation
Exponent in raindrop soil detachment equation
Concentration of wash load
Rainfall rate
Time increment in computation
Mean wash load sediment size
Bed-material load transport rate
Maximum depth to which a raindrop.can penetrate
the soil layer
Wash load transport rate
Total number of time increment for computation
Total number of time increment at the end of an
event
Number of space increment
Number of storm for computation
Percent of bed-material load size in a typical
soil sample
Percent of washload size in a typical soil
sample
Symbol
'1,2
w
V
At
257
-------
FORTRAN Label
QOUT
SB
SKI
SLEN
SLOPE
SNU
. STORM
SW
TITLE
ZE
LIST OF FORTRAN LABELS (Cont'd)
Definition
Water discharge
Depth of loose soil for bed-material load size
Constant representing grain resistance without
rainfall for Nr < 900
Length of an overland flow plot
Bedslope
Kinematic viscosity of water
Alphabetical or numerical identification of an
event
Depth of loose soil for wash load size
Alphabetical or numerical identification of the
problem
Net depth of loose soil
Symbol
L
S.
w
258
-------
PROGRAM SEDIM
SEDIM Ti>AO: «-OC 6<.OU FTN VJ.O-PJtib OPT=I
PWlKiflAM btOIM (INf-ur.OUlHUTl
MAU FKUM OVtMLANt) FLOW
. IULNTIMCATION oe THE
NO' = KilM,;1!. ff ir^Cr. |l»C^L"t •:! S
NSlOM = NU'*!»F •' CU tVrUT I- OK CO''IJUT AT I OM
ni - Tl"«- !-r«.r''t-T r'nw iiU'it^lCiL COMMUTATION IN MINUTES
ViO^^TlC VIM'UMTY u> i«rj> IN bUUAKF MErtWb PtH SLCONO »E5
.
C SLfN - l.^'KiTM Of AH (jvt>L«'N.) fLu« PLOT IN
C FVW*'" F»tLL'?t:L:-'UT» Vf MtAM PAKTICLt SIZk FO^ V.ASH LOAO SLUIMtNT
-- OK 4«bn LOAD l'< M 1 LL
-..
P>4 = Pf«-CrNt t F nh 0-^'A fr i I AL LO/VO aUt IN bl-OlLS
' to
8 MM.BlS^V^^^Hir^cUSSKbCKldlNo S^OIL OUACHHtNT Ur «A1NOROP
C Ar,^l.A^^J.lJ^X.SE-.'=r'^A'-^KA•^t TL«^i UtSCHI'llMo SLtMMtNT THANSPOHUNO
r *DF = nr f»O"F'«T cof r " 1C I tNT MI^ SJKKACt ^UwOM EHOSION
I |Vo^5 Dr.r^Al.tTiC^L 0- .iL-icMtAL I r«t N Tl f I CA I I OJ* OJ fcVh M b
r ?T''ftx = TOTAL -.u^-c- ov H'.-. ini.nt >'tfiT «T int eNii Vl',,,ft'N tvtN'
5 U.,,M : n it rU..'., o^ Il^t lMC.«ficNT Kil- CUM-'l.MIlON
$ lTY«'t = I'- SI"ULCTJuli OF i)0xi tVtNl W I I nUU T »» J ^ ALL
r 1TYPF = I. SI "01 a t I<»(| u>- S !()•"•' hVtNT wjln M4INHALL
.5 DM « PA INF ML 0* sio*»-tLl rfUN^r- IN^UT IN MHLIMtltHS P£R
C lNCMr:"£Hl
N TIH.I. I loJ .UTUT JrtOU)
s..i.u.«!"Js.I «i4*,
0 • ' t nun ' • <->•• If'UU I » t .K < oUif » i rV 1 r U « , !•"•. . " « . f V »
.CO1"1 Or- X'tlN/ Juf i«l"t I I
C IMPUT AM) OUTPUT .TIUE
Pf.AO "iO'it TllLt.
PKINT J) Ot T I ILt
C INPUT A\"J OUTPUT TtKt OK
RF.AO :M-« IT
C INHU1 A';U OUTPUT OLNcKAL
PCMJ
c INPUT AMU OUTPUT ^tOMLiM* UMA
PRINT J^oi S.LL««
C INPUT A^'i) OUTPUT SOIL DATA
PUJMT
C 1WPUT A"'0 OUTI'UT
pr AD tll« SI:L
C INPUT A-D OUTPUT bLOlMtM HOUTINC,
Rf40 ?Jfl« «JIM.«.'.u. A-l'tMt'^ «t«tw«bvA
A iMt A0^
C ESTAHl IS* SOMt IMVA^UNT I NKOMMA I 1 ON
I " . TI ""i^s.'ju/ i ooouu t
= <»"f ft/fcll n>. ^3
IF (ITTPt ,tO« I 1 uO TU <«nif
259
-------
A 1 Mill.
r.n TO i, o<»
402 »| ••=«!»' /I OOUOU.
40* f»T =1)1 /' <• •
f)TS = l>l u "'"0.
r)fNrljTV'> UJwt I -MM)
MCl = l<'.".»M,0.
Pw=l .-«•*"
ji"-(/.i|j* .a
00
C INPUT *r.'l) OUTPUT J.-J.MUKF fcvtNl CM*P»CT t'HlST ICb DATA
c INPUT AMD
nu|»iT »S(J
PI' An- * NO. ('•••>( IT) ,1T»1 • ITMA*)
no 1 1» u = i . i IM.\A
nn< I I)rf).»(lTJ/OT
ION DATA
Oil I '.il '. I U
prin ••>"• ('(IMHi
00 12^ I fsl , I F'li
C
123 p|Jii{Mbl»a. (IT.NU (in .11 = 1. UMAX)
i n
cwnn=nuF
127
bSO, «lT.t-*(lTl.lT«liITHAX»
DO 1 •"• ' M • I • 1 I ••«« *
J2*
00 !?'•
ElKlU
126 CONTI»
C 0=0
, ,
00 1JO J=1»NDX
SW(J)-0.
SM1J)=U.
130 COMT1NUI
C KOUTIMJ PC* tACn Tl'-lt INCHLMtNl
ooui tiii = ').'•" r.-'
-------
c
c
180 CON I
i) > >i;-«Si'
OF UlXtCT
ACT«<: f j. jv/
Pf/I'M
190
270
2PO
300
F IIW-U1
FO..H.AJ
320
330
3*0
310
340
390
410
420
430
440
4*0
460
470
b)0
(I.UUU1 I J J
.^ /i<«,t.(i I^'t Mf
/« I • « {••j-ifn-ciuN I 0* *> i MUt aft-.t)
/ JJ* . J .'lAMiiuMf Or vfUL«»tt\l
'jl n s i-'iu. A tti> -»«SM t.-jntj s
II) «G«uUTU) ttiTQT (I) »l = l«ITCOMI
\
*«MU.i»
« ) OArtl
s».2 »-ii /Pt Of "UUlIN'i EVKNT= til)
(/ 1 1"« J1- in. .t)
{//«. «.< < f /i 10.M
{//-> J« • I -^-^<^JNO^ f "iVtNl tXCESS)
(••«<.I) O.'.XtrJU.n)
CNO
SUBROUTINE ROUT
c
c
c
,NSTOM
lio
120
COMMON /hi";/ AliKiAlMi I
THIS SO^UOUI Ihf. ^oUItS THt FLOW DOWN TH£ OVF>LANU FLOw PLOT
COMPOIf AT TlML IT |T«1)!>
IF (IT.RT.ITMAX) GO TO 1JU
DKF =.)(?< I J I
00 TO I?U
o.
si u i r T i •
c
c
IF IIT.M . I T'-lAfc) OU 10
fill 10 /?'•">
rue.to.n i>o TO
l
r.o fi si •
"S': = I>"t«^A.lt.?»tt'
OLATsl < PM/'M Jt'OU.
I ••',.'.»( SCHF,~«t. F0'<
DO
'
,|-I«NOA
J.I)
?>« / J
r.t. A ) « A TfM«r..T A«-;UF»
If (ASl/"i.LL. t.Ol-IOl PO TO
261
-------
C M:i UP A-O HtLAT l
OMr/»bli"/')T X
At f'r (PS* »S'IU/
ftf_ r = i ./i.
»lf MrHf I - I .
fU«,:»l'«-l.
ALH»" ' r t\.'"'»'t. I
• L'il Msfll ""!!( 1
01 Xt-(>1 X 'AC L.
C LINC«H scum? 10 f itiU tnt MHSI APPROXIMATION
"?'»
fx I Alt«-'//M.>'» «« I 1 ,/l-tf)
oo 10 i:<»o
op>«UU •«»'£. T
. .
If ( Af M /VUt V) .L^ .t"«UtV) 00 TO JOO
ir IIT*.'',U r. i--j«i to 10 «:ou
I I • J
?60 i m »=
U »rs/
=f OF f/snf.u
IF (STL'-'.'jl'.O. ) 00 TO
or=»;t« Antv/Kt/iK
00 10 «"-U
270 SirM=«;:jt.r (Sfr.MJ
If (ACLtf.'jt t'J. > o() TO
.
IF an ,t;T .C. ) oo 10 o TO
GO TO i"<0
2VO
Of** I
if « AM .«.-?.*&<•> UIF = J<<;
CO TO irsu1
C
C OttlP"lNATI'jN Of FLO* CUN'JI TlOflS
300 nr°i»'=.M Ue"f 0«"1. T
•
t Ffi
WNsi»L/S";
Vsyfc/MPH
C Dt"1F.Kv'lMa riO-4 UF Int. AMfiUNl 01- SOIL OKTACntO HY «AIN[)HOP IMPACT
IF (P'J.«T.'»H£A
W(HJ) =(»"C/t.-r
C Cf«CCK MJt AVAlLAiflLHr OK HtO MATtHlAL LOAD
.
IF (Lt'M M.OL.O. ) (>0 TO JoU
C
C OETt»**'lMATION uF IMC AMOUNT Of SOIL OETACHCD ur FLOW '
GO 10 J°U
262
-------
C Of UWMIM-UIOM uF ui.K)Sli|ON uut lo LACK OF CAHIIYINIJ CAPACITY
370 PfUjl
C niTFWMUAllON OF AVrtLAblLIfY Of WASH LOAD
380 Z'rfTfMs^v
C WASH IOAU KOUTlNti
IF (TfO.Lf.l.OL-IO) C'O TO WO
C OFTC««IMA1 ION OP *MT1O OF SUSPtNUCU WASH LOAD
r i A')/ i .
* V/SV
ru . «->sv )
IF M>r)t
IF HI i.K'j.n .iiw . /Tbui-'.rj I . u*"/VA) 00 TO
SOOs^lf" ( I .-
420 »«(j>sO.
(-^ 'i"-o | r »k If H
t)nF.M*tlNATl(i«j uF AVAJLAHlLlTt > 0»J Nt«T STtH
440 COM
C DEUrfM»'M ION OF *ATtM. «A^M LOAOf AHO btl)-MAT£H I A 1. LOAO OIbCHA"'r-
460 fOMl'AT (3nx.«.lnOO NOT CONVtUOt F0f« Tut COMPUTATION *
two
263
-------
SUBROUTINE POWER
C
C THIS SU^OUT Kt f,VALUMKS Jl ANO J
c HOI AT tu'i«»
C XJl a WftL'lt Of Jl INUOWAL
C AJ2 = V*I.UL or Jt' [vfMj-'iL
C N a uWUr » UF 4->-K')* In.i I (us • 1
C CONV s
XJl=0.
.
l 00
.
fN=l.
At'XsA
CiO 1U I
110 M-N> I
c*c»o/rv
0=6
r=u» i.
FNI-fLf"*1
lifO P (A»-S(M .LE.O.l'OI ) GU TO JJU
sXJl «O' 1 1 . -«r A J /r.
(,n TO
130
140 IF (N,
CJl-sAI.SH ,-f J
CJ^SAKSI 1 .-
IT (Cjl .LE.CUNV. JND.CJtT.tt .CONV)
ISO
00 TO 110
264
-------
APPENDIX III
UNSATURATED FLOW AND CHEMISTRY MODEL
PROGRAM MOISTRE
"OrSTPE
I < INPUT ,ouT»yT,PUNCH.T»prc, »pirl0)
4P»1 ?0,
• SftS«tf2
rr> rn*rnwM TO Dor,1 MHUJON, JAN. n»i*rs one
miu (i, i t/lndi i «nN i s i . ,IA f r < 10 i (AMtl^O) « TM^ ( SO ) •'if
.» f »l { ft <1 I « T N < MM • f N I F> 0 ) . A N t I Ml , ,!« .M)l • O t Mil ,S(60) ,E(60l ,F(60)
JUH (J«sft) .HP1IM » THIMJI , AUFNT :>»0 >
[MO
',MONTH,UH.KP <
:.ICPOP
CCMMON/AJST/O
ISTAPT
PEAL K.IP.KP3,KSATD«K54T MQQ 1
C-_.—R£»O P9[NT OPTIONS, WON e>4PlMETEPS. «ATE» APPLICATIONS, PROFILE DATA
PEAO 9156, IPUNCH.MESTfl, J \VE, ITENTM«lNFlL5,LLSTRr.HMSTOP, I STOP,
i IOEF,FCAP,JPOPT,ICSOP
HEAD 100.AA.Qfl.CC»LL'*M»PBC •TBC«VEAP, CSOP,M,
C--— OEAO SOILS OAT* FO« CQMPuTI ': HfC»AULlC PftOPEPTICS, COMPUTE HOD
€**•• CONSTANTS MOO
IF( IPOPT.EO.11 CALL CHAH MOO
IF....COMPUTATION OF TIME OF wATrj APPLICATIONS
ST*PT«fl
00 29 L=1»AP»S
AHTIL'**UTIL)•?.5* •«•••••!
29 TME(LI*OAY(04T£(L) .START,M.',,N(L) >
ESTARLISH MOISTURE oiSTRiet'irN AND CONSTANTS
OELTM«l./M
CALL PPOP(TS.TS,TO.KSATO.OS*TO
OELT«OELTM
G«OELX
C tF 8«C EOUALS 0 THFN PRC IS TWANSIFNT
i0.2)i
:O.HTBCaTM
IF»OIF»CON5T»0,0
CL3»0.0
TMli«TRC
SFIl)«0.0
00 43 J«2.0
265
-------
»•»
t j i *o . n
(w I "S"
( «;.«.vi.
„ • 1 1 B^n
[Mir
fttlttlV
ttttttts
IFirc%f».li
i TOI
•<)!
•' t -.Ed. I > TN< j) «rOiji
Th f ji «TN ij>
CONTINUE
00 20 J»l»0
/i J>»0.0
JF« ,.EO.O) TNU«1)»TN(J)
>>TN( J>
20
16
CO 1*
CCM$T»CONST-0.5«(TN(1)
IF
wE
»e*n
GO TO
C4ROS
t CL- CHECK
SStSlSlO
ooc
•••«*•« I
tssjsssa
••••••«!
•••*»•»
'•lt-3 91P2. (ANT ( Jl t J»l
(Z( Jl
---- «£tn «JEST»flT 0»T* FRiOn S*VEO RESTART TAPE 10
8401 -F.fcIND 10
-f A 7 (10) t^EAP.NONTH. IDTE»I I»CL»CHECKiI»tC«HEOtCONST»CI«ETS»£T,
1 C>
sstssssa
SSSISflO
tsssttio
SilltSlO
R-FAOIlO) I »NT ( J) i J«l «Q>
••#0(10) (Z(J)tj«l«0)
8*92 "PINT 9155
PRINT 9tSl.L*E»ft«*0*TH.IOTE«II«CL»CM£CK,IRtLiHEOtCONSTtCltETS«ETt
SSSfSflO
sttstsio
sttsssio
stsstsio
SSfSSf tl
(TM J) • J=l »Q)
9153, (FN(J» tJ»l»0>
PfilNT 915»» (ANT GO TO 8700
1FUBESTR.E0.2.AMD. JSAVE.EQ.U GO TO 8700
I5C TO 8710
S
C— —
8700
I GO TO «73A
00 8733 tal.lMl
CiLL SKIPI5)
CONTINUE
IFtLL.EO.LLSTBT) GO TO »710
00 8T05 I»1.10
PEAOI5) IDUMYtKOAY
IF (EOF(S)iaeOO.8705
870S CONTINUE
IF (LLSTOT-K-OAT-1 H802.tf7lO.a70l
c-— SET INDICES OF YEARLY LOOP
8710 Itf
IF if A.'.'O.fSTOP) IMJ»MM$TOP
C— -— OATS ulTi-iN TOT*L "UN
7700 00 32 n«!L">«[M
»••••• i 2
SSSStf 11
• ••••• 12
••••••12
sststsi
SSSSttl
SlJJJtl
SlftSSl
sstsssi
873<.
8701
• ••••• J 2
sssssss*
••••••12
••••••12
SSSSISSl
>S1SffS7
ittsttt?
(ttttSfl
266
-------
t«0
DPF»M»O.O
OEFA«o»o.
ISTCT-O
titssstt)
•••••••4
•••••••9
CALL 'HPo«TE(STAwT . II I
C-—--NOTF TM»T TH[S PBOGWAM CAN ONLr «E WESTAHTgO ON FIBST OB
IFftOTF.FQ.l .<>». trm .EO.HS.OU.KF|.AG.EO.O.OB.IC»OP.F.0.3) 11.8500
U n,| 1 j.j.o
CALL ^ONUM> iffcOP.Df L*«.I.IM.D • 11» ILO, IMII
I
IFCKFLAG.EQ.O) GO TO
<>15l. VEAP,MONTH, IPTE, II.CL.CHECK, M.L.HEO. CONST,CIi
••••••16
••••••I
( TM( J) ,j«l .3)
9153. (FNIJI .«!•! .0)
POINT Q15*. (»NT ( J) ,J»1 »O)
sttssssa
•••••••1
•••••••1
•••••••1
esoo iF(ii.co.iLO) no TO 10
HEAD INITIAL O»TA FOB M£W HAV FPOM TAPP *. C»C£PT ON DAY ONE
IF (I I.EO.l.ANO.TEiR.EQ.l )GO TO 1
C
C
HEAO(*>(TNij),FN(j>,CL.CHECK,IR,L«HEO«ANT(ji ,j*i,o>
:- -INITIALIZE Mf.0 AT STAwT OF EACH OAY
1 IF(L.GT.APPS) GO TO 6330
!F(ii.F.O.TMC(L>) 30.3*
6330 IFtHED.LT.0.01.AND.CNA.LE.0.05) GO TO 3*
*CNA
CCCCCCCl
••••••15
CCCCCCCl
CCCCCCCl
GC TO 6332
30
6332
C IF(HEO.LE.O.O)HED»AMT(L)
t»L*l
IP>100.
IR'1000.
102.II.HEO
••••••13
••••••15
••••••10
C—... TIME INTERVALS WITHIN EACH OAY LOOP
3* 00 21 J«ltQ
21 TOIJ)»TH(j)
C- COMPUTE SIZE OF TIME INTF*VAL« OELT
1*1*1
0€LTO«OELT
(DELX«0.03S/IS.OELTM)
IF(HEO.GT.0.01.ANO.KSATO*OELT.GT.hEO> OELT«HEO/KSATO
!F(X.OELT.LE!O.1)00 TO *
OELT«0.1-X
* CONTINUE
X«X*DELT
XT««T»DEL7»10.»«t-10.)
Y.0.7
C"——EXAMINE UPPPB SOUNOARY CONDITIONS
IF
CALL pBOP (ANT(2).TN(2),TD«K«2),0(2))
K(1)«(K(]).K(2)1/2.
CALL AOIF IANTH ) .ANT 12* .o< n )
60 TO (Sfi6
670 CALL APBQP1ANTU I .Z«
-------
|F <0 I 1 I .Lf.O .
6ft* t ( 1 ) «0.0
GO TO IP
667 f ( 1 ) «-•< f 1 I •Or|_»/0< 1 I
<3C TP l«
IT TIM ( 1 1 «Tor ooc
e d)»o.o
Tni t t »Tni£> i /;.
IF < tPOPT ,NF . \ i nn TO */l
f»U **TIP t *»jf i 1 1 « TNI I i • TO«I« ( 1 1 »ni 1 1 i
r'»(.| t»l»pr ( ANT ( / I . f* I f' > « III jH I •'» f |) 1/1 I
« I I I • I n I I I • r ( Jt I I / / .
C*U •0|l> I »Nf I 1 I • *HT < »D *P(J)
TN( J«1)»TN( J>
ANTej»«TN(J)«y»(TN(J>-TOlJM
IF(ANT«J).GT.TS»AMT(J)=TS
tF(ANT(JI .LT.TO»ANT(J)»ro
ANTU»H*ANT(J)
IF(BBC.EO.TM) rN(OI»TM
d«0-l
268
-------
*• TNIJI «F (ji »TN(J» I i »F< ji
»NT(J)«TMJI »r» ( TN< JI-TO(J) I
IF i ANT ( j> .<;r.Tsi*NTiji»rs
IF < ANT «|K(jl-(OlJl»(TN(J.ll.TO(J»l)-TN(JI-rO(jll /'(^.•OEL* I I » »nfl T
*O«FSI ji /CELT
i [
SF I j)"SF I J) «F>t< J)
tF»TNCj».LT.O.O.OR.TNt JI.ST.TSl ISTCT- t STCT* 1 • ••••••!»
J»J-1
IF IJ.GT.OIGO TO 48
IF«J.GF.7) CL3«CLJ»FN<7)
£T$«ETS*FN< 11
IF(FN( 1 ) .Lf .O.O.ANO.HCO.LE.O.Ol) 00 TO 8793
ct-ci»FX(n
HEO»MEO-FN( 1)
If (HEO.Lf ,0.0)MED»0.0 ••••••10
SO TO ?3 ••••••13
CM»CNI-FNI1) ««»«»«13
C^^-CHi-FMl) ••••••13
(JO TO ?3 ••••••13
CONTINUE
C— — »BIT? ON T*P£ s 0* PSINT THCTA AND FLU* AT 0.1 OAr INTERVALS
IFCX.LT.0.1 JfiO TO 2
•BITE (5) YEAff, t I iXT.CI •CL'HEOtCTStOEFAHC. ( J. TN ( J) « Z ( Jl i SF ( J) i
1 U(JI.J«ltQ>
IF(ITENTH.NE.l) GO TO 9713 OOC
PRINT 103 •••••••«
PRINT 121,re»». t l«XT.MONTM.IOT£.CL«CH€CKi£TStET,01F,CNA.CNl tMEO. OOC
1 CI«OEF4MCt[ OOC
PBIKT IDS. (TN(J) , J«1.Q) •••••••*
PRINT91M.CL3
«161 FORMATtkX.'CL AT3.5 FEET THIS TENTH OAV«.F10.3)
IFIISTCT.NE.OI o«INT <»166.ISTCT •••••••«
9U6 FOHVAT (10X«»UNSTA«LE SOLUTION SITUATION ENCOUNTERED •« I8t «
1ES THIS TENTM o*v»>
IF(ICTOF.NE.O) PRINT 9170. ICTDF.OEFAH
9170 FORMAT dox. ^DEFICIT MOISTUPE SITUATION ENCOUNTERED «t I8t •
IS THIS TENTH DAY. AMOUNT is •» F6.2t • CH»> -
•71? CONTINUE ••••••1Z
ISTCTO«ISTCTO*ISTCT
ISTCT-0 _
ICTOrO«ICTOFO»ICTOF .
ICTOF.O
OEFAMO*OEFAMO»0£FAM
DEFAM.0.0
DO 6 J«1.Q
6 SF(J)»0.0
2 IFtXT.LT. 1.0)60 TO 3*
c. ---- COMPUTE -CHECK" TO VERIFY CL
CONST1«0.0
DO 19 J«1.0
19 CONSTI»CONST1»TN(JI
CONST l»CONST 1-0.5* (TNI I )*TN«1»
OIF» ( CONS Tl -CONS T)»OEUX
CHECK-ETS-OIF-ET
C-— -*rR[TE FIN-L VALUES FOR LAST (I) IN DAY (II) AS INPUTS FOR DAY «I|»1)
C t»PITE (*>( TN(J) .FN(J) .CL. CHECK. IP, L. MED. ANT(J) tj»l ,0)
O... -PRJNT ONE OF TtaO OPTIONS FOR DAILY OUTPUT
IF
-------
. inU . a .rMrrn,£Ts.r».oir.CNA.CNI
I 0 I %TCTO
91*7 FQPMAT ( tO«.«"JNSTA«Lf: SOLUTION SITUATION ENCOUNTERED •• 18,
1ES THIS PAY*!
IFdCl
33
ooc
OOC
•••••••9
TIM•••••••9
•••••••9
tttttttn
1ES
. AMOUNT (S ».
• CM«)
••••••12
n»tr-i
C—- POINT HFSTAPT DATA »T END OF
irdnTE.EQ.l.OB.IOTE.£0.161 GO TO 4721
PPINT 9151, vEAR,MONTH.IOTE«ll.CL'CHECK.IR.U.HED.CONST,CI
9152. (TMJJI.j«1.0i
P»lN.| 915*. ( ANTtJ),J»1«0>
PRINT 9172* (ZtJ).J*l»0)
917? FC»/»AT ( X. • Z (J) t J» 1'O4/. 7(9£U,6«/>>
C-... INCBEHENT YEAW. SET DAY NUMBERS. CALCULATE MOISTURE DEFICIT, RE»0
C-—- laPfGATION APPLICATIONS. »E5£T
c—- COMPUTE MOISTU»E DEFICIT
IFdOEF.NE.l) "50 TO 8738
SUHDEFaO.O
SSSSffS3
SSSS1SS3
SIIS1SI3
ISSSSS13
SSSSSSS3
ff ISSSS8
SffSfSSS
00 8736 J«2>0
IF L»1'APPS
AMT(L)=AMT(L)»2.5*
873«)
KFLAG=0
GO TO T700
9300 If (tPUNrH.E'3,2) 00 TO
IFtJPUNCH.NE.1) GO TO 99
C—— PUNCH PESTAPT DATA AT END OF RUN
PUNCH 91*1. YEAR.MONTH. tOTf.n.Cl..CHECK-
PUNCH •J«l«0)
PUNCH 9l«j, |ANTu>,J»l.0)
PUNCH 91<»2» (Z«J)«J«l«0>
GO TP 99
9181 roflmTUS. 13. 13* !»• 3F13.6. 13. £13.6)
9162 FOPMAT (
-------
idi i TN iji . j« I .(ji tftttllO
iiii i > -I ( ,j i . j« I .'J i tlttttlO
lOi i »N' i .;) t J- I .iji
ini i i < ji . j» i ,iji
I.F 10
STOP
n«oi roo««r(s«, • END 'ILF FOUND ON TAPF 5 MCFOPE OAY NO, LLSTRT-I FOUNOOC
10.- f «EfijTioN TrwurwATEn •> ooc
r,o TO QO ••••••12
a*o2 PRINT ««c>3 ••••••u
»803 FgH»AT(5x. • OAY »t»0 F»OM T4PE 5 Eau*tS OB tS SPEATFB THAN ST*PTI»»«»»»12
ING C»Y, EXECUTION T£OMlNAT.EO *l ••••••12
GO TO <»g ••••••12
C--.--POIST QUN PAPAMCTEAS AND INITIAL CONDITIONS.
10 KFLAfial , •••••••!
IF(AA.FQ.1)<5.12 ....... 1
9 PPI*T 10Q
110
91«Sfl. I PUNCH. I RES TR, I SAVE t I TENTH, INF II.5.H.STRT ,MMSTOP« I STOP. tl$J«««l
IPOPT "00 1
PRINT 120, ( T"E(JI .MQN(J) ,OATE{J> .AMT(J) .AOENT(j) ,J»1 ,APPS>
PRINT 10S, (TM( JI ,J»1.Q)
127
107
PRINT Hl.AA,8P.CC,LL»»"*»dBC.T8C.YF»R. CPOP.-. APPS.PELX «TS. TM.TO.
IS*
PBINT 113
PRINT 11*. I IPENT.MQfi (N) ,Nsl ,0f
IFICOOP.NE.3.00.ICflOP.Eii.3)GO TO 12
PRINT 115
POINT llA
P«INT 11?
PRINT 108
60 TO 12
100 FC«**T(5I5.2F5. 0.315, 5<»5F5.0> - OOC
101 FORHAT(2X, US.1F10.0)
102 FCRM4T(/,2X.«W4TEB APPLIED. DAY NUM8EP •.!*,». AMOUNT « »F7.2,
!• C".»>
103 FORMAT
io« FORMAT i inn
109 FORMAT ( iwi ,3x .•PAPAMETESS, CONSTANTS. AND INITIAL CONDITIONS USED
UN THIS REPORT.*)
110 FORMAT(/,3X,»- ---- NOTE ----- DIFFUSIVITV AND CONDUCTIVITY RELATIONS
IHIPS uusT BE IVSEBTEO INTO SOU»CE DECK.*./)
Ill FORMAT |9X.5I5«2FS.2,*I5.5F5.2t
112 FORMAT (/,7x.»tNtTi AL sou MOISTURE CONDITIONS, •»
113 FORMAT!/. 7X. 'SOIL IDENTIFICATION AND HORIZON DEPTHS. •!
1U FORMAT (<)X.*IDENTIFICATION« «.A8,». DEPTH" »,F5.1I
US FORMAT «/,7x .»CONSUMPT IVE use OATA,»J
H6 FCRMAT(
119 FORMAT I/,7X.»WATER APPLICATION DAYS. DATES. AND AMOUNTS. •>
120 FORMAT (Qx,»DAY NUMBER*. I*. 7X, "DATE ••t2»»/*«I2«7x»«AMOUNT«».F6.2«»
I CM. SOURCE » ».A1>
121 FORMAT (X. 2[«. Fft. 3, 13. l«. X. 10F10.*. 16) OOC
271
-------
i|10.F-».A.llO,F20.4l
l?T FQRMt T I / , ?> • •>*VlN °»U»MFTFi»S. »N() H') IM TMW« ( i / . >. *i»^TAuf n*f4 > OM "f *" •« I "S. • MONTH • 12, * OAV • OOC
^ •» rii.*« • L« •« IJ« • NFH« •« fit.A/. t. • rnNST» • (li.is. • ctitiitii
II* •• Ml * i • *TS« *iri1.>i, • f • •» Ml,*/ ». • CNA« •» EM.*. Mtlttll
* • CNt> • Ml.«i« • 0''*"C" •• C l.A) tttVttf.4
OIHJ rpnuati » • TNI j> ..)• I >u*/. M>Jr ).*,/ii •••••••I
rruwtfi i • PN i j i • j« i • ij» / • 'c*f l.n«/>i •••••••I
rct)MAri i »4N r iji . j • i ,ij»/. ri^i I.^./M •••••••)
FQHM»T|/,«, .uf5T«Wt l)«r« «S HEAD F«QM CAPOS*) •••••••!
fORMAr i Mis«>^.o•is1131
FOMMtTl/, x.* HCPLETEU MOISTURE FOP VEAP«*. IS, • MONTH«*tI3< • OADOC
IT£««. 13. • IOAY NUMftgQ»»,14, •) is "«.F5.2. • INCHES*) fSSSStSl
IPESTP ISAVE ITENTH INFIL5 LLSTPtitJJ$$l
IT MMSTOP ISTOP IDEF FCAP IPOPT»/t 18. 8110. *QO I
2 F12.3. IH> MOO 1
StBHOUTJNE CONUSE(CPOP.OELX•J»U«K.LL»MM)
C——CONSUMPTIVE USE SURWOUTtNE
COMMON/ XYZ/ I DT£ ,M0^4TM.UH, Kf
COMMON/CHCK/I CHECK,ICP^P
DIMENSION wEANT1112) ,MEANT2(12>tKAl(12l ,KA2(12) .KB 1(12) ,K82(12)«
1P1(12),P2(12)«"P(>|),KPI(6).KP2(6.3).KP3I6).U1 (36«) .UH(366)
DIMENSION FACTOP(24) •••••••2
WEAL MEANT! ,uEANT2,KAl.<(A2.K81.KB2>KP«KPl tKP2iKP3
C——THIS ROUTINE 9ETUPNS CONSUMPTIVE USE IN CM WATER/DEPTH SEGMENT/
C DELTA TIME
C--— ICROP*3 USED KMEN DAILY CONSUMPTIVE USE VALUES ARE READ
C ——MEAN TFMPERATUBES (F) FOP FIRST HALf OF MONTH
DATA! (MEANT1 (I) «I»1 • 12>*
f«2.4t47.0.46.3,55.a,66.6«71.9.79.1,81."5.75.0,65.0.55.0,45.0),
C——MEAN TEMPEPATUBES (F) FOP SECOND HALF OF MONTH .
1((MJ4NT2(I) .[ = 1.121*
t47.0,»fl.l,5*.8.62.5.73.9,73.8.82.6.77.6.75.0,65.0,55.0.41.6).
C— CONSUMOTIVE USE CONSTANTS FOR 8ARLEV (FIRST HALF OF MONTH)
f< IKAl( I) .1*1•12>»
$0.50.1. 58.1.17.3.00.1.22.0.4*.0.00,0.00,0.00.O.Ofl.0.00,0.12),
C CONSUMPTIVE USE CONSTANTS FOR BABLEY (SECOND HALF OF MONTH)
f
$0.00,0.00,0.00,0.00.0.00.0.1*,1.29,2.20,2.20,0.10.0.10.0.10).
C---—CONSUMPTIVE USE CONSTANTS FOR MlLO (SECOND HALF OF MONTH)
SO.00.0.00,0.00.0.00.0.00,0.74,1.67.2.20,2.20,0.10.0.10,0.10),
c—-—PERCENT OF O»VLIGHT HOURS FOR FIRST HALF OF MONTH
t UP1(I)«I»1.12I«
C——PERCENT OF DAYLIGHT HOURS FOR SECOND HALF OF MONTH
SIIP2111«I*1• 12)»
C—.—ROOT OISTPIBUTION WITH DEPTH FOR BARLEY
DATA((KP(I).1*1.61>0.«0,0.24.0.19,0.13,0.04,0.00)
C——ROOT DISTRIBUTION *ITH DEPTH FQR M[LO
DATAIIKP1 i I), 1*1,6)«0.31,0.22.0.14.0.09,0.08.0.0*1
DATA ( ICHFCK»0)
DATA (F*CIOH*15.,16.,tS.,13..15.,16.,15..15..l5..l6.,15..15..l5.. •••••••2
U6..15..16.«1S..15..15.,16..15..15.,15.,16.1 •••••••2
C
C
COMPUTE DEPTH IN FEET
J-1 1/30.48
MOO
HEAD U OFF DATA CAPOS IF CROP* 3.
IF (CHOP. NE. 3)00 TO 11
IF(J.NE.2IGO TO 9
ICOUNT»MONTM»2-0.5
272
-------
.o .s
If I l
IF (-C"ECi'.ro. 1 ) GO TO 7000
READ 100.
00 11 I»l.6
JCHCCK-1 fSSSSSSl
CALL AJST(KP1.0ELJ.AOJUSTI MQO 2
7000 CCNTINUF
IF(ICROP.EU.3.AND.CtfQP.EQ.lt GO TO 201
»f AO lo?. i CODE. iYEAH.(ui •••••••!
FQRM4T (?IS.12FS.OI •••••••!
PPJNT 103.ICOOE.ttEAP ••«••••!
103 FCRM4T(7X .»CnNSUMPTIvE USE CONSTANTS RFAO IN FROM OATA CARDS IN IN»»«»»»»1
1CHES °EB SEMIMONTHLY PERIOD. ICODE**. 15. • IYCAR» •» 15) •••••••!
PPINT 10*.(Ul I I> . I»l .24) •••••••!
10* FORMAT (H*.?4F3.2> ••«•«••!
DC !•» 1*1.2»" ••«••••!
UllP'Ul (I)»2.54 •••••••!
tiH(i)*ut(n ••«••••!
19 CONTINUE •••••••1
20 U«Ul(ICOUNTi
GO TO 21
' 201 CONTINUE
C READ NU^HER OF DAYS OF ET
READ 106.ICOOE. IVEAP
1A6 FORMAT(2I5)
READ ios«(ui»i«i<
105 FORMAT<1AF5.3)
PRINT 107.
107 FORMAT<7X.»CONSUMPTIVE USE CONSTANTS READ IN FROM DATA CARDS IN IN
1CHES PER DAY. ICODE« *.I5. • IYEAR- »,I5)
PRINT 10«.(Ul(I).I«1.LMI
10R FOP.M4T (1 IX, 10FS.3)
7001 00 7002 I*1>LM
Ul Ct»«Ul (M«2.5*>
700? CONTINUE
21 CONTINUE
IFIICROP.NE.3I U»U1(ICOUNT)
|F(ICBOP.EQ.3.ANO.CROO.£0.3) U»U
GO TO 7
11 00 12 1*1.6
KP2(I«1) * «P(H
12 KP2(I.2) r P2(ll/100.)»2.5»
C«—ADJUST CONSUMPTIVE USE FOB LENGTH OF TIME INTERVAL
7 CCNTINUF
IFdCROP.EO.3 .AND.CROP.FO.31 GO TO 200 -
273
-------
IM inir ,iF . i«• i •'
« T * | 5
HO TO
• A « |S
.—.--•rji'ST r
oti,i«,
200 CCNT|NI)F
tvi
or
«OOT
if ( IFOOT.GT.M u»o.o
IF ( IFOOT.Lf .*> U«U»hP?( IFOOT.CBOPI ••OJUST
uf TURN
too ro«««*T (f>r 10,0)
CKO
• ••••*•;>
•••••••»'
•••• •••£
*oo
MQO
INTEGER FUNCTION DAY
2
^
4
5
6
7
10
II
12
INTEGER ^JNCTION OAT IK,L»M'
GO TO «I .2«3«*.5i6t7.8.9.10»1 I.12) M
BETUPN
RETURN
RETURN
RETURN
RETURN
RETURN
RETURN
RETURN
RETURN
END
7
A
q
10
11
1?
IF1M.GT.31.
SUBROUTINE THEDATE
SUBROUTINE T"EOATE (K.D
CCl*MON/XYZ/ir>T£tMONTH.UH.KP3
M»K«L
'.LE.311 GO TO
.M.LE.59) GO TO
... I.AND.M.LE.90) GO TO
IFIM.GT.90.ANO.M.LE.120) GO TO
IFJM.GT.120.AND.M.LE.151)GO TO
IF(M.GT.151.AND.M.LE.1«1IGO TO
IFCM.GT.191.AND.M.LE.212)00 TO
!F(M.GT.212.ANO.M.LE.2*3)GO TO
IF(M.GT.2*3.ANO.M.LE.273)00 TO
IF(M.GT.273.AND.M.LE.30*)GO TO 10
IFIM.GT.304.ANO.M.LE.334IGO TO 11
IF
-------
SUBROUTINE CHAR
CMAU CMAR 10
r---- CHAN *0
C — — P*or.RAM TO «JK»n INOUTS *no COM*"
-------
11 ««0.0
GO TO 100
10 IPlZx.fiF. TSI 00 TO 20
IF if* .LT.S»> I/.L * .%HI l»0 T.T 111
0»C/» < I/-SW) ''C*!
nn TO IUOD
120 o-or>«i*T
\r toTS.ro. TSI
IOOQ RETURN
END
loo
PWOP ITO
PROP ito
PROP 220
PROP
"HOP 250
SUBROUTINE AJST
SUBROUTINE AJST(KP3.0ELX.ADJUST)
.... PROGRAM TO COMPUTE ADJUSTMENT FACTOR FOR CONSUMPTIVE US€
DIMENSION KP3I6J
COM»«ON/AJ5T/0
INTEGER Q
REAL «R3
U»U-KP3(1)»DELX
.—. CONVERT NODE SPACING FROM CM TO FT
OEL»OELX/30.*a
UxO.O
ENTER LOOP FOR ALL NODES
00 100 I»2.Q
COMPUTE DEPTH OF NODE IN FT.ET
D«OEt»
COMPUTE SUBSCRIPT FOR PROPER 0£PTn ZONf
IFOOT«0»1
TEST IF BEYOND ROOT ZONE
IFHFOOT.GT.M GO TO 110
COMPUTE RELATIVE CONUSE REMOVED
U»U*KP3( IFOOT)»DEL
CONTINUE
ADJUST CONUSE FOP REMOVAL AT TOP NODE
IFQOT«OEL»1
U»U«*P3UFOOT>»OEL*0.5
ADJUST»1./U
RETURN
END
C
C
C
C ----
100
110
10
20
30
*0
SO
60
TO
ao
90
AJsr
»JST
»JST
AJST
AJST
AJST
AJST
AJST
AJST
«jST 100
AJST HO
AJST 120
»JST 130
AJST 1*0
4JST 150
AJST 160
AJST 170
AJST ISO
AJST 190
»JST 200
AJST 210
AJST 220
AJST 230
AJST 2*0
AJST 250
AJST 260
AJST 270
AJST aso
AJST 290
276
-------
SUBROUTINE AD/F
iOlF< TZ,Tx,0>
T. Al_"iAlHINF •
If ( T7.nT.SH.AND.??.I T . rS.ANO. *• .1 ! . %«l > HO TO »
ir i T» .lit .su.ANfi. r* .t r .TS.ANO. f /.i ? .SM> no TO t
if c T/.KO. TS.ANU. r» .c>:.:•>«> 10 ro «.
ir (T«.r(j.rs.iwn.r/.Lt.swi ao TO ft
ic TO ip
S T7»TS- .0001
r,o TO *
* TK-TS-O.OOOI
r,c TO T
» C4LLSlMPrrZ.SH.AVO)
0«»VO/(T7-TK)
GO TO 100
7 CALl.SIMplTK.SR,ftVD>
0«»VO/(TK-TZ) ;
CO Tn 100
10 ir(T2.GT.SP.ANO.T?.LT.TS.ANO.TK.GT.5B.ANO.TK.LT.TS) 60 TO 12
GO TO 32
12 IFCTZ.GT.T") GO TO 15
IF
-------
SUBROUTINE SIMP
111 I I Sf si "V ( W . s
• •Hi. /(T5-SB !)••(<»l.M»KSAT«f.MP»AI<»INF
H«(l./Ml./lAi.e«rS)i««£M«>H« > •• < BE T •£•«?/ ALP > I
SECTION COMPUTES INTEGRAL VALUES USING SIMPSON
M«(SZ-SK )/2.o
ThO'O.G
FCUB«FU
OLOI»'T»
C EVALUATION LOOP
00 2<« I»1.N
C CHECK FOB COMVFRr,£NCF 0» EXCESSIVE NU.-9ER OF ITERATIONS
IFt»PS(OLOl^T-INTEG) .LT.l .OE-6.0W .N.GT.10000) 60 TO 10
OLOINT«INTEG
00 TO 25
10 SVO«CONST»IHTF.G
RETUPN
PROGRAM USCHEM
PROGPAH USCHf.M
C 1
C»-.»ONSATURATED CHE>«ISTBV PROGRAM uSBH VERSION 1.2.0—NOV 197*
C 3
4
DIMENSION XI7.25I 5
6
7
COMHON/erP4S/NPYP»S.IDYSTB.IOYSTP.ILOtIHItINFILL 1CONT1iJPAS 9
COMMQN/APLE/TITLE(10).SMQNTH.MM.O.IPPINT.JPHINTtINK,lPUNCHtlSTOPt 9
1ITEST.IPEAOP.IMASS.IAOO(2S)tIORNAP(25) .M0rt(9) ,TOTN(9-»I, YEA* , 10
ZAIRRC91 ,[OR(25>,TT(f>0)
3ITOT.JTOT.IOTOT.NT
CCMMQN/XX2/A1tA2.*3>x 13
COMKON/AFG/FNH3.II.LLL«IOP.ANETLIM(2S» 14
co**ON/rYY/START. IOTE.MO»»TW. j. ILAP
COMKON/XXr/KHECK.ICOUNT.CONV.PK.P*! • CROP, f ACT I ft
COMMON/xxx/OELX.OELT,MS.wTAfir,fiD<25 I.TEN(25 ).CHECK(25 I.MOISIN 17
KZS )»CMM20U2I> i.MOISOUT<25 ».AN03(25 >.ANM3<25 I.UHEAI25 I ,OPN 18
2(25 I»CA(25 I.ANAI25 ).AMQi25 1.MC03I25 ltCU(2S )>C03(2S ».SO*l25 19
3)iE5(2S ).CS(2b l,SA5(25 ).XX5(25 I.CASO(25 l.AQso(25 ),BNM*«25 ). 20
*EC«25 >,CNl(?-s ).SAMT(25 ),»N(25 i.PC<25 )tTEM(25 ).C*L(25 J.OtSPO 21
1P.XTBACT.SUMN03.THOP. («> •TO
-------
I so«ouT . amur* . ! ] n i ?•> > . A/F i ."11 «ri> 1 1 o »
«> > « 4K c« < i** >
COMMON «TPN[f'U(^^l . AC TfA ( /7)
31
32
REWIND 3 SHEwlND A trttv«lMD 9 JREWJNO 10 33
3*
SET INJTUC V4LUES 35
ICMECK » 0 %CMM201(U « 1.0 36
JPAS»0 37
OC 6^3 fj=1.7 3«
00
c ----- BEAD TITLE CASO **
BEAD 8 GO TO 9012 *2
00 «OU !»1
SKIPI2)
CCNTIMUF
9012 IFIICONT2.E0.1) GO TO 901*
SSKIP*lLO-tDYSTR
IFINSKIP.LE.OI GO TO 901* 6'
00 <5013 tsl.^SKlP 6*
READ (2) U • **
^013 CONTINUE I,
901* CCNTIMUF - i*
C—— COMPUTE NO. OF TIME INTERVALS PER DAY ~\
LLL» l.'/DELT .0.5 "
IF(IP«INTI.NF,OI CALL P»NT ( IPR INTI . IPRINT J) '.
IF (NPVPAS.EO.l ) GO TO 9000 15
--— -READ TEMPERATURE HORIZON DEPTHS II
BEAD 107. (TnOP (J) .J*l . T0> I,
9000 CONTINUE f?
-— — BEAD COMPONEMT HORIZON OFPTHJ . •*
READ 100. (HOR(J) «J»1>0) ?
IF JNBVPAS.EQ.l > fiO TO "001
c ----- STORE TEMPERATURE PROFILE DATA ON TAPE e
00 800 J»1.NT
PEAD aoi.iTTin .1*1. TOI
279
.**
-------
• H[TP (A) ( TT ( I ) . I« I .TO I
r-----Ufifl
ur»n ion. AMM H I i .ANOII i i .r» i | i . ANA t l 1 , ANII< | i .MCO »( 1 1 tCL« I I tC <^
-srnur ruANS»i»«M»a |ow|.»»r ION **rnj ANAI VMS <»r
fiumi | i ^«MT i i > .u . 0 >»M
0*11 t'N| TS| I I > •»«
t I I ••NW » I 1 • «* I >JU I .» I «AN
TOT*L MUMOEP or COMPONENT HORIZONS io»
(0) /OEL<« I . 1
Sf 4DHOO. (SEBAT IO(N9) ,U(M9) . ACTCA ( N9 ) »N9«2i'J)
not. (ISETNIN?) «N9«2.oi
^92.783 10ft
"EAO 78*« (C«Haoi < J) »MOISlNtJ> .MQISOUTlj) «TE»*(J) »U(J) »J»l«0) 107
IOB
C — .,-PBlNT HEADING - 109
793 IF(ieERUN.EQ.O) POINT 201 110
111
C— -- -SET COUNTERS • 112
11*
C— --- CALL OUTPT TO ZEBO INITIAL VALUES US
CALL OUTPT(Kl) 116
1F(IO£RUN.EO.O»22.7D1 ' ' 117
US
C ----- (JEAO INITIAL SOIL ANALYSES 119
22 PEAO 100<*NH3<1) .AN03I 1) .UREA! 1 ) .CA< 1 t .ANA(1 ) ,AMG(1) .HC03I 1 120
H.CLtl»»C03(H tSO*(l) .EC(1) tXXSCl) .CAL(l) »BO ( I ) .SAMT ( 1) .CNl (1) 121
123
C— .-PRINT INITIAL SOIL ANALYSES 12*
PRINT 200.L.*NH3(1 ) .AN03 ( 1 ) «(JP.£AU I . SAMT ( 1 > ,CA ( I ) tANA(l) .AMG(l) . 125
1HC03(1) tCL»C03 .S0*(l) 126
BEAD lol.xTPCTd) tPC02il) .AKCSJD .AKCMID
127
C—— COMPUTE SEGMENT NUMBER OF COMPONENT HORIZON - 12H
130
C --- —STORE INITIAL SOIL ANALYSES IN PROPER COMPONENT ARRAYS 131
00 23 J«N,KK - 132
A*M3(J)«ANn3(l> fAN03(J)«AN03(l) $UBEA(J)«UHEA( I) 133
CA(J)*CA(1) »ANA( Jl «ANA( 1 ) $ AMQ { J ) « AMQ ( 1 ) 134
HC03(J)*HC03(U JCL IC03 « Jl »C03 ( I ) 135
S04(J)aSO»(l) SECIJ)»EC(D 1XX5
-------
C-----FCP « apauN. "FAO *»0* fAPFl QP F°OM CAPOS
IF I IRFACP.Fd.Ul IS?
| .4N03I-JI .URR4 « J)
REWIND 10
IF{NTEMPIN.E«.OI GO TO 523
IF (NRVP»S.£J. I ) GO TO 900*
C— ..-SPACE TAPE8 FOPEWASO TM£ PROPER NO. Of RECORDS
00 510 I.I.NTEHPIN
510 READ (8)
900* CONTINUE * .
ir i tor APP.AIF .0 i
|nr»n *(>*, i COMNT .NFJ u r tsi,*ouiiiN.*rr*»|N. i »NM 1 1 n i *NOI < ,11
i > ANA i M < 4*n i ,ii ,H« •)! i .ji i.rMii 1 1 .(IM«I i ii.«N(ji.«rijit»>i ji.c^iji.SASut.CASii I1**1
> ( j l . An-.n i ii tHiMM* i j> . i t we T i j i , *Nf n. IM u i ,»;» ( j) . 1 1 * (j) , 16*
4f>( O/l II . »H(. N( Jl f AKCX I.JI » J«?tUI
C.....SFI INITTAI. VALUES
703 00 1 J-?.Q
IF ( IHERUN.EO.O) 780.T81
780 OCN(J) -we (J) »«JN(J) "CHECK (J) «0.0
$•5 ( J) »RNM* ( j) «o .0
CKM201 (j) »*TOCT ( j) •«0( J) »OELX
IT*
CC(J)"€C(J)/l.E5
C— --- CAUL U*!T CO^JVEMSION SU8POUTINE
CALL UNITSHJ) ' I78
179
O— -POINT TOANSPORMEO DATA I"
POINT 200. J.ANH3 .UREA tj) .SAMT (j) tCA (j) .ANA (j) .AM<;cj) . lfll
1«C03( J) «CL( J) «C03( J) .50*1 J) la2
781 IF (10EPUN.E0.1) CH£CK(J)«l.O 1«3
1 CONTINUE lfl*
185
c- ---- OEAO FEPTILIZEO APPLICATION DATES !••
BEAD 10*. ITOT, ( UOO(K) .K«l . ITOT) 1^'
188
IF(NRYPA5.EQ.l) GO TO 9002 I**
C ----- BEAD 0»GANIC-N APPLICATION DATES . 1*°
P-EAO 10*.JTOT. ( IQPNAP(K) ,K«1 , JTQT) I'1
192
900? CONTINUF I*3
C— — P-EAO IBPIGATION WATER APPLICATION DATES I9*
READ 104.IRTQT. ( IPR(K) tKml , IRTOT) l"
196
C—— STOPE FEPTILIZE9 APPLICATIONS ON TAPE 9 I*7
DO «02 I«I.ITOT lq*
READ 100.
«03 CONTINUE
9003 CONTINUE
C- ---- SET SEGHFNT ONE VALUES EOUAL TO 2E»0
16 *NH3a)«AN03(l)«C*a)»4NA«0.0
IFdPEOUN.NE. 0)508, 720
508 REWIND S 215
281
-------
522 IF (NrE"»TTN.EO.O» GO TO S10
221
C-— ~SP»CF TAPF.Q FOOE«AHO Tut PBOPF.B NO. OF RF.COHOS ?2«
DC 511 I»l.^rCortN 229
511 "EAO (<» 330
550 If (NDBGIN.EO.OI GO TO S1J 232
2J3
IF CMqYcis.Fr). n rfo TO <»005 23*
C.-.—SPACF T»PF.IO FO«FWAI»O TMF PPOPTH NO. or RECORDS
DO SI? I»1.NO»OIN
*00«» CDNMNUr 2J9
00 Tn M ] 240
« 242
10 2*3
* NO»OTN « NTE»«PIN « 0 2«*
513 CONTINUE
tSWO » 1 2*6
If {(PBlNTj.Nf .01 CALL PPNTI i IPBINTI.IPHINTJJ 2*T
c*— — CALL suflpouriNE TO EXECUTE PPOG»»« FOP. EACH OAILV TIME INTERVAL 2*9
CALL EXFCUTE 250
251
'C— — CHECK fQ° END OF »UN 252-
CNOFILE 2 253
IFIMOOUDAVtlDYSTP) .EQ.0 1726.721 25*
T26 IFlYCAR.ea.ISTOP) SO TO 721 255
256
c — — PESET COUNTERS 257
ICOUNT a o $YE*R * *EA* • 1 *LL » i 25a
ILO-IOYSTB 259
ILAP " ILO
IHI.IOVSTP 260
IFCVEAP.rO. ISTOP) IHI»HM 261
IF( ICONT1 .EQ.O) GO TO 720
PEMND 10
C— - READ IPP.I&ATION WATEH APPLICATION DATES FOP NEXT ?EAP
«EAO 10».I»TOT, ( IBS(K> .K-l . IBTOT)
C— — REAC LAST OH<;ANIC-N APPLICATION FOP NEXT YEAH
READ 100. (OFEWTUI .J«1.31
IOCOU « JTOT - I
JPAS • 0
00 1331 I=1.IOCOO
1321 RCAO (10)
. (101 (OFEPT(J) .J«1.3)
10
GO TO 720 264
721 CONTINUF 265
C 721 £».OFILE 2 266
C EKOFILE 15 267
NTEwPIN»*TEMP[N-l 268
1COUNT»ICOUNT-1 269
270
C— — -EITWE« PUNCH A flEPUN DECK OR wPITE REPUN (PESTAHT) DATA ON TAPC3 271
IF (lOAY.EO.IOYSTP) ICOUNT « NFERTIN « NOR6IN * NTE"PIN • 0 272
IFUPUNCN.EO.O) 502.503 2/3
502 PE»INO 3 274
•RITE (3) ICOUNT. NFEPTfN, NORGIN. NTE^PIN, (ANH3IJ) .AN»03(J) .UREA(J) 275
l.CA(J) .AN AC J) .AMG(J) .MC03I J) .CL(J) «C03(J) .S04 «AKCS(J>
GO TO 561 279
503 PUNCH 505. ICOUNT, NFE»TIN, NOHGIN. NTEH»IN, (ANH3IJ) .AN031J) .UPEA(J) 290
1»CAIJI .ANA(J» .A- .XX5(J) .CAL( 281
- 2J) t80(J) .SA"T(j) .CNl (J> «ORNjj>,»N(j» .«C(J> .E5( J) *C5( jl , SA5IJ) .CASO 282
3(J» .AGSO(J) ,ONH4( J) . XTBCT(J) , ANE TL IM ( J| . A2E ( J, . I IK { J) , 283
4Pco2u> (AKcsi ji . AKCHIJI ,j«z.ai
282
-------
C
C <>e«l"0 1 285
r---- POINT «TST»OT DATA 2*6
S»| PU[M UlOOi IOAY»YCAH. ICOUNT .NFFOT IN.NOBfllN.NTrnOIm.O 2BT
PHfNT U|Ql, (A*H 1 ( Jl .AN03 IJ» .U«t* U>
I«CA«C*L<
?JI*art .OUNIJI .rfNl.ll .or IJI «eS< Jl .C* 302
100 FO«MAT(16r5.0)
10* FCPMAT (1615) 30*
105 *"OR»-4TmFS. 0/1615/1615) 305
106 FO<"<«T<11I5I 306
107 FC»M4T(5F5.0l 30T
200 FCBM»T(l5.liri0.3l • 30«
201 FOPwtTt ////IX'IMTtAL SOIL *NAL ^SES IH€0/L OF SOIL EXTPACT1 — (0»G» 309
OF SOtLl •//2X»H2N« 310
311
312
202 FOH*«T(//U»TP*NSFOSMED SOIL ANALYSES (uG/SEGMfNT OF SOIL ) «^/2X»SFa 313
«HC03*8X»CL
505 FOBM»T(»I5/. (6E 1 3. 5/6E1 3.5/6E 13.5/6E1 3.5/«E 13.5 . I5/3E 13.5) I
101 FORMAT (F5.0.F10.0.2FS.O)
7(J» FCfixiT (5F10.0)
601 FCRM*T(2X.Ffl.0.7F10.0» 318
1100 FORM»T|3F5.0)
1101 FORM»T(80I1)
1*0
SUBROUTINE COMBINE
SUBROUTINE COMBINE!
COMBINE
COMBINE
COMMON/SALT/S£RATIO{25>.S8YPAS
COMMON/3/IFLBVPA.ISEGST •
COMMON/BYPAS/NBYPAS.IOYST«tIDYSTPtlLOtlhlt2DUM(3)
COMMON/SABLE/SUMS(3) •
COMMQN/EEE/PSUM.OIFNH4.0IFN03.TPLANT *
COMMON/XXY/ICHECKtICOUNT.ZOUMl(5)
COHMON/YYY/STAPT.IOTE««OMTH,III.LL °
COMMON/AFG/ENM3«II»LLL«I OP»ANETLI«(25) \
COMBINE
C——fhis SU8POUTINE C*LLS THE COMPUTATIONAL SUBROUTINES AND ASSEMBLES
C——THEIP OFLTA VALUES
COMMON/XXX/DFLX.OELT.MM.wTAPT.BO(25 ).TENI25 ).CnECK<25 I.M01SIN
1125 I«CMH201(25 I . NO [SOCL(25 I.C03I25 ».SO*(25
3)»E5(25 ).CS(2b ).SAS(2S ).xx5(25 liCASOl^S )>AGSOI2S )
4CCI25 l*CNli?5 i.SAnT(25 J.RNI25 >.RC(?5 ).TEM(25 )iCAL(25
1 P.SPACE(36).ISWCH.CUMSUM.SUMQuTfPEOoCE
COMMON/GIPL/UHEA1.UREA2.r>NH31«ONHJ2C031.S041,KKK.PPPP(«)
COMMON/C02/PC02I25).IPC02
283
-------
I*
IT
. of LMH J i £M .nf i o«»'*N( ?*> i .OELuor. A l<"> i tCXHCf>3 i «C «CO 1 1 ** > !•
. f «ci. < /* i , r <«M->» i ?s > .FI NOI i«»* > . FL *M 1 1 i* > «'I»JHF. A i ?
' '">' .'LA«iiif M «FI. HdM i JM »FLCL i?si «Fi.Cf)3H*.i .FI so*
iM.NM* (,»*.> .nr.LriNH* «,»*i , iNf M I;M . ANF 1? 1^*1 . ANCT)!/*,) « 41)0 I
liMtrr.ro r
• i.o
NO* • /
IF I l«*Fn*T.Mj. | I NOW • I
IFACT • PtnuCE 28
ISCT • IFACT . i IF « i.o 29
IFIII.CO.LLL) K»2 30
31
c-.—COMPUTE DELTA VALUES FOR EACH SOIL SEGMENT 32
SO 00 1 I*fcO*i3 33
34
C CALL SHUT-OFF SUBROUTINE 35
C CALL CHKitl tL2*L3< I .EXNH3 IT I .EXCA(H ,EXANA( I) • F.XAMG(I) tOELN03( I) . 36
C IOELNH33.4 . . 39
3 Ll»L2«L3«0 39
4 CONTINUE 40
C——SET A UNIT CONVERSION CONSTANT 42
CONVERT(I) » OELX'HOCII 43
C 44
C——IF so INDICATED 45
C ENTER SECTION TO COMPUTE AMOUNT OF LIME WHICH HAS PRECIPITATED 46
IF(IflP.EQ.O) GO TO 205 47
C 48
C——COMPUTE AMOUNT OF LlME IN SYSTEM EXCLUSIVE OF SOLID STATE 49
C——UMTS ARE UG OF CAC03 P£» SEGMENT OF SOIL 50
ASUMl - CA(I)*2.497 « CASO*BNH4 £X8NM*(I> » 0.0
134 IF(IOP.FO.O) GO TO 206
C 59
C——AGAIN COMPUTE LIME IN SYSTE* EXCLUSIVE OF SOLID STATE SO
ASUM2 « CA ANETLIM(I) * ASUMl - ASUM2 65
C 66
C- COMPUTE POROSITY OF SOIL SEGMENT* ASSUME PARTICLE DENSITY IS 2.65 67
C 6<»
C——COMPUTE UG OF CAC03 WHICH CAN PRECIPITATE IN PORE SPACE 70
APOR « DFLX«POR»2.82B£ft 71
C 72
C——COMPARE UG OF LIME PRECIPITATED WITH UG OF CAC03 NECESSARY TO 73
c FILL THIS SPACE 74
C-..—ASSUME DENSITY OF CAC03 (CALCITE) IS 2.828 75
C 76
C---.-IF PORE SPACE HAS BEEN EXCEEDED* PPINT DAY, SEGMENT* MASS OF CAC03 77
C— WHICH CAN PRECIPITATE IN PORE SPACE* AND MASS OF CAC03 WHICH HAS 78
C——PP£CIPIT*TED 79
IFUNETllMiIt.GE.APORI PRINT 201, 111,1, APOR * 80
i Hi
284
-------
20* CONUNUF •**
if ( L i .Nf . o> F «NM3< n-Fxc* • ^3
IE ISC* I I ) •(••»«»NH» iIl»e*CLiti«0.0 8*
85
IF(NRYO»s.F(j. 1 i (JO TO 9008 <* tOEL«N< t) .OEL»C( t) •!!) *•»
40
900* CCNTJNUF 41
lFtiFL«vP4.fo.i) r,o TO aooo «*
C ----- CALL T*F; *XO* SU««OUTtNE • «
C*LL FL ( I «FLNOit t > .FLNMJJ I > tftu»CA 1 1 ) «FLCAI I) .FLAN* in •»*
•*oon CONTINUF 96
IFJII.NP.H r,o TO 20 «T
ir«IS€T.LE.IFACT) GO TO 20 9«
•»«
IFCKBYOAS. E0.lt GO TO <*00*
63 AOOITlTt « PLN03UI • l*r
C-- — TEST FOP LOw NM4 COHCENTPATION
6* IF(COM.LT.0.2)6S.66
65 ' AOOITKI) » 0.0 12*
00 TO 67 122
66 AOoiTim « PLNH*U» i23
67 CONTINUE l2*
C —— COMPUTE NET CHANGES FOB NK4. UREA, AND N03 12*
•NETKIJ * DFLNH3(I> • FLXH3U) * EXNH3U) • AOOITllI) If!
ANET2III* OELUPEAd) » FLUPEA(I) - 1*'
«NET3(I)« OELN03(t) • FLN03U* • AOOZT(I) 12*
130
CCNT1NUF 13*
TEST TO DETERMINE IF SEGMENT ONE is BEING CONSIDERED |3?
..
TT SNH3l»ONH3l SSN03l«ON03l tSBEA I =U«EA I $SAl«CAl t?NAl«AN*l * '*
SMQlaAMfil «SC03l=HC03l iSLl'CLl SS03l«C03l *BO*l»Sn%l I3?
I CONTINUE l~Z
138
139
C—— TEST TO DETEPMINE IF ADDITIONAL TI"E STEPS ARE 8EING USED *?
IF(ISET.LE.tFACT) GO TO 16
IFCNflYPAS.EO.l) GO TO 9010
C., --- TEST TO DETERMINE IF MftSS IN SYSTEM *ILL 6E EXCEEDED e
00 5 1*2.0 } 6
IFJANhSd) • ANET1 (I) .LT.0.0) (JO TO 14 }"T
IFJUPEAU) . ANET2H) .LT.0.0) GO TO 1* } „
IFI4N03(I) • ANET3I I I .LT.0.0) SO TO I* }4Q
5 CONTINUE J-fl
6C TO 1* Jj,
c— — USE SM*LLE» TIME STEPS IF NECESSARY J*.
I* TSET « 1 *F « tFACT
9010 CONTINUF
C- ---- UPDATE THE MASSES IN STORAGE
16 00 6 I*NO*«U i5«
ANH3(I) a ANH3(t) » AN£Tl(t)/F SUPEA ( I I a UP£A(H '
285
-------
niJTPT in I
r---.-cALL *m H»I ANCF unuMN* row N
tr iNHVpft.t.f -JOM
ir ( (4MCM.FD. I .»"(>. r | .tU. J*>«MN' I C*| L
r— ---ufTuuN TO su"HouriMt e«cruTf
100
201 FOBM*r(lX»THE SOIL POPOSITY EQUALED ZEPO OU£ To P«ECIPIT»TEO LI"*
10*. 0«Y NO.«.r5»/lX»OEPTM SEGMENT NO. •« ISi / 1 OX-OQPOSITY *LLOWS»3X,
2C10.3»?X«UG OP LIME TO PPECIP I T*TE» « 5x tE 10. 3«2X*UG OF L!M£ M»V€ 2*3
END
SUBROUTINE XCHANGE
XCHANGE
XCHANGE
O_...T»IS IS Ti-F. EXCHANGE SUBROUTINE XCHANGE
•CHANGE
XCHANGF
COMMON/* ION/11
COMMON/TPNIT/ISTPO5) tACTCA(2S».IOPN,ISETN<25)
CCMWON/XXX/PEL»«DELTtMH,STA«TtRO(25 J.TEN125 ).CHECK(25 ).MQISIN XCHANGE
1(25 ) •C«H20?(25 >«MOISOUT«ANOZ(25 liANHZtas ).U»EA(25 tiOPN
2125 I.C2I25 I.ANZC25 l.AMZ(25 ).HCOZ(2« ).CT(25 1.COZ125 )tSOZ(25 XCHANG1
3>«CZ(?5 ).CX(2^ ).SAZ(25 it**Z<25 ).CASZ(25 )tAGS2(2S ).flNHZ<25 I.XCHANGJ
4CYI25 )«CNl(?S )iSA«T(25 )t»N(25 )iBC(25 ) t TE" ( 25 ) tCAZ (25 ) tQ tCPOXCHANOl
lP.XTPACT.SUMN03tTMOP(*) tTO. IOATiU3(^5) .CH.CH1 .JPEPUN. SPC (» ) « I IK (25
tAKCS(25) «AKCM(2SI
COMMON/C02/PC02 (251 . I°C02
XCHANG1
OlMgNSION C»H201(25) XCHANG1
XCMANG1
OATA(TES»1.E-100)
XCHANGI
C— — SET EXCHANGE CONSTANTS
OA " AKCS(J)
0 « »KCM(J>
•SET SEGMENT VOLUMES . XCHANGI
CMH201 ( J)»CMW202(J) . XCHANG2
-COMPUTE MOISTURE CONTENT ON « PEPCENT BASIS XCHANGZ
Bl • CMH201 (J) / IRO(J) »nELX) XCHANG2
81 » Bl-100. XCHANG2
XCHANG2
C— --- COMPUTE SEOMPNT VOLUMES RASED ON INITIAL SOIL ANALYSES XCHANQ2
IFICHECKljI.I-Q.O.O) CVH201 (JI»XTPCT(J»«OELX«PO(JI •«•« 2
XCHANG2
C ----- CONVF9T UNITS F»OM UG/SEGMfNT TO MOUCS/LtT^P . XCHANG2
XCMANG2
<; ----- PESET STOPAGI: LOCATIONS ro" USE IN THIS ROUTINE XCHANGJ
1005 ANH4 a ANhZU) /TMH201 (Jl/UOOO. XCHANG3
A m r.lt J1/CMM201 (Jl /*0080. . XCHANG3
S • ANZ (J)/C"H201 ( J) /2?9<)0. ' XCHANQ3
F • AMZ1J1/CMH201 ( Jl /2*320. XCHANG3
HC03 » HCOZI JI/CMM201 ( Jl/61000. XCHANG3
C03 • COZ
-------
AN03IP « AN03IP • ANET3IP/F tCA(p • CA«P • FLCA ( P/F • F. XCA ( 159
IP 160
ANA(i) x ANA < p • FLANAIP/F » FXANAIP tA"G(p * AMQIP <
1 I /F »fXAMfi I I )
HC03IP • *C033
HP/' » FXCL ( P
CC3(P • C03(P » n CO IIP/* « F*rn3lp ISO* IP • SO* I P •
I P /r . risn* i I i
• QNM4ip . FIMNfetfp « nH_HNH*»P/F |nMN - SLl/''
« C03(l) • S03UFf SO* I I ) > SO* ( p - RQM/F 185
37 CONTINUF 196
IF/F
SOMS(l) = SUMS(l) • ON032/F
SUMS (2) » SUMS(2) • ONM32/F 175
SuHS(3) « SUMS 1 3) » U»EA2/r 17^
177
9011 CONTINUE
!FIO«»N(t). LT.0.0) OPNIP « 0.0
IFICAIP .LT.0.0) CA(P a 0.0 194
IF (ANA ( I ) .LT.0.0) ANAIP a 0.0 - 19$
IF (AMQ( p .LT.0.0! AMQil) a o.o |96
IF(HC03(P. LT.0.0) HCOJ(I) « 0.0 197
' IFICLU). LT.0.0) CL(P « 0.0 196
IFIC03IP .LT.0.0) C03(P » 0.0- 199
IF (SO* ( P .LT.0.0) SO*(P « 0.0 200
201
IF(NRYPAS.EO.P GO TO <»012 203
C ----- KEEP TOACT OF PLANT UPTAKE OF N 203
IFUSET.Le.IFACT) 17tl8 204
17 PL1 » AOOITl P/IFACT %PL2 ' AODI Tl ( P /IFACT 205
TPLANT = TPLSNT « PLNH4( p /IFACT » PLN03 ( P /IFACT 204
Ll*L2»L3=0 207
60 TO 25 208
18 PL1 » APOIT(I) JPL2 » AOOITKP 209
TPLANT = TPLANT • PLNH4IP * PLN03U) • 210
25 PSU" = P5UM • PL1 • PL? 2.11
IF ( ANM3I P .EO.0.0) OIFNH* = OIFNH* » PL2 212
IF4AN03I P .CO.0.0) OIFN03 » OIFN03 » PL 1 213
9012 CONTINUE 21*
IFUSET.LT. IF«CT) GO TO 6 215
tF
-------
c* (.11 /rx
<••
A
F*
\ft»tc*
« • * - r«sn
^ • f •
no
ff « rrrji « CAL •
|F (CHECH (JI .FU. 0.0 1200.^01
c-—-CALL TNF eout_ipi>tuM EXCHANGE SUBROUTINE IF
C--— TIME INTPBVAL
IS TM£ FIRST
Zn* CALL EOfxCM8
JJMOO.
1001*602>«03
IIK(J) a \
AZEIJ) e 3.46737E-5
GO TO *>04
IIK(J) = 2
ZE ar (A»MCn3«»2»EXP(-7.033»U/(l.«Ul I I
602
603
60*
CC3 i
HC03 « HC03«OATtO
29P B » 1.ES/B1
IFtIPC02.£Q,ll GO TO 300
2f » AZF.
.301 If *
2* AlBA
*»0.0
U»SO» T (
96-A-G
.0«(A*F*GI *0.5*(S-M»HC03«ANM*»AN03)
(-9.366»U/(UO«U) J
S-A»
-------
«*«r»p i -u. i*4»u/ i \ . .n
f Mi>*-f».A*AA»<. I
Cf»4A»*«r. -4 ,-»f - i«r 4 V1
I » 1 1 •HM*aq.«, 0 • »A«CC
!» i I » « I i «v. \*> . ,l
IF(AfiS(ZZ) ,LT.TES.OR.»9S(ZZZ) .(.T.TES) 00 TO 515
zz«zz/zzz
tF(*RS(ZZ) .LT.TES.OR.AOS(Z) .I.T.TCS) r,o To 515"
zzz»zz/z
z«z»zz
IF(A8S(ZZZ)-.001)83.83.ai
A»A.=t«Z -
IFA
AA • R» (
.0-r>SlH4)
(CHANlQ
ICMAN|i)
•CHAN)"
»CMAN id
(CHANlO
•CMANlO
XCHANlO
XCHANll
XCHANll
XCHANll
XCHANll
XCHANll
XCHANll
XCHANll
XCHAN12
XCHAN12
XCHAN12
XCHAN12
XCHAN12
XCHAN12
XCHAN12
XCHAN12
XCHAN12
XCHAN12
XCHAN13
XCHAN13
XCHAN13
XCHAN13
XCHAN13
XCHAN13
XCHAN13
XCHAN13
XCHAN13
XCHANl*
XCHANl*
XCHANl*
XCHANl*
XCHANl*
XCHANl*
XCHANl*
XCHANl*
XCHANl*
XCHANl*
XCHANl 5
XCHAN15
XCHANlS
ICHAN15
XCHANlS
XCHANlS
XCHANlS
XCHANlS
XCHANlS
XCHANlS
XCHANl*
XCHAN16
XCHANl*
XCHAN16
XCHAN16
AMH4
XCHANl*
KCHANl*
XCHANl*
XCHANl*
XCMANl?
289
-------
cc *
'•SCUT (pq»qfl-4. J»44»CC 1
YT , .qo.O | / <2.J*A4 I
r}\H4 s PM*» - T
S»T « S4T . *
ANN* *
S . S
tr ini 790 i r^o . TI\
ro\ \r ir i r<*n « »HM-« . i)»»*»rC
U (XI *» I 7SO
650 Z1«Z
63 ZZ=-< ( <
Z2Z»( (
ir«*«S(Z2) .LT.TES.OB.»flS(ZZZ).LT.TES) GO TO 600
zz«zz/zzz
IF(»«S(ZZ) .LT.TES.OW.ARS(Z) .LT.TES) (30 TO 600
ZZZ«ZZ/Z
44
ir«AfiS(ZZZ)-.001)6*.6*t63
Z«-Z1
r,0 TO "S.I
651 I^U) 7«i2
7S3 CAL«C*L-7
600 IFII<.FO,2I GO TO 606
60S
A06
»fl
49
50
51
1K«3
AZE(J) » (fll ••1.6«)«Zx
OEL»*-»1
IF (OFL-CHl) »q.49.24
IF
IMDEL-OD51 «5l«2*
IF (DEC «CHl ) ^* .52.52
IF«DEL-CM1I8.8.24
OEL«»-A*
66
1000
A7
CONTINUE
|F(CH€CKij) .F'J.O.OI CMH201IJI « CMM2021J)
4NH4 « 4NM4»CMH201 (J) »l*000.
ACTCAIJI * A»io.»»(-2.03ft«u/r i .»ui )
A • (4 • C*SO) «CMM201 ( J) «*0080.
S • S«CMW201 (J>
XCMAN17
XCHANl 7
XCMAN17
XCHAN17
XCHAN17
XCMAN1?
XCHAN17
7
ICHANl*
ICMANIM
XCHANlfl
XCHAN18
XCMANlri
XCHANlA
XCMAN1H
X CM AN 19
XCMAN1Q
XCHAN19
XCHAN19
XCHAN19
ICHAN19
XCHAN19
XCMAN20
XCWAN20
XCHAN20
XCHAN20
XCHAN20
XCHAN20
XCHAN20
XCHAN20
XCHAN20
XCHAN20
XCHAN21
XCHAN21
XCHAN21
XCHAN21
XCHAN21
XCHAN21
XCMAN21
XCHAN21
XCHAN21
XCHAN2?
XCHAN22
XCHAN22
XCHAN22
XCHAN22
XCHAN22
XCHAN22
XCHAN22
XCHAN22
XCHAN22
XCHAN23
XCHAN23
XCHAN23
XCHAN23
XCHAN23
XCHAN23
XCHAN23
XCHAN23
XCHAN23
290
-------
i sii ••> I pun .
>»
rc\ « r,M"' wni>n i i j i ••>!)!> i' 11 .
,-. , ,«f o ii' . f . j . .) . o i « t> ii . • ill
• Aft •C.H/I.II • »N»'« »< / ' Jl • »
*N/ , Jl . «. »*M/ IJ.I - »
MCI!/ ui • i-ro i *t')/ ( ji • c« »
r* iji • M \^ni ui »it
j) , DNH* XCHAN25
CMrc««.M«c*Kri«»i. XCMAN«
401 CONTINUE XCMAN2S
C COMPUTE DELTA VALUES ro* COMPONENTS . XCHAN25
- * 8N** - 4N«Z(J) *EXCA = A - CZ(J) «rWAN2S
e...-,. > S - »NZ SEXAMG = F - AMZ(J) XCMAN25
EXHC03 . HC.13 - MCOZCJl SE*C03 • C03 - COZlJ. JSlSl
«»CU « - - Cv.JI ,«*SO* - G - SOZtJ. ,CMAN25
.. XCHAN26
EZ«JI • FT icxut • CT
SAZIJ»SAT *XXZ1J)»XXT XCMAN26
~""~ 1AG^Z(J)**Q50 -«*^A4.odL
5EYU1.EC -CHAN26
ISTPCJI • U»»?
c*scji«c*szij) »CMMaoi i ji • 13ft 180.
XCHAN26
XCHAN26
TO SUBROUTINE COMHINt .j-uAUJh
. Jlt™^ ™
1001 STOP XCH4N26
END
SUBROUTINE EQEXCH
SUBROUTINE FO£XCH(CA.AMG,SOS.CL»SO.HC03tC03tEC.ANM4,E5.C5.SA5
j. AN03.01tO)
C THIS SUPPOUTINf COMPUTES THE AMOUNTS OF IONS CONTAINED ON THE F.*-
C CHANGE COMPLEX (BASED ON INITIAL SOIL ANALYSIS! eoeXCM
OA m 1.414/01 EOEXCHl
" ?2 EOEXCHl
EOEXCHl
AGSO»O.O FQExcHi
42 ACT2»EXP<-9.3
-------
TJ?/ I4.<
CX-C»-C»S*
IF I»PS (uu/u-i . i -i .OF-»
*1 U«UU
SO«SriT
r'0 T0 »2
4ft CASO«C»SX
»oso«
c*«r«
0.5»(SOS»MCr>3»CL»ANH**AN03)
*CTl«S«JhT(»rri»
*CT««iOBT<»CTl I
C»"C*»2.
«AMfi
^C/ '
"ACT
CS-EC-E5-S45
i »cri
IO«*CTl»AMf5/(ACTl»CAI I )
looo
EQEXCH3
EQEXCH3
EOEXCH3
EQEXCH3
EQEXCM4
EOEXCH4
F.QEXCM4
EOEXO4
COEXCM4
COCXCH4
COCXCM5
COEJICM5
EOEXCHS
EOEXCHS
COEXCM5
SUBROUTINE EXECUTE
SUBOOUTIME EXECUTE
C-— -- SUBOOUTINE TO EXECUTE P»OG«»AM FOP EACH
TIME INTERVAL
c- —
EXECUTE
EXECUTE
EXECUTE
EXECUTE
>IOrSTR,lOYSTP,lLO«tm»lNFILl«ICONTi,jPAS JtJS A
:tloj.SMO*TM.«M,O«IPRINT,jpsiNT.iNKtiPUNCH*isToPt EXECUTE
1ITEST • IPE*00. IMASS t lAOO< 25 > » lORNAp (25) tHOPl1*) »TOTN|99» • yf.»B . $11$
2AIRO(9).IP»(25)tTT<60».FF.RT<7)iOFEPT(3l»NOPGIN.NFEPTINtNTEMPIN. EXECUTE
3ITOT.JTOT.IPTOT . EXECUTE
CO***'ON/XX2/A 1 < AS.A 3.x EXECUTl
COMKON/ryr/STAPTiIOTE.MONTH.I.LL CXECUTl
COM*ON/XXV/ I C^ECK • ICOUNT. CON VtPKtPKl tCRQP
COMMPN/AFG/ENH.3. I I .LLL
.CELT.MS.wTAST,80(25 ).TEN(25 >«CHECK<25 I.
)«MOISOUT(25 >tAN03(2S )«ANH3(2S )iU«EA<25
2C25 )«C*r25 )»ANA(25 >.AHG(25 L»HC03(2S )tCL(2«5 t.C03(25 )
3l.ES(2S ).C5(25 )fS*5(25 ).x«5(25 >.CASO(2S
4£C«25 ).rNl(?5 ).SAMT<25 I.HNI2? )t»C(25 ).T£M(25 I«CAL<25
IO«XTPACT«SIJMN03 GO TO 8002
JPAS«1
REMIND i
SSSSAC2
SSSS 02
SSSS E2
SSSS F2
SSSS 02
1SSS H2
SSSS J2
SSSS K2
SSSS L2
SSSS *<2
SSSS N2
SSSS 02
SSSS P2
292
-------
.n i *tO *o •toij'j
IF
OC
CALL
8001 CONTINUE
GC TO *009
800? RE AD i nil
IF (FOF ( 1 I )«007.8003
T2
SStS U2
0003 Bojsr HOT*.
H004 FO»"«T(/, 5«. • FUUOP- E\.J> OF MLE MOT FOUND ON TAPE I AT START
1 *EAD NO. •• I?'. S*. » F1FCUTJON TERMJNATFO •)
CALL EXIT
C- — -LL » STAPTlNf? HAY. MM s TEOMINATION
or * I»ILO.IHI
If (».RYP«< .PQ.l) tE5TM.EO.O.)600.60t _
600 CCC » CONV
6o« rs « i SIOEPTH « i
00 TO 602
ft01 IOEPTH > OEPTh/OELX • 1
IF(IQEPTM,LT.2) IDEPTh » 2
IS « Z
CCC » 0£LX/0£PTH»CONV
602 SAVF1 a JANH3»CCC»O.T7T7
SAVE2 » AAN03»CCC»0.225d
SAVE3 a »UBEA»CCC«0,*466
SAVE*»AC»»CCC
SAVEIO*ASO*»CCC
DO 302 j « IS.IOEPTH
-AOO THE FE&TILIZEB TO THF. PPOP£R APflAVS
A».M3(J| m ANH3IJ) • 5AVE1
J» * »N03(J> » SAVE2
J) * UB£A(J) • SAVE3
CAtj)»CA (j)-SAVE*
C03(Jl»C03IJ)»SAVE9
SO*
-------
302 CU»'«1<.-CU»4Sn*»^AVEl 0
« 'f 0 OW.ANlC-N APPLICATION HATE
IF I I ,r«J. TOWNAPf* M T.I*
T TPNTJNUF
PFAO iioi IIF«TM. ATNI
• Ntim. I* » i
r ....
*M>
OC 303 j«>ttr£PTM
STQOE ORGAMfON APPLICATION INTO PBQPEP
I • SSA*T«CCC
C — — -
C-— -STQPE ACCU"-1 AMOUNT OF ORGAMC-N AOOEO
C SAVE « SSAMT»CCC»O.*/ACNI
SAVE • SSAMT'CCC
CUMSUM « CU«SUM * SAVE
303 CHlij) * ACNJ
GO TO 17
•t CONTINUE
t
IT
COMPUTF Tf.wpf OATUPE ^EAO-IN DATE
IF t»oo < i .n .FO.O.OP. i .eo. ILOJ
C-—
590
CALL T£MDEOATUPE INPUT SUBOOUTINE
CALL TfMP tNTE^PlN » NTFMPIN • 1
CONTINUE
900A CONTINUE
.EQ.O) K * 2
LOOP TO EXECUTE PROGRAM FQP EACH TIME INTERVAL
LLL IS THE NO. OF TIME INTERVALS PER OAT
C-..--
790
;-_..-EKTF.P aPUTINF. TO AOO IRRIGATION WATER COMPONENTS
EXECUTE
H$tt AA
FXECUTH
CIECUTB
CKfCUTH
tifCUTO
eircurq
etecuTo
e*ecuT«»
CXCCUT9
CXECUT9
CXECIJIO
EXECUIO
CXECUIO
CXECUIO
EXECUIO
EXECUIO
EXECUIO
F.XECUIO
EXECUIO
-THE P°OG°AM MAT OR MAY NOT CALL ALL OF THE COMPUTATIONAL SUB-
•ROUTINES FOR EACH INTERVAL
-ALL CRITICAL ROUTINES AWE CALLED AT LEAST ONCE PE» DAY
00 10 11*1.LLL
READ INPUT OATA ON TAPE1 FROM MOISTURE FLO* PROGRAM
IF(ITEST.EO.O) READ!1)(II.12.13.I3.CMH201.MOISIN(j).MOISOUTlJ)
1 TEN! J>tU(Jl•J»l.01
IF(II.E0.1.*NO.CMM20l(l).GT.O.O)790.79S
•CHECK TO SEE If ™IS IS AN IRRIGATION HAY
00 7<>2 L*. » l«l
EXECU11
CXECUII
CXECUH
EXECU11
EXECU11
SS1SAU
EXECU11
EXECU11
exEcuii
EXECUI2
EXECU12
EXECU12
EXECU12
EXECU12
CXECU12
EXECUI2
EXECU12
EXECU12
CXECU12
EXECU13
EXECU13
.« SSSA13
793
SAVE1"AIPB<1)*CMH201(1)
SAVE**AIRH<3)*CMH201 (1 I
n
fSAVe2*AIRR(2> »CMH201 <1)
1SAVF.53*tRR(4)«CMH201 (1)
1SAVE7»AIRB(6)»C"M?Ql(I)
1 (8) •OMSOl < 1 '
AKH3U»ANH3>
CA(l)»CA(ll
CLU)*CL 111
SAN03(il>AN03(l)*SAVE2
tANA(l)*ANA I 1) *SAVF.S
SHC03M )«MC03(1)»SAVP7
«C03(l»C03(ll*SAve9
*SAVFIO
j—...STOP-IT ACCU" AMOUNTS OF COMPONENTS
CXECU13
EXECU13
EXECU1*
EXECU1*
EXECUl*
EXECUU
CXECUU
EXECUU
EXECUU
EXECUU
EXECUU
CXECUU
EXECU15
exccuis
EXECU1S
EXECUIS
EXECU15
EXECUIS
294
-------
SAVE,'1 »ClJMr»«CUMf A«S« WE 4
»riiM4
7 tC;i"CL
01 .••co.
POINT ?0*.C**21l (1 I .1
CONTINUE
79S CCNT IMIf
I' (M00< ir>AY. t"0[NT I .EO.O.AUO. II .EO.JPRI NT I 400.401
400 PRINT ?Q*>«VEAP. 1111
PP[»iT ?0*
*0I . CCNTINUF
f-~ —CALl COMPfNF Ml8»OUTINF
CALL COMBINF(IDAY,IPB[Nf•JPHINT)
10 CONTINUF
4 CCNTtNUf
WfiTIJRN
EXECUl*
IX'PPEOICTEO AMOUNTSCUG/SEGMENT OF SOIL*—FGVCL»'SX»ES«>M»»RC02 ( ATM) »i
20*i FORMATI//IX»VEAR« ••I».IO*»OAY« ». 14, iox«r IMF. INTERVAL* •
1.14)
20T FCHMAT»///10X»AN IRRIGATION OF*.F6.1.«C« -AS APPLIED ON OAT
SUBROUTINE OUTPT
OUTPT\<\
c—-—THIS SUBROUTINE WRITES PBEIMCTEO TOTAL AND DELTA AMOUNTS FOR THE
C —--COMPONENTS »NO VOLUMES ON TAPE2 (UNITS APE EXPRESSED IN UO/UNIT OUTPT
C—--APEA ANP ML/UNIT AREA). OUTPT
OUTPT
OUTPT
DIMENSION A-T(IO),AMTIIIO).OEL(IO) - OUTPT
OUTPT 1
INTEGER Q.O.STAPT.CPQP.TO * OUTPT i
INTEGER YEAR
"I[*L MOISOUT — OUTPT
OUTPT
COMMON/SA8LE/SU»*S (3) OUTPT
K25 1.CMH20K2? ) .-OISOUT (25 ).AN03(25 1.ANH3I25 >.UREA<25 (.C^N OUTPT
2(25 I«CA(25 >.ANA<25 l.AMG(.SO*;?5 OUTPT
3I.E5I25 1.C5I25 ).SA5(25 ).«X5(25 )«CASO(25 ).AGSO(2S )»BNH*(25 I OUTPT
»EC(2S ).CNl(25 ).SAMT(25 ).«N(25 ).»C«25 ).TEM(25 I.CALC2S ).0 CR')OUTPT
1P.XTRACT.SUMN03. THOP-(4) ,TO.IOAY,Ut25) .CH.CHl. IRERUN OUTPT
COMMON/ARtE/TITLE(10) .SMQNTM,MM.O.I PR INT.JPRINT,IHX,I PUNCH.I STOP' 2222
UTEST.IRFAOP. I^ASS.uoo(2S). IORNAPI25I »HORc») .TOTN<9<» » YEAR » 2222
2AIPRI9),IRPI251.TT(60).FFRT(7).OFERT(3).NOROIN.NFERTIN.NTEMPIN. 2222
3ITOT.JTOT.IRTOT.NT 222
COMMON/JP/CASI25).AMG5I25) 2
OUTPT 2
C —--ESTAPLISH STATEMENT FunCTION OUTPT t
SOHA(X.Y) > X»Y OUTPT ^
IF(K.EO.l) 1.2 OU 2
c--—ZERO INITIAL VALUES OUTPT *
1 SUMOUT « SUMQUTl s 0.0 OUTP" ^
00 3 I" 1.10 OUTPT '
3 AMTd) « AMTld) • 0.0 OOTPT
00 TO 5 ' OUTPT :
2 Y • CMM20KO) OUT!!T 1
Z • -OISOUT(O) ~'
295
-------
V
ir
7'Y
O-
-SLM TMf COM(>ONPNTS
AMT ( I » • SUMS I I I
1) « SUM5<3>
A • SURA(CAS(01iY)
H«SUBAiAMOS(Qi•Y)
(4)
1 (9)
AMT(10)
r*i • sup*ICA toi.vi
(51 • SUM* IAMA « AMTCII - AMTl | I I
.—COMPUTE DELTA VALUE FOP VOLUME ou'
OFV.N » 5UMQUT - SUMOUTl
..-•RITE SUMMATIONS ANO OELTA VALUES IN
• RITE (21 Yf AB, IOAY .SUMOUT.OELN. (A1- - ! I) ,OEL 11) i I«l«10)
C
C
c
T
OESET VALUES FO» DELTA
OC 7 I«i.lO
tl » AMT(T)
— — -HETUON TO MAIN
RETURN
100
END
OuTPT 3
OUrPT
OUTPT 3
OUTPT 3
OUTPT 3
OUTPT 3
OUTPT 4
OUTPT *
OUTPT 4
OUTPT 4
OUTPT 4
OUTPT 4
OUTPT 4
OUTPT 5
OUTPT
OUTPT
OUTPT
OUTPT
OUTPT
OUTPT
OUTPT
OUTPT
OUTPT S
OUTPT 6
2222 6
OUTPT 6
OUTPT 6
OUTPT 6
OUTPT b
OUTPT 6
OUTPT 6
OUTPT 6
OUTPT 6
OUTPT T
OUTPT 7
OUTPT 7
OUTPT 7
5
5
5
5
5
5
5
5
INTEGER FUNCTION DAY
12
1
2
3
A
•>
A
7
10
U
13
FUNCTION OAYETJRN
RETURN
RETURN
RETURN
RETURN
DAY.
0»Y
OAV
OAV
DAY
OAV
DAY
DAY
OAV
OAV
DAY
DAY
OAV
OAV
DAY
OAV
DAY
296
-------
SUBROUTINE IDAY
I.'1** '
t MONTH, TOTE> J0»" •«»
,OAY
JJD«V « 0*' ( IDTE ."ONT" I
JOAY ' J.;CAT - JO»Y • K
If I „•£•-». L£. Oi 1.2
JCAV JPAY « 3'>5 • «
IDAV
IOAV
13**
IDA*
IOAV
SUBROUTINE THEDATE
12
I
2
-Me
tKtL»SMONTH«it i j
M .« ..«• nv . t - •fAN03(25 I«ANM3(25 ).UR£A(25 ).ORN UN'leil
2(25 )»CA(25 I.ANAI25 )t*Mf,(25 )tHC03(25 ltCL(2«> >.C03<25 )tSO"(25 UNl,!ij
'2*- ).C5(2S 1.SA5I25 ).xx5(25 ).CASO(2S """ "~ ' ""
4ECI2S I«CN1
• AMG(J)«rKH201(j)»12.16
UNITS}'
UNITS}
UNITS
• UNITS}
UNITS
297
-------
CC3(J>
CL(J)
SC*U>
-rn*vr»T mi>«4 ^
MlM 1 I II * AMU 1
• *'1 1111 • «'in l
UbP * I it • uw» «
f.Hjt • C* i (i '
i J> •C*n«>01 i J> «M .
J' «C"M,?oi ( J> «10.0
("L IJI »CMH?0 I Ijl »35.»6
- SO* (Jl •C"«?0l I J) •*H.3
»of.u<
r N T TO
| | / | r
j I / i r MM/<) |l (I • I «. . 0 I
j i / I (.'MM /I) | I 1 1 »/n .u I
Jl •('U .0* )
• ••GUI • »MI', (Jl / !CMM?01 ( J) • I? . 1M
HC03U) » HC«U< j) / (CMH?oi ( Jl »61 .0 I
C031JI • Cnil j)/ / .CASO( 25) .AGS0125) .aNM4(25) •
25).SA»T(25)»PN(25).»C{25l.TEM(?5».CAL(25) .Q.CBO
5P.XTOACT.SUMM03.THnP (» ) .TO. I DA Y.US (25) . CH.CH1 , [«CPUN. I SWC
6S(JHOUT.P$PA (60)
»ONH32.0N03l .ON03^.C»U ANAl,
DIMENSION ANH3 (25) . AN03(25) tU&£A (25) .CA (25) .AN* (25) iAM6(25> .1-003(2
15> »CL«25) »C03 (25) .SO* (25)
INTEC-F3 0
IFIJ.NE.?' GO TO 1
00 1" 1*1.0
ANH3(I> * BNM3(I> SA*03 SSO*(I) * BO*(I)
CONTINUE
OPMOIS(0*1> * OPuOIS(Q)
ANM3(0*1> : 4NN3(0) SAN03(Q*l) * AN03(0)
URE*<0»1> * URfAIQ) fCA(0«l) > CA(Q)
ANA(0»1) * ANAlO) $A«G(Q»1) * AMG(Q)
HC03»0»1) = HC0310) tCUO*l» « CL(O)
C03(Q*U * C03(0) SSO*(0»D » S0*(0)
CONTINUE
IF(00*OtS(1 I.LT.0.0) OPMOISIU » 0.0
IF(MOIS!N(J).LT.0.0) 2.3
CCEFIN * MOISINfJ)/ORMOISIJ)
GO TO *
IF(0»MOIS(J-l).GT.0.0) GO TO 1*
CCE^IN « 0.0
GC TO 15
COEFIN r «OISJN(J)/OPMOISIJ-I)
CONTINUE
IF(fOlSOUT(J),LT.O.OI5.6
COEFOUT m MOISOUT(j)/OPMOI')«6L
-------
• « • J
'if! TO 10
« * • J-l
10 IF iCnEFr>nT.LT.O.O> I I » 1?
II (.».)•!
N «Q TO 105
UOEAl > COEFIN»U»E* (K)
UO TO 106
105 U«C«l • UPF»3
106 'JSE»2 » COF.FOUT»UPE» CD
•NA1 » COEF[M»4NA(K ) SANA? a COEFOUT • AN A
HC031 • COEFIN*WC03("') IHC032
CLI » COFFI»I»CL ix> set? .
CC31 » COEFlN«C03c<) NC033 » COEFOUT»C03 CL )
50*1 • COEFI^«SO*IK) ISO*? « COEFOUT»SO*«LI
FLNOI « ON031 - ON032
FLNH3 « OMHJl - ONM32 .
« UBEAl - UPFA2
« c»i - CA2
ANAl - ANA2
AMQI - AMG2
FUMC03 * HC031 - HC032
FLCL * CLI - CL2
FLC03 « C031 - C032
FLSO* " S0*l - S0»2
LSET1 « LSET2 * 0
SUBROUTINE PRNT
PRIMTI IPRINTI t IPPlNTjl P9NT
PPNT 3
C-— -THIS SUPROUTI»4E PRINTS CONTROL AND INPUT DATA PRNT
PPNT *
COM»«ON/APLE/TITL£(10) .SMQNTH.HM ,0. JPRIMT . JPRINT t INK t IPUNCH, I STOP t PRwT -,
IOPNAPI25) .MOR<<»
t TT«60> »FFRT (7) tOFFP.T(3) .NORQIN.NFERT INtNTEHPIN. °PNT
3ITOT.jTOT,IBTOT.Nr PRNT
.a2. A3.X P«NT .
COMUON/YYY/STAHT . [OTEi"ONTH. I .I.L
CCMuON/XXY/ICMfCK. ICOUNT .CONV.PK.Pirl .CROP PRNT
CCMMON/X»X/nF.LX.OELT«MS.kiT4PT, 80(25 ) .TEN(25 I .CnrCK(25 (.MQISIN
1125 I.CHH?Ol(2S ).MQt50UT(?5 I.AN03C2S ).A».TEni2!> I.CALI25 ).O.SPflPRNT g
1P.XTPACT .SUVN03. THOP. <*i .TO. I DAY ,u (25) .CH.CHl » I RERUN. ISwCM.CUMSUM.PRNT g
ISUMQUT.RF.OUCE PPNT in
_ PRNT
INTCGCP TITLF.SMONTM, START, o. TO. YEAR ......
299
-------
PONT
22
23
-----PO (MI
i r-o ro l
«•«•
no mt l-<> rruifH«»i
in?. |uu i«*r . i
, CR<1P.(.L,P« t"M t*»K I .OH.*
nurr
. t rr M . INK . i"»ss .
,IP«INTI,
e— —
S*IP PARF
POINT 103
POINT 104. (THO«( J) « J»l «TO)
109
00 10 J»l»'iT
BC*0 ««) (TT(I) t!»l.TOI
10 PPJNT ins. j, (TT ( I) . I«l »TQ»
8
P»6E
103
c— — PC:*T KATES ANALYSIS
107
C--— -PPJMT jooifiJTiON WAT^P ANALYSIS
PRINT 10*. (AIPBU) .I»ltO)
C. ---- PPINT IRPI6ATION 4pPLlCATION DATES
PP.INT 110.
-PRINT FERTILIZER APPLICATION DATES
PRINT 111* (IADOU) tl»lt
PRINT 112
j»i.7>
OC 2 I»UITOT
BEAD ««»)
( 3 I
FORE* * FFRT<») »CONV«.»<>66 $PCA m fEPT(5)»CONtf
FSO* « FFBTI6) »CONV SFC03 » FERT(7)«CONV
C ----- PPINT FPBTILIZEB APPLICATIONS
2 POINT H3.I .FE(»T< 1) «FNH*.FN03.FUREAtFCAiFSO*«FC03
PEKING Q
REWIND 10
PRINT 100
C- — -PRINT ORGANIC APPLICATION DATES
PRINT 1I». ( IORNAPI Jl .J»l .jTQTl
POINT 115
00 3 I«1»JTOT
READ <10l IOFERTU>»J«1«3»
FCfiN * OFE"T(3l«CONV
C-— -POINT ORGANIC APPLICATIONS
3 PRINT U3t I.OFERTil ) fOFf«T(2» .TORN
REMIND 10
C— .-PRINT COMPONENT MQWIZON DEPTHS
PRINT 106< IHOP2
33
34
35
36
37
41
42
43
44
45
49
50
51
52
53
54
64
65
66
73
74
75
76
77
78
79
80
81
A2
83
84
AS
PRNT 47
PRNT 88
PRNT 89
PRNT
PRNT
90
91
PRNT 92
PRNT 93
300
-------
PBfNT 103 PPNT 94
PBNT 9*
100 rrjwuAT I 1 «1//. )A» « 1 QAH// I PBNT 97
101 FOWw/t T 1 t(,j( •CONTOOL CAMO ^U^mW Y •/S 7 I • ( B AS I C PAflAMETFPSI•//35» PBNT 94
{•ST*aT|»i(; unsirn *•« |5*10**(TRACT ••.F5.lt/35x PUNT 99
|»ST»PTjMr. ftAV ••• (5. 10X»C*>OP •••I*/JSX PBNT 100
«*«n|l *rr,ttfnt Vf/f •••."S.I." (•*"• If , »CONVf »0l •• i* ^»rf « / J*1 • PPNT I"'
4»t|wr |MTF"V nif/Ni •• • (••• • I o i »r»«tf.* / ••«»•»,l/Jlx PMNT |0*
•• • !^ t t 0 « "T» AH •••f^/l'iX PBNT 10T
•«.»H.O////| PBNT IOH
l»IPPINT ••.!•}. 10x«IPEAOP »».I5/35X PBNT 110
•».I5.10X«ITEST »».I5/35X PRNT 111
-•.15/35* PBNT 112
. ... . «»,I5////) PBNT 11*
103 FCR"AT(lMD PRNT H'
104 FCRM4TI//1SX»*E?KLY TEupfRATuRE OATA»13X»HO^IZON OEPTHJCH)* PBNT H6
1/46X.6(3X.F6.1)) ' PHNT 11^
105 FCP»'4T(2CiX.I3.2X»TEMPERATUBEiDEG C)»»2x.6F9.1l PRNT 110
116 FORMAT(//IOt»COupONENT HORIZON OEPTHSICMI* . 6x.6l3XF6.il) PflNT 11'
107 FORM AT(1Ox»I°BIGATION «ATE° ANALYSIS(P°"l•/10X»NH4»7x«N03*7X»C4 «7PRNT 12"
IX'NA «7)i«Mf, •*X»HCQ3'»7X»CL *7X»C03»7x«SO**) PPNT 121
108 FCBMAT(3X»9F10.2//I PBNT 122
_ 1 'i 1
104 FORMATi//> PBNT iej
no FORMAT 9999 12^
115 FORMAT(//10X»ORGANIC-N APPL ICATIONS(UG)•/10X»OEPTH»5x•C/N»5X»OBN*)PRNT 13"
PBNT 131
SUBROUTINE CHK
SUBROUTINE CHK(Ll>L2>L3tj»EXNH3.EXCA.£XANA.EXAMGtOELN03«OELNH3.0ELCHK ^
10BGN.OELUBE4) CHK ^
CHK *
C THIS SUBROUTINE DETERMINES IF THE NITROGEN TRANSFORMATION ANO/OR CHK *
C. ION r«TEN(25 )»CHECK(25 )»*"OISIN CHK
1(25 )*CMM201(25 ).WOISOUT(25 )«AN03(2S ).AMH3(25 I.UREAI25 ).OBN CHK ||
2(25 ».CA(25 ).ANA(25 I.AMGI25 ).HC03(25 )«CL<25 ).C03<25 ).SO*(25 CHK J*
3UE5(25 ).C5<25 ).SA5(25 ).XX5(25 ).CASO<2'i )«AGSO«25 ).BNH4(25 ).CHK {^
4EC«2^ ).CNl(25 J.SA-TI25 I.RNI25 I.RCI25 ).T£M(25 ).CAL(25 I.O.SROCHK j
1P»XTRACT.SUMNO.T.TMOH(4) ,TO» lOAYiU (25) .CH.CH1. IRERUN CHK J
COMMON/XX2/A1.A2.A3.X . CHK '?
CHK J!
BEAU MOISIN. MOISOUT CHK }"
CMK Jo
DIMENSION XI7.25) . CHK «"
CHK |»
LI • L2 « L3 « 0 CHK I*
lF(ABS(fXNM3).LT.Al.ANO.*8S(EXCA).LT.AII1.2 CHK e*
1 IFIAPSIEXANA) .LT.Al.AND.A3SIEKAMG) .LT.41)3.2 CHK e_,
3 IFIARS(X(2.J) - 4Nn3(J)).LT.il)*.2 CHK c.
4 IF(*flS(X(5.J) - CA(j) I .LT.AD5.2 CHK ^
5 IF(ARSIXI6.J) - ANA(J)).LT.A1)6.2 CHK f~
t> IF UflSur <7, j) - -""•••• • T «•»••-•» *•** *•
301
-------
11
12
13
»4
16
ll » 1
\c (sP»P»«;.ro. 1 1 *0 TO 0007
IF IAPSIX < i • ;> • *M>3 1 „ M ,LT.A2! 1 1 •<»
IF I4MS < X (Z> Jl - ANH3 ( jl I ,LT.A2» l£.
IF I4MM uiltMNI Jl I ,t»r.A J,Oi».AM«»«««OHOUTl Jl I .(§T . A J I 1 S t 1 6
LJ • I
Ml.jl • **03U> *«(<*• Jl • ANMJ(J)
X(3t.j) • O^F*U> IX(t.J) » 0«N(J>
115. J) » C»IJ) »X<6fj) « 4NA1JI
CHK
CH«C
?•*
CHIC
31
32
33
3*
».»• 436
CH< 37
CMK Id
CHK 40
CHK 41
CMK *2
CHK 43
CHK 44
SUBROUTINE SKIP
SL8ROUTINE
c—
c—
c-
c-
c-
19
74-28
.— PBOGPAM TO SKIP F»OM PRESENT LOGICAL FILE TO NEXT LOGICAL
— IUMT*LOGICAL UNIT
IF
120.10
20
ESO
SKIP
SKIP
SKIP
SKIP 40
SKIP 50
SKIP 60
10
20
30
SKIP 80
SKIP <»0
SKIP 100
SUBROUTINE BACK
SUBPOUTINE
c—
c
c—
c
c
c
10
74-2**
TO ^ACK f»0" PRESENT LOGIC»L FILE TO END OF PREVIOUS
- LOGICAL cILf HE.JUST BEFORE ENO-OF-FILE MARK)
- IUMT»LOGICAL UNIT
IF (EOF ( IIINIT) I 30i20
20 PACKSPftCE [UNIT
r,C TO 10
30 «4CKSP*CF IUNIT
END
BACK 10
BACK 20
BACK 30
BACK 40
SACK 41
BACK SO
BACK 60
BACK 40
BACK 90
SACK 100
BACK 110
BACK 120
BACK 130
BACK 140
302
-------
APPENDIX IV
WTQUAL1 MODEL
Program WTQUAL1 is a numerical model utilizing finite difference
techniques to predict transient, two-dimensional areal groundwater level
or piezometric head fluctuations, the coffesponding flows, and convective
transport of conservative ions. The program is capable of modeling either
confined or unconfined aquifers (but not both simultaneously) with or
without leaky conditions present and streams or lakes that are hydraul-
Ically connected to the aquifer. The hydrologic and geologic parameters
that define a particular study area are incorporated into the model and
each parameter may vary in both space and time. Variables considered
include water applied as irrigation; recharge due to precipitation, lakes,
ponds and recharge areas; withdrawals from the aquifer by pumps and
phreatophytes; and geologic parameters of permeability, storage coefficient
and bedrock, ground surface and initial water table elevations and source
of contaminants either as a slug source or continuous source. The source
1s considered to have initial concentration of 1.0 and values computed
are relative concentration (relative to 1.0). The purpose of this Appendix
1s to describe program WTQUAL1, the procedures required for its use and
the required data.
, APPENDIX A-I. PROCEDURE FOR ANALYSIS
Program WTQUAL1 consists of a main controlling program and several
subprograms. The main program's primary function is to control the
execution of subroutines for all time steps at which calculations are
desired. The basic sequence of events is shown in Figure 1. A detailed
flow chart and program listing is contained in the Appendices. The sub-
programs are designed for specific tasks, such as physical parameter
Input, solving a set of simultaneous equations, and mass balance
303
-------
(START)
\ INPUT/
\DATA /
IS NUMBER OF \
TIME STEPS )—YES
COMPLETE /
NO
CALCULATE
HYDROLOGIC
PARAMETERS
CALCULATE
HEADS
CALCULATE
FLOWS
CALCULATE
•%
'o
\ OUTPUT /
\RESULTS /
Figure 1. Basic sequence of computations.
304
-------
computations. A description of each subprogram is contained in A-VI.
Problem control parameters and certain variables are organized into
labeled common blocks by type of usage. All matrices are packed in blank
common so that variable dimensions can be utilized. Section A-IV contains
an alphabetical listing of the primary variables used by WTQUAL1.
The program is written in Fortran IV utilizing several CDC system
subprograms. These particular subprograms could be easily adapted to
other computer systems if desired. WTQUAL1 does not require changes
within the program for most problems analyzed. The exception to this is
1f the ET subprogram (evapotranspiration) is used. Function ET permits
calculation of phreatophyte use based upon the depth to water table.
The algorithm used would vary from one study area to another, therefore,
necessitating a change in Function ET.
The area to be studied is overlain with a grid system such as that
shown in Figure 2. The selection of the space dimensions (DX.DY) for the
grids is dependent upon the availability of geologic and hydrologic data,
and the desired accuracy and detail of analysis. The accuracy of the
solution is enhanced as values for DX and DY are decreased, providing
the availability of geologic and hydrologic data justifies the additional
computation time. The space dimensions are also dependent on the storage
capacity of the computer being used.
The grid system selected should be oriented to allow for easy boundary
approximation, provide for easy adaptation of hydrologic and geologic
data, and provide the necessary model accuracy. Space dimensions should
be small enough that the geologic and hydrologic conditions may be reason-
ably assumed uniform over the entire grid. In areas where detailed values
of water level or piezometric head is desired, smaller values of DX
305
-------
J-direction
direction
Aquifer Boundary
Grid (3,2) is identified as grid (NR+3) where numbering
1s done by columns.
Figure 2. Sketch of grid system.
306
-------
DY should be used. Grids located outside the boundary of the study area,
such as grid (1,NC) must be defined as impermeable grids. Hydraulically
connected lakes and rivers must be specified as constant head grids.
Boundary conditions due to geologic and hydrologic influences include
(1) Impermeable or no flow boundaries, (2) constant head or hydraulic
boundaries, and (3) underflow or gradient boundaries. Program WTQUAL1
uses an initial water level coding to distinguish the type of boundary.
H(I,J) is the initial water level or piezometric head in grid (I,J) and
the coding used is:
0 < H(I,J) < 10,000 - actual water level elevation.
10,000 <.H(I,J) < 20,000 - impermeable grid.
20,000 <. H(I,J) < 30,000 - underflow grid.
30,000 <. H(I,J) < 40,000 - constant head grid.
A buffer zone is used so that irregularities in flow boundary conditions
at the physical boundary are damped out. Buffer zones are automatically
set in the program at three grids.
The relative concentration-(C/C0) of the conservative ion is speci-
fied for each grid as an initial boundary value and is identified as C5.
The coding used for initial values is:
0 < CS < 1 - slug injection of concentration C/C .
— 'o
2 <_ CS <_ 3 - constant source of concentration for all time
(C/CQ-2).
The maximum size of time increment, DT, which will provide adequate
accuracy should be used to conserve computer time. The optimum DT may be
determined by performing short period analyses with varying DT values
for the selected grid dimensions. Smaller grid dimensions may require
shorter time increments. The number of rows (NR) should always be less
307
-------
than the number of columns (NC) to conserve computer time during the
solution of the set of simultaneous equations,
In addition to space and time dimensions, (DX.DY.DT), the following
average or representative parameters must be determined for each grid:
1. G - Ground surface or top of confined aquifer elevation (feet)
2. Z - Bedrock elevation (feet),
3. PHI - Specific yield or storage coefficient (dimensionless
fraction).
4. FK - Permeability (feet/day),
5. H - Initial water level or piezometric head elevation (feet).
6. CA - Part of grid irrigated (decimal).
7. PHR - Phreatophyte use (ac-ft/yr) or phreatophytes present -
use to be calculated (unconfined case only).
8. IWELL - Well number if grid represents a well (integer) - may
vary per year.
9. IPIT - Recharge pit or line if grid*represents a pit or line
(integer) - may vary per year.
For leaky confined aquifer conditions the following additional
parameters must also be determined for each grid:
1. HL - Unconfined water level elevation causing leakage (feet).
2. TL - Thickness of leaky layer (feet).
3. FKL - Vertical permeability of leaky layer (feet/day).
To complete the model, the following hydrologic parameters must be
determined for every year of analysis:
1. PPT r< Precipitation (inches/year) - assumed uniform over
the entire study area.
308
-------
2. CPT - Relative amount of precipitation that reaches ground-
water (decimal fraction).
3. YPT - Distribution of precipitation for each DT during one
year (decimal fraction).
4. APW - Applied water as irrigation (feet/year) - assumed uni-
form over study area.
5. CAW - Relative amount of applied water that reaches ground-
water (decimal fraction).
6. YAW - Distribution of applied water for each DT during one
year (decimal fraction).
7. NW - Number of wells in study area (integer) - may vary
per year.
8. RPUM - Amount each well pumps per year (ac-ft/yr).
9. CPM - Relative amount of water removed from groundwater due
to pumping for each well (decimal fraction).
10. YPM - Distribution of pumping for each well for each DT
during one year (decimal fraction).
11. NP - Number of recharge pits or lines in study area (integer)
may vary per year.
12. RCHR - Amount each pit or line recharges per year (feet).
13. YRC - Distribution of recharge for each-pit or line for
each DT during one year (decimal fraction).
APPENDIX A-II. DATA INPUT
Input is by computer cards as outlined below. Detailed coding forms
used for data input are contained in APPENDIX A-YIII. Note that zero is
distinguished from a blank in data input and that decimal override
be used for all non-integer values,
309
-------
Card 1: 8A10 Format.
TITLE - Title for the particular run.
Card 2: 6I5,3F10.1 Format,
NR = number of rows in the grid system (NR should be less than NC)
NC = number of columns in grid system.
NW = maximum number of wells used in analysis.
NP = maximum number of pits used in analysis.
ICFAQ - 1 for confined aquifer analysis, otherwise blank.
ILKAQ = 1 for leaky aquifer conditions, otherwise blank.
DT • time increment (days).
ST = total time of analysis (days) (integer multiple of DT).
FWTOP = desired interval of times for printed output (integer
multiple, of DT).
Card 3: 8F10.1 Format.
DLX = value of uniform DX (blank otherwise) - x dimension
(feet).
DLY = value of uniform-DY (blank otherwise) - y dimension
(feet).
FFK = value of uniform KF (blank otherwise) - permeability
(feet/day).
ZZ = value of uniform Z (blank otherwise) - bedrock ele-
vation (feet).
GG = value of uniform G (blank otherwise) - ground surface
or top of confined aquifer ele-
vation (feet).
PPHI = value of uniform PHI (blank otherwise) - specific yield
or storage coefficient (decimal).
310
-------
CCA = value of uniform CA (blank otherwise) - part of grid
irrigated (decimal).
PPHIC = value of uniform PHIC (blank otherwise) - confined aqui-
fer (blank for unconfined aquifer
(analysis) storage coefficient
(decimal).
Card 3A: 8F10.1 Format - if DLX is blank, otherwise omit.
DX J-l.NC'
Card 3B: 8F10.1 Format - if DLY is blank, otherwise omit.
DY 1=1,NR
Card 3C: 8F10.1 Format - if FFK is blank, otherwise omit.
FK I=1,NR*NC by columns (example, the permeability of grid (4,3)
for NR«11,NC>11 is the second entry on the fourth
card).
Card 3D: 8F10.1 Format - if ZZ is blank, otherwise omit.
Z I=1,NR*NC by columns (see card 3C)~.
Card 3E: 8F10.1 Format - if 6G is blank, otherwise omit.
G I=1,NR*NC by columns (see card 3C).
Card 3F: 8F10.1 Format - if PPHI is blank, otherwise omit.
PHI I=1,NR*NC by columns (see card 3C).
Card 3G: 8F10.1 Format - if CCA is blank, otherwise omit.
CA I=1,NR*NC by columns (see card 3C).
Card 3H: 8F10.1 Format - if PPHIC is blank, otherwise omit, (for confined
PHIC I=1,NR*NC by columns (see card 3C). aquifer only)
Card 33: 1F10.1 Format.
CSS = uniform value of relative concentration CS (blank otherwise)
311
-------
Card 33A: 8F10.1 Format - if CSS is blank, otherwise omit.
CS = initial value of relative concentration.
CS I=1,NR*NC by columns.
Note:If CS(I,J) is less than 1.0, then that grid is treated as a slug
Injection of contaminants. If a grid (I,J) is to be treated as
a source of constant relative concentration, add 2.0 to the
value of CS(I,J).
Card 4: 5F10.1 Format.
HW * horizontal H - constant value of initial water level or
piezometric head elevation (feet).
LBC = left boundary code (single value along left boundary).
RBC = right boundary code (single value along right boundary).
TBC = top boundary code (single value along top boundary).
BBC = bottom boundary code (single value along bottom boundary).
NoterCorner grids are not critical. Use 10,000.00 for impermeable
boundary, 20,000.00 for underflow boundary, and 30,000.00 for
constant head boundary,, to denote outer boundary conditions.
Card 4A: 8F10.1 Format - if HW is blank, otherwise omit.
H I=1,NR*NC by columns (coded H values for boundary conditions).
Card 44: 3F10.1 Format - only if ILKAQ=1 (leaky confined aquifer condition),
otherwise omit.
»
HHL = value of horizontal HL (otherwise blank) - water level
elevation causing leak (feet).
TTL = value of uniform TL (otherwise blank) - thickness of
leaky layer (feet).
FFKL = value of uniform FKL (otherwise blank) - vertical perme-
ability of leaky layer (feet/day).
312
-------
Card 44A: 8F10.1 Format - if HHL is blank, otherwise omit.
HL I=1,NR*NC by columns.
Card 44B: 8F10.1 Format - if TTL is blank, otherwise omit.
TL I=1,NR*NC by columns.
Card 44C: 8F10.1 Format - if FFKL is blank, otherwise omit,
FKL I=1,NR*NC by columns.
Card 5: 2F10.1 Format.
PPT = precipitation (inches/year).
CPT = coefficient of effective precipitation to groundwater
(decimal fraction).
Card 5A: 16F5.1 Format - if PPT is greater than zero, otherwise omit.
YPT = distribution of precipitation for each DT during a one
year period (decimal fraction).
Card 6: 2F10.1 Format.
APW = applied water as a result of surface irrigation (feet/year)
CAW = coefficient of deep percolation of applied water (decimal
fraction).
Card 6A: 16F5.1 Format - if APW is greater than zero, otherwise omit.
YAW = distribution of applied water for each DT during a one
year period (decimal fraction).
Card 7: 15 Format.
NGPU = number of grids with phreatophyte use (zero if none).
Blank indicates PHR array to be read in.
Card 7A: 2I5,1F10.1 Format - if NGPU is greater than zero, otherwise omit.
I = row number of grid.
J = column number of grid.
PHR = phreatophyte use in specified grid (I,J) (acre-feet/year).
313
-------
If PHR is any negative value, phreatophyte use is calcu-
lated from ET subprogram, (see listing),
Note:0ne card for each grid containing phreatophyte is needed.
Total number of cards equal NGPU.
Cord 7AA: 8F10.1 Format - if NGPU is blank, otherwise omit.
PHR I=1,NR*NC by columns (see card 3C).
Card 7B: 16F5.1 Format - only if PHR is greater than zero, otherwise omit.
YPR » distribution of phreatophyte use for each DT during a one
year period (decimal fraction).
Card 77: 15 Format - if miximum number of wells is greater than zero,
otherwise omit.
NW - number of wells.
Card 77A: 3I5.2F10.1 Format - if NW is greater than zero, otherwise omit.
IWNO = well number code (integer).
I = row number of grid with well.
J = column number of grid with well.
RPUM = amount each well pumps per year (acre-feet/year).
CPM = coefficient of groundwater removed by each well (decimal
fraction).
Card 77B: 16F5.1 Format - if NW is greater than zero,.otherwise omit.
YPM = distribution of pumping for each DT.during one year
(decimal fraction).
Note:Cards 77A and 77B are to be input as a set for each grid.
Total number of sets equal NW,
Card 777: 15 Format - if maximum number of pits is greater than zero,
otherwise omit.
NP = number of recharge pits or lines,
314
-------
Card 777A: 3I5,1F10,1 Format - if NP is greater than zero, otherwise omit.
IPNO = pit number code (integer),
I = row number of grid.
J = column number of grid.
RCHR = amount each pit recharges per year (feet/year),
Card 777B: 16F5.1 Format - if NP is greater than zero, otherwise omit.
YRC = distribution of pit recharge for each DT during one year
(decimal fraction).
Note:Cards 777A and 777B are to be input as a set for each grid.
Total number of sets equal NP.
Card 7777: 1F10.1 Format - for every year of analysis greater than one,
otherwise omit.
REPEAT = data input code for multiple year analysis. If blank,
repeat cards 5 through 777B, otherwise program uses
same data as previous year.
APPENDIX A-III. OUTPUT OF RESULTS
Output is in a tabular form. All input data are reproduced for veri-
fication. Calculated results are printed only at the desired time steps
which is controlled by FWTOP in card 2. The exception to this is if a
grid changes from confined to unconfined and vice versa or experiences
flooding or overdraft, a message is printed at the appropriate time
indicating such conditions. Following is a list of output results.
1. Amount of core required to run program.
2. Title of run.
3. Problem control parameters,
a. Number of rows,
315
-------
b. Number of columns,
c. Time increment,
d. Total time of analysis.
4. Input data.
a. Delta-X map, spacing in J-direction.
b. Delta-Y map, spacing in I-direction.
c. Surface elevation map or top of confined aquifer elevation map.
d. Bedrock elevation map.
e. Specific yield map (confined and/or unconfined).
f. Permeability map.
g. Coefficient for part of grid irrigated.
h. Yearly applied water.
1. Yearly distribution of applied water.
j. Matrix of phreatophyte use.
k. Yearly distribution of phreatophyte use.
1. Yearly precipitation.
m. Yearly distribution of precipitation.
n. Well number map - if wells present.
o. Well table - if wells present.
p. Recharge pit number map - if pits or lines are present.
q. Pit table - if pits or lines are present.
r. Initial water table or piezometric head elevation map.
s. Initial storage.
5. Results - only printed at desired time steps by using the input vari-
able FWTOP. The exception to this is if a grid changes from confined
to unconfined and vise versa or a grid is overdrawn or flooded, a
message 1s printed at the time step of occurrence,
316
-------
a. Matrix of phreatophyte use (acre-feet of water used/grid/DT) -
only if calculated from ET subprogram.
b. Matrix of Q (acre-feet/day) - net vertical withdrawal of water
from each grid. This includes precipitation, applied water,
pumping, recharge, phreatophyte use and leakage.
c. List of grids, if any, at time steps where the main aquifer
changes from confined to unconfined or unconfined to confined.
This results in negligible errors of the resulting water table
or piezometric head elevation.
Example: GRID 5 6 CONFINED TO UNCONFINED
d. List of overdrawn or flooded grids at all time steps of occur-
rence.
Example:
ROW COL NR.
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
2
3
4
5
6
7
8
9
10
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
17
17
17
17
17
17
17
17
17
OVERDRAW
FT AC/ FT
162.73
160.72
159.33
158.20
157.25
156.36
155.43
154.45
153.31
151.92
149.94
145.87
136.86
119.12
83.43
12.19
12.19
12.18
12.18
12.18
12,18
12.18
12.17
12,16
12,11
0.00
0.00
0.00
0.00
- o.oo
0.00
0.00
0.00
0.01
0.01
0.01
0.03
0.03
0.04
0.05
0.01
0.01
0.02
0.03
0.04
0.06
0.08
0.14
0.22
0.42
FLOODED AREA OD. AREA FD
FT AC/ FT AC/ FT AC/ FT,
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00 •
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
317
-------
ROW COL MR. OVERDRAW FLOODED AREA OD. AREA FD.
FT AC/FT FT AC/FT AC/FT AC/FT
TOTALS
11
12
13
14
15
17
17
17
17
17
11.98
11,62
10.58
7.87
2.05
0.69
1.07
1.82
1.81
0.70
7.3
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00 6,1
0.00
The first four columns refer to the total area and the second two
columns refer to between stations or buffer zones. If overdraw
occurs, the piezometric head or water table elevation of the partic-
ular grid is set to the bedrock elevation and the program proceeds.
If flooding occurs, the piezometric head or water table elevation of
the particular grid is set to the ground surface elevation and the
program proceeds. Results will be in error if this message is
printed and should be checked.
e. Matrix of discharge (acre-feet/DT) between grids in the I-
direction. Flow down is positive and flow up is negative.
Discharge in the first row of the matrix is the flow between
grids in row #1 and #2,-and so on for the remainder of grids.
therefore, the last row is always zero.
f. Matrix of discharge (acre-feet/DT) between grids in the J-
dlrection. Flow right is positive and flow left is negative.
Discharge in the first column of the matrix is the flow between
grids in column #1 and #2, and so on for the remainder of grids.
The last column is, therefore, always zero.
g. Matrix of net flow (acre-feet/DT) from constant head grids.
The heading for this matrix is:
RIVER FLOW IN EACH GRID MINUS MEANS FLOW FROM AQUIFER (AC/FT/DT)
INCREMENT NUMBER 5
318
-------
h. Table of water balance computations (all entries in acre-feet/DT)
Example:
OUTPUT FOR SIMULATED TIME 6 - 9, AS AT THE
END OF'PERIOD 3 FOR CONFINED AQUIFER
APPLIED WATER (BETWEEN STATIONS -
TOTAL AREA) -974.97 -1282.89
INFLOW FROM RIVER (BETWEEN STATIONS -
TOTAL AREA) 0,00 0.00
BOUNDARY INFLOW (BETWEEN STATIONS -
TOTAL AREA) 138,18 88.48
TOTAL AREA STORAGE AND DECREASE IN
STORAGE 1332.68
BETWEEN STATIONS STORAGE AND DECREASE
OF STORAGE 62.42 637.08
STORAGE OF OVERLAP AREAS 1270.26
ILLEGALLY WITHDRAWN (BETWEEN STATIONS -
TOTAL AREA) 199.71 200.10
TOTALS (BETWEEN STATIONS - TOTAL AREA) 0.00 0.00
The first column is storage quantities for the total area,
area between stations or buffer zones and the overlap or
buffer zone area. The second column refers to the area
between stations or buffer zones and the third column refers
to the total area. Mass balance between stations (second
column totals) must always be satisfied otherwise results will
be in error. The only time that mass balance for the total
area (third column totals) is satisfied is when all boundaries
are impermeable.
i. Matrix of water table or piezometric head elevations (feet above
datum).
If data changes for a year of analysis, a message is printed
Indicating the year number and output 4g through 4q is printed.
j. Matrix of relative concentrations C/C
319
-------
APPENDIX A-IV. LIST OF IMPORTANT VARIABLES IN WTQUAL1
Symbol Description
APW Applied water ,as a result of surface irrigation (feet/year).
BBC Bottom boundary code.
CA Fraction of grid to which water is applied (decimal).
CAW Coefficient of deep percolation of applied water (decimal).
CCA Uniform fraction of each grid irrigated for all grids (decimal).
CCS Uniform value of relative concentration.
CMATRX Coefficient matrix.
CPM Coefficient of groundwater removed by pumping (decimal).
CRT Coefficient of effective precipitation to groundwater (decimal).
CR Right hand side matrix.
CS Relative value of concentration.
DLX Uniform X-dimension of all grids (feet).
DLY Uniform Y-dimension of all grids (feet).
DT Time increment (days).
DTWT Depth to water table from ground surface.
DX X-dimension of grid (feet).
DY Y-dimension of grid (feet).
ET Evapotranspiration (feet) - calculated from ET program.
FACFTA Amount flooded between buffer zone boundaries (acre-feet).
FACFTT Total amount flooded (acre-feet).
FFK Uniform permeability of all grids (feet/day).
FK Permeability (feet/day).
FKL Permeability of leaky layer (feet/day).
FVA Amount flooded between buffer zone boundaries (acre-feet).
FVT Total amount flooded (acre-feet).
FWTOP Desired time of output (multiple of DT).
320
-------
G Ground surface or top of confined aquifer elevation (feet).
GG Uniform ground surface or top of confined aquifer elevation
of all grids (feet),
H Initial water table elevation or piezometric head of aquifer
being modeled (feet),
HL Constant head value causing leaky aquifer (feet).
HP Water table elevation or piezometric head at previous time
level (feet).
HT Present water table elevation or piezometric-head (feet).
HW Uniform water table elevation or piezometric head (feet).
ICFAQ Equals 1 for confined aquifer analysis, otherwise blank.
ILFAQ Equals 1 for leaky aquifer analysis, otherwise blank.
INYR Number of time increments per year (360./DT).
LBC Left boundary code.
LCIE Right (J=NC-2) buffer zone.
LCIW Left (J=3) buffer zone.
LCJE Bottom (I=NR-2) buffer-zone.
LCJW Top (1=3) buffer zone.
NA Number of rows in reduced band matrix.
NB Number of columns in reduced band matrix.
NC Number of columns in model (integer).
NGPU Number of grids with phreatophyte use (integer).
NP Number of recharge pits or lines (integer).
NR Number of rows in model (integer).
NW Number of wells (integer).
OACFTA Overdraw between buffer zone boundaries (acre-feet).
OACFTT Total overdraw (acre-feet).
321
-------
OVA Overdraw between buffer zone boundaries (acre-feet).
OVT Total overdraw (acre-feet).
PHI Specific yield or storage coefficient (decimal).
PHIC Storage coefficient of confined aquifer (decimal),
PHR Phreatophyte use (acre-feet/year) or phreatophytes present -
use calculated in program ET.
PIT Recharge pit or line number code,
PPHI Uniform specific yield or storage coefficient of all grids
(decimal).
•PPHIC Uniform storage coefficient of confined aquifer for all grids
(decimal).
PPT Precipitation (inches/year).
Q Net value of hydrologic and artificial withdrawal per grid
(acre-feet/day).
RBC Right boundary code.
RCHR Amount each pit or line recharges per year (feet/year).
REPEAT Data input code for multiple year analysis (blank indicates to
read in data).
RPUM Amount each well pumps per year (acre-feet/year).
SQA Total Q per DT between stations (acre-feet).
SQBA Total inflow through buffer zone boundaries (acre-feet).
SQBT Total inflow through boundaries (acre-feet).
SQGGI Flow between grids in I-direction (acre-feet).
SQGGJ Flow between grids in J-direction (acre-feet).
SQR Inflow from constant head grids (acre-feet).
SQRA Inflow from constant head grids within buffer zone boundaries
(acre-feet),
322
-------
SQT Total Q per DT (acre-feet).
SQRT Total inflow from constant head grids (acre-feet).
ST Total time of analysis (days).
STA Between stations storage (acre-feet).
STATEM Decrease of storage between stations (acre-feet),
STOL Overlap area storage (acre-feet).
STT Total area storage (acre-feet).
STTTEM Total decrease of storage (acre-feet).
TBC Top boundary code.
TL Thickness of leaky layer (feet).
WELL Well number code.
YAW Distribution of applied water for each DT of one year (decimal).
YPM Distribution of pumping for each DT of one year (decimal).
YPR Distribution of phreatophyte use for each DT of one year (decimal)
YPT Distribution of precipitation for each DT of one year (decimal).
YRC Distribution of pit recharge for each DT of one year (decimal).
Z Bedrock elevation (feet).
ZZ Uniform bedrock elevation of all grids (feet).
323
-------
APPENDIX A-V. FLOW CHARTS
(START)
READ
NR, NC, NW, NP
ICFAQ, ILKAQ
,DT, ST, FWTOJV
CALCULATE FIRST
WORD ADDRESS FOR
ARRAYS TO BE PACKED
IN BLANK COMMON
CALCULATE LAST
WORD ADDRESS
PRINT
TOTAL AMOUNT
vOF CORE USED/
SET BUFFER
ZONE BOUNDARIES
PRINT
vTITLE/
t
PRINT"7
NR, NC. DT, ST/
324
-------
IS ILKAQ
LESS THAN
OR EQUAL
TO ZERO
SUBROUTINE READPH
Reads and Writes
Physical Data
Describing System
and Concentration
SUBROUTINE MATROP
Organizes Data into
Suitable Form for
Printing and Prints
SUBROUTINE READH
Reads in Initial
Water Level or
Piezometric Head
Elevation
SUBROUTINE LEKAQF
Reads and Writes
Leaky Aquifer
Parameters
SUBROUTINE MATROP
Organizes Data i«to
Suitable Form for
Printing and Printf.
SUBROUTINE STORAG
Computes Increase
or Decrease in
Storage
LOOPYL=ST/DT
INDX = 1
325
-------
fcropV—
NO
/
4
V
' 1=1 \
Is I.LE. \
LOOPUL /—YES-]
1=1+1 /
1 '•
c
r
IDT=DT
J1=1DT*(I-1)
J2=IDT*1
FI=I, TI=JI, T2=J2
SUBROUTINE QFIX
Reads and Writes
Hydrologic Param-
eters and Computes
Hydrologic and
Artificial Inputs
for Each Grid
SUBROUTINE MATSOL
Sets up Coef-
ficients Matrix
and Right Hand
Side Vector
Matrix
FUNCTION ET
Computes the
Phreatophyte Use
Using Water Table
Elevations
SUBROUTINE MATROP
Organizes Data into
Suitable Form for
Printing 5 Prints
FUNCTION PARAM :
Computes Coefficients
For Finite Difference,
Equations ;
SUBROUTINE NSCONT
Checks for Known
Grid Values, Such as
Boundaries, and
Transfers them to
Right Side Vector
Matrix
SUBROUTINE BSOLVE
Solves Banded Matrix!
Set Up in MATSOL
by Gauss Elimination
326
-------
CALL
BJUST
SUBROUTINE BJUST
Adjusts Under Flow
Grids Water Level
or Piezometric Head
Elevation .
CALL
ODFLOD
SUBROUTINE ODFLOD
Checks for Flooded
or Overdrawn Grids
and Writes Results
if Any
•YES-
INDX=INDX+1
STTTEM=STT
STATEM=STA
CALL
STORAG
SUBROUTINE STORAG
Computes Increase
or Decrease in
Storage _
©
327
-------
SUBROUTINE BYFLOW
Computes and Writes
Flows for Each Grid,
Thru Boundaries and
from Constant Head
Grids. Computes
Relative Concentra-
tion Based on Mass
Balance
^ rAKAM p
FUNCTION
PARAM
Computes Coefficients
for Finite Difference ,
Equations :
CALL
MATROP
SUBROUTINE MATROP
Organizes Data into
Suitable Form for
Printing § Prints
1
CALL
BALCOP
SUBROUTINE BALCOP
Writes Out Balance
Computations for
a Time Increment
Set
HP = HT
328
-------
APPENDIX A-VI. DESCRIPTION OF SUBPROGRAMS
Subroutine READPH
This subroutine reads and writes the physical data describing the
study area. The following variables are read and printed: DX, DY, FK, Z,
CS, G, PHI, and PHIC. CA is also read but printed later. Coded values
of CS are printed. Only one data card is required if all variables are
uniform for each grid, otherwise each parameter that is variable must be
read in matrix form. Variables DX and DY require only NC and NR values
respectively.
Called From: Main Program
Subprograms Used: MATROP
Important Variables: DX DY FK Z G PHI PHIC CA CS
Subroutine READH
This subroutine reads the initial coded water level or piezometric
head elevations. H is decoded and set equal to HT and HP. One data card
1s required if the initial water level is horizontal, otherwise the entire
H-matrix must be read.
Called From: Main Program
Subprograms Used: None
Important Variables: H HT HP
Subroutine LEKAQF
This subroutine reads and writes the leaky aquifer parameters. The
following variables are read and printed: HL, TL, and FKL. One data card
1s required if these variables are uniform, otherwise each matrix that is
variable must be read.
Called From: Main Program
Subprograms Used: MATROP
Important Variables: HL TL FKL
329
-------
Subroutine CSET
This subroutine initializes the relative concentration throughout
the aquifer.
Called From: Main Program
Subprograms Used: None
Important Variables: CO CT H G CS
Subroutine STORAG
This subroutine computes the initial storage and increase or decrease
of storage. Total area and between station (between buffer zone boundaries),
'storage is calculated. Also storage of overlap areas is computed.
Called From: Main Program
Subroutines Used: None
Important Variables: STA STT STOL H HT Z
Subroutine QFIX
This subroutine reads and writes the hydrologic parameters. The
hydrologic and artificial inputs are then calculated for each grid. A
value of zero on the input card Indicates a particular parameter is not
used (see listing). The exception to this is the number of grids with
phyratophyte use, NGPU. If NGPU is blank, the entire PHR matrix must be
read, otherwise the number of grids specified is read. NGPU equal to
zero indicates no phreatophyte use. -
Coding PHR less than one indicates that phreatophyte use should be
calculated every time increment from the previous time step water level
elevation. The ET subprogram is used for this.
The factors considered in QFIX are (1) precipitation, (2) applied
water as irrigation, (3) phreatophyte use, (4) wells, (5) recharge areas
or lines, and (6) leaky aquifer conditions.
330
-------
Called From: Main Program
Subprograms Used: ET
Important Variables: PPT CPT YPT APW CAW YAW NGPU PHR YPR WELL
RPUM YPM PIT RCHR YRC Q SQT SQA REPEAT CPM
Function ET
This subprogram computes the phreatophyte use for each grid using
the water level elevations from the previous time step. If the depth
of water table DTWT is negative, an error message is printed. It is
anticipated this program, if used, will change with each study area.
Called From: QFIX
Subprograms Used: None
Important Variables: ET DTWT
Subroutine MATSOL
This subroutine sets up the coefficient matrix, CMATRX, and the
right hand side vector matrix, CR. CMATRX is a reduced matrix containing
only the band of known values in the left side of the difference equations
and is written vertically rather than diagonally. Its dimensions are
(NR-2)*(NC-2) by 2*NR-3. The coefficients are computed using Function
PARAM and checked for adjacent boundary values of H in subroutine NSCONT.
MATSOL treats known grid values of H. BSOLVE is used to solve the matrix
equation set up.
Called From: Main Program
Subprograms Used: PARAM NSCONT BSOLVE
Important Variables: CMATRX CR
Function PARAM
This subprogram computes the coefficients in the left side of the
finite difference equation. For confined aquifer analysis, saturated
331
-------
thickness is compared to aquifer thickness and the smallest of the two
is used to calculate the coefficient.
Called From: MATSOL BYFLOW
Subprograms Used: None
Important Variables: PARAM
Subroutine NSCONT
This subroutine transfers the coefficients, in CMATRX, multiplied by
their respective H-value, to the right hand side vector matrix in case
of adjacent head or known boundary conditions. It also- sets coefficients
equal to zero in case of adjacent impermeable grids.
Called From: MATSOL
Subprograms Used: None
Important Variables: None
Subroutine BSOLVE
This subroutine solves the matrix equation set up in MATSOL by Gauss
Elimination. BSOLVE is designed specifically for a diagonal matrix that
results from analysis of groundwater systems.
Called From: MATSOL
Subprograms Used: None
Important Variables: None
Subroutine BJUST
This subroutine adjusts the underflow boundary water level elevations.
Gradients are calculated three grids in from the exterior boundary grids
and the gradients are projected back to the exterior boundary grids to
obtain new water level elevations, This calculation is performed at
even time steps. At odd time steps the water level elevations are held
constant and the exterior boundary grids are treated as constant head
grids.
332
-------
Called From: Main Program
Subprograms Used: None
Important Variables: H HT
Subroutine ODFLOD
This subroutine checks for overdrawn or flooded grids. If either
should occur, a message is printed indicating such. For confined aquifer
analysis the flooded grid computations are bypassed. Total flooded and
overdraw amounts are computed for the total area and between stations.
Called From: Main Program
Subprograms Used: None
Important Variables: OACFTT=OVT OACFTA=OVA FACFTT=FVT FACFTA=FVA
Subroutine BYFLOW
This subroutine computes flows for each grid. Total flow through
model boundaries and buffer zone boundaries is calculated as well as flow
Into the system from constant head grids. The flow equation used is
developed from the finite difference equations and uses particular values
of the CMATRX. These values are" transferred from MATSOL except for
boundary values which are calculated in BYFLOW using Function PARAM.
Flow is not allowed to or from an impermeable grid and between any two
adjacent underflow grids. I-direction and J-direction flows are printed
and flows from constant head grids are interpreted atid printed as flow
from river grids. Relative concentration calculations are made using
the flow between grids.
Called From: Main Program
Subprograms Used: PARAM MATROP
Important Variables: SQGGI SQGGJ SQBT SQBA SQR SQRT SQRA CS
333
-------
Subroutine BALCOP
This subroutine writes the balance computations at the desired time
steps specified by FWTOP. Mass balance for the entire area cannot always
be obtained, due to accounting procedures used to compute mass flow at
exterior boundary grids. However, for between stations, which refers to
the area between the buffer zone boundaries, mass balance must always be
satisfied except for the case when a confined grid becomes unconfined.
This error should be small and is indicated by the "TOTALS" in the mass
balance output being different than zero, to reduce this error, decrease
the value of AT. For confined aquifer analysis, a message is printed
indicating if a grid becomes unconfined.
Called From: Main Program
Subprograms Used: None
Important Variables: SQA SQT SQRA SQRT SQBA SQBT STT STTTEM STA
STATEM STOL OVA OVT
Subroutine MATROP
This subroutine organizes data or results into a suitable form for
printing and then prints.
Called From: READPH LEKAQF QFIX BYFLOW
Subprograms Used: None
Important Variables: NR=NOROW NC=NOCOL
334
-------
APPENDIX A-VII. PROGRAM LISTING
PROGRAM WTOUAL 1
..,V,-i1 WTOUAL.1 i:f.f.iI,-.iir:....T.|i;-.-,,T-,f^1.T<'.-:.'J-...
r • 3J
0 ,-cr.-;M u c.i'-iu /.MT; IC.,L .•;.* HMWL'.-MC x-.fi.uE-.C-- «•< v«ou'O u'.
C K-1TC*. Lr-VcL. f, ;
c
r
c
C scr.
I -ES:."?S j?",:,,"::^,^;.?^^:^1^^. gj
K.IS s;-fs;M!K'/?ss.srs3;^«s:l'» ;i; ;s f «ff Sj.
C PFi' -'-C F" ^FCSIItfE .JSL . NK.CT.7 A»3 JU..GI.7 . iv
c
c
C U'-JhsTOf Itl nurr1: ; 22C
c. LC,JC«ooTTOM t:» turrEj. /ONC !3,3
C FOR rxPuI^TJOr. OF VAP.KBU3 «'0 S'JdWOTWfS, SEE SUMOUIINES. «6 '
C 2ft.
C COMT30L VAPIi.TL •
C NfjEtiuHyf.p or
C NC=M).1'T:i< PF .
C N^ SIICULO 4L-M3 « LESS THltl 0». £ IUAL TO NC .
C »a=H4XI.-!U.-l f.U''.i-ER C- 'J£LL'> Jtfl
0 MPrHAXIMI/.l NU1UC.R. OF •'LC'n?,'JE PU> _Tll-D,.T..r ,,.,.„ 320
C ICFtQ-1 F:j? C3N-INCH iOJIF-'I S-.SLYSI3. OTHiR'yist 3L-h< «B
C IL«0-.l FJ5 Lt;4
JOC*NK»NC ' . 59^
T.FK-1 . 6M
JPHI=ID',<1 • 61i
JCA*fc«ICCH 6*1
1M = 7«IOC»: (.70
JHT-1'IX'l , ' 05i
jM«'si»I?CU (,.};
io«:3*ncu .. ;,-;
:»i • 7i!
l u;
. 73i
ir.3
-------
I-H.MS-M I'.* !'.*•• i* 1 'C *10
ir*L=i U< iu •'/•'. 'I!",'; '''I
ura,i» ?*• ::i-»K' 0'icci I L<•:
lftM»:,-vin J*'«
It PIJMS T, *<•"••;>( *'••
TMuHssl *fcC«.-i"*!'•''-' '*•*'
,E!.J14:.I..»HfM«» IYP»."r!WlN.>'I«t^».irii»-«EMl *'«
e''J'"'"M'"*''JI 99;
lHA*LCCf ICtin.'X'U 1DOC
CSLL COS-. UK: i "20
CO 13C LI*ltIfUC3
CIL1 IsUNPEf
ico CONTIKU-:
c io>o
I 1HJ
ic-e M"
LCJW'3 1060
IP.ITE «6,21fl TITLE »"*
IF IILKAC.Lc.3) 60 TO 113 11"
WHITE U,25i» Hi3
60 TO 13; Ji?l
lid IF (KF40.L-:. ci GO TO i?: 113^
W TC 13C ' if.;;
tz: W..ITE c6,2&:i . u..%
133 KKITE CStl'K) NR,NC,BT,5T lir.-
C • in;
CALL p.fto^H er'»,i«c,ct:?<» .c (iPHD ,: tCt IPKD) 1«:
C
1*C CALL STC'iG llt(?,MC,C(lHI |!7« IHT ), C «I Z» tC IIOXI ,C (I3Y ),C(IPHI) ,(
lCtl"HIC» _ . 12SO
C»LL CS£r 131C
• FI«I . ' 1380
• IMC
HO 3
c u:.
CALL QFU inr;t'ic,c(nx».c( i?»<,;tic;.> .:.-« I7> ..:t»Ht i .i.ct:: i«.?:
l I.IJK.'.C.C i ;pt"*i ,r iirft LLI .c <:-n;i) ,c dri'T) ,n (lYAwi ,r.i
3 tlIC>,C <:*>!.3TK.>) )
c
l x> tcd'jvi ,r.iK;I ,c,ciio 1333
i TJ.cii'-mcti t3<«o
C ISGi
STTT-HrSTr
CALL ST09J.O
-------
t. i '>''
POM-ITU Ifcl!
IF i (r r • n i.'[-. (t C'.f *F V*TT>I i oo TO i'JJ u>C6
i-, u
C . '."'•(•
,M«.L "/FLO'* t .:v, • i:..i., j i,.: «;; •• >, ; i: u ,•• i r «ri.: 11, i, .i; j* i .r i: i» :-.;.',
i i.«. u ;T.'.I i n, <. v .M M .; i: :•). ".c ;••••-•.r '<),:.<; •.) ,i «i i), ;i I.HM ,i •; i-mu
in):
c 1/1:
H.'ITC «ti,l«J> fi 1/10
CALL MU •<•>••• <>;•<,.MI.CIIHTM iJ?i
C IJ'J'j
• 1 A
1 So "C i - * * ' *"
00 Iti L=1.X . 175C
00 165 K-t.'W l/bv
itc ccMirue !?">«
rro? i»i«.'
c i*2r-
tT AdCVt UiTJH>»
t90 ro=irtAT ii5H-^c« nin;i.r;ofi*T'ti2iH COLUIN ijlfcNSior.si^.'iiM TIM ii5c
1C INC«£."i:'4T20 'JY -j'i.2, = '< OA/3.27H IOI4L TIC£ OF £t;«LV'• ij> G9.2.5 U
8H nevsi • H'M
U8.
230 FO^AI (IH ,///, 3cn FIELO Lifir-Tn = ,0231
Z?C FOT1AI CJH-,!,B<, 22HL£-(
f.-^t • , t. .... » >- .-•- ,o . ,*„»
C'^'IC"' /3LTtJT>I!^""A'),[L
9°
«iP
ilP
«p
HP
QP
JY KP
np
.{P
tfp
yp
s«'
fp
ro
JP
4"
;p
Ri-
i(P
HP
HP
HP
'P
ftP
2:
JC
••C
SO
60
73
di
9J
too
114
12 C
llj
IOC
i50
164
t'l
is:
Hu
20«
216
ZtQ
23S
^'iC
2'<0
Z7C
510
343
J35
i: j
J2C
3K
i«i.
?%3
jre
jij
330
<«C i
<.!•;
W2C.
337
-------
I-'." '" U,J» if. X
13 TO f.J
no --.in lOt'i^di loxttijii j=tfNci
09 me I**,**
no itc j-i ,NC
ifcc 'j*ti,j) = o< (i ,j>
tSJ COUTIhUE
IF (Qir.L £.:..} I 1,3 10 170
00 Irt J-l.MC
no i%u i«i ,;,y
ICO fJfCI.JI --OI.Y
co ro ma
170 «t*o (t, 1.1. '.i (ovii.ii. i n.-ai
13 MC J--.VIC
,n
RP U30
t<»C COI.TIf.U3
216 •'.H
If IFFK.Lf.l.O) GO TO
on 2JC j-i.?.^
00 2~r : «flO
PP
A-3
A.
833
8*:
450
««• 6JJ
»P SCO
°.p 91 :
RP <)25
»P 933
RP 3% 5
RP 97C
SP *»C
SP <»'JO
«!P 1030
*.P ICU
RP 102U
RP
RP
RP
f.T TO 3J
ruiO l!i,
r.o;iriMj£
HRITE
CAUt
H'.IT'I
C»LL
HM.NC.CSI
TO
; (c.1.70)
GO TO i.2i;
P.P 10)0
P.P nac
«>.F 1UO
*P 1120
«F :iJ5
*P ir»o
*:• I15u
*P 117C
PP I I 4 11
338
-------
CALL >'ATi-Ot- tf.w. NC.FHil
TF 4ICi'»O.L.-.:. C> 00 TO <«SO
CAUL HAf-OP IM.I.C.PHICI
'.x.
A-X IUP,3°'.CII1G -iCHO^S IN J-OIEECTION
-t .lif, SPftCrHC DCKS IS I-OIKiCTION
»
N M-P (FE£T «riOV
51? FO'L'UT o
-- j.:ic j;t:_ ir* .-HP (FttT/oayi
C= OF C.irif i:i£-l AC'JIF'.K ELi
f.f
P.P
/: JvOO .< 1
CO T? U;>0,17C
NT iiiJi -"i i .j»
CC TO 20J
HI iiiJ»"i -icc:o.
CO 10 2Cw
HT«IiJI-"M,JI-iOOC3.
ci n t'ju
HT(I,JlrH( [, J)-.UCCO.
Hn(iiJI = M II. .M
»ir lit j> »n (It Ji
.
-------
C . . HH 55J
J20 FOVIAT ieFll.il • ' AH »<•,(,
C RH 5/0
SUBROUTINE LEKAOF
LA U
C LA 3C
c rt'THic< )L cf OF Lfix» L/w-.f: (>"V-,TI " ie 60
CUM fg
LA SO
IS» LA SO
LA 13 i
00 110 1 = 1, 'if. % U 11C
00 lit! J = t,:iG LA 12)
rtlUiJI'J.g L1 jj-,
TLd.JMJ.C Li ua
FXLII.JMC.O tt 156
r.£/0 <5,2tC) HML, f TL.FFirOLITI»E C^MTI,TF.$ f'<; IN^ ?i;Ai'E OR DCCRfASi IM jTOR-Cc. ST VJ
C STT = T'IAL'\s-I *F' ST 3j
C STJUiCVf«L«I% JRfA IIC-'Ai£ (AFI |j *J
s sr .*,
O'.nenSKN KIN; ,HCI , ».MNP,HCI, ZI :j*, fie i, c»lf.P,N:», i)M'l?,nj», -;<.j ^r -ij
,, ' ' " '••''' • it»>-s-t'.CI . S| jj.
ST UC
ST 123
>-Jvr sr ui
KC1>HC-1 ' 51
IK t,» .HCl , • sl
1 If Kr2,-J?.« . -
IF »<»;,L ».r.r,ucio. t •••) ro uo 3*
IF Ut.F
-------
Ir li*TK,U.Lr..i,«,LI) 'iJ TO 110 ST f»j
STPs I (•'.(>; i LI- I (-KiLH* j(K,L I I lfirC(«D«G(KtL»*$t I <,L > ) «J< «<. L >"' ST {1.9
X YdC,L)/'.3t.ftJ. . 2T Z3L
GO TO 1?0 ST 26'
110 SIPMMI C<,UI-/(v,LII-U««,t)'jrU.L»»S«K,L)/t «s^
IF iL.r.T.LCh ) co rn ua *t .,ft
IF C1.LT.LC.MJ GU 10 131 ST 31?
if «.cr.LCJc> GO ro iii ... ,t.'
St/UMAtiTF „ — 3
t.io c "*
»• • $T
ST ,u
EWO ST U3
..'.... sr j*c
SUBROUTINE: CSET
ZVltKOVllt!* CSeT I"P ,f C.CO.CT.H.IJ.C-n
01 «f N£ ICt. C.IKViNOJ . CTIfi't-JC) .Mf'io ,'
OJ m jsi.::c
30 t3C fslt.'i!?
CMI,J)'O.C
IF tCSU ,Jl .CT.2.C» GO T3 11
COd.JI -CC([, J)
GO TO 130
ic tr,«i,j)=cs(i,ji-z.o
ice cc'irifuit
KtlUKN
£«»
"SUBROUTINE QFIX
pT, OF ic
OF 20
C OF IS
C . OF <,fl
C THIS SU9ROUTtHE COMPUTES f(-:C I°I T JT IOU TO G^OONJWITE^ IDtCIHlLI OF 7(1
C YPTsOISMiaunON CF PStCIPt T^TIOM FO^ EACH OT OF CNE V£i^. OF :
C APHsAF?LICn rl-T£R iS 1 ?i'iULT OF SUfFiC? I?.t;i31 TIC't (F-riT/YtSR) OF 102
C CAHsCCiF. Of CE?P 0-r^coLiTtos OF a»PLI-:0 MAfcx (OCCIMiLI OF ItC
C YAHsOIS'FMUTlON CF APPLltJ HATE'* FCK EACH OT OF Cf*£ Y£AR OF 122
C UieCI.-IALI Qf 130
c N6Po»».u«. THIS MiY a£ OF l5C
C CALCULA'^O FP;C^ rrt£ £T S:;8P»..')G'?i.1 •!/ COOIUS ?tlr? VALUES Qc l%ii
C . LESS THAN ZFKC. OF 17 C
C Y"R*OISTS'.t5UTION' CF PHRtATOPMYfE U?E *0i». EACH OT OF OM£ Y£A.» OF HC
C- - CCtCI*ALI ~ QF Wj
C. v HELL»WCUL HUrtfE'i CCO^ . Of iOC
C «PUf'--'irn;:ir .-.ACH WELL P-J^OJ °ER YEAR CAF/Y«A^) • or jic
C CPH'COtF. OF GRC1UM">t«eTe f ^ilOVE'l .IV PJ-fPINO IDHCI^iL) OF Z2J
C »PM=OI57C.I5UTICN CF auMP(.-& PCS tACH 3T OF C.'Ji TEiK. «OfiC I!1AL > QF 33*
C PITsC;CHA-£3 OT (AF» » OF 29C
C 30»*T3TAi. n Ptft Or PiT«rf.;< STATION'S I.'.F) gr 33;
C »;0EM*OArA INFUt CO.)!; F0< 'fllLTI'ti ft.^ AfULfSIS IF 313
c "IF ^^';
nn OF no
3). HCHr<(i;^). t-tCK.^CI. TL (.•;'{. \CI i r,5C
»«?ic> OF jr;
C - OF Jflt
Lc:r:.ii:m.LCj?.tLOJii.F'tTG9 OF 39i
-jrij'KMisrjKi, jjf»l OF !.<»:
itc no izt
-------
HULL -C .1
PCTU, (.)-}. 3
lid COCTlNUi
oo I.K lu
TPTUI--3.1
m CO-0,fl
IK COtJTII.u'J
Ff r*Mt < t, ji-r. GO TO ISO
• «SAO (5, -.1C) t Y-u«) ,< = i , I
15C CONTI.MJi
ft tO ISiH1) .V.PU
IF OI&PU.EU.JI r.O.TO 1/0
OJ 1'iC Af-stt'^.-'U
stsAO «5,«'n) J,K,P
03 rc isc
170 IF CSIC.'tll.Q .MfilMJI .ST. 0.0 I 3J TO
1« -.0 t-90 «=ltN*
03 130 LU-li'lC
IF IP»^«KiLL) iLT.J.]) r.o TO Z10
IF {PH^KK.LO.&T.C.C I 03 TO Z6J
190 CONTINUE
ZtO IF rtH.lf.OI GO TO ZIC
na
00
C3MTIHUE
00 2<»C J=l ,'IH
I5,30C) !'.VNC,<,L,i?
,CPM(J»
15.M3)
IF IVP' I .0-:. 5.3) 50 TO
oo ZJS <»i,ir.r*
Ztfl CONTINUE
ZSC IF (SF.Lf.iTI GO TO 293
fi£« I!
00 ZftC <*!.
03 Z« 1=1,
ZftO CONTIHUr.
S.??1)) IPHC.Kt L.RSHP(J)
15. -«.1C • l»SCIl°'li)1i».<(•:.J.G> 50 r(
oo zu <--i. tf.r-? (
2 ,<
K <£,-•'..'> Pf'T
W'Ift, (
C/.LL «A
1.. J> C.O TCI JC
fl.9,!lCi 'Xf LLI
j=t.
295 'I1 IT- (6,fv--,OI J.-JPU-H J) . IVifJ.O ,Kzl,I'
3CC I' «:;i'.i.1:. 3) r.o t) Jiu
If ITt «^,-5J3l
C'.LL t".4fi"ip ti,.j,f::,pn.)
5U
IF
OF
OF BZt
OF 51J
OF 51,1
OF 530
si:
OF
Of
OF
OF •,)•)
GF TlC
QF 60
OF 77-]
OF 7*0
OF 71J
OF ajj
OF 410
QF »2Q
QF 810
OF ft'. Q
OF 454
QF 863
OF 870
QP C
OF 1470
OF 134C
OF US)
QF UOJ
QF ItlJ
OF xizu
'JF
QF
OF 1153
GF 1166
OF 1171,
QF H«:
OF mj
OF 12CC
OF
3F
OF
OF
IF
OF
I MO
342
-------
sir i* (l.i.£. vi r.r TO 320
3 H'-iti t *,
f.LL liftt;':? (!N,'i;fH»
U.hJCI SK.iiTCl.fSTT
10
I2U KTslLTU
;r UCr.LF.lM*'
K.T«ICT-rit.<-l
Kt»0 IS,"./:! rtKftT
IF |?EPf.*T.'tE.O. J» HO TO 321
• c.'.lul
CO TO 110
J33 P^nT=
. F;*I
S(! 4=0.0
.
IF IMK.U.CT. :.£<.) CO TO V.3
IF |»-T«<,L).GT.ZCK,Ln '0 10
A»C.a
IF JKCT.GT.il) fiO TO 3ltO
K9ITC I6,5i:i I
H = 3.3
CCfiTINUF
CTCHT*!.
C3 1C ^1
CONTINUE
jJ=W£LUKtLI
IF (JJ.LE.U CO TO <«1U .
0»'Ort=S.PU-1 oo TO v.3
IF IK.LT.LCJKI GO T) <|V)
IF «<.r.T.LCJi£> uO TO ^!.J
BOTO
IF C£TCf-T.L£.C.» C-C TO <.53
«,oc K-
%7iJ FOCMAT (•*r10.1l
I..1C FO^HJt tl'TS.H
l/'!5.7>" 1C .It
52- F'.;t-At (ini.5/<.
««j Forur !i,it,«.j<,«r.i^Ki.
^r r.v..i;.i (ii-t M: mo>A«Ai.
555 FI'MAT (l'»l,30t, .'JH'tA
1C I
FOR
run
OF U»C
OF
OF
OF
OF I MC
OF IJ'K
,JF ivj:j
ii,c()
OF
OF
OF
OF
OF 1510
OF 152C
OF isic
OF 15"»C
OF 1550
OF I56i
OF 157)
qp 15,C
QF
Cf 16;J
QF 163C
op i6i»(,
"
QF
Q(. I6.j.,
OF 171',
QF 1720
OF 17.I&
OF 177C
OF 17s;
flr
OF
OF
OF 18JZ
OF 1*30
OF 1^1^0
OF ift^c
OF IJ55
^
QF 199:
«»0 OF «CC
OF 1120
Of 195C
cr
Qr l!)f>c
or j^/j
OF
Z052
QF ,
OF j
QF 7
CF ^
OF 213C
OF 2113
OF 2123
OF
OF
343
-------
•.••:•-.i;ui « OF
. 56C F37M1J :*,6Cii"*f°,ix nr "HKc Fif'MA'i (:•<:.///iJ/x. ci>itt'.i\Lf f;'.rci''iT«riCNa,e'i;.3»5x, IHINCMCJ/ OF >uj
iv,7A*> • - of 2)(>o
•i*. 1QHMELL TA>1LEI OF 2J50
,r,-x. 4MH-L!. no..''<,iio ,//,5x, ifitkaTt UC-FT/Y-:A») or 216C
ll,*)X.F:C.i.,l / ,5x, K-'-i.'! >T"I A.TIr.'J. JX.liF/.tl I 'IF
67J rc^H-tT (t^i./////,5>', iuH».tniii-:".c fiii.xi»iOi»owiNCcri.i , o=LV«. t t'S uc
20(MU»OU.I.JCJI.> i C^.Tr 41 tP,!'.! . C'ttf'l. i(!.0-Oll.-
-------
™ "" MC 1T"1
C«Mllll=0.3 KS 38J •
IF IMU.JI .r.c.ict :c.c» ca TO 151 ««s iJ>!
J4*l "S -.OC
40*1 M* *ii
c lEfiui MS die
C ' 'IS V«3
li. J-l) iUtL<(ItJI irtUM Ji, J-ll ,r.rLt IIiJJ .G( J'. i J r«5 ViJ
HS H7j
C ' HS I.9J
C • TOPIRI MS .viilji t KLX< t-:»ji ,j;.< t:i.ii iCti-iiji ,'> 'i*. 3j;
2 (IiJM . 1i 53;
C IS 51.3
C BCITO'KO 1S 550
C >IS SSJ
CMATK*cir.rc> =p*f ^vir "s s?"i
C HS 63fi
I»H iF<(It Jltf TC J3> J«ll trtTtlt Jl »Z I lit "IS £''«!
ILK U, Jl iBCLT (JO,J»1I ,li£Lf .GC,J)> MS 661
9111JI = C~*4'-'X CJT ,IOt IS 670
HS 64 Z
hSCOHT (^( Jt,J-l l,'*T(Ji,J-11 ,'IT(I.JI.7.IJ1,J-U ,i(I, JI.C^iT 15 age
kM> US 7 .MfC n, J) ,HTCI .Jl i2J)l/OEi.ri MS ^7C
-IHT (I,JI •STC3EF«0-i.Xll , J)»0£i.r«I . J It /0:LT—J 11. J > • 'IS 6S^
HS 903
CO TC :<>C . HS 91C
HS 920
HS 930
1*5 CONTINUE MS 4<»3
BEMIN3 7 US 05*
tMITE (71 C'lAT^X.CK . "S 96C
CSLt CSCLtfi (CXATiiX,I!>,IR.C«l MS 97J
MTrJ . «S
-------
CST ,N"i) K ».ii d , ji -Pci:; i; , ji I'Oc.t U .Ji'OS "1 i:i«
i LYII.J>/>;:I.T us t^ac
£ lf-,?'.5l I ,J • MS 1211
co T: 2i~. is :?j
100 l9VQ>r'.ii »<*>iici itJi-o'icitji Moau 1 1 ji • x -0
PO ^<«C J=2,»Ct • «S 1!53
00 223 I = :.f.Rl MS tJfitl
NTsMTU . • US U7:
Hri:.,ji=cv(! MS 13*5
MFCI.JI -C'-fC.!l ... »t3 1 J95
??0 CO'irifJUE « 11.3;
23! 3u'ru>?h • ItS tdlj
c MS mso
?*0 F3?M4T CH .li.HXi «.HMIO,»I3,SX, :>lf(IMCONr IhiO TO CG;6c
FUNCTION PARAM
FUtiCTiOrl P«=i'>l',iri,ira,A;i,4G2l PR ". 13
J PR ZO
c THIS suuPsooRiri co.ipurcs THE coEFFicunn USEO IN MAT-JOL AND PP HO
c " ofFiou. :T ii APFL:cii:L; TO CAS,:S OF MRIVIL- nx, CY. F<, PR 53
C ArO SATU^Tt.0 THICKNESS. . PR 60
C PR 73
COM.tON /9LK1/ OT|SI,ICF40,ILKAO,LCIE.LCIM,LCJC»t:jW,FKTOP " PR 80
c PR 90
IF ucF&a.L£.('i r,n ro uc • PR uo
03 10 120 PP J4, j
110 Sf TMCKO'liXl (iHTl, AHT?)-A1i NC 1C
C UC 20
C UC JO
c THIS suapouri'ie niN$f£RS Ttie cotrricicNrs, MULTIPLICO JY THEIR NC CC TC 1 U ' NC ISC
12C IK" (IHT:<-Z>t|.l£. 1.0.4i-!O.Hr 1.0T...r:.MTi'l) OT TO l'«3 . UC 170
13! CP*5 = C.J HC Hi
C . • HC VK
l«tfl Ktruiill • - NC 2)0
C »:C 210
CNO . HC i2fl
346
-------
SUBROUTINE BSOLVE
C . w» 1 t
C ns iv>
C THIS i'tH?OIJTU? SELVES MIE rt»TkIX, 3ET U? IK CftTSCL, :JY GiUSS «S *)
c rtisxwrioM. „ ...
• '
oo m i-i.L* - - .2! .?.c
irui.M.M « ":
DO 120 t--=l.IK . 'I "i
00 JK J^.M "
-
00 19C 1 = 1, t'l «
WXM
00 «1B
IF
13C CONTlNJi
• IF u.p:v.',r..ti oo TO 150
03 XnC J--t,,l
SUBROUTINE B JUST
OII,JI • 21 .,«
00 .17(1 J*Z,« • "
IF IUP.LT.M URalR+1 •' 5 *^
I9C iOMTIMUS . f* J«
v(Ni=vtM/o(n, ii _ . . • r* 7*7
Jfts2 ' ' **w*>
CO 216 1*1, 1 H ' *f J'J
l*N-I ' ™* *'3
00 289 J=2,J« - ff i?^
KM=L»J ' . I! ??i
23C Vll>sVfLl-na,J) »7(«-i> 21 I.i
IF
-------
Htti,i?>!"tt;,ifi»n:TP 9J "J
itc if 51 '•O TO t'O "J J&9
I.F I ill i, NC I < L T. i ."C ? 1. 2 > T.O TO tf»(J *JJ •""
— HJ 'Qfl
HTlii.N.^lyHTiIlffilill^SAXi *•*«
!%0 CC'ITIhuC ^J *lfl
c «J 1,35
RETURN 8J WC
e 3J «,",5
EMS 8J H50
SUBROUTINE OOFLOD
5U3 FO1? OVE?0 > ZIN'i,MCl| CtN'i,^;), SC:i*»UC), OXCNS OF 100
^ I,NCI, oru, ^uu.R,.iC) op ^^
COMMON /8LK1/ OT,ST,ICFiO,IUKAQ,LOt;.UCIH.UCje,l.CJ«,F1'ITCP OF 1JO
'COfHON /3LKa/ StA,STCu,^TT,SnS,S:lT,3a^i,SO-3fi,sqiCF.TA*C.C • *" 1"
Stf01.0*0.a . •' - 2r S?S
SVCLF*0.0 OF ZiC
KCT'C . OF 2ZC
00 HC J*l,MC ' OF 2JO
00 18C 1*1, Hi* _ - OF 2*0
IF CK;,J» .GE.HCCC.C) CO TO ISO Of 250
IF mtt.ji.<;e.z(i.j)> oa 10 in ' <>F ?*«
OOFtsZCl,JI-l'T tl,j> OF 27C
rr.Tsa.i* ' OF 2^0
SS=SU,J) OF JOC
IF IICF40.LE.GI GO TO Hi) OF Jll
SS'SCII.JI OF 320
IF tt-TU. J» .LS.GI1.J) ) SSsSa.JI OF JJ3
. it «H(i, n .r,c.2ico;.5) GO TO 1*3 or jsc
SVCliOfUCFT . Or 36C
1?C MTU .Jl =ZtI,Jl • ' 'OF 37C
f.O TO li'i ' OF JSC
uj JP i:ceiD.r,T.?) co TO 133 OF 19;
I* IHTCI.J) .LE.MI. J» I GJ TO 140 0" «iCO
OOrTsC.C CF <»JC
svoco-a.J CF m,c
IF (fi; ,j).r,T.2occ;.o oj TO luo . . OF «,so
SVOLF-^iCFI OF 1,63
i;o TO i?i • ' OF i»7ii
!<.; HI < i,j»-;, (i,,)) OF <»3t'
,3VCtf*0.1 ' OF VJO
IfO OACf t f-Q^CP TT I'JVOLO ' OF. 5tlU
FlCF1tsFACrTT»'Vll.r OF !10
i^ ij.i.r.ici»i ».'J TIJ iii or 52C
IF U.M.lClr:) OCT) X(0 OF 5,«fl
if ii.it.LCxn r-c TO i*.i OF $*u
IF ll.r'l .LCJi I C-9 '0 1«,3 • OF 551
" " ' .or 5*.;
348 :
-------
Hi If IKCT.r.T.OI CO TO 171
H9ITC I6.2CO)
t>3 :i«ir- I*., 2^3) I, J.OOKT, J-UCFT.FrTtFiCFT
.-ra.r.j.it.o.riCKTr. ro.s.ii r,o TO tio
N*IT<- «6,2i:> O'CFTTiFACf TTtQSO TAtFACFTA
t'••
Of *•"
OF •>?'!
3p **C
OF «i»j
0* 45Jj
M,2X,.11HF LOO OF 7iC
OF 71C
f .(iur:Frx» OF 72;
OF 75C
SUBROUTINE BYFLOVV
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
T"IS
r.»IO«.
THOUGH
; ,CC,CI.°'t;iPHIC.CS>
r FIO/; .JJl'Ur/HJScJ.
60 T3 1/3
9F
3^
IF
3F
BF
dF
IF
3F
3F
BF
OF
OF
OF
9F
OF
BF
JF
tr
IF
3F
3F
OF
8F
BF
PF
-JF
-»t-'
9F
3F
iiF
3F
HF
f.F
-IF
rf?
BF
JF
f»F
OF
1M
60 TO 171
SQOGKI.J): ].
. J)«(HF(I,J) -
•OT/'»J'!60.
30
to
5u
«j
73
s:
9d
100
IK
120
1JC
160
19C
20C
i'16
2Z5
23'.
2n3
£50
280
201
3JO
11)
aea
3S3
3<*i
J51
J6C
J7C
J^iC
H:
«.0d
<»r>'i
i,6i
Vlj
500
511
szi
53;
5<<0
55C
563
I7C
59;
OF 6JC
349
-------
i /:
if ii.'i-. i j GU re isu
If IKK. En. >.l Go 1C 105
acr.i «t. j)
SUM OF t-FLOU
CO TO 193
1AJ tf (I.tit.N'll f.O TO
IF i«.f.n,<>) GO TO
SU»1 OF I-F
TIHUUGM
c
c
c
IES
IF It.Nr.N!) GO FC Jt C
IF I J.LT.!.Cr-l» GC T'J 210
IF I J.r.T.v.CIr:) GC TO ti'1 .
IF IKK.EU.I.I r,c ic «?ie
T H»i--s^r.:i(i,ji
fu»! CF i-FLf!H rue JUGN 9UFF-t», ZONE 9 C'JM
GO TO 210
285 IF CI.'I£.NSI GO TC 210
IF IJ.UT.LCISI GO 10 dl3
IF .IJ.GT.LCIr.1 GO TO ZU
IF KV.E;.-.) GD TO 210
c
c
c
SOH OF I-^LOH TMf«OUC-H SUFFER ZO'lc HOIIHOA^IES
i ir;n£
<»!R,f.:,SOr,GI)
211) CONTINUE
CALL r-
0-3 «.3 LCCP CO'"UT£S FLOUS n THE J-OI«tCCTION.
DO 320 l:l.Nf>
0* 32C J=l.NCl
CO TO
CO TO
CO TO
GO TO
: »/icccc.» »i
tt »7!»,r3C,3i.C» t
GO TO 360
soco m ».i>:PAr«AM{Fici; ,j*:i ,F
(I, JI.OXII, j»n ,CX(t, j» ,jv (I,
23P
2 ".0
250
GO TC 2H!
265 Srj&GJ « (HF Ui J) -HF 11, J»l» »[>?/'« 356B ,
CO TO 2d3
270 SOG1J(I,J)=0.0
, 2*3 IF ( J.'lE.ll GO TO iO
IF (kK.Fu.'4| GO TO 330
c
c
c
c
c
c
SU!* OF J-
THfCUCH 'JCUH34-'. If5
CO TO J01
IF IJ.:J£.:JI:D GO TO 3CJ
IF
SUM OF J-rLOW PV.CUC.h uOLimAKIES
i^ u. iF..i.'.t »'.j fc nc
IF u.Lf.tcj.n co n 3^•<
IF U.r,T.L.:jC) GC T>) J3:
IF (KK.ru.iti no TO J20
S«H OF
i* TMft'JUOH
?0'/li
C<3 TC J2i
310 IF |J,'I=:.L?.I GO TO JJO
•XF II.LT.I.CJ4I CO TO 3ZJ
•OF 6'.1
IF 60,
8F b/C
OF 610
OF i>9i
OF ?CG
OF 710
OF 72:
0
rtF 160
BF »70
BF 44]
SF 8)0
OF 930
BF 32-3
BF -9JO
BF 950
»F 970
aF Mi
ec uic
R-% • ive
BF 1010
OF ijti£
3F 135C
BF 1370
BF 108d
BF 159C
3F 11JC
OF i:ic
9F 1120
'IF 1130
dF ll-JC
OF H6i
BF 1170
3F
8F
BF 1Z3Q
BF
3F
UF
8F
OF 1250.
3F 1240
•JF tZ75
OF
HF
•IF
IF uii
•IF I U5
JF lli»1
•IF 1 !f,0
OF IJ7C
HF 1JJC
IF 1 IK
OF f.00
OF J!.JJ
9F
BF
BF
350
-------
.- »....-• !•• i '•<• • •> •'-' ,jf.
IF (<".. rl'l.'. I 'iO TC -J^i • jj.
SCO !F(
CT I
o:c IF isoGijc .j-t).ot. i-.i-i ••-• >" ••"«
cTfi,jj=crii,j)*-i-:<-.ii,ji|-'-j:i.J-i»»»sQOGjii,J-i
CHIiJt=i>"J" * J'IG'»J(I »J>
Cfl TO 560
55C CTCI.J»=CSIIi
56C COHTIliUE
C VJM Of J-FLO-J T.IP3ur.H WF» /ONt 10'n'jA»:< S jj J«»
C ( . luuu or 13'!l'^»u***i*** •* •»» -• _ . 7i»t
50PL*£n«U,J) QF 175i
C BF i760
C TOTAL UFLOvl F3CH CCNSTiNT rll&O GftlOS 9f t/,.
C . ...... 5F 1780
J55 CTlTII.Ut gf
C ,^r ' ' "F
C TOT4L IfFUOH Ti^Cof.H &OUKO-44IE5 . HF
SOBT=5Q3J4»3r.£lIA OF IH,,J
;r isessT.i.!:.?.oi GO it J6l aF l(,5l,
V'.IT£ i£,«.<>:i ITIME •!,- 154tl
CALL HAT.'iOP IIRiNC ,SC«<) (|r il7;
C . 6F
C , -OF
309 Stl?JiC«0 . . • $f
00 37C J-1OIH >LCIE , sip
saoGic.'i.-n -j.c • OF iiit
S.l&CI I'iS. .1)-C . 0 (if 193:;
37? CONTINUE 'OF U^l
M JrtC lsLC,J'!i(.CJ£ ' fjf iw.
SOC.fi JII .L'll =C .u gp ; j^(,
* ' • _ii . . • l*r I s» j
00 310 i-LC.IKiLCJE . t,,: ^^
UO J9C .J--LCiM,LCIF . or ^j.ji,
sor.ii.Ji-o.o . nr' iii:
351
-------
C TOTAL ICFLw* F-C.I'. CCl.STillT f--IUi WITHIN 'OFF-:"? ZOI.f. PCUr.
C
MC CGtlTII.UE
C TOTAL IfcFLCil THOUGH tVJFFrR ZCTIE 30UHOA3IES
DO .ru j-i.-ir.
1.3 :n/ :-I,NK
codtJi -cr d,ji
397 COHIIMH
fciu
«,»3
IAOUIFER
(!'(;, KX,
OISCMi'ilE IN I-OlKtCTICN eA>FT/U7>
ISCHSP'ie IN J-'J t-'SCTION CA>FT/OT>
HF
JF Z11C
V 212J
3F 2110
OF ^ll.5
6F i'i-;.J
nr
1:1 ;i':i« Gem MI-US M-:tf.» CLC.< F?.OM .o
tK-i-'i/nUi/ti" .'.?<•
tMO
3F
THIS SUSfOUTI.-Jr: W.MTti OUT Tin: OAL4llCi COMPUTATIONS FCR lACH
TIME inco.iiiE.'iT. ALL UNIT; ATE ;;l AC-FT PSR TIM*
SUBROUTINE BALCOP
SUBROUTINE 3ALCOP J Jl, J2 11,3TTTE*,STATErt»
C
C
c
e
c
c
c
c
c
c
c
c
c
c
/3LK2/ STA,3TOL,STT,;>:ia,S'JTtS-}f;A,S1BA,5«HT.S1iJT,OVA,OVT
J1.J2.I
3'JA.SCT . • "
SORi.SCRT
SOA.SOT - Jpr»LltO WAT^P, 0£TH£t».
SU»*iSO>T - INFLOW F
-------
SUBROUTINE MATROP '
SUHKJUI '.s.e. MAK'CH (M^iJH.NOCJLi'll *f 1C
C ** ft
C MF »0
c rnii suvi'wTpic G-t,.:M.!-s turA o* ^u'utr; TMO Q suiti-iLE FORM
-------
IBM
APPENDIX A-VIII. DATA CODING FORMS
1/4 FORTRAN Coding Form
GX2B-7327-6 U/M 050*
Printed in U.SA.
*.-OG.«
"C='"™'t' n. K. Stmada
i
s.
;-•
l; Vil»
N R
I 5
D
F 1 0
D
•
D
f
,
Z
G
P H
C
P H
Z
L
,
X
y
K
z
G
I
A
I
N C
I 5
X
1
•
C
N W
I 5
D
F 10
0
0
• F
I
G
P H
; C
; PH
N P
I 5
L Y ' '• ••
. 1
X
l
y '
• ; •
K ' '
Z
G
' ' l
I ' . :
A
1
1C
14 17 It n »
1C FAQ
I 5
1 - ' F
F 1 0
-/-*• (
• •
-f-*- (
-+-* (
-f-* i
• ,
-f-» - (
j
-f-* ' (
I .. (
71 a a it a
Mn Junp 1976
SSSS-.
G*49HIC
PUNCH
FORTRAN STATEMENT
11 KA
I
Q
5
FK : .
. 1
MR* N
C
NR* N
C
1
NR* N
NR* M
C
C
NR* N
C
| ;
NR* N
N R * M
c
C
H R * N C
» n a v x
; Til
1 8 A 1
i ,
: D
. i F 1
, . •• '
: : i Z
: :F10
* t 1
)' *'( Fl
1 ' 1
);* ('Fl
! i ' '
) * ( Fl
\
)''* ( Fl
• ;
)'* ('Fl
' j :
) * t'Fl
) * ( F 1
) !* ( F 1
L E '
0 '
; I
T
0. 2
Z ! : 1
.1 '
; : !
o. i );
• j
0.1)-
; i
0. 1 );
. i
0.1 ) i
; i
0. l )
.' i i !
0. 1 )
0. 1 )
0.1)
,
:
F
s
1
• i
, i i IG
' F 1
0
! ; '.
' ' !
' ! i
| I
; j j
READ
READ'
RE A D
READ
RE A D
i
R E'A'D
I '
1 - '
; i •
T ': ' ':
0 .2 '
= ! l !
G1' : ' i
. 1'
• i . ;
i '
• ( 1 !
i ] ! ;
I N: !BY
1 ; ;
IN' BY
. ! i
IN BY
; . !
i N; 'BY
IN :B Y
: • j '
IN! !B'Y
'. ! ;
! i
|
!
•
i
t
1 F
1 F
u
i
i i • :
i
! P
Fl
p
0
1
; i ; •
!
1
1
i .
! •
i
C;OL U
: ! . •
COL
n
i i •
C OL
D
' i ' ;
',
;
r;ni
i
C OL
i
. i
n
'u
i
51 a S3 S4 55
i ?
: : !
T 0
0
H1
T
p
2
'. \
: \
. 1
1 • ^N
i i :_
• •' ' /
M
H
N
N
S
s
I
\
j ' j
M
N
S
: • 'T
M
M
M
56
N
N
n
17
^
<;
\
\
J,
. ' ••
• ' ;
i '
• i
: 1
i 1
i i
F 1
t
1
C
1 o
: i i
1 i :
; I
0 HI
I N;
1
1
T
r
•: i ;
! 1 ;
; l i •
t i
.
i
—
i ;
a a i
•4 49
MCI O>
CA1O IUCI1O nvvjii-
' ; "
: ; ; f
c
! i
A '
! :
t
1
. ' !
' >
; ,
: ' . i ! !
1 • ' • !
A
'IF
'R n
5 P F
3
1 .
i . ,
'
!
41
4> M M W 71 71
(OlNflflCATlON
P P H I C
F 1 0. 1
l
1 •• i ; ; ,
p I F IF n
i
n ?« TJ '4 ?* ?i T* ac
-------
GX28-7327-8U/M050*
IBM 2'4
33
33A
4
4A
44
44A
44B
44C
5
FORTRAN Codlnt Turn
«-:•:•«
tcfwvFt D. K. Sunada
.™«r
•34:
c
F 1
c
F l
H
H
T
F
P
F 1
?
<;
D
n
M
H
n
i
i
K
P
n
i > K
1
1
L
1
L
T
. 1
{ UN I
C
l
F 1
T
F 1
H
T
F
C
F 1
11 K .6 '» X
r 0 R M
S
B C
0 . 1
A
TL
0.1
L
L
K L
P T
0 . 1
A L U E
-/-* (
U
R
F 1
i
/v (
1 U
F F
F 1
/ * (
• '
:'•<-' (
• . • :
I .. (
i I F N
• • 1
0 F
wn June 1976
PUNCHING
NS1RUC1ION*
CU4WC
fUHCM
FORTRAN STATEMENT
R
NR * N
SE C
C
0
BC '
0. 1
NR * N
C
SE C 0
KL
0 . 1
NR * N C
i
NR * N C
: i
N R * N C
ONE,
EL AT I
: . ; .
) * ( Fl
DE Q V
1 ;
! ; i T
i : .
r* ( FI
D E D V
! 1
{0 M IT
) * ('Fl
; • j
) * ( Fl
) * ( Fl
BLANK
,
i ; i •
VE CO
0 . 1)
A L U E S
B C: ' '
0.1
i
0 .:1 )
A L U E S
•
U N L E
0.11
1 • i
0.1)
•'•'.'
0 .1 )
CARD
M C E N T
READ
1 , !
: ' ; B
' F 1
• • : '
R EA D!
i ! i
S S L E
: , i
R EA D;
. i
R EA D
' ' i i
R EA 0
: '' . ,'
! i l 1
AT 10
N
BY
i '• '
B C ' :
0 .
1
j .
IN
B Y
\ i !
A K
Y A
i i
I N
B Y
1
IN
I N
B Y
BY
i ,
: i i
:
i
j i
'
0
T
C 0 L
II
i
• .
!
C 0 L U
. * i
Q U I F E
,
CO
I
. ' i
r. o
1
r n
1
i t ^
!
'• . i
j i
H_
E R W I
M
N
1 i i
1
i j |
! ;
•1 N; i
; t : :
R
CON
•: '•
: '
M.
N
i i :
M
N :
M
N; -)
i
: i
:
SE B L
0 M I'
I H
T
CA
t
1
0 M I
I N' ;
T
C A
' ' ! '
D I T
I 0
•
0 M I
I N
T
r A
• •
i
A N K
Printed in U.S A.
C A.R
F
RD1
SPEC
3 3
!
I F
R D
SPEC
4
1
N S
D E S I
I F
RD
SPEC
4 4
'°— -°"
n
I FI E D
IF I E D
RE D
IF I E D
-------
IBM
3/4
FORTRAN bdtaf Fm
GX2S-73274 U/MOSO"
Printed in U.S.A.
5A
6
6A
7
7A
7AA
78
« JSMM
«oc.«w« D. K. Sunada |O»K June 1976
nJNCHIMG
*
",r;r
YP T
A
F 1
YAW
HGP U
I 5
•
J
I S
D
F 1
YPR
|
P
0
ft
0
U
. 1
YAM
K
I S
\
t
H GPU
R
- 1
YPR
C
F 1
YAW
P
F 1
P
Fl
YPR
:
i
_^
AW
0. 1
-y-».
{ZERO
HR 1
0.1
' :' i
i
i
1
- :.
HR
o. r -
._/»•
. - . :
I N YR *
{IF N
I N YR *
1 •
IF :N
•
FOR
, j
I , (
: 1 ' ;
• i ' !
I N YR *
f t
(F5 . 1
ON E .
(F 5.1
OM E -.'
- . ! ;
, , •; .
' t
• ' i
: ' ; i
\\ G P U
I ' -
• ,
' ! ! '
NR * N C
! 1 1
' , ! 1
{ F5 . 1
', '. • '.
• 1 i
j ;
fORTRAN STATEMENT
{
i
, i
B L A
NK
! !
i ! •
f
: • : '
BL:AN K
i ' ! '
i ' '
|
•
I'll
i • ! 1
. GTVO
i '
1
j ;
)'*'{
f ,
1 1
i ;
V •
;•
. . .
Fl
•
;
f
: i !
i j
!
0 MIT
CAR
0
OH IT
i i
IN D
I
i : '•
: • •• •
i : '•.
• ' 1
; 1 i
. o i.
; . J
' 1
0 Jl )
' ! :
1 * !
ONLY
i '.ij>
'' '• \'
""<
IF P P
t
I F: A P
C'AT E s
: ; i :
i • ! l
: i ' !
• ! ; "I
t i ' i
• i ' 1
O-TH'E'R
i i :
i , i
. ! : I
READ
• i i !
,iii
N E ED E
;. i i !
-\ i i i
OWMC
IUNCH
T .;L'E'.
• ; • :
1 i i '
' •
. : .
w
. • L E .
: i !
E N T'l
i : ' '
• : i
! ! !
1 t ! i
• i ! i
w iis!E'
l : J !
, j !
i i :
I
N BiY
; i ;
: ! ! i
D JI'F)
i i i ;•
Mh
• ! ! 1
0. 0
i ' '.
1 i ; !
' ' i i
0.
0
• ; : ;
RE: P
H
• • ' '
• ! i !
1 j i
.
I
1 — ' —
3 MIT
1 t
: : ,
! ! i
^CO!LU
T
i
1 I
< \
P H R:.
r,
i : I'i
i ; 1
!
' 1
; « i
j
t
i !!
i ' i
; ' i
R'
M
A
T
1 * i •
i
i
1 •
i
i i
i !
i 1 i
i , ,
\
M'N
i
i
j I
l
i i :
T'
i
0
n
,
i ' ! •
i
I
: ' i
!
R
J
T X'
1
•
— ! —
T
I
; i ;
i 1
i j
',
!
1
I
' ; i '
i : ! '
; 1
F
OR1
N
1
f
' : i
i '! . 1
' ! ' "
i :
'
MCI Of
CAC3 UtCIftO KUMKf
'' l '. '
->,-.'•.
i i ; ; ' •
I.I
0
BE RE
J '• :
1 ' '. ' •
: . ; I '
i
i i
i I ':
i i J ! ] ,
' i , i
i ; i • ' -
1
GP U '
-BL A
.
: - . :
•
IO£N1JFiCAflON
1
AD IN;''
• .
l;|.
. • • : i . j
-' i . ' • ,
. •
. '..-',.
NK
<*>
U1
0>
-------
TRiyr 4/4 FORTRAN Coding Turn
»CO'«M
«CG'«VMH D. K. sunada h" June 1976 1
fUNOONG
msnucnoM
CHAfWC
fUNCM
OX28-7327-«U/M050"
Printed in US A
CJUQ ELECTRO Mj***t<*
?
r.i.*. •
).''{!
• ' i
T
i . 1 :
i ' ;
1 ;
i '
R'E AjT E
' :
l
r , '
i ' '
1 ;
1
) ! Ul
i :
T
: , [
C A
! i
R n
! :
F BL A
R U B U T
. ! : i
;
I ,
R 'T HA
• : i '
' ' r (
- B LA
? I BUT
i : '
• 1 i
•
S 5 1
. ; . ;
1 : .
i
.
N K AS
ION A
! i
: ; •
i • i
i :
n; ;z E R
' • i ;
! '
. . I
N K A S
I 0,N A
'i:l
• i '
H R 0 UG
! i i |
; . '
SUM E'S
S P RE
:
i
1 ' i
O'.i ;OT
; ; i ;
1 i i i
. ; . i
1 ; 1
S U MEi
S PR E
' : i !
t "•
H' :7!77
: ' 1 i
: i • .
: •
< '•• . '.
• ;
!
S AM E
V I 0
US
: 1
1
HE RW.I
; ! ; '•
: ', • :
, '
i ' . '
S A
VI 0
H E
US
! ' '
! ! i •
• •
B F
1
O'R
!
•
• : j .
; . • ;
i •
i ,
D I
' LI C
1 Fl t
i
Si-
L L
i
1 i i
1 ! •
! ! ,
S E
0 M
•• 1 ' ;
1 ! : :
i 1
I :
! i
D I S!-
p:iT'
i ! i t
t :
i 1 .
N EX T
1 ;
i :
l '
I
1 ' !
, N
N
C
1 i
: i ! i
i ;T ! i
' ^ ' ;
1 : ! i
•
•
i
N
t
C
i
YEAR
j ! '. •
1 •
i ! : i i •
OT NEED
W * 0 IN
A R D
' • •
! ,
•;;•!!
i i : i ,
: • ' . -
i • ! !
!
0 T NEED
W = 0 IN
A R D'
OF ANA
WCUiNCt
1
E D IF
CONTROL
• •
' • • ' i '
! i : '
• '
ED IF
CO NT R 0 L
LYSIS
77
77A
77B
777
777
7771
777
-------