468
A.S. Anderson and R.C. Hanson
880
NN=NN+1
890
78 CONTINUE
900
S=H-SUM
910
WRITE(2) TIME,S,(CON(L),L=2,Kl ,3)
920
IF(.996*SSAVE.GT.S) GO TO 67
930
IF(.998*SSAVE.LT.S.AND J.GT.6) GO TO 400
940
DT=DT+DT1
950
67 SSAVE=S
960
100 CONTINUE
970
400 CONTINUE
980
NC=I-1
990
WRITE(6,600) NC
1000
600 FORMATCNO OF ROWS =' ,18)
1010
END
APPENDIX B
100
REAL K, VTR(501),VEL(501),VN(501),V(501),F(501),P(501),CON(501), %
110
U(501),FOLD(501),VR(501),PS(501)
120
DIMENSION 11(501)
130
REALLTH
140
QIN=750.
150
QOUT=200.
160
DR=5610.
170
DT=80.
180
SDT=8.
190
M=501
200
VI=.00055
210
VIN=QIN*231./(1728.*1440.*60)
220
LTH=VIN*DT*12.
230
CO=QIN*DR/QOUT
240
PBW=2.41
250
R=3.
260
SG=1.41
270
WW=DR*(SG-1 .)/(SG* .061*453.6*1 .E06)
280
DEPTH=DR*LTH/CO
290
PI=WW*LTH
300
FI=1 .-DR/SG*1 .E06)-DR*PBW*1 .E-06
310
AK=PI*FI**3/(VI*(1.-FI)**2)
320
FU=.90761
330
VLT=LTH/12.
340
DO300L=l,M
350
I1(L)=4*L
360
VTR(L)=0.
370
F(L)=FI
380
FOLD(L)=FI
390
V(L)=VLT
400
CON(L)=DR
410
U(L)=1.