-------
(FLECS VERSION 22, 46) 13-MAR.81 13l20l2U PAGE OOOOJ
1 30,00, 1.E2/
DATA CF /9.00, l.E«, 5,E«, 1.E5, 1.5E5/
DATA P X0.60, 0.90, 1,0, 1.0, 0.83, 0,6o, 0,«U, 0,25, 0,15,
1 0.09, 0,05/
OATA OP /0.10, 0,15, 0,20, 0,30, 0,40, 0,50, 0,60, 0,70,
1 0,80, 0,90, t.OO/
DATA DC /0.10, 1,00. 10,0, 100, /
DATA Vfin.O, 1,5, 2,0, J.O, 0,0, 6,0, 8.0, 10. /
OATA 050G/0.10, 0.20, 0,10, 0,40, 0,60, O.»0/
DATA TEMP/32,, «0,, 50,, 70,, 8».f 90,, 100, /
OB50 * l>50 * 1000,0
FHRAO a HHAD * 3,2BOfl3J
TMPH m tEMPR * 1,8 * 32,0
V 3 (UTOT * ALEN/VOL) * 3.7975E-5
H s VOL/(HRAI)**UEN) * 5,280833
»«* FSL....FINE SEDIMENT (i.e. COHESIVE SEDIMENT OR WASH) LOAD
IN MICRO GRAMS/LITER **»
FSL * 0.0
00 (iXsl.NELEM)
, FSL * FSL » o,s*(C(ix,2) + cdxti.2) » C(ix,3j » c,dx + i,3))
...FIN
FSL » FSL/NELEM * 1000,0
IF((OB50 ,LT. 050G(1)) ,OR, (0*50 .GT, D50G(6)))
. FEKRUR a .TRUE,
, WRITt(b,l)
, FOHMAT 1//10X, < ***** FATAL EMROH «> SU8HOUTINE COLBY *****•)
, NRlTt(b,2)
, FORMA! l 10X, I ***** DOUT *****')
...FIN
IF((FH«AO ,LT, UG(I)
. FEHRUR a ,T«*JE.
, WRlTt(6,l)
, HRITt(b,3)
, FDRMAH ) ox, '»»***
...FIN
IF((V ,UT, VG(D) .OR, (V ,GT,
. FERRUH » .TRUE.
, MRITb(b,l)
. HHJTEtb.t) V
, FORMATl 10X, ******
...FIN
001 10
OOlt 1
001 12
001 13
001JH
00115
00116
00117
001 IB
001 |9
00120
00121
00122
00123
00120
00125
00116
0012T
00129
00129
OOljO
00131
00132
00133
00130
00135
001 36
00137
0013S
00139
00140
00141
00142
00)43
00140
00145
00146
001 47
001 uft
00149
00150
00151
00153
00)53
00155
00156
00157
00)58
00159
00160
00161
00162
00163
00164
00165
C
C
C
C
C
C
C
C
C
C
1
2
3
4
, (FHRAO ,GT,
HUUT
»***«')
*****', F I 0 ,5 )
IFCTMPK .LT, 32.0 .OR. TMPR ,GT. 100.0)
, TMPH s 32,0
1D1 a U
102 « 0
00 (1*1,3)
72
-------
(FUCS VERSION 22.«6)
13-MAR-at I3l20|24 PACE 00004
00166
00167
00168
00169
00170
00171
00172 114
00173
00174
OOJ75
00176
00177
00178
00179
00180
00)81
00182 1(8
001fl3
00184
00185
00166
00187
00188
00189
00190
00191
00192 122
00193
00194
00195
00196
00197
OOJ9*
00199
00200
00201
00202
00203
00204
00205
00206
00207
0020H
00209 123
00210
00211
00212 12"
0051J
00214 126
00215
00216 1
00217 2
00219
00219 127
00220
00221 128
IF <(FHR*0 ,Gt. 06(1)) ,*NO. JJ(J)
. . mj) « ALOGtO(VG(JI))
. . 00 (Kzl.2)
. . . Kl » KK(K)
. , , ZZ(M « »LOr.lO(OSUG(KD)
. . , If (G(I1,J1,K1)-0.) 123,123,127
. , . CONTINUE
... DO (J)*J1,7)
.... CONTINUE
. . . CONTINUE
. . . X(J,K) , *L()G10(G(I1.J3,K1))UALOG10(VG(J1 )/VG(J3)))»
• • • (ALOG10(G(Il,J3f 1,K1 )/G(llf J3,K) ) ) ) / ( ALOG1 0 ( VG ( J3* I) /
... VG(J3)))
. . . GO TO 12*
. . . CONTINUE
. . . X(J,K) » »LOGIO(G(H,JI,K1))
. . . CONTlwije
73
-------
(FLECS VERSION 22,«6)
13-MAH-81 13l20|2« PAGE 00005
00?22
0022J
0022«
00225
00226
00227
00226
002*9
00?30
00231
002J2
00231
00235
00236
00237
00258
00239
00210
002U1
002'J2
002«3
002«5
002«6
0021?
0 0 ? '> ft
002«9
00250
00251
00252
0025J
00250
00255
00256
00257
0025?
00259
00260
002bl
00262
00263
00?fa7
002b«
00269
00270
00?71
00272
00273
0027«
00275
00276
00277
C
C
C
C
C
1
C
C
C
. ...FIN
...FIN
XD >
XM
XM2
ALOG10(OB50) - ZZ(l)
= X(l|2) -
' X(2,2) -
x ZZ(2) -
XA(1) a x(l,l )
a X(2. t)
XA(2) -
XU, 1)
X(2,l)
ZZ(U
t XN1*XO/XOEN
* XN2«XO/XOEN
XA( 1)
WHEN (TMPR
ELSE
•
t
•
•
t
f
TM
IT2
DO
.
.
,
a
a
u*
IF
.
,
0
0
1
(
1
I
.EQ.
,6)
(TMPR
Tl a I
T2 a J
60.)
,GE.
»)
XM2)
XNH 3
XOY a YY(2) • YYU )
XG(1) a XA(l) » XNM*XV/XOY
...FIN
XNM a XG(2) - XG(1)
XO s HOGIO(FHHAD) • XX(1)
XOEN » XX(2) • XX(l)
GTUC » XGU) » *NM*XO/XOEN
GTUC s 10,**GTUC
*** ttTUC IS UNCORRECTED GT IN L8/SEC/FT **»
*** '-EXT APPLY FIME SEDIMENT LQAQ ANg TEMPERATURE COnRECTlONg •**
CFT « 1,0
TEMP(I)) .AND, (TMPR ,LE, TEMP(Itl)))
. . . GO 10 IJb
FIN
. ...FIN
. CONTINUE
. xTii.n * ALOGIOUUTI»IPI »
, xr(2,l) a AunGlO(TClT2,IDl))
. xr(l,2J * ALOC,10fT(ITl,IL>2))
'. XNT a ALOG10(TMPH/TEMP(iTl))/ALOGlO(TEMPcIT2)/T£MP(lTl))
. XCT(1) a XT(1,I) t XNT*(XT(2,1) - XT(1,1»
. XC!<2) a XT(I,2) * KNT*(XT(2,2) - XT(I,2))
. CFT a XCUt) + (XCTC2) - XC T ( t ) ) "XU/XOEN
. CFT a 10,**CFT
...FIN
*«* FINt SEDIMENT L0*0 COHHF.CTION ***
WHEN (FSL ,Lt. 10.) CFFal.o
ELSt
. 101 a 0
. IU2 s o
, DU (131,9)
. . IF((FHK»0 ,GE. DF(I>) .AND, FHRAO ,LE, OF(IfD)
. . . 101 * I
. . . 102 r U|
. GO TO 102
74
-------
(FLECS VERSION 22.46) 1J.MAH-8I 13«aOl2'4 PAGE OQ006
00278
00279
00280 M?
00?8)
00282
00283 5
00?8«
00285
00286
00287
00288
00289
00290
00291
00292
00293
0029«
00295
00?9b
00297 M8
00298
002V9
00300
00301
0030?
00303
0030«
O'SOS
C. jOb
00307
00308
00309
00310
0031 1
00312
00313
0 0 3 1 t
00315
00316
00317
00318
00319
00320
00321
00322 153
00323
0032"
00325
00326
00327
00328
00329
00330
00331
00332 C
00333 C
fl
•
t
•
•
t
t
t
t
•
t
•
t
t
f
•
•
t
t
*
t
f
t
V
t
•
•
«
t
»
t
•
•
•
t
t
•
•
f
•
t
t
•
•
•
•
•
•
f
f
t
•
•
. . ...FIN
. ...flN
, CONTINUE
. WHEN (FSL ,GT, t.E»5)
, , 1RI TE (t>, 5)
, FORMAT(//IOX, »***•* SUBROUTINE COLBY »• FSL WENT > l.Et5«)
. . in m 9
. . IF2 * 5
, , , ,f IN
, ELSE
. . in s o
. . 1F2 » 0
, . oo d»i ,2 ) )
. XF(2,1) a A|_OG10(F (102, 101 ) )
. XNT t (FSL - CF(IFI))/(CF
-------
(FUCS VERSION £?.'J6) ij.MAK-fll t3l20i?« PARE 00007
003J5 ...FIN
OOlib RE.TUHN
003J7 ENO
(FLECS VEUS10N
76
-------
(FLECS VEHSION 2?,U
Iiil9i2b PAGE ooool
00001
00002
00003
0000<(
00005
00006
00007
00008
00009
00010
00011
0001?
0001 3
DOOM
00015
00016
0001 7
00018
000 1 9
00020
00021
00022
00023
0 0 0 2 «
00025
00027
0002P
00029
00030
00031
00032
00031
00031
00035
00036
00037
0003S
00039
00011
000«2
00043
OOOU5
00016
000«7
OOOUB
00050
00051
00052
00053
SUBROUTINE COLLAP(ALfN,AREA,C,OELZ»El»lELM,K,NELEM,
1 CELCr*,TMASS,AwlD,VSET,OFZ,CNOUE,XSAREA)
c INTERPRETS INITIAL CONDITIONS
C CROSS-SECTION
C SEDIMENT CONCENTRATIONS KQ/M3
C PARTICIPATE CONCfNTRA 1 IONS PC/KG
C OJSbOLVED CONCFNTHATIONS PC/M3
C ARRIVES At BOTH
C NODAL CONCENTMATIONS
C ELFMENf AVERAUE VALUES
C NOTE (I,-U£LZ**S/EZ) MUST BE » 0
C
REAL MSAfl,MSH8,MSAT,MSBT
C
INCLUDE 'SYJELMSIZ.PRMI
c
DIMENSION AH£A(MX£LEM),C(MXELEM,MAXCON),IELM(MXtLEM),
EL(MXELEM),CELCTH(MX£LEM,MAXCON),TMASS(MAXCON),
A«IO(MXELEM), VS£T(3),DFZ('»)f CN()DE(MXELEM,MAXCnN),
XSAREA(MXELEM)
BLENst,/ALEN
ELTOPsnfcLZ
wisAREAtI)*BLEN
Wj*AMtAt2)*f5L£N
D£L6VeFL(2)-ELU)
C1«C(1|h)
CJ»C(2,«)
NELMHTsJ
MSAHsDELEV«(CJ*«J/3,
MSSBsO.
IFfK .Nh, 7)
GIsCO,K»;
PM6AB«0£L£V»C,25*(WJ*CJ*GJ»WI*cI*GI.
"I*CI
OU( ;
ELTP=£L(NELMTP*l
ELBTsEL(N£LMTP)
CI«C (NELMTP,K)
FACls(eLTOP-ElBT)/(ELTP-EL8T)
')/(ELTP«tL8T)
77
-------
VERSION
13-MAK«81
I3ii9i26 PAGE
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
OOOfaU
00065
00066
00067
00068
000&9
00070
00071
00072
0007S
0007«
00075
00076
00077
00078
00079
00080
00081
00082
0008J
0008U
00085
00086
00087
00088
00089
00090
0 0 0 9 1
00092
0009}
0009U
0009S
00096
00097
0009H
00099
00100 C
00101 C
00102 C
001 0 J
0010"
00105 C
00106 C
00107 C
00108
00109
, CTUPsCJ*F»Cl+CI*FAC2
, NTOP*'*J»FACItWi*FAC2
. MSAT=(ELTP-EUOP)»CCJ*«J/3.+CJ*WTOP/6.»CTOP»wj/6.
1, »CTOP»vUOP/3.)
, MSBT = !/
. . . 1F(K ,Nt, 7)
. . . . GI=C(J,K»3)
.... GJ = C ( J+ 1 ,K + 3)
. PMASSsPMASS»OELEV»(.25*(Wt *CI*GUWJ*CJ*G
1. . . . CON*(«i*Cj*GJ»'<'J»Cl*GJ + ''J*CJ»6l-»
2, ... wj»ci*GI»rfI »CJ*GI + rfl*CI*GJ) )
FIN
PIN
. ...FIN
, 1MASSrH)rTMASS(K)»CM4SS*ALE'J
, IF(K .NE, 7) TMAS5 ChLCTHf I,KM>=PMASS/X3AP£A( I )
•
, OEVFLOP NODAL VALUES OF CONCENTH AT I ON
f
. HHtMU .Ef). 1)
. . wHENfK ,LE. 3)
78
-------
(FLECS VERSION 22.«6) i3.MAP.fll 13»)9j26 PAGE OOOU3
00110 C ...
00111 C ... NOFEi CMASS IN KG/M
00112 C ...
00115 , . . COEFsO.
00114 , . , «S*VSET(K)*AR£AU)/(AWID(1)*AL£N)
00115 , . . EZsOFZ(K)
00116 , . . CNOOEf1,K)S(2.*CMASS/XSAR£A(1}-COEF*DELZ/EZ)/(2,-n£LZ*wS/EZ)
00117 , . . CNOD£(2,K)*2.»CMASS/XSAH£A(I)-CN00E(1,K)
00)18 , . . CNOOE( 1 ,K»3) = (2,*PMASS/XSAREA(1)-COEF*OELZ/EZ)/(2,-WS*OELZ/EZ)
00119 , , . CNOOE(2,K*3)»2,*PMASS/X9AREA(I).CNOOE(1,K»3)
00)20 , . ,..MN
00121 , . ELSt
00122 . . . CNOOE(1iK)aCHASS/XS*REA( I)
00123 , , . CNOOEC2,K)«CNOOE(1,K)
00124 FIN
00125 . ...FIN
00)26 , ELSE
00)27 . . CNOUEf Itl,K)«2.*CMA3S/XSAREA(I)-CNOOE(I,K)
00|2fl . . IF(>\ ,LE, 3)CNOO£(Ifl,Kt3)a2,*PMAss/XS*REA(J)-CNOOE(I,K*3)
00)29 . ...FIN
00130 C
00131 C . OVERWRITE BOTTOM ELEMENTAL NODE INFORMATION
00132 C
00133
00131
00)35 , EL10PrELTOP»DELZ
00)36 , NELMHT»N£LMTP
00137 . IF(K ,N£, 7)
0013fl , . PMS*B»PMSAT
00139 , , PMSB8*PMSBT
oouo , .. .FIN
0014) ...FIN
00142 RETURN
001«3 END
(FLECS VEHSION 22.46)
79
-------
(FLECS VERSION 22,«6)
13,MAR.81 t3l22lO<* PAGE OQOuI
00001
oooo?
00003
oooou
00005
00006
00007
OOOU8
00009
00010
0001 I
00012
0001)
OOOIU
00015
00016
00017
00018
00019
C
C
C
C
C
C
C
C
SUPHOUTlNt
THIS SUfoPUU
BY THE KNOW
COMB(M,
TINE MlJL
S, 1, R)
TIPLIES THE UNSYMME.TRIC BAND MATRIX (
LOAD VtCTOR AND ADOS TH£ RESULT TO ,
CALLED BY TRANSP.
INCLUDE 'iL
DIMENSION S
R ( ) ) rH ( 1 ) tS
H(M)»«(M) t$
DO ( I s2 j M« J )
, Y(1)»S(I
. H ( 1 ) SR ( 1
...FIN
RETURN
END
(MXEIEM,
( 1 r 2 ) *Z (
t
31, Rt^XtLEM), Yt"xELEH), Z(MxELEM)
l)fS(l,3)«Z(2)
(M, 1)*Z(M-1)+S(M,2)*Z(M)
,!)*Z(I-
) > Y ( 1 )
l)tS(I,2)*Z(ntS(I,3)*Z(I + l)
(FLECS VERSION 22,06)
80
-------
(FLECS VERSION 22,16)
iyio9i3« PAGE ooooi
OOOOI
0000?
00003
00000
00005
00006
00007
00008
00009
oooio
00011
00012
00013
OOOH
00015
00016
00017
00018
00019
00020
00021
00022
00023
0002'J
00025
00026
0002?
0002R
00029
00030
00031
00032
00033
SUBROUTINE OIAG(ECH02, ECM03( ECHO'S ECH05* ECH06, ECHO/,
1 ECHQ8, ECHU9, ECHOIO, ISEG, JSEG, SAVECH)
LOGICAL*! ECH02, ECH03, ECHO'4, ECH05, ECH06, ECHO?, ECH08,
i ECH09, ECHOIO, SAVECH, WRTSEG
DIMENSION JSEG(5), SAVECH(IO)
ECHOJ», FALSE.
ECH04*. FALSE.
ECH05*. FALSE.
ECH06=. FALSE,
ECH07sSAV£CH(6)
ECHOES. FALSE.
ECH09a. FALSE,
ECHOlOs. FALSE,
WRTSEGs, FALSE,
WHEN(JSEG(1 ) .EQ.O) .iRTSEGs , TRUE ,
ELSE
DO (J=l,5)
, IMJSEG(J).EO.ISEG) WR TSEG* , TRUE,
...FIN
..FIN
F(WRISEG)
ECH02»SAVECH( 1 )
ECHOisSAV£CH(2)
£CH()«sSAVECH(3)
tCH05agAvECH(U)
tCH06sSA\/ECH(5)
ECHOB=SAvECH(7)
£CH09sSAveCH<8)
EChOlO«SAV£CH(9)
..FIN
RETURN
END
(FLECS VERSION 22,
81
-------
(FLECS VERSION 22.46)
10-APH.ei
PAGE oooot
00001
o o o o i
0000 J
oooou
00005
00006
00007
OOOOB
00009
00010
oooi i
0001 2
0001 J
0001 0
0001 5
00016
0001 7
0001 8
000 1 9
00020
00021
00023
0002 J
ooo2«
0000.0)SHtAR STRESS COMPUTED USING HOTTO"
82
-------
(FLECS VFHS10N 22,16)
iji<*oiD2 PAGE oooo?
0005«
00055
00056
00057
00058
00059
00060
00061
00062
00063
OOOo'l
00065
00066
00067
00068
00069
00070
00071
00072
0007J
0007«
00075
00076
00077
00078
00079
0 0 0 fl 0
00081
00082
0008J
0009«
000«5
00086
00087
OOOflfl
00069
0009Q
00091
00092
0 0 0 9 J
0009(1
00095
00096
00097
ooo9e
00099
001 oo
noioi
0010?
0 0 1 0 i
001 ou
001 05
001 06
001 07
00108
001 09
c
c
c
c
c
C
C
C
C
C
C
SLOP, HYDRAULIC RADIUS AND SPECIFIC
"EIGHT OF KAlER (FREE FLOWING
"I V£R)
WHEN (ISEG ,EU, I )
, READd,!) NELEM,NRED,DLZSAV,B01V, BED, ALEN, tLEV,POR, PELEV
. WHEN (PELEV ,EU, 0,0) HIVER a .FALSE,
, ELSE RIVER e .TRUF.
...FIN
ELSE
, REAndii) NELEM, NBEO,DLZSAV,HOIV, BED, ALEN, ELEV, PUR
...FIN
»«. XYSU .- THICKNESS OF THE TOP BED LAYER
XCSO * BEO - (NBED-1) * BOIy
IF (ECHO)
. WRI TE (6, s ) ISEG
. WRITE(6,2) NELEM,N8ED,DLZSAV,BOIV,8ED, ALEN, ELEV, POR, XYSO
, IF (IStG ,EO, I)
. WRIT£(6,6) PELEV
, . WHEN (PELEV ,IQ, 0,0) W«ITE(6,7)
. ELSE *RI1E(6,8)
, , , ,F I N
...FIN
IF (NELL* ,LT, 0 .DM. NEL£*M ,GT. MX£LEM)
, CALL PUTERR( 1 J,NUME»R,HLDERH)
...FIN
IF (NI^EO ,LE, 0) CALL PUTERR(6,NUMERF»,HLO£RR)
IF (NH£U ,GT, MAXLEV) CALL PU I fcHK ( a , NUMERR , HLDEHP )
IF (OL7SAV ,LE, 0,0) CALL PUTERS ( 7 , NUMERR , HUDERR )
jF (HOjv ,LE, 0,0) CALL PUTERRfP, NUMERR, HLOERR)
IF («£0 ,GT. NBEO«BOIV ,OR, BED ,LE. ( ^BEO- 1 ) *BC> I V )
. CALL PUT£SR(9,NUMEHR,HLDE"R)
...FIN
IF (ALfN ,LE, 0.0) CALL PUTEHWtJOf NUME»H , HLDEHR J
IF (ELEV ,LE. 0,0) C«CL PUTERR< i I,NUMERH,HLPEHR)
IF (POR ,GT, 1,0) CALL PU f ERR ( 1 2, NUMERR , HLDERR )
CARD 2 AR£A OF EACH ELE^E^r
00 (Is|,MXELEM)
. AREA( I )sO.O
. EL( I )«0,0
...FIN
RtAD(l,9) NAREA, VPAREA, OELEV
WHEN (f4AREA ,EO. U)
ELt so.
D0( Is) , ^XELEM)
. AkEA ( l JSVPAREA
. EL(I)«FLE
. ELEsELEfOELEV
...FIN
..FIN
ELSE
83
-------
(FLECS VERSION 2?,«ti>
1U-APH-8I I3»«0|02 PAGE 00003
00) 0
001 1
00| 2
001 3
OOt U
001 5
001 6
00! 7
001 8
001 9
00120
00121
00122
0012}
0012«
00125
00126
0012?
0012*
00129
OOliO
001 Jl
00132
001 JJ
001 ju
001 iS
001 Sh
00137
00158
001 39
001 «0
00141
00142
00143
001 ID
001 «s
001 Ub
00147
C
C
1
2
1
i
3
a
s
b
7
8
S
4
5
6
7
1
8
1
9
10
C
KEAD (1,3) (AREMJ), I=l,NARtA)
«FAn(l,3) (ELCJ), I*l,NAH£A)
If (NAkEA.LT .MXELFM)
, Ou (I a NARCA* I ,MXELEM)
. . A»tA(I) « VPAHEA
. . tUl) « EKl-M » OELEV
. ...FIN
IF (ECHO)
wR!Tt
ViP]Tt
...FIN
RETURN
fj, »REAfl),
FORMAT (IHO.IJX, 15, '...NUMPER OF- VEHTICAL ELEMENTS'/
l«x, 15, • ...NUMBER nF RED L*vER8'/
7X.1PE12.5,
7X.1PE12.5,
7X,1PE12.5,
7X,lPt'2,5,
7x,|PE12.5,
7X.1PEU.5,
7x, 1PE12.5,
FURMAT(«FIO,0)
FORMAT(1H0.5PX,
.STANDARD ELEMENT THICKNESS (METERS)'/
.STANPAHD BED L*YEH THICKNESS CMETEHS)'/
.INITIAL BED THICKNESS (METERS)'/
.LENGTH OF THE SEGMENT (METE«5)'/
.SEGMENT ELEVATION (METEPS)i/
.POROSITY' /
.THICKNESS OF THE TOP BED LAYER (CALCULATED)1)
ELEMENT AREAS'/«(27x,5(I3,lpE12.5)/)>
INPUT DATA FOH SEGMENT 1,13)
FORMAT(7x,1PE12.5,'.. .UPSTREAM ELEVATION (METERS)')
FORMAT(19",'...SHEAR STRESS VALUES COMPUTED USING METHOD1
FOR RESERVOIR1 )
FORMAT( 1SX,»,..SHEAR STRESS VALUES COMPUTED USING METHOD'
FOR FREE FLOWING RIVERS')
FORMATtlS, 7F10.0)
FORMAT(1HO,S8X,'NOOAL ELEVATIONS'/4(27X,5(I3,lPEI2,5)/»
END
(FLECS V6HSION 22.U6)
84
-------
(FLECS VtHSION 22.46)
10123106 PAfJE 00001
00001
00002
00003
oooo'i
00005
00006
00007
OOOOPI
00009
00010
00011
00012
00013
00011
00015
00016
00017
00019
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00031
0003S
00036
00037
00038
00039
00010
0004J
00042
00043
OU044
OOOU5
00046
00047
00048
00049
00050
00051
00052
00053
SUBROUTINE. OISOLV(ABA«, 8, BU1V, C, COIN, COLO, OECAY, DELZ,
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
1 DELTO, OEMS, r>IAM, I, KAY1, KAY2,
2 NEtEM» NBEO, PCOEF, POK, WHIN, UHOUT, ^v, SQRBK,
3 SETA, veu, vet.2, HETAI, BETA2,
4 OEPO, SCOUR, HEDSD, XYSO, AREA, 090P3)
THI3 SUHHOiJTINE CALCULATES COEPFICIENT3 OF CCHVECTIVE, DECAY
ALFA,
ANO SOURCE/SINK TER^S OF DISSOLVED POLLUTANT CONVECTION-OIFFUSION
EQUATION
INPUT PARAMtTt»Sl
ABAH - AVERAGE AREA
AREA - VERTICAL PROJECTION AREAS ("2)
B . bEO CONCENTRATIONS
BOIV - STANDARD BED LAYER THICKNESS
C • NilDAL CONCENTRATION
COIN • CONCENTRATION OF INFLOW
COLO - CELL-CENTERED CONCENTRATION
'JECAY - OtCAY VALUES
OELTO - TIME STEP IN DAYS
OEMS • DENSITY
OEPO - DEPOSITION RATE ( KG /M2/OA Y )
OIAM « PARTICLE OIAMETFRS
OELZ • THICKNESS OF THE ELEMENT
I • ELEMENT INDEX
KAYl . LIGHT EXTINCTION COEFFICIENT OF HATER
KAY2 . LIGHT EXTINCTION COEFFICIENT OF SUSPENDED
SEDIMENT IN WATER
NELE*. . NUMBER OF VERTICAL ELEMENTS
PCOEF - 1ST TERM OF THE PHOTOLYSIS RATE EUlJATION, COMPUTED
IN SUBROUTINE PHQINP.
POP • POROSITY
dHIN « INFLOW DISCHARGE
UHOUT • OUTFLOW DISCHARGE
UV • VkRTICAL DISCHARGE
3COU» - SCHUR BATE (K6<>»C) /M2/OAY)
SORrtK • AOSOR8TION ON SEOTMENT
OSO«« - DESO^PTION FROM SEOI^ENT
XYSO . TOP HEU LAYER THICKNESS
OUTPUT PAflA"EfERSl
ALFA . DECAY TERM
8EDSO • SCnU« OR OEP03ITION OF AHSORPED CONTAHINANT
(PC/^2/DAY) WHEN NO SCOUR IS TAKING PLACE
BETA • Si'UHCE OR SINK TERM
8ETA1 « INFLUENT SOURCE TERM FOR I-TH NODE
BETAS • INFLUENT SOURCE TE*M FOH i+t TH NODE
VELI - FIRST CUNVECUVE TER"
VEL2 • SECOND CONVCCTIVE TERM
CALLED BY TRAIJSP.
INCLUDE 'ELM3l/,PR'i'
85
-------
(FLEC3 VERSION 22.46)
10|23»06 PAGE 00002
00050
00055
00056
00057
00058
00059
00060
00061
00062
00063
0006«
00065
00066
00067
00068
00069
00070
00071
00072
00073
0 0 0 7 '1
00075
00076
00077
00078
00079
00080
00081
ooo«2
00083
00084
000«5
00066
00087
00088
00089
00090
00091
00093
00095
0009a
00095
00096
00097
00098
00099
00100
0010J
00102
00103
ooioa
00105
00106
00107
00106
00109
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
REAL KAY,KAY1,KAY2
DIMENSION ABAR(MXELE"), B ( H A XLE V, *A XCUN- 1 ) , BE03D(3>,
1 CCIki(MXELEM,MAXCON),COLlJ( 1XELEM,HAXCON), DECAYC6)*
2 OHIN(MXELEM), QHUUT(MXF.LEM), UV***************t**ll******»*****l>k»*l>*»k*******l>***
* *
* wAWNINGi THE VALUE OF C3AR SHOULD HE UPDATED BY ITEHATIYELY *
* SOLVI'iG FOR C AT THE ADVANCED TI«E, AND APPROXIMATING *
* C8*» \S THE NEW AVERAGE CQNCENTHAT ION OVEH THE TIME *
* STEP. *
* *
it*********************************************************************
DO (IE « l,MXeUEH)
DO (1C a 1,MAXCON)
1
2
3
1
2
3
1
2
. CHAH{IE,IC) • CCIN(IE,IC)
...FIN
..FIN
0 (J » 1,3)
JP3 « J + 3
IPl « I + 1
IF(CBAF»(!,J).GT,O.O.ANO.CHAH(IPlfJ).ST.O,0)
AOOS1 a (SURHK(J)*30RHK(JPJ)/12.*(3.*CBAR(I, J)*CBAR(J,7)
*CBAR(I,J)*CHAR(IPl,n»CBAH(IPlfJ:i*CI4Arl(U7) +
CBAR(IPt,J)*CHAR(lP! , 7) ) - SORRK f JP3) X6 . * ( 2, *
CBARU.JP3) + CHAR(IP1,JP3)J)
OSA01 a (DSO«B(J)*nsOHH(JPS)/12,*(3,*CBAPi:i,J)*CHAH(l,7)
»CB»R(I,J)*CHAmiPl»7)+CHAH(IPl,J'J*CBAB(I.7)»
CBAH(IP1,J)*CBAR(IP1.7)) - 030«iHJP3)/6.*(2,*
CHAR(I,JP3) + CBAH(IP1, JP3)))
AODS2 a (SORHK(J)*Sn»BK(JP3)/12,*(C(*AH(I, J}*C8AR(I»7)+
CH*H(I»J)*CU»«(ICl»7)+CaAf»(IPl»J)*CH»H(If7)+
3,*C8AR(IP1, J)*CKAR(IP] ,7) ) - 30PWK(JP3)/6.*
86
-------
(FLECS VERSION 22.<
2Q.jAN.e2 f«f23l06 PAGE 00003
00110
00111
00112
00113
001H
00115
00116
00117
001 18
00119
00120
00121
00122
00123
00121
.00125
00126
00127
00128
00129
00130
00131
00132
00133
00131
00135
10136
00137
00136
00139
00110
001«1
00142
00113
00111
00115
001«6
00117
St .liT.O.O.UK.ADDSl.EQ.nSAOl )HETA 1»bETA 1-A0031
IF(U3A01.LT,0.0)9ETAl»HtU1-OSAD1
IF(*noS2.GT.O.O.OR,ADI}S2.Eg.OSAD2)RETA2aaETA2-APOS2
IF(l)SAD2.LT,0.0)dETA2*eETA2-l'SA02
..FIN
..FIN
C****************************************************************4
c TRANSFER HETWFE-I DISSOLVED STREAM CONTAMINANT ANN ARSORPEO
c BED CONTAMINANT is INCLUDED WHE^EVEW NU SCOURING OCCURS FOR
C » PARTICULAR 3EPIMEMT SIZE (EG SAMO, SILT, "« CLAY)
C**« **•!
8ETA a 0,0
If (I.E'J.l)
00 (J»l,3) BEnSO(J) • 0.0
BEr>SO(J) •
IFINPED.GT.O)
00 (J«l,3)
«HEN(SCOU«(JJ,GT,
, HfD30(J)aO.O
...FIN
ELSE
O.O.OR.H(NMEL), J) .LE.ZERO)
. 0 • OIAM(J)
. IFf!),GT.XYSO) 0«XYSO
. RATE a SOReK(Jt6)*(SORHK(J)*fCOAH(l,7)+CB*R(2,7))/2,
. -B(NBED,J+3))* D * RHOJ
. HETA a bETA . RATE /DEL?
. PEDSO(J)a-RATE
...FIN
..FIN
..FIN
..FIN
RETIIHN
00118
00149
00150
00151
00152
00153
00151
00155
00156
00157
0015?
00159
001*0
00161
00162
10 COMPUTE-PHnTOLYSIS-RATE-FOR-ELFHEHT-I
AVGSEU a o.o
tiiHEN (i ,E(J. NELE*U
. 00 (IJI«I,3) AVGSED » AVGSEO t C(Ifl.IJI)
...FIN
ELSE
. 00 ClKaI+l,NELE*+l)
. . 'JO (tJt«l,3) 4VGSEO « AVRSED * C(IK,IJI)
. ...FIN
. AVtSSEO a AVG3EO / (NELE^-f 1-1)
...FlN
KAY • KAYl + K4Y2 * AVGSEO
"HEN (I ,£iJ, NEI.FM) TE^"1 a 1.0
ELSE TERM1 a EXP (-KAY*(MEL£M-I)*OFLZ)
» (1.0 - EXP(-KAY*OELZ) ) / f*»Y*r>ELZ)
87
-------
(FLECS VEHSION 22.<»6) 2fl-JA*-P2 10»?3:06 P»GE 0000«
0016S . PHOTO
CROSS-REFERENCE
OOIA8 CttMPUTt
00072
(FLECS VEWS10N 22,
88
-------
IFLIC3 vEKSION 22,06)
20-MAR-8I
10122125 PAGE 0000!
00001
00002 t
0000) C
00000 C
00005 c (
00006 c r
00007 C
00008
00009 1
00010 C
00011 (
00012 1
0001} C
oooia '
00015
00016 t
00017 I
00018 (
00019
00020 C
00021 I
00022
00023
00020
00025
00026
00027
00028
00029
00030
00031
00032
0003J
00030
00035
0003*
00037
00038
00039
00040
0 0 0 ') 1
oon«2
00003
ooooa
00005
00006
00007
ooooe
oooso
00051
00052
00053
SUBROUTINE EQUPCSCPXSAR, UWIO, UEL, XSAREA, NELEM, MELEM,
RATIO, IELP, HEQXS, FERHOR, DELTA)
THIS SUBROUTINE FINDS CROSS-SECTIONAL AREAS AND HEIGHTS WITHIN THE
JPSTREAM CROSS-SECTION WHICH CORRESPONDS TO THE SEGMENT IMMEDIATELY
DOWNSTREAM,
INCLUDE 'ELMSIZ.PRMi
.OGICALM FERROH
DIMENSION pxsAR(MXELEH), uw ID ( MXELEM >, UELCMXELEM),
XSAREA(MXELtM), IEIP(MXELEM)» H£QXS (MXELEM)
•£RROR», FALSE,
[P«l
'XS3PX5AR ( 1 )
I8TMS(JW1D( 1 )
:LBT*sUtL(l )
'EMpxjso .
)0(I»1, NELEM)
XS*»AUO*XSAHEA( I )
UNTILCXS .(.I. PXS ,OR. IP ,EO. MELEM)
. 1P*1P»I
. TtMPxs»PX3
, U8TM«U«(ID(IP)
. ELUTMsUEL(IP)
. PXS*PxStPXSAR(iP)
t t t F J N
IELP(1)«IP
WHEN(X3 ,EO, PXS)
, HtOXS(I)sUELflPtl)
. IFU.EO. NELEM) HEQXSU)* UELUP) * DELTA
ELPTMruELl 1P+1 )
, UBTM«UWID(IP» 1 )
, PxS«0
, TLMPXSrQ
, , ,F IN
ELSE
t *s(UwiD(lPfl).UBTM)/(2t»fUEL(lPfl)»EL8TM))
. RsUBlM
. C«XS-TEMPX5
, WHEN (A ,tu. 0.) HEUXS< n*C/8tELBTM
. ELSE
. . HSUUiC = 8*Bf a . *A*C
, . IF (8SUOAC ,LT. 0.) GO TO 200
, . MEUXSU )«(SORT (BSUOAC )-B)/2,/A+£L8TM
. ...^IN
. ELFHMsHE'JXSU)
, Pxs=°xs-xs
. TtMPXSsO.
, , ,F JN
89
-------
(FLECS VERSION
ooosa
00055
00056
00057
0005?
00059
00060
00061
200
PAGE
...FIN
Rt TURN
CONT]NU£
FERNOhs,IKUE.
WR11E(b»1)
FORMAT ( 1 OX, 'FATAL ERROR
RETURN
END
- 8SD«AC IN EIJUPCS SUbROUTINC IS < 0')
(FLECS VERSION
90
-------
(FLECS VERSION
iji2«i3o PAGE ooooi
00001
0000? 1
OOOOi C
0000«
00005 C
00006
00007 1
OOOOP C
00009
0001 0
0001 1
00012
0001}
OOOH C
00015
00016
00017
00018
0001 9
00020
00021
00022
00021
00025
00026
00027
0 0 0 1? fl
00029
00030
00031
OOOJ2
OOOJJ
0005'J
OOOJ5
00036
00037
0003S
00039 F
00040 1
3UBHOUTJNE £QUPXS(PX3AR, PW ID, POELZ , XSARE A , NELEM , M£LE M ,
RATIO, IELP, HEUXS)
INCLUDE ' SY lELMSIZ.PRMi
DIMENSION PxSAP(MxELEM), Pw I D( MxELEM) , xSAR£ A ( MxELEM) ,
IELP(MXELEH) , HE(JXS(MXEL£M)
;Pa I
'XSsPXSAHCl )
,JBT, M = pMD( 1 )
•LBTMio,
TEMPXS«0 .
J0( let ,N£UEM)
XS=KATIO*XSAHEA(I )
UNTUtXS ,LE. PXS .OR, IP .EO, M£LEM)
. IPeJPf 1
, T£MPxgsPXS
, wBTMxpMlO(lP)
, EL8 TMsPDELZ» ( IP-1 )
. P*SSPXS»PXSAR < IP)
lELP(I)«IP
wH£N(xs ,ew, PXS)
. HEQXS(I)aPDELZ • IP
, PXS=0,
, rEMPxsso.
, EL8 fMa (P*PO£L I
, wBiM»PwID{ IP* 1 )
, , ,F IN
ELSE
, HtUXS(I)s(XS.IFMPXS)/W8TM+ELBTM
, Pf 3*PX3-X3
TEMpXSsO,
. Et6TMSHEOX3( I )
...FIN
..FIN
!ETURN
NO
VERSION ^^,Ht>)
91
-------
(FLECS VERSION
13!24|
-------
(FLECS VERSION 22.46)
PMiE
0005«
00055
00056
00057
00058
00059
00060
0006 t
0006?
0006J
0006
-------
(FLECS VERSION 22.46) 13-MAR-ei IJ|2'J|«1 PAGE OOOOi
00)10 C
oont c •** iNstwi NULL CHARACTER ***
00112 FNAM£( ICA«)s()
0011J RETURN
0011« END
(FLECS VERSION 22, (16)
94
-------
(FLECS
J.MAR.
PAGE
00001
00002
00003
o o o o a
00005
00006
00007
00008
00009
00010
0001 1
0001?
00013
0001 U
00015
000 \ 6
0001 7
00019
00019
00020
00021
0002?
00023
0002«
00025
00026
00027
0002^
00029
00030
00031
0003?
00033
ooo 51
00035
00036
00037
00039
00039
00040
000«1
0001?
00013
Q 0 Q 4 11
00015
00016
OOOU7
oooan
00019
00050
00051
0005?
00053
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
SUBROUTINE FDC ODE (F N»ME, 8ASE, NHRTP.F TYPE rDEV, UIC 1, UIC?)
THIS ROUTINE SPEHATES FNAM£ INTO 6 COMPONENTS
RASE - 5 CHARACTER HASE FILE NAME (BYTE ARRAY)
WHRfP . TIM£ PLANE NUH8ER THAT IS THE LAST « CHARACTERS OF THE
9 CHARACTER FILE NAME (INTEGER)
FTYPE . FILE EXIENSION (BrTE ARBAY)
OEV . PHICAL DEVICE SPECIFICATION (BYTE ARRAY)
UIC1 - 1ST UIC
U1C2 • 2ND UIC
THE OPTIONAL PARAMETERS OEV, UtCl, AND ulC2 WILL BE SET
TO BLANKS IF NOT PRESENT IN THE ORIGINAL FILE SPECIFICATION.
CALLED BYl STHTUP
BYTE F NAME (27,, BASE (5) ,FTYPE<3) , OEV (3), UIC 1 (3) ,UIC2(3) ,
1 LHRAK,RBRAK, COMMA, PERIOD, COL ON, BLANK
DATA LBRAH/ i ( t /
DATA RBKAK/ * ] ' '
DATA COMMA/' , >/
DATA PEMUO/','/
DA TA C'JLON/ ' : » /
DATA BLANK,/ i 1 /
*** FILE SPECIFICATION HAVE FOUR POSSIBLE FORMS ***
(1 ) FILENAME. EXT
(2) DEViFILENAMt.EXT
C3> (UIC1,UIC2IFILENAME,£XT
(«) DEvi (UIC1,UIC2)FUENAME. ExT
TMF FORM CAN BE DETERMINED B* COUNTING THE FOUR SPECIAL
CHARACTERS [l.i
Nso
DU ( Ts| ,27)
. SELECT (FNAME(l))
, . (COLON) NsMt)
(LBRAK) NaNH
. .
-------
(FLECS VERSION
mzsiio PAGE 00002
OU05«
00055
00056
00057
OOOS8
00059 C
(2) DECOOE-FORM?
{}) OECODE-FORM3
(U) UECODE-FORMil
RETURN
00060
OOObl
00062
0006J C
TO DECOOE-FORMl
. DECODE-nUNAME-EXTENSlON
OOOb"
00065
00066
00067
00068
TO OPCODE-FORM?
, DECODE-OEVICE
...FIN
00069
00070
00071
00072
0007J
TO DECOUE-FOBH3
. OFCOUE-UIC
EXTENSION
...FIN
0007a
00075
00076
00077
00078
00079 C
TO DECODE-*ORM
-------
(FLECS VERSION 22.06)
13.MAR-8i i3i25tio PAGE 00003
00095
00096
00097
00098
00099
00100
00101
00102
00103 C
ootou
0010S
00106
00107
00108
00109
001 10
00111 C
00112
00113
001 Id C
001 15
00116 C
00117
001 IS
00119
00120
00121
00122
00123
0012« C
00125
00126
00127
00128
00129
00130
00131
00132 C
001 33
0013«
00135 E
t ]C*R*IC*R*1
, NBRTP»0
...FIN
DO (I«l,3)
, FTrP£(I)«FNAME(ICAR)
IC»R»ICARH
• • .FIN
..FIN
TO OECOUE-OEVICE
1*1
REPEAT *HILE (FN»ME(ICARJ ,NE, COLONJ
, l)EV(I )sFNAM£( 1CAR)
, ICA«*ICAR»1
, I»I » 1
...FIN
»*• SKIP OVER COLON ***
jCAR«ICARfl
..FIN
ro DECOUE-UIC
*** S«IP OVER LEFT BRACKET ***
ICAPBlCAH+l
1 = 1
REPEAT rtHILE (FNAME(ICAR) ,NE. COMMA)
. U1C1 ( l)*FNAME(ICAR)
. 1C AR»ICAR+1
, I * I » 1
...FIN
»** S"«IP OVER COMMA **•
I CAP* iC AR» i
1:1
REPF.AT *HILE (FN»ME(JCARJ ,NE. RBRAK)
, U1C2( I)t'FNAME( 1CAR)
. 1C A&«KAflf 1
. 1*1 t 1
...FIN
«** SKIP OVER RIGHT BRACKET **«
ICARsICARt 1
. , F I N
.NO
PROCEDUHE
00060 DECOO
0005J
no<>6« DECClOt-F
TABLE
97
-------
(FLECS VERSION 22,16) 13.MAN.81 I3I25HO PAGE 00001
0005«
00115 OECODE-UIC
00070 00076
OOOfjP DECODE-FOKM3
00055
0007« PECOOi
00056
00080 OECOOE-f
00061 00066 00071 00077
OOlOU OECOOE-UEVICE
00065 00075
(FLECS VERSION 25.«6)
98
-------
(FLECS VEKSJON
13132:1,? PAG£ OOOUl
00001
ooon?
0 0 0 0 J
ooooo
00005
00006
000(17
00008
00009
00010
0001 1
00012
0001 3
00014
00015
00016
0001 7
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00050
00031
00032
00033
00031
00035
00036
00037
0003B
00039
oooao
000
-------
(FLECS VERSION 22,
. . . tNDTIM«(ENDTIM/IOELTH*l)*IOELTH
FIN
If (fcNOTIM .NE. -9999) PflEUM *
UNLESS (ENOTIM ,EO. -9999}
DEPTH FLAG TO DIVERT FROM FURTHER CALCULATION
I.
VUL«0.
VtL'O.
D0( 1»1
. AHAR(I)sO.
. l£LM(I)aO.
, XSAREA(I)»0,
...FIN
IF (OEPTH ,GT, OEPMIN)
', *«* COMPUTE! NELEM,,..NUMBEH OF ELEMENTS CONTAINED WITHIN DEPTH
ABAR(I),.AVERAGE AREA OF ELEMENTS I AND 1+1
. AWID(I ) ..WIDTH OF F.LEMENT
, VOL TOTAL VOLUME OF TH£ SEGMENT
. VUL=0.
. NELEMrDEPTH/DLZSAV
. uELZ«OEPTH/NEL£M
. WHEN(NELEM ,LE. I .OR. NELEM+1 ,GT. MXELEM)
. CALL PUTERR U«, NUMERR, HLDEHR)
. WRIIE(6,6)
. ...FIN
. fcLSE
, CALL TRNPOStABAR,AREA,AWID.ALEN,8WIO,OELZ,EL»IELM,
. , XSAPEA,VOL)
. CALL RADIUS (ALEN, AREA, C»OSEC, DEPTH, EL, HRAO)
. ...FIN
* IT IS IMPLICITLY ASSUMED THAT A DOWNSTREAM COURANT
. NUMBER AT OR NEAR UNITY HAS BEEN EMPLOYED IN THIS
, ANALYSIS
t «o)/2. * ALEN/VOL
TIH,NEL£MrOELZ,oI,gC
DEPTH,H*ID,AHAR, HRAO, CHuSEC
...FIN
WHl 1E(«
X5AHEA,
IF (ECHO)
««IT£(6,2) ENOTIM.NELEM,TEMPH.yI,QO,DEPTH
', 00 (|sl,NEL£M)
. , WRnE(6,y)I,A«IO(I),ABAH(I),AKEA{I),EL(I)
. ...FIN
100
-------
CFLECS VERSION 2?.«6)
18-MAH-81 13U2I12 PAGE 00003
00110
001 1 I
00112
001 13
001M
001 IS
00116
00117
001 IS
00119
00120
00121
00122
00185
0012"
00125
00126
00127
0012ft
00129
00130
00131
00132
FIN
. ...FIN
...FIN
IF (PREIIM .LT,
REMIND «
RETURN
CALL PUTERH<25. NUM£RR, HLDERH)
FURMATJ MO,IF 10,0)
FORMAT!!Ho,5eX,'HYDROLOGY DATA'/
9x,no, '...DATA su ENDING
MX,15, ' ...NUMBER OF ELEMENTS WITHIN THE FLOW DEPTH'/
7x,IPf12.5,'...WATER TEMPERATURE'/
7X,1PE12.5,'...TOTAL DISCHAHGE OF THIS SEGMENT'/
7x,lP£12.5,'...TOTAL DISCHARGE OUT OF THIS SEGMENT1/
7x,lPE12,5,!,,,FLOW OEPTHi)
3 FORMAT(I HO,'ELEMENT',«*,'SEGMENT',1«X,'AVERAGE'/12X,'NODE VP-AR£A' ,
1 I1X,'NOUE ELEV'/IX,'NUMBER',6X,iWIOTHl,12X,'ELEMENT AREA')
a F(jRMAr(3x, !2,'(X,lPEl?,5,l?X,lPei2,5,10x,lPE12.5,10x,lPEI2.«)
6 FORMAT(//50x,'DEPTH TOO GREAT FOR THE MAXIMUM NUMBER OF'
1 ' ELEMENTS'//)
C
END
(FLECS VERSION 22,16)
101
-------
(FLECS VERSION 22.U6)
28-MAK.Hl 03127108 PAGE OOOOt
00001
0000?
00003
ooooa
00005
00006
00007
0 0 0 0 B
00009
00010
0001 t
00012
00015
OOOt <4
00015
00016
00017
00018
00019
00020
00021
00022
0002i
0 0 0 2 U
00025
00026
00027
0002B
00029
OOOJO
00031
OOOJ2
00033
00030
00035
00 Oib
00037
0003K
00039
oooon
0 0 0 1 1
000^2
OOnu j
0001411
000(45
00050,
1 ELFV, ENOHYO, ETIME, FERROH, HRAD, Nfc'LEM,
2 NF.WUI, PELEV, QHIN, QHOUT, ov, HTVEH. SLOPE,
3 STRESS, TE*PH, VEL, VOL, OEPMlN,
1 XSAREA, Hwlf), ABAR, 01, CROSEC)
THIS suppouriNE rs CALLED EACH TIME STEP TO READ ANY NEW HYDROLOGY
DATA THAT WAS WRITTEN TO LUN « a* SUBROUTINE HYDDAT,
FORMAI PAKAMflEHS:
ALEN - LENGTH OF TH£ SEGMENT
A)7fcA - AREA OF fcACH ELEMENT
AW10 - tLErt£NT wJOTHS
"ELZ - STANOAHO FLEw£NT THICKNESS
DEP1H • FLOW DEPTH
DSO - MEDIAN BCD SEDIMENT DIAMETER CMETER)
ELEV • SEGMENT ELEVATION
ENOHYD - ENDING TIME OF THE CURRENT HYDROLOGY DATA u*«)
ETIME - ELAPSED TIME OF THE SIMULATION (IM)
FERHOR • FATAL ERROR FLAG CL*I)
HHAD • HYDHAUuIC RADIUS
NELEM - MUMSER OF ELEMENTS
NE^QT - N£« 01 DATA FLAG
-------
(FLECS VFHSJON 22,«6)
0005U F£RHl)R
00055
OOOS6
00057
00058
00059
00060 1
OOObt
00062
00063
00060
00065
00066
00067
00068
00069
00070
00071 1
00072
00073 1
0007« C
00075
00076
00077
00078
00079 C
00080 C
00081
00082
00083
0008«
00065
00086
OOOfl?
0008B C
0008"? 200
00090
00091
00092 1
00093
ooo9a 3oo
00095
00096
00097
0009B 3
00099 2
00100
00101 C
00102
28.MAH-SI
.FALSE.
03l27»0fl PAGE 00002
IF (ETI*E ,GT, ENOHYO)
NEWQJ a .THUE,
REPEAT UNTIL (ETIME ,LE, ENDHYD)
ENDHYO,NELEM,DELZ,QI,UOrV(JL,VEL»i
IF (DEPTH ,GE. DEPMIN)
*H£N (RIVER)
. C*LL SHEARSfALEN,ELEV,HHAD,PElEV,SLOPE,STRESS,USTAR)
F.LSt.
. CALL SH£ARH(DgPTH,050,STRESS,U3TAR,VEL)
' (RHO*HRAD)
, CALL PROFIL(ALF.N, AWID, DELZ, DEPTH, NELEM, QI, UST*RI
, VULf QHlN,DUMMY,o,OELZ)
. CALL PROFIL(*LEN, AWID, DELZ, DEPTH, NELEM, 00, USfAR,
, VOL, OHOUT,DUMMY,0,OEuZ)
. «** CONVERT UNITS TU M**3/OAY **•
. DO (Js1,NtLEM)
QHIN(J) « (JHIN(J) * SECDAY
. . UHOUl(J) e QHQUT(J) » SECDAY
. ...FIN
[ *.* COMPUTE VERTICAL FLOWS ***
, (3V ( 1) « 0.0
, 00 (J*l,NELEM)
. OV(JtJ) a UHINCJ) - QHOUT(J) » QV(J)
. ...FIN
...FIN
..FIN
CONTINUE
FERROft * ,TRUfc,
WHITE(6» 1 )
FO«M4f(]OX, 'FATAL ERROfl » HYDROLOGY DATA EXHAUSTED')
RETURN
CON!INUE
FERROR«,TRUE.
FORMAT(10X,'FATAL ERROR » VERTICAL FLUX COMPUTATION')
FORMAT(15X,I5,1PE]2.'J)
RETURN
END
(FCECS VEKSION 22,«6)
103
-------
(FLECS VERSION 22,46)
2P.JAN.R2 10121121 PAGE 00001
00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00015
00014
00015
00016
00017
00018
OOOTJ
00020
00021
00022
0002J
00024
00025
00026
00027
0002«
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
000<40
ooo«t
00042
f>00«3
00044
00035
00046
OOOU7
00048
00049
00050
00051
00052
00053
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
SUBROUTINE ICFLU(CCIN,OEPTHf DELZ,r>50,KNOIC,ETIME,FERFl'nK,
1 TNFLO, ISEG, NELEM, NEWIC,
2 UHIN,OI,OEPMIN,ALEN,
3 HEL, UWID, XSAHEA, AR£», A«(IO, i)FZ, V3ET,
4 ELfELEV.PELEVfRIvER^twuifNEwTRB)
THIS ROUTINE 13 CALLED EACH TIME STEP TO READ THE INITIAL
CONDITIONS TO THE FIRST SEGMENT OR THE INFLOWS
FRO* THE PREVIOUS SEGMENT.
FORMAL PARAMETERSI
CCIN CONCENTRATION OF INFLOWS- CELL. CENTERED
DEPTH FLO* OEPTH OF THE CURHENT SEGMF.NT
OELZ STAMOARO ELEMENT THICKNESS OF THE CURRENT SEGMENT
ENOIC ENDING TIME OF THE CURRENT INITIAL CONDITIONS DATA
ETIME ELAPSEO TIME OF THE SIMULATION (1*4)
FEHHOH FATAL ERROR FLAG (L*l)
INFLO LOGICAL UNIT NUMBER FOR DATA FROM PREVIOUS
ISEG CURRENT SEGMENT NUMBER
NELEM NUMBER OF ELEMENTS IN THE CURRENT SEGMENT
NELEMH NUM&EH °p ELEMENTS IN THE PREVIOUS SEGMF.NT
NEWIC INITIAL CONDITIONS FLAG a*i)
POELZ STANDARD ELEMENT THICH^ESS OF THE PREVIOUS
POEPTH FLO« OEPTH OF THE PREVIOUS SEKMfcNT
UHl^t INFLOW DISCHARGE
SEGMENT
SEGMENT
UHQLD DISCHARGE INTO THE SEGMENT FROM THE PREVIOUS ONE
CALLED fln SERATRA
CALLS! EOUPCS, fcQUPXS, PROFIL, RADIUS
INCLUDE 'ELMSIZ.PRH'
INTEGER*a EMOIC,ETIME
LOGICAL*! NEWIC,Nfc»QI,NEwTRa,FERROH,RI VtH
DIMENSION CCI'-t'^XELEM.MAXCON), CNOOE(MXEL£M, MAXCON) ,
1 QHINMXELF H)» (4HOLD(M)(ELEM)» XSAREA(MXtLEM) ,
2 PxSAt!(MxELliMj,ieLP(MXELEM),Ui«ir)(MX£LEM},uxOAVG(MXELEM),
3 HEUXS(MXFLE^) , H£L(MXELeM),TMA3S(^AXCON),Pi»ID(MxeLEM)f
U OCNDOE(MxELEM,MAXCON) ,
5 ARE»(^XELEM), AWID(MXELEM), MFZ(a), VSET(3)»
t> FL(MXELEM), UAHEA(MXELEM)
DATA RHO/1000.X
DATA SECoAY/BbUOO./
FEPROR » .FALSE.
NE"IC * .FALSF..
iniM3Mj(ELEM
JOIMaMAKC"M
104
-------
(TIECS VERSION 22.«6)
10I2H21 PACE
0005"
00055 C
00056 C
00057 C
00058 C
00059 C
00060 C
00061
00062
00063
00061 C
00065
00066
00067
00068
00069
00070
00071
00072
00073
ooo/a
00075
00076
00077
00078
00079
00080
00081
00082
00083
0008'!
00085
00086 10
00087
00088
00089
0009Q 1
00091
00092
00093
0009'4
00095
00096
00097 C
00098
00099
00100
00101
P0102
00103
00100 C
00105
00106
00107
00108
00109 1
»HE
t (ISEG , E'v, 1)
DISTRIBUTES INITIAL CONDITIONS UPSTHEAM OF INITIAL SEGMENT IY
CONSERVING RELATIVE CROSS-SECTIONAL AREAS AM[) QISCM&kuES---
ASSUMES LINEAR PROFILE OF MOTH 3EUIMEWT AND PAWTIC'/LATE, AND
NI|iIO{I) + l.il<(IO(Itl))/2.
U'
-------
(FLECS VERSION
28-JAN-H8 10I2H.J1 P4GE 00003
00110
om 11
00112
00113
001 1'4
00115
00116
00117
0011B
00119
00120
00121
00122
00123
00124
00125
oni2b
00127
0012H
00129
00130
001 H
00132
00153
00134
00135
00136
OM37
00138
00139
00140
00141
001«2
001 43
OOlaa
001'45
00146
00147
00148
00149
00150
00151
00152
00153
00154
00155
0015*>
00157
00156
00159
00160
00161
00162
00163
00164
00165
UU
iJHOL'>(l)»«ttOLO(I)*SECDAV/{UwO*VGCB*WTOP/6,+C8*wB/3. )
1 ..... *UPS'JO
..... IF(K.NE,7)
...... GTsCNOOECNELMTPtUK + S)
...... GB«CMUOE(NELMTP,K+3)
...... r,TOP»GT»F»Cl»GH*F»C2
..... . PMS*l3(ET-HEL)*( .25*(*T*CT*GT+WTOP*rTnp»GTaP)
1 ..... , +CON»(rtT*CTPP*GTOP + ^ IOP*C r*6Tr)P + «iTOP*CTOP*GT +
2. ..... xT"P*CT*GT + wT*CTO"*r,T + rtT*CT*GTOH) )*UI'SOO
TOP + i'<8*
*GTOP +
2. ..... wH*CTOP*GT00+«TnP*C»*fiTuP*«TOP*CTOP*Gn) )*t)PSQO
........ fin
..... IMDIC«NE
..... CM4SS»0,
106
-------
(FLECS VERSION
10i?H21
00166 PMASSaO.
00167 IMTNOIC ,£'J. 0)
001^6 ...... CMASSaCMS3T-
00169 IF(K.NE,7)
00170 FIN
00171 JFdNOK.GE.U
00172 C"t*SSsCM9HT
00173 ...... IF(K.NE.7)
00170 FIN
00175 ..... IFUNDIC.GE.2)
00176 (>r)(JsNELM>
00177 OELauEL(J*n-UtL(J)
00178 CTaCNOOECJ-H.K)
00179 CtJaCMOr>E( J,K)
00180 ....... UP3QOar3HOLO(J)
001*1
00182
00193 .......
00181 IFCK.NE.7)
•90M5 STaCNODf (J+l,Kf 3}
00186
00187
00188 1
001P9 £ *UP500
00190 FIN
00191 .FIN
00192 FIN
00193 ...
00194
00195 CM3AB«CM3»T
00196 CMSHHaCMSBT
00197 ..... NELM8TSNEIMTP
00198 IF(K,NE.7)
00199 ...... P1SABaPM9*T
00200 ...
00201
00202 ...
00203 FIN
00504 C .....
"0205 C NQTtj CMASS IS IN (KG/M**3)
00206 C
00207 ..... COMPUTE-P^nFILE-VALUF.S
00208 FIN
00209 . . . ...FIN
00210 FIN
00211 . ...FIM
00212 ...Fin
00213 C DISTRIBUTES INITIAL CI'NOITIHNS UPSTWEAM OF SUWSF.UUENT SEGMENTS HY
00210 C CONSE^VlMG MASS FLUX — ASSUMES LINEAR UPSTHfAM 01S TRI *»KT IO'
-------
(FLECS
10I2H21 PAGF 00005
00222
00223
00224
00225
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236
00237
00238
002J9
00240
002U1
00202
002«3
002«5
002«6
002«7
002<»8
00250
00251
00252
002SS
00250
00255
00256
00257
00258
00259
00260
00261
00262
00263
0026U
00265
00266
00267
00268
00269
00270
00271
00272
00273
00274
00275
00276
00277
XS»0,
PX3«0.
00 {1»1,NELEM) X3*<3+XS4REA{I)
00 (JaljMELErt) FXSaPXS + PXSARd)
RATIUaPxS/XS
CALL
ALLOCATES HASS
no
CONSERVING RELATIVE CROSS-SECTKU AL AREAS
upsuo»uHOLO(:
CT»(C*tE(I'I8,K) + OCNODE(N8,K))/2.
HEL«Ht'JX9(I)
F*C2«(E1-MEL)/PDELZ
CTOPoCT*FAClfCB*FAC2
OSAT«FAC2*(CT+CTOP)*UPSJi>/2.
CMShT»FACl*(CTOPfCB) /2.* UPS 'JO
GT»(OQOE(NT,K3) * OC'TJOE (MT, K3 ) 1/2.
G8»(CNOOE(H8,K3) + OCNOOECMR, K3) ) /2.
GTtiPaGT*FACl*G»*FAC?
IF(K.,ME,7)
108
-------
(FLECS VERSION
2fl-JAN.«2 10I2H21 PAGfc 00006
00276
00279
00290
00261
00202
002*5
00284
002«5
00296
00287
00268
002P9
00290
00291
00292
00293
00294
00295
00296
00297
0029P
00299
00300
00301
00302
00303
00304
00305
00306
00307
00300
00309
00310
00311
00312
00313
00314
00315
00316
00317
00316
00319
00320
00321
. , . . ifciNoic.Gfc.2)
,K) * >JCNOO£(J»l,K)}/2,
C8«
-------
(FLECS VERSION 22.46)
28. JAN, 82 10J2U21 PARE 00007
00331
00332
00333
0033<*
00335
00336
00337
00338
00339
00340
00341
00342
00343
00344
00345
00346 I
. "«(SsVSET(K)*AREA( J ) /( AHJO(
. EZ«DpZ(K)
. CCIN( 1 ,K )«(2,*CH*S8»COEF*
CCIN(2,K)B2,*CMASS-CCIN(1
. CCINO rKK)»(2,*PM*S8»COEF
. CCIN(2.KK)«2.*PMASS-CCrU
..FIN
...MN
ELSE
. CCIN(I+1,K)«2.*CMASS-CCIN(1,
. IF (K.NE.7)
. . CCIN( 1+1 t K+ 3)"2,*PHA35-CC
. ...FIN
...FIN
..FIN
:NO
i)*4L
DELZ/
»K)
»ntLZ
1,KK)
K)
n(I>
PROCEDURE CHOSS-HEFERENCE TABLE
00322 COMPUTE-PHOFIIE-VALUES
00207 00304
(FLECS VERSION 22,46)
110
-------
(FLECS VERSION SS,06)
IJ.MAH-61
I3J27M2 PAGE 00001
00001
00002
00003
0000')
00005
00006
00007
0000*
00009
00010
00011
•00012
00013
0001U
00015
00016
0001 7
OOOlfl
00019
00020
00021
00022
0002J
00n2<"
00(125
0002b
00027
0002P
00029
00030
00031
OOOJ?
OOOJ3
00034
00035
00036
00037
0003R
00039
00040
OOOU1
00042
OOOU J
00001
OOOM5
OOOUb
000«7
OOOafl
00019
00050
00051
00052
00053
C
C
C
C
C
C
C
c
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
SUBROUTINE IN1DAT (ANALMT, ANALYS. OEITH, ECHOr HLDERR. ITPRT,
1 NSEG, N3TEPS, NUMEKR, SlMlEN, DEPMIN)
THIS ROUTINE READS THE INITIAL DATA COMMON TO ALL SEGMENTS AND IS
ONl.r CALLtO FOR SEGMENT NUMBER 1,
FORMAL PARAMETERS!
ANALMT ANALYSIS CONCENTRATION LIMIT
ANALYS TIME SERIES ANALYSIS CONTROL VARIABLE (L*l)
UELTH TIME STEP LENGTH (SECONDS)
ECHO LINE PRINTER ECHO CONTROL VARIABLE (L*l)
HLDEHR HOLDING ARRAY FOR ERROR NUMBERS (BYTE)
ITPRI PRINT FREQUENCY
NSEG NUMBEH OF SEGMENTS
NSTEPS NUMBER OF TIME STEPS TO BE TAKEN (I*«)
NUMERR NUMBER OF INPUT ERRORS DETECTED
SIMLEN SIMULATION LENGTH (SECONDS - I»«)
DfcPMIN • MINIMUM (CUTOFF) FLOW DEPTH (METE«S)
CALLED BYI SERATRA
CALLESl PUTERR
BYTE HLDENR(IOO)
INTtGEfl*« SIMLEN, NSTEPS
LOGICAL*! ECHQ.ANALYS
DIMENSION TjrLE(«0)
DATA M»xsE(i /35/
If (ECHU)
ff
, *»» PRINT HEADING •**
. WRIU(6,1)
, »(RlTE(b|2)
,..F1N
CARDS 1 AND 2. ....SIMULATION IDENTIFICATION TITLE
READ(1,3) (TITLE (I ) , ls| ,i|0)
IF (ECHO) WRITE16.U) ( Tl TLE( I) , I* 1 , -45. .ANALMT, . ,.LO«ER Ll*iT OF AVERAGE DISSOLVED
CONCENTRATION, BFFOKE THE RESULTS OF A
TIME STEP AHE SAVED, THE AVERAGE
111
-------
(FLECS VERSION 22.06)
1 3,MAR.61
I3i27ji2 PAGE 00002
00051
00055
00056
00057
00058
00059
00060
OOflbl
00062
00063
0006U
00065
00066
00067
00066
00069
00070
00071
00072
00073
0007
00077
00078
00079
00080
00061
00082
00081
0008U
00065
00066
00067
00068
DISSOLVED CONC, MUST BE >
06-55. ..DEPM1N... .MINIMUM (CUTOFF) FLOW DEPTH BELOW WHICH
THE CHANNEL is CONSIDERED DRIED,
READ (1,5) NST EPS, NSHG, I TPRT, ANAL YS, DEL TH, ANALMT, DEPM IN
ttt COMPUTE SIMULATION LENGTH (SECONDS) »**
SIMLEN « DELTH
S1MLEN « SIMLEN » N3TEPS
IF (ECHU)
, WRITE (6, 6) NSTEPS,NSEG, ITPRT,ANALY3»DELTH,ANALMT,!3IMLEN,DEPMIN
...FIN
IF (NSTtPS ,LE. 0) CALL PUTERR ( 1 , NUMERR, HLDERR)
IF (NSEU.LE.O .OR. NSEG , G T , MAXSEG) CALL PUTERR (2, NUM|=RR, HLOERR)
IF UTPHT ,LE, 0) C*LL PUTERR(3, NUMERR, HLDERR)
IF (OtUfH ,LE, 0.0) CALL PUTERR (5, NuM£HR, HLOtRR )
RETURN
FORMAT dHOrjox, 'SEDIMENT AND CONTAMINANT TRANSPORT SIMULATION',
| I PROGRAM . StRATRA')
FORMAT ( lHo,5qx, iPROKLEM SPECIFICATIONS')
)
25X>20AU/?6X,20AU)
2I5,L5,2F10.o,E10.3)
6X, I 10, '. ..NijMBEH OF TIME STEPS TO BE
..NUMBER OF SEGMENTS'/
. , PRINT FREQUENCY <* OF TIME STEPS)'/
TIME SERIES ANALYSIS CONTROL1/
,,. TIME STEP LENGTH (SECONDS)"/
FORMATf IMO
FORMAT(110
FORMAT f IHO
1 1«X, Ib, »
i «x, is, '
I 8X,L 1, '
7x, IPt 12
7X, JPE12
9X, I 10, '
9x,E10.i
END
..
5, '
5, » ...TIME SERIES CONCENTRATION LIMIT'/
. .COMPUTED SIMULATION LENGTH (SECONDS)'/
' ...MINIMUM (CUTOFF) FLOW DEPTH (
(FLECS VERSION 22, U6)
112
-------
(FLECS VERSION 22,«6)
PARE opoo)
OOQOi
00002
00001
oooou
00005
00006
00007
0 0 0 0 8
00009
00010
0001 1
00012
00013
00014
ooois
00016
0001 7
oooie
00019
00020
00021
00022
00023
0002«
00026
00027
00028
00029
00030
00031
00032
00033
00030
00035
00036
00037
OOOjfl
00039
000«0
000«l
00002
00003
oooaa
c
c
c
c
c
c
c
c
c
c
c
c
c
1
c
c
200
(106,IJMSGENT.FLX
THIS UTILITY PROGRAM IS USED TO ENT£R TEXT INTO THE ERROR MESSAGE
FILE USED BY SEKATRA, EACH ERROK MESSAGE IS Two RECORDS (160 BYTES)
LONG, MAOfc UP Of J PAHfS.
(I) FATAL OR WARNING TAG (BYTE 1)
(2) ERROR IDENTIFICATION NUMBER (BYTES 2-«)
(3) MESSAGE TEXT (BYTES 5-t60>
THE oiwtcr ACCESS FILE WILL CURRENTLY HOLD 100 SEPARATE ERRORS.
BYTE «tiFF
1 MAxRtC*200,RECORDS IZE»80,FORMS'FQRMAT TED',
2 ASSOCIATEVARlABLEslII,INITlALSlZ£s25)
( 1 , 1 ) , I s?, 77 } ,
(BUFFU,
WHITE(2'NREC,3)(BUFF(1,2),1*\,60)
GO TO 1UO
CONT INIIE
00 (INO«1«NM5G)
, NRtCa( IND-1 ) *2»-l
. REA()(2lNREC»3)CBUFF(l,l),Ial,Bo)
, NRECsNRECtl
, l'EAO(2'NRECt3) {RUrF{I,e), 1*1 ,60)
, hlRI TE(3,a)8UFF
...FIN
STOP
I FoRM«TC>l , I3.76AI/80A1 )
2 FORMATfA),I3.76A))
3 FORMAT(80A1)
« FORMAT( IX,80AJ/lX,flOAl)
100 CONTINUE
READU , I ,ENDs200)BUFF (1,1),
NHEC»( IND-1 ) *2»1
(FLECS VERSION 22.U6)
113
-------
(FLEC3 VERSION 23.
SB-JAM-82 10118131 PAGE 00001
OOOOL
00002
00003
00004
00005
00006
00007
OOOOfl
00009
00010
0001 1
00012
00013
00014
00015
00016
00017
OOOlfl
00019
00020
00021
00022
00023
00024
0002S
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
0003P
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
SUBROUTINE PAMTICtABAH, ALE*, 8, C, CCIN, COLD, DECAY, OCLTD,
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
1 DFZ, lt J» N8ED, NELEM, UHOUT, C3HIN, QV, 80«BK,
2 30, SH, ALFA, BETA, VEL1» VEL2, VSET,
3 OELZ. DEPTH, 8*10, AHIO, BETA1, BETA2, DSO*B)
THIS SUHROHTINE CALCULATES COEFFICIENTS OF CONVECTIVE, DECAY
AND
SOURCE TERMS FOH TRANSPORT OF POLLUTANT ATTACHED TO SEDIMENT.
INPUT PAHAMETF.RSI
ABAR - AVERAGE ELEMENT AREAS
H • BED CONDITIONS
C - WATER CONDITIONS
CCIN . CONCENTRATION OF INFLOW
COLO - Cf LL-CENTERFO CONCENTRATION
OfCAY . FIRST UP.DER DECAY
OELTO - TIME STEP (DAYS)
DFZ . DIFFUSION COEFFICIENT
OELZ - ELEMENT THICKNESS
I . ELEMENT INPEX
j . PARAMETER INDEX
NBED • NUMWER OF BED LAYERS
NELEM - NU^ER OF ELEMENTS
IJHIN - I*oFLO^ DISCHARGE
UV - VtRTICAL DISCHARGE
SORBK . AOSORHTION ON SEDIMENT, (1-3) ***J/KG, (M-91 1
/UAY
DSORB - OESOHPTIUN FRU"« SEDI^E^T. (1-3) M*»3/KR, (4-9) I/DAY
5H - EROSION RATE, KG
-------
(FLFCS
10116131 PAGE 00002
0005U
00055
00056
00057
00059
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00071
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
10087
00088
00089
00090
00091
00092
00093
0009«
00095
00096
00097
00098
00099
00100
D0101
00102
00103
00104
00105
00106
00107
00108
00109
C
AIJ a iJV(I)
VF.LI »(A>4»VSET(JMJ)»HwIO(I)*ALENl/ABAH(I)
C
A (J • T V (I P 1 )
C
VeL?a(AiJ-VSET(JM3)*ewir>(I-H)*ALEN)/AHAR(
I )
C *** DECAY TERM •** (EFFLUENT)
C
C
»LFAadHUUT(I)/(ABAH(I)*OELZ) + OECAr(l)
C «** SOURCE UR SINK TEBM *** (INFLUENT)
C
HETAl«UHlN(i)/(ABAR(l)*nELZ)*(CCIN(I,J)/3,tCCIN(I*l,J)/6,)
BErA2«aHIN(i)/(ABAW(I)*OELZ1*(CCIN(I,J)/6,+CCKUI*l,J)/3.)
C ****
C
C «
C
C
C
C
C
C
C *
****************************************
A*******************
IfMINUi CCIN IS w«tTTEM IMTO CHAR AS A FIRST
APPROXIMATION TO THE EVENTUAL AVERAGE
CONCENTRATION. THE ONLY MEANS 9Y WHICH
TO ASSURE THE ACCURACY OF THIS
ASSUMPTION IS TO ITERATE TO THE
CORRECT
SOLUTION, AND USE NEW ITERATES TO BETTER
APPROXIMATE CHAR,
C A***************************************************************
C
c *** ABSORPTION / DESORPTION ***
c
DO (IE « l,MXF.LEM>
, DO (1C • 1,MAXCON)
. . CBAKUE.IC) « CCIN(IE,IC)
. ...FIN
C
C
...FIN
IF(CBAR(I,J*3).GT,0.0,ANO.CBAR(I+1,JM3).
I
2
3
j
2
3
1
2
3
1
?
3
AOnSl» SOHHK(J)*SORBK(JM3)/ia.*(
+ CfA»?(I»JM3)*CHAR(I + tj7)fCRAR(I*l/JM3
CBAR(H.l,Jl3)*CRARU + l,7))-SOR8K(J)/6
CBAUdtl, J))
GT.0,0)
3.*C8A«(I,JM3)*CBARd,7)
)*C8AI?(I»7)t
.*(2.»C8ARd,J)»
OSA01» DSOHB(J)*DSOR8(JM3)/12,*(3,*CBAR(I, JM3)*CBAR(I/7)
»CBAH(l>JM3)*CHA»(Itl«7)+CHAH(I + l»JMJ5
CBAW(!-M,JM3)*CHAR(Itl,7))»OSOR8(J)/6
t BAR (I+l»J))
AOi'52* 30r(Hn(J)*30RBK(J13)/l2,*(
+CBAP(I,JMJ)*CH4W(I+1,7)*CHA»(I+1»JM3
*}.*CHAR(I + 1,JM3)»CB43(J»1,7))-S()RHK(
*2. "C^AR ( I* 1 , J ) )
)*CBAR( 1,7)*
.*(2.*CBAR(I,,n +
C3AR(I,JM3)*C»Ai?(l,7)
) *CSAf<( 1,7)
J) /6,* (CUArtd, J)
OSA02" OSn»8(J)*DSOP8CJ«3)/12,*(C"AR( I, JMJ) *C8AR( I, 7)
*CHAR(I»JM3)*CHAU(I+l»7)*C"AH(I+lfJHJ
) *CHAK ( [ . 7 )
+ 3.*CrtAR(I*l,Jr*3)*CBAR(I»l,7))«O.SORLUJ)/6,*(CaAR(I,.))
*2.*C8AR(i+i,J) )
IF(AOOSl,GT.O,O.Oa.AODSl.EJ,riS1iitl^ET
IF(US*LH.LT.O.O)8ET»l»nET»l*oSAl)i
IF(ADOS2,GT,0,O.OW.ADOS2,E'3,PSA02)'iET
A 1«HET 4 1 + 400S1
A2SHE TA?tAOr)$2
115
-------
tFLECS VERSION ?3.(I6) 28-.UN-83 lOliejJl PAGE OOOOJ
OOtlO . IF(OSA02,LT.O.O)BETA2«8ETA;?tDSAD2
OOlil ...FIN
00112 C
0011J C
00110 C *** SCOUR OR DEV031TinN *«*
00115 C
00116 BETA • SR(J) • SO(I,J)
00117 RETURN
00119 END
(FLECS VE«3ION 22.06)
116
-------
(FLECS VERSION 22,
1j.MAR.81
m27i«a PAGE oooui
00001
00002
0000)
00004
00005
00006
00007
OOOOB
00009
00010
00011
00012
00013
00011
00015
00016
00017
00019
00019
00020
00021
00022
0002J
00021
00025
00026
00027
00026
00029
000)0
000)1
000)2
000))
000) | s o.O
KAY? « 0.0
IF (ECHO) WRITEf6,2)
,.FIN
ELSE
*
117
-------
(FLECS VERSION 22.46)
13.MAR.81 I3I27«42 PAGE 00002
00054
00055
00056
00057
00058
00059
00060
OU061
00062
00063
00064
OOObS
00066
00067
0006B
00069
00070
0007 J
00072
0007^
00074
00075
00076
00077
00078
00079
oooeo
00061
000«2
00063
0008U
00085
00066
00087
OOOflft
00069
00090
0009»
00092
00093
0009"
00095
00096
011097
00098
00099
00100
00101
001 0?
0 0 10 S
ool ya
ootos
001 06
00107
00108
00109
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
£
c
c
r
c
c
c
c
c
c
c
c
, ...SECOND DATA SET - ADSORHTION COEFFICIENTS TABU:
THIS IS THE. TABLE OF AOSURTION COEFFICIENTS FOR |8
DIFFERENT WAVELENGTHS THAT ARE A MEASURE OF THE CHEMICAL'S
ABILITY TO ADSORB LIGHT AT THE DIFFERENT WAVELENGTHS,
WAVELENGTH UNITS ARE NANO METERS.
... .RECORD 1
COL, 1-10. ,.E( 1J...COEF, FOR WAVELENGTH OF 300.00
1I-?«,..E( 2). ,.CO£F. FOR WAVELENGTH OF JO}, 75
21-30, ..EC 3),,,COEF, FOR WAVELENGTH OF 306.75
31-40. ..EC i).,,COEF. FOR WAVELENGTH OF 313,75
01-50, ..£( 5). ..COEF, FOR WAVELENGTH OF 318.75
51-60, ..£( 6),,.COEF. FOR WAVELENGTH OF 323.10
61-70. ..£( 7),.,CUEF. FOR WAVELENGTH OF 346,00
71-«0.,,E( 8),,,COEF. FUR WAVELENGTH OF 370.00
....RECORD 2
COL. 1-10. ,.£( 9)...COEF, FOR WAVELENGTH OF 400.00
1 1-20. ,,f( 10) .. .COEF. FOR WAVELENGTH OF «30,00
21-30. ,.E( 11) .. .COEF. FOR WAVELENGTH OF «60,00
31-40. ,,£(12).., COEF. FOR WAVELENGTH OF 490,00
"1-50, ..£(13), ..COEF. FOR WAVELENGTH OF 536,25
51-60. ,,£(!«). ..COEF. FOR WAVELENGTH OF 537.50
61-70, .,£(15) ,, .COEF, FOR WAVELENGTH OF 637.50
71-80.. .£( 16) ,,,C06F. FOR WAVELENGTH OF 687.50
....RECORD 3
tUL, 1-10, . ,E(17> ...COEF, FOR WAVELENGTH OF 756,00
11-20, ,.E(16),,, COEF. FOR WAVELENGTH OF 800,00
REAOil,4) (E(L),Lel,|8j
....THIRD DATA SET - SoLAR INTENSITY TABLE
THIS TABLE CONSISTS OF FOUR SETS OF is VALUES. THE
FOUR SETS CORRESPOND TO SPRING, SUMMER, FALL* AND WINTER,
RESPECTIVELY, THE l *J VALUES CORRESPOND TO THE 18 WAVELENGTHS
AS DESCRIBED ABOVE IN TH£ ADSORPTION COEFFICIENT TABLE,
THE INCLUSIVE DATES FOR EACH SEASON ARE GIVEN BELOWI
CALENDER DATES JULIAN DATES
SPRING MARCH 1 . MAY 31 60-151
SUMMER JUNE 1 - AUG. Jl 152-243
FALL SEP. 1 - NOv, JO 244. JS4
"INTER U€C. 1 - FEB. 28 ii5"365> 1-59
HE AD (1,4) ((SI(L,I),Lsl,18),Isl,4)
IF (fcCH(j)
N K I T E ( 6 f 3 )
«HITE(6»7) JULIAN,PHI,KAn.KAY2
«HI 1E(6,5)
HO (Let, 1«)
, wfiiTE<6,&) WL(L),ECL),(SI(L/D, I = l,«)
...MN
..FIN
*»* THE JULIAN DATE is ADJUSTED TO MAKE THE FIRST DAY OF
118
-------
(FLECS VERSION 22,06)
ij.M4H.et i3i27i'J2 PAGE 00003
00110
001 II
001 12
00113
001 1 U
00115
001 16
001 1 1
00118
00119
00120
00121
00122
00123
0012V
00125
00126
00127
00128
00129
00130
00131
00132
00133
0013"
00135
00136
00137
00138
00139
SPRING JULIAN DAY 1 ***
(JULIAN ,GE. 60) JULIAN
ELSE JULIAN « JULIAN + 30«>
JULIAN - 59
FIRST TERM OF THE RATE.
THE FOUR SEASONS ***
OF CHANGE EQUATION
*»* COMPUTE THE
FOK EACH OF
00 (Isl.U)
. PCOEF(I) » 0.0
. 00 (Lst,ia> PCOEF(I) « PCOEF(I) * E(L) * Sl(L>I)
PC,Otf(I) a PHI * PCOEF(I) * (2. 303/6,02E2U) * SECDAY
...FIN
..FIN
RETURN
FORMAT(15,3FIO.O)
FURMATfIHu,13X,'NO PHOTOLYSIS DEGRADATION WJLL BE COMPUTED')
FORMAI{iHO»i8Xi'PHOTOLYSIS TABLES AND COEFFICIENTS')
FORMAT(PF10tO/3FlU,Q/2FlO.O)
FORMATdHU,56X,'PHOTOLYSIS TABLES•/39x»'ADSQRBTION',
1 23x,'SOLAR INTENSITIES'/23x, H.AMDA CENTER',
2 2X,'COtFFlCIENTS',5X,'SPRlNG',8X,iSUMMERS9X,'FALL',
3 9X,'WINTERI/23X,6(12(1-'),2X1)
FORMATf26X,F6,2,3X,5(2X,lPEl2,«))
..JULIAN STARTING DATE'/
.REACTION QUANTUM YIELD'/
.LIGHT EXTINCTION COEFFICIENT OF HATER'/
.LIGHT EXTINCTION COEFFICIENT OF SUSPENDED',
FORMA] (26X, 15,
1 |«X, 1PE12.5, >.
2 1UX, 1PE12.5, '.
3 Mx,)PEI2,5, ',
t>
SOLIDS IN WATER')
END
(FLECS VERSION 22,16)
119
-------
(FLECS VEHSION ?2.U6)
2fl.jAN.92 10H7I41 PAGF. 00001
00001
00002
00003
00004
00005
00006
00007
00008
00009
U0010
00011
00012
ooois
ooo t a
00015
00016
00017
00010
00019
00020
00021
00022
00023
00024
00025
00026
10027
00028
00029
00030
00031
OOOJ2
00033
00034
00035
00036
00037
00038
00039
oooao
oooai
00042
00043
00044
00005
00006
00047
oooae
000tt9
00050
00051
00052
00053
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
SUBROUTINE PRQFILCALEN, AWID, DE.LZ, DEPTH, NELEM, Q,
1 USTAR, VOL* OH, EL, IFLAG, DELTA)
THIS ROUTIME ASSIGNS A LOGARITHMIC OR UNIFORM PROFILE FOR THE
BULK VOLUMETRIC FLO*S.
INPUT PARAMtTERSi
ALEN SEGMENT LENGTH
AHIO SEGMtNT rtlUTH
DELZ STANDARD ELEHENT THICKNESS
OEPTH FLO* DEPTH
NELErf NUMHER OF VERTICAL ELEMENTS IN THE SEGMENT
'i FLOW TO BE OIST»IHUTEO
USTAR SHEAR VELOCITY
VOL SEGMENT VOLUME
OUTPUT PARAMETERS!
Urt . DISTRIBUTED FLOW
CALLED BYI HYDFLO, ICFLO
INCLUDE "ELMSiz.pRM"
DIMENSION HH(MXELEH), Z(«XEL£M), Zl(MXELEM), A^IO(rtXELEM) ,
1 tL(MXELEH)
DATA XK/O.U/
DATA G/9.S1/
CO (I • l.MXELE*!)
. OH(I) • 0.0
...FIN
uBA««fj / VOL * ALEN
Zn»OEPTH/(10.**(URAR*XK/(2.3*USTAR) t 1./2.3))
ICOUNTaO
Z0«0.001
Cl«uaAR*XK/USTARtl.O
C ON«OEPTH*(C 1-2. 303 *ALOG10 (OEPTH) )
REPEAT UNTIL (»BS( EPS ).LT,O.OI. OR. ICOUNT.GT.IO)
ICfJUNT»ICOUNT+l
FZO«CON-ZO«Cl+2,303*nEPTH*ALOG10(ZO)
FPZO«DEPTH/ZU"C)
ZP«ZO«FZO/FPZO
EP5"(*P«ZO)//P
;.C»ZP
IF(ZO.LT.O.O)
. ZO»i)EPTM/(lo.«*(U8AR»XK/(2.30i*USTA9)+l./2.303) )
, ICOUNT«tl
...FIM
..FIN
WHEN (ZO ,GT. DELZ/'l,O.OR.NELEM,EQ.l)
*
120
-------
LtCij VtKSJON 2
00051 C
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068 C
00069
00070
00071
00072
00073 C
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
OOU9J C
0009« C
00095
00096
00097 C
00098 C
00099 C
00100 C
00101 C
00102 C
00103 C
00104 C
00105 C
00106 C
00107
00108
00109
>,46) ?«-.)AlJ.S£ 10U7l«l PAiJE 00002
*«* DISTRl^UTF VELOCITY UNIFORMLY ***
XH£N( if\. AG.EU.O)
, DO (J«l ,NELEM)
. . >JH(I)silHAK*nELZ*A*iIO( I )
. . . *f IN
. . «F I"1
ELSE
, 00 ( 1*1 , ^ELE*- 1 )
. . ;JH ( I } «UH *** (EL ( I + 1 ) "EL ( 1 ) ) *Ax Jo ( 1 )
. . . .Flfo
, ljlH(N£LE1>O»IJ»AH*OfLT**AWJD(NELEM)
...FIN
..FIN
LSE
*** niSTtfllUTE VEL'ICITV INTO A LOGARITHMIC PROFILE ***
MUM»»'JELEhtl
Z( 1 )*Z'>
Z(NtJM8)BO£MTH
Zl (N|JMH)«0.
JX a 1
JaNUMB
REPEAT UNTIL (JX ,FU. NgLEH ,0«. Z(J) ,LE, ZO)
. Jt m JX + t
, JsMELE* • JX * 2
. *HEN (IFLAG.EU.O) zi (.n«zi t J+D+OELZ
. ELSE
. . "HEN (J.EU.NELEM) Z1(J)»Z1(J+1)+OELTA
. . ELSE zi(j)»zifjti)*£L(j+i)-eL(j)
. ...FIN
. Z(J)«OE»TH . ZHJ)
. . . F I rt
"
-------
(FltCS VERSION
00110
00111
00112
00113
ooim c
00115
00116
00117
00118
00119 f
00120 f
. T2«
. 9H(
» SUM
...FIN
DO (I
•
• t *
..FIN
RETURN
INO
28-JAN.H2 10il7l'41 PAGE 00003
T2«A*(Z(IPl)*ALOG10(7(lPl))-Z(n*ALOG10(Z(I)n
UH(I)BIJH(I)/SUM * Q
VERSION tz.nt,i
122
-------
(FLECS VEHSION 22.46) IJ.MAH-ei 13«28ll<4 PAGE OQ001
00001 SUBROUTINE PUTEKR(If), NUMERR,
00002 C
00003 C WHEN AN EKROR IS DETECTED IN THE INPUT STREAM ,
OOOOU C THIS SUBROUTINE IS CALLED TO PLACE THE E»ROR IDENTIFICATION
00005 C CODE (ID) INTO THE HOLDING ARRAY (HLDE«H) AND INCREMENTS
00006 C TH£ NUM8EK UF EKRORS
-------
(FIECS VERSION 22,
19.MAR-81 13141102 PACE 00001
ooooi
00002
00003
00004
00005
00006
00007
00008
00009
00010
0001 1
00012
00013
oooi u
OOOI 5
00016
00017
OOOlfl
00019
00020
00021
00022
00023
00024
00025
00026
00027
0002R
00029
00030
00031
00032
00033
00034
00035
0 0 0 S 6
OU037
00038
00039
OOOuO
OOOUl
00042
00043
00044
00045
00046
00047
C
C
C
C
C
C
c
c
c
c
c
c
c
c
c
c
I
10
20
SUBROUTINE RADIUS ( ALEN, AREA, CROSEC, DEPTH, EL, HRAD)
THIS ROUTINE CALCULATES THE HYDRAULIC RADIUS OF A CROSS. SECT ION.
INPUT PARAMETERSI
ALIN . SEGMENT LENGTH
AREA - SURFACE AREA AT NODAL DEPTHS
DEPTH . DEPTH QF CROSS-SECTION
EL . ELEVATION OF NODAL AREAS
OUTPUT PARAMETERS!
CRUStC - TOTAL CROSS-SECTIONAL AREA
HRAD - HYDRAULIC RADIUS
CALLED bYt HYODAT, KFLO
INCLUDE 'ELMSIZ.PR*'
DIMENSION AREA(MXELEM) , EL(MXELE^)
CROSEC » U,
BLEN x 1,/ALEN
WETPER s AR£A(U*f»LEN
ELBTM = EL( 1)
WBTM s AREA ( 1 ) «bLEN
DO ERRO« IN SUBROUTINE RADIUS'/
1 "0 CROSEC s ', E12,«,'M**3'/
2 '0 wETTF.O PERIMETER a i, E 1 2, 1 , ' M * *2 ' )
GO TO 20
MOP » WHTM f (wrrjp - WHTM) * (DEPTH - ELBTM)/ (ELTOP - ELBTM)
ELTOP * DEPTH
CROSEC * CROSEC » (EL^OP - ELBTM) t (WTuP f WBTM)/2,
WtTPER « SORT ((ELTOP - ELBTM)**2 » «WTOP -«BTM) /3 , ) * *2) * 2,0
t t WETPER
HR*D = CROSEC / WETPFR
RETURN
END
(FLECS VERSION 22,
124
-------
(FLECS VERSION 22.46)
6.APR.
17121109 PAGE 00001
00001
00002
00003
00004
00005
00006
00007
OOOOfl
00009
00010
0001 1
00012
0001 3
00014
00015
00016
0001 7
0001P
00019
00020
00021
00022
00023
00024
00025
00026
00027
00026
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
oooao
00041
00042
00043
00044
OOOU5
OOOU6
00047
00048
00050
00051
00052
0005S
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
SUBROUTINE RDSFLOIALEN, AREA, AWID,C, OELZ, OFZ, NELEM, NELEMP,
1 PDELZ, P"1D, PXSAR, VSET, XSAREA)
THIS SUBROUTINE IS CALLED EACH TIME STEP (EXCEPT THE FIRST) WHEN
DEPTH WITHIN THE S>E«MENT HAS CHANGED. ITS TASK is To REDISTRIBUTE
THE CONCENTRATIONS.
FORMAL PAHAMETERSI
A*IO • "IOTH P CURRENT TIME STEP
C • THE NELMpTti NODAL CONCENTRATIONS THAT ARC TO BE
REDISTRIBUTED
THE
OELZ • THE STANDARD ELEMENT THICKNESS FOR THE CURRENT TIMf ST£P
DFZ • DIFFUSION. DISPERSION COEFFICIENT
NELEM • THE NUMBER OF ELEMENTS FOR THE CURRENT TIME STEP
NELEMP - THE NUMBER OF ELEMENTS DURING THE PREVIOUS TIME STEP
PDELZ • THE STANDARD ELEM£NT THlCKNESS USED DURING THE PREVIOUS
TIME STEP
PXSAR • PRESS SECTION f PREVIOUS TIME STEP
PWIO • WIDTH » PREVIOUS TIME STEP
XSAREA • CROSS SECTION # CURRENT TIME STEP
VSET " SETTLING VELOCITY OF SEDIMENT
CALLED BVI SERATRA
INCLUDE 'ELMSIZ.PRM'
DIMENSION AREA(MxELEM), AwIO (MXEL£M) , C (MXEL£Mf MAxCON) ,
KP(MXlLtM,MAXCON), OFZ((f), PW ID ( MX ELEM ) , P XSAR ( MXELEM) ,
2VSEH3), XSARE*(MxELEM), JELP (MxELEM) , HEQXS (MXELEM) „ TMASS(7)
PXS « 0,
XS * 0,
00 (Isj, NELEM) XS » XS t XSAHEA(I)
00 U*i«NtLEMP) PXS 3 PXS » PxSAR(I)
DO ( ISI,NCLEMP+D
. DO (J«1,MAXCON)
, . CP(1,J) - C(I,J)
, , , ,F IN
...FIN
RATIO * PXS/XS
CALL t(iupxs(pxSAR, PWIO, POELZ, XSAREA, NELEM, NELEMP, RATIO,
HELP, HtUXS)
00 (Ksl,MAxCON)
T"ASS(MSO,
CT«CP(2,K)
CHsCPt 1 .K)
NELMBTsI
CMAS^QeO ,
DO ( i * i , NELEM i
. N£LMTP*lELP(I)
, NbsNELMTP
125
-------
(FLECS VERSION 32,06)
6-APn-oi
ooosa
00055
00056
00057
00058
00059
00060
00061
00062
0006J
0 0 0 fc »
00065
00066
00067
00068
00069
00070
00071
00072
0007}
00071
00075
00076
0 0 (1 7 7 C
00078
000/9
00080 C
00081 C
0003? C
OOOH J
ooonu
OOOP5
OOOflfr
00087
00088 C
00089 C
00090 C
00091
0009?
0009J
00094
00095
00096 C
00097 C
00098 C
00099
001 00
00101
00102 C
0010J C
ooioa c
00105
00106
00107
0010M
00109
. . N 1 * N (U 1
. . F 1 5N8«PDELZ
, , FH=t!«POELZ
, , C I *CP ( NT , K )
, . CBsCP (NB,K )
. , HfcL = HEf)XS( I)
, , f AC 1 » (HEL-EH) /POELZ
, , F AC?,s ( E T »HEL > /POELZ
, , c TOP=C T *FAC t « CP»F AC
-------
(FLECS VERSION 22.«6)
00110 C
00111
00112
0011}
00114
00115
00116 .,FIN
00117 RETURN
00118 END
6-APH-ei I7l21|09 PAGE 00003
NUOAl, VALUES ABOVE BOTTOM ELEMENT
ELSE CfI*l,K) * 2,*CMA3S'XS»REA(I).CU,K)
(FLECS VERSION 21,Ub)
127
-------
(FlECb VERSION 22.U6)
Onoyl
00001
00002
00003
QQ004
00005
00006
00007
OOOOfl
00009
o o o i n
0001 I
00012
0001 3
00010
00015
oooi*
0001 7
00018
000 19
00020
00021
00022
00023
00021
00025
00026
00027
0002S
00029
00010
00031
QOOS2
0003$
0003U
00015
00016
00017
OOOjfl
00019
00040
C
C
C
c
c
c
c
c
c
c
c
c
c
SUBROUTINE tfPTE*R(NUMERH, HLOERR, FERKOfO
THIS KouTiNt is RESPONSIBLE POR REPORTING ANY INPUT ERRORS
UNCOVEHEO BY THE INPUT HOUTIN£S AND DETERMINING THEIR SEVERITY,
CALLED bYt SEHATRA, STRTUP
BYTE HLDEKROOO) ,BiiFF(so,2)
LOGICAL*! FERHOH
CALLED BY: SERATRA
OPENfUNl T«)0, NAME »<5EH*TE'"'.^S6I, TYPE"' OLD',ACCESS='PI"ECT',
1 FONMt'FOHMATTEO' , ^AXREC«200 , RECORDS UEsftO ,
2 ASSOCIATE VARIABLES 1 1 1, READONLY)
NUHflsO
NUMfsO
WHI T£<6. a)
DO ( 1 s 1 i NU^ERR)
lNOs(HLOtRR( l)-l )*2»1
KEAD( 1 0' IND, 1 ) (HIJFFCK, 1 ),Kel ,80)
IN|")slNL>+l
WEAnilU'lNl>,l)(RUFF(K,2)fK3l,80)
IFfRUFF ( 1 , 1 ) ,EQ, 'W ' )NUMW»N(JMH+t
IF(«UFF ( 1 , 1 ) ,EQ, 'F ' )NUMF=NUMF*1
wRITt(0,2)((BUFF(K,J),K*l,80),J»l,2)
..FIN
WKITE(6,3)NUMW,NUMF
CLOSE(UNIT«10)
If (HWf .ST. 0) FERROR » ,T"UE.
RETURN
1 FURMA1(80A1)
2 FQRMAT(lHg.Al,2X,lA\,Jx,7bA|/loX,BOAl)
3 FORM*T
-------
(FLECS VfRSION
lO-APR-81 I3I5H37 PAGE 00001
oonoi
00002
00003
00000
00005
00006
00007
00006
00009
00010
0001 1
00012
00013
0001"
00015
00016
00017
00018
00019
00020
00021
00022
00023
00020
00025
00026
00027
00028
00029
00030
OOOJI
00032
00033
00031
00035
00036
oooj;
00039
00039
00000
00041
00002
00043
00040
OOooS
00006
00007
00008
00009
00050
00051
00052
00053
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
SUBROUTINE SANDUBAR, ALEN, AREA, 8, BOIV, CCIN, OELTO, Of.L
1 DENS, D50, HRAD, N0ED, NELEM, POR, QHIN,
2 GHOUT, SCSHR, SLOPE, SM£TM, STRESS, T£MPR,
3 VSET, vou xrsn, OEPQ, ILAYR, su, SR,
« XNT, COLO, C, CROSEC, BwIO, ECH07, SCOUR)
THIS SUBROUTINE COMPUTES THE SOURCE/SINK TERMS REQUIRED FOR
SCOUH/PEPOSITION OF SAND. TRANSPORT CAPACITY IS CALCULATED
ONCE PF> SEGMENT,
INPUT PARAMETERS:
ABA« . AVERAGE VERTICAL PROJECTION AREA
ALEN . SEGMENT LENGTH
ARtA - ELEMENT (REAL) VERTICAL PROJECTION AREA » DATA
B • 6tO CONDITIONS
BDIV - STANDARD REO L*YER THICKNESS
BWJD . Rt*L HIOTH AT CROSS-SECTION BREAK POINTS
C . WATER CONDITIONS
CC1N . CONCENTRATION OF INFLOW
CROSfcC - TOTAL CROSS-SECTIONAL AREA, M*»2
OEL'L) - TIME STEP (DAIS)
DELZ - STANDARD ELEMENT THICKNESS
DfLZT - THICKNESS OF THE TOP ELEMENT
OENS . DtNgiiY
050 . MEDIAN BED SEDIMENT DIAMETER (METER)
HRAD - HYDRAULIC RADIUS
NBED . NUMBER OF BED LAYERS
NELEM - NUMBER OF ELEMENTS
POH - POROSITY
OHIN . INFLOW DISCHARGE
OHOUI « OUTFLOW DISCHARGE
SCSHH . CRITICAL SHEAR STRESS FOR SCOUR
SI UPE - ENEHfiY OR RIV£R 8EO SLOPE
SMETH . METHOD TO HE USED WHEN COMPUTING SAND CAPACITY
sTl TOFFALETTI'3 METHOD
sc» COLBY'S METHOD
STRESS - BED SHEAR STRESS
TEMPR . WATER TEMPERATURE
VSET . PARTICLE SETTLING VELOCI^y
VOL - VOLUME
XYJO . THICKNESS "r TOP BED LAYER
OUTFUI PARAMETERS!
ILAYR . NO, OF BED LAYERS AFFECTED BY SED, DEPOSITION
SO - DEPOSITION RATE, (KGtPC)/M**3/r>AY)
SR - EROSloN f»ATE» (KG (PC! /M**3/UAy)
xNt . WEIGHT OF TOP BED SEDIMENT LAYER, (KG/M**a)
OEPO • BED DEPOSITION RATE ( KG (PC J /M2/DA Y )
SCOUR - BED SCOUR RATE (KG (PC 1 /Ma/pA Y )
CALLED BY: TRANSP,
CALLS: fOFML, COLBY
INCLUDE 'ELMSIZ,PRM<
l,
NOOES
(BYTE
AND ER
129
-------
(FLECS VERSION
0005«
00055
00056
00057
00058
00059
00060
00061
00062
00063
0006«
00065
00066
00067
00069
00069
00070
0007 1
00072
00073
0007«
00075
00076
00077
00078
00079
00080
OOOS1
00082
0006 J
0009"
00095
00096
00087
OOOBS
00089
00090
00091
00092
00093
0009U
00095
00096
00097
00098
00099
00100
00101
00102
'00103
001 0«
001 05
00106
00107
00)06
00109
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
22. «6)
BYTE !
LOGIC/
D1HEN!
1
2
3
u
5
i N r T 1 1
ILAYH
RSsO.i
CSsO.l
OEPO (
DEPOO
SHC s
G5I *
QTOT '
SR(1 ):
SR(« ) :
scou«
SCOUR
RATE '•
VPCROi
VOLU^I
DO (K
. 80
, SO
...Ml
XNTO
IF (NP1
1 ,
...M
IF (E
,
. CA
t
9
•
ff
f
. oo
* •
1. .
• •
• 1 ff
§
. IT
, AT
t
1Q.APR-81 13I5UJ7 PAGE 00002
LOGICAL*! FERROR, ECHO?
DIMENSION ABAR(MXELfM), AREA(MXELEM), 0(MAXLEV,MAXCON-l),
COLD(MXELEM,MAXCON),CC1N(MXELEM,MAXCON), OEMS(3),
OHIN(MXELKM), OHOUT(HXELEM), SCSHRO),
6), SR(6), VSET(3)» XNT(3).
C(MXELEM,MAXCON),
OEPO(6), SCOUP(6)
INITIALIZE SCALAR AND ARRAY VARIABLES
« o
VPCROS * ALEN * ewiod)
CROSEC * ALEN
NELEM)
80(K. 1 )*0,0
SD(K,4)30,0
0.
Nt.O)
)»(! ,0-POH)/(8(N8£0,1)/OENS(1)•
(ECH07)
CALCULATE ACTUAL SANO IRANSPOHT WITHIN THE RIVER
SRC (KG/DAY)
NOTEi cci'J is A TI^E AVERAGED QUANTITY, COLO is NOT, AN
ITERATIVE LOOP IS CALLED FOR ''HERE SRC IS UPDATED UNTIL
RESULTS ARE UNCHANGED.
00 fix*!,NELEM)
3MC*SHC*(JHIN(IX)*(CCIN
-------
(FLECS VERSION 22.
to-ApR-81 msiiJT PAGE oooos
ooj jo
00)11
00112
001)1
001 1U
00115
00|16
001 17
00118
00119
00120
00121
00122
00123
00121
00125
00126
00127
00)28
0012"»
00130
00)31
00132
00133
0 0 1 3 'i
00)35
00116
00137
0013"
001 39
001 10
/% n t it f
o o i y i
001«2
00113
001 14
00105
00116
00|U7
00 1 18
00119
00150
00151
00152
00153
00151
00155
00)56
00157
001 5fi
00159
00) 60
00161
00162
00163
00161
00)65
C
C
C
c
c
c
1
1
1
c
c
c
c
e
c
50
1
1
1
f
,
\
<
1
1
«
1
c
c
c
100 ,
c
c
c
c
c
c .
150
•
t
I
, c
, s
1 •
k •
I •
* •
i «
k •
k •
k •
k t
k t
k t
4
k •
i •
\ «
•
t •
k t
G
0
0
I
s
0
.
I
D
*
*
•
R
C
n
c
D
I
I
I
CALCULATIONS OF STREAM CAPACITY FOR 3ANO
SELECT (SMETH)
( • T « )
* THE TOFFALETI TECHNIQUE MAr NOT BE w£LU SUITED FOR CHANNEL
. CROSS-SECTIONS WHICH DIFFER MARKEDLY FROM RECTANGULAR,
'. CALL TOFFALfALEN, 050, GSI, HRAD, QTOT, SLOPE, T£MpR, VOL, VSET,
GSU, GSM, GSL, GSB, YU, YM, YL)
{ 'C' )
. CALL COLBY(ALEN, C, DELZ, 050, HHAD, NELtM, OTQT, TEMPR,
. VOL, SSI, FERROR)
. IF(FERROH)
. , CALL TOFFAL(*LEN, D50, GSI, HRAD, QTOT, SLOPE, TEHPR, VOL, VSET,
. GSU, GSM, G5L, 558, YU, YM, YL)
. , FERHOH s .FALSE.
. ...FIN
...FIN
..FIN
GSI » liSI » 8wiO(NELEM t 1)
OETF:HHINE JF DEPOSITION, SCOUR, OR NEITHER OCCURS,
f G3I • 3PC
IF (OIF) 50, 100, 150
SAND HITHIN THE WATER COLUMN EXCEEDS CAPACITY . DEPOSITION OCCURS
OEPO(l) « -OIF/ AREA(l)
RA1E » -DIF/VOLUME
hf t) x -1
DO fK*l,NEL£M)
3D(K,1 ) s RATE
VULK a A8AR(K)*DELZ
StOsQHIN(K)»(cCIN(K,l)fCCIN(K*l, 1 ) ) /« . t CJHOUT (K ) »COLD ( K , l>/2.
CONIsQHIN(K),(CCIN(K,a)»CCIN(K*l,«))/q.tQHOUT(K)*COLD(K,a)/2,
RATtK » RAT£ * VOLK
SU(K,a) * RATEK * CONT/SED/VOLK
DtPO(i() « DEPO(
-------
VERSION 22,«6)
io-APn-81 instil? PAGE oooo<»
00166
00167
00168
00169
001 70
001 71
00172
00173
00171
00175
001 76
00177
001 78
001 7<>
00180
00181
00162
OOtBJ
0018'4
001B5
00)86
00187
001*9
001B9
00190
001V1
00192
00193
Q019«
00195
00 (K«I,NBEO)
NB * NBEO - K » 1
luENS s (1,0 . POR) / (B(NB, 1 )/DENS(l ) + B (NBf ?. ) /OEN8 (Z J
» 8(NB,3)/OEN3(3))
DtL s BOIV
IF (NB.EO.NBED) DEL « XYSO
TJNBED « TOENS * DEL » B(NB,1) * VPCROS
WHEN (TRJOSP.GI: .TINBED)
. ftS « RS «• tINBEO
, C3 « CS + TIN8EO * 8(NB,a)
. f«SUSP « TPSUSP • UN8ED
. IINBEO » 0,0
GO TO 175
ELSE
, ft$
. C3
» RS t TPSUSP
» CS * TRSUSP * 8(NB,«I)
HNBEO « UNUEO - TRSUSP
1 GO TO 175
...FIN
..FIN
175 SCOUH(l) » »3 / DELTD/ AREAC1)
SCOIJR(t) * CS / OElTD/ AREA(l)
XNT(l) « TINBEO / VPCR03
SR(1) » RS /OELTO / VOLUME
SRt«) * cs / DELTO / VOLUME
..MM
200 RtTURN
END
(FLECS VEHS10N 22.U6)
132
-------
(FLECS VERSION 22,ut.)
2fe.MAH.81 12l27|03 PAGE 00001
00001
00002
00003
00004
00005
00006
00007
00009
00009
oonio
0001 1
00012
000)3
oooiu
00015
00016
0001 7
00018
00019
00020
00021
00022
0002J
00021
00025
00026
00027
00029
00029
00030
000}1
OOOJ2
OOOJ3
00034
OOOJ5
00016
OOOJ7
OOOJfl
00039
00010
00041
00092
0004)
00044
00045
00046
00047
00048
00049
00050
00051
00052
0005)
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
SUBROUTINE SAVEIT(R, BDIV, BED, ELEV, C, OELZ, NBED,
1 NtLEM, NX£<5, R£sELN» STRESS, XYSO,OLDC,
2 ALEN,QHIN,QHOUT,CCIN,QV,AwiD,BWlD,VSET,OENs,D£LTO,DFi,
i PUR, TBEO)
THIS ROUTINE WHITES THE SIMULATION RESUUS TO THE RESULT FILE
(IOH 5) THAI HAS BEEN OPENED BY SERATHA
FORMAL PARAMETERS:
8 BED CONCENTRATIONS
BUIV ST*NOAH(> BED LAYER THICKNESS
BED BED THICKNESS
C WATER CONCENTRATIONS
OELZ STANDARD ELEMENT THICKNESS
ELEV SEGMENT ELEVATION (DATUM ELEVATION)
NBED NUMBER Of BED L*YERS
NELtM NUM9ER OF ELCMENTS
NXEO CURRENT TIME STEP (I*'i)
REStLN WATER SURFACE ELEVATION
STRESb SHEAR STHES3
XYSO THICKNESS OF THE TOP BED L**ER
CALLED 8Y| SER*TRA.
INCLUDE 'ELMslZ.PRMi
INTEGERS NXEQ
DIMENSION 8(MAXLEV,MAXCON-1 j ,8AVG(MAXLEV) , BEL(MAXLEV),
1 C(MXELEM,MAXCON), CMASS(MXELEM), CVOLM(MXELEM),
2 CTOTL(MXELEM), wELEV (MXELEM) , D (MXELEM, MAXCON)
3 ,QHIN(MXtLEM),OHOUT fMXtLEM),CCIN (MXELEM, MAXC UN),
« UV (MXELEM), AW 10 (MXELEM) , tJi*ID (MXELEM) , VSE t ( 3) ,
5 OE NS (3J,'JClN(MxELEM,MAxC ON ),QCUlJT (MXELEM, MAXCON),
6 CEL*V(MXEL£M, MAXCOff),OVC£L( MXELEM, MAXCON),
7 ULDC(MXELEM,M*xCON), DF?(4), OVD I F (MXELEn, MAXCON) ,
8 0 TBED ( MAXCON ), TREO ( MA XCON),BAL( MAXCON), OIF (MAXCON),
9 OCELAV(MAXCON), TO f 02 F {MAXCON )
DATA EPS1/1 .OE-30/
NELMPlzNELEM»i
NBEDP1 * NHED » 1
*** WATER CONCENTRATIONS ***
J * NELtM » l
WRITE(5) NX En,J,Nbeu,ELEV,OtLZr BO IV, XYSO, STRESS
REPEAT UNTIL (J .EH. U)
CMASb(J) s C(J,«) + C(J,S) » C(J,6j
SUM = C(J, 1 ) » C(J,2) » C( J,3)
WHEN (SUM ,GT. 0.0) CVOLM(J) s CMASS(J) / 3IJM
ELSE CVOLM(J) s 0,0
CTOTL(J) a CMASS(J) + C(J,7)
SELECT (J)
133
-------
(FLEC3 VERSION 2?.06)
iai27:oj PAGE 00002
0005U
00055
00056
00057
00058
00059
00060
00061
, , (NELEM+l) WELE.V (NELEMt 1 ) s HES£LN
. . (UMtRWISfc) WELEV(J) a «ELEV(J+1) • UtL2
,
, ,r IN
, J 3 J - 1
...F
IN
DO(KTl , HAXCON)
. ood*if NEUEM+I )
. . WH£N (K .LE. J ,OH. K ,EO, 7} D(I,K)SCd,K)
00062 . . ELSE
0006}
0006«
00065
00066
00067
00068
00069
00070
00071 C
00072 C
0007} C
0007«
00075
00076
00077
00078
00079
00080
00081
00082
0008J
0008U
00085
00086 C
00087
00088
00089
00090
00091
00092
0009 J
00090
00095
00096
00097
00098
00099
00100
001 01
001 0?
0010 J
001 0«
OU1 05
001 06
00107
00108
00109
f
•
, 0(I,K)=0.
. IF (CU,X-S) ,GT, l.OE-10) 0(I,K)sc(I,K)/c(I,K.JJ
FIN
,
. , ,F
UN I
1
wttl
* * *
..FIN
IN
Et5) (WELEV(J), (C(JfK)rK*i,MAXCON),cMASS(J),CVOLk'(JJf
CTOTL(J)>J*NELEMtl,l,«l)
E(5) dD(J,K),KsU,6),J«NELEM»l,l,«l)
BED CONCENTHAT IONS ***
BELEV s (NBED-1) * BOIV » XYSO + ELEV
J s
RtPl
NHEO
.AT UNTIL (J .E'J. 0)
BAvO(j) * B(J,l)*fl(j,«) t 8(J,2)*B(J,5) + 8 ( J, } ) *R (J , 6)
SELEC1 (J)
(NH&D) BEL(NHto) a BELEV
, (N9ED-1) SEL(NBEO-I) x BELEv - XYSO
fUTH£RWlSE) BEL(J) » 8EL(J»l) - BDIV
...FIN
J * J-l
..FIN
WW J
fE(5) (BEU(J)»(B(J,K),K«l,MAXCON.i),t)AVG(J),JsNHEO/l,-l
*** ELEMENT MASS »N0 rONvECTEO MASS ***
00 (J»1,MAXCON)
. OC£l_AV(J) s 0,
• '
,EL A v ( NELMH 1 i J ) = 0
;C I N (NtLMP 1 • J ) sO
. WCOUI (NELMP1, J)«o
. 00(1 al,N£LE«)
f
.
.
•
.
.
.
.
.
.
.
•
.
.
.
.
.
vULs*wlO(l)*riE.LZ*ALEN
XS=AW jo( I ) *OELZ
VPX&*AwiO( I ) * ALEN
ULMtAN s (OLOCd,J) + OLOC ( I + 1 , J ) ) /2 ,
CMEANa(CdiJ)tC(I + l»J))/2,
CEL*^(IfJ)»VOU*CMEAfj
'3CIN(I,J)aQHIn(l)»(CClN(IfJ)+CCIN(I»l,J))/2t*OELTD
WCOUT(T(J)=QH()ilT(i)»OELTf)*(CMEANt(OLfC(IfJ)»OLOC(I*l
CtLAV(NfcLMPl,J)sCEUAV(lMtLMP|,J)+CELAV(l,J)
ULlN(NELMPl»J)srJClN(NELMPl,J)+QClN(l»J)
f)CUUl(NELMPl,J)BQcOUT(NELHPt,J)fOcOUr(l,J)
ULEIAV(J) r OCELAV(J) » OCM£4N*VOL
K Ij
IMJ ,GT, 3) KsJ.J
WHEN (K ,£Q, H 1 w$30.
ELSE wgeVSE! (K)
w"EN(I ,EU, t) QVQBTMrQV ( 1 ) • (C f 1 i J) ^OU^C ( 1 , J ) ) /2,
134
-------
(FLECS VEKSION 22.46) 26-MAH.8] 12127103 PAGE OflOOS
001 10
00111
00112
OOt 1 J
ooi ia
001 15
00116
001 17
00118
00119
00120
00121
00122
00123
0012"
00125
00)26
00127
00)28
00129
001 30
00131
00132
00133
OOj 34
00135
00] 36
00137
00138
00139
001 UO
ooi/2,
. WHEN (I ,EO.
, . . ClVCTOP»OV(NfLEM»l )*(C(NELEM*l»J)+0(.DC(NfLEMtlf)J))/2.
.....
, UVCTOPa(UV(I+l)-wS*8Wll)(I+l)*»UEN)*(C(I+l»J)+ULDC(l»l
FIN
, J)a(OVC8TM-UVCTOP)*DELTD
1. . (Cf U
. ...PIN
...FIN
»«4**t»««»«*
CHANGE 3/26/81
, J)tOLOC(I»t,J)-OLOC(I,J) J/(Z,»OELZ)
00 { Jsl
WHITE (5) VOLUME
DO (J=t,MAXCON)
. SALfJ) s GC1N(NELMP| ,J)-(1CUUT(NELMPI, J)+OCEUAV(J)
. DIF(J) « CELAV(NELMPJ,J)-BAL(J)
...FIN
TBALsO.
TOIF»0.
DO (jB4iMAXCON)
. TH*U « TBAL t BALfJ)
. TOIF « TOIF t OIF(J)
...FIN
,) (BAL(J),J»l,M4XCON),T8AL
t» BED SEDIMENT AND CONTAMINANT (KG, PC)/EL£MENr *»»
DO
.
,
...
VP
DU
(Jst,M»xCON)
oTBtU(J) * TBEO(J)
IBEP(J) a 0,
FIN
t | )»AUF.N
00 (l=l,N8EO)
WHEN (J ,UI, » )
. UNSlTY'd ,(I-POR)/(B( I, 1 )/OEN3d) »B ( I , 2) /DENS r 2) *B ( l.3)/D£N3(3) )
. OELxBDIV
IFtI.ElJ. NBFO) DELsXYSO
VOLsOEL*VPXS
. CELAV(I,J)»P(I, J)*VOL*ON8ITY
...FIN
ELSE CEUAVU, J)»CELAV(I,J-3)*8(1.J)
135
-------
(FLECS VERSION 22. U6)
00166
00tb7
00168
00170
00171
0017Z
00173
0017"
00175
00176
00177
00178
00179
00180
, , TBED(J) s T8EO{J) » CEUAV(1,J)
, ...FIN
...FIN
no
-------
(FLEC3 VERSION 22.46)
2B-JAN.f)2 10124130 PAGE 00001
00001
00002
00003
0000/4
00005
00006
00007
00006
00009
00010
00011
00012
00013
00014
00015
00016
00017
00019
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
000<47
OOOifl
00049
00050
00051
00052
00053
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
SUBROUTINE 8E'JDAT(UECAY, DENS, DFZ, OIAM, OSHR, 050, ECHO,
1 ErtODE, HLDEUK, MUMEHft, SCSHR, 30RBK, VSET,
2 OSORB)
THIS SUHROUTME 13 RESPONSIBLE FOB REAOING ANU PROCESSING THE
SEDIMENT CHARACTERISTICS.
FORMAL PAMAMETE^Sl
1>EC»Y ntCAY PARAM£!E«S
OENS SPECIFIC WEIGHT
DFZ VERTICAL DIFFUSION COEFFICIENTS
DIAM PARTICLE DIAMETERS
OSHR CRITICAL SHEAR STRESS VALUE FOR DEPOSITION
D50 MEDIAN »EO SEDIMENT DIAMETER
ECHO LHE PRINTER ECHO CONTROL VARIABLE (1*1)
ERODE ERuDABILITY
HLDEHR HOLDING ARRAY FOR ERROR NUMBERS (BYTE)
NUhERH NU»'HER OF INPUT ERRORS
SCSHK CRITICAL SHEAR STRESS VALUE FOR SCOUR
SORBK AOSOHBTION VALUES
DSORH • I;ESO»PTJOM VALUES
VSET - VERTICAL SETTLING VELOCITIES
CALLED BYlSERATRA
CALLESl PUffcRH
BYTE MLOERR(IOO)
LOGICAL*! ECHO
DIMENSION DECAY(6),I>FNS(3),OFZ(«),DIAM(3),03H«(3),EROPE(3),
t SCSHH(3),SOR6K(9),DSORR(9),V8ET(3)
....PARTICLE SETTLINE VELOCITY (M/SEC)
COL. 1-10.. ..VSET(l)... .SANQ SETTLING VELOCITY
n-20....V3ET(2)....SlLT SETTLIMP, VELOCITY
21-30. ...VSET(3),,. .CLAY SETTLING VELOCITY
f?EAO(l,l) (VSET(I),Ial,3)
IF(ECHO)
. WRITE(6,2)
. V.RITE(6,12) (VSET(I),I»1,3)
...FIN
....DENSITY (KS/M**3)
COL. 1-10, ...D£NS(t J... .DENSITY OF SAWD
n-20....ieNS<2)... ."Easily OF SILT
21-30. ...r-ENSO).,. .DENSITY OF CLAY
9EAD(1,1) (D£NS(I),I«1,3)
IF(ECHO)
137
-------
CFLECS VERSION 22,a<>)
10»2«l3« PAGE 00002
ooo5u
00055
0005»>
00057
00058
00059
00060
OOObl
00062
00063
00064
00065
00066
00067
00060
00069
00070
00071
00072
00073
00070
00075
00076
00077
00078
00079
00080
000«t
00092
00063
0006*1
00095
00086
ouoer
000«9
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103
0010U
0010S
00106
00107
00106
00109
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
. »RITE(6,3)
. *RITE{6,12) (OENS(I),I*1,3)
...FIN
IF(0£NS(1) ,LE. 0.0) CAUL PUTERRt \ 7,NUMERR,HLDERH)
IFCOEN8(2) ,LE. 0.0) CALL PUT6 RR(ie,NIJMER*, HLOE«R)
IF(I)EN3(3) ,LE. 0.0) CALL PUTE««( 19, MJMEPR, HLOERR)
....DIAMETEH (METERS)
COL. 1-10. ...OlAM(l)... .DIAMETER OF SAND
11-20. ...DIAM(2)..,,01AMETE* OF 3ILT
21-30. ...OIAiMJ).., .DIAMETER OF CLAY
31-«O....D50.,.,...,MEr>IAM OEO SEDIMENT DIAMETER
READ(l.l) (DIAM(I),I«1,3),D50
IF(ECHO)
. «RlTE«>i'J)
. ^BItE«>,125 (r»IAMH),I»lr3)
. *«ITE(6,13) D50
...FIN
IF(DIAM(U t|.E, o.O) CALL PUTEHe(20,NUMEHR,HLDE«R)
IF(D1AM(2) ,LE. 0.0) CALL PUTESR (21 , NUMERR,HLDERH)
IF(OIA«(3) .Lfc. 0.0) CALL PUTERR(22,NUME«B,Hl_OERR)
....CRITICAL SHEAR STRESS FOR SCOUR (KG/M**2)
COL. 1-10. ., .SCSHRd) ,.,. CRITICAL SHEAR STRESS FOR SAND
11-20. ...scsn«(2) ,ti. CRITICAL SHEAR STRESS FOR SILT
21-30.. ,,SCSHK(3) ....CRITICAL SHEAR STRESS FOR CLAY
REAO(1,1) (3CSHH(I),I»1.3>
IF(ECMO)
. WHITE (0,5)
. ^RITE(6,12) (SC3HR(I),I«1,3)
...FIN
IF(SCSHRd) .LE. 0.0) CALL PJTERR (26. NUME«R,HLDE*K)
IF(SCSH«(2) .LE. 0.0) CALL "UTERR (27,NUMER»,HLOE«R)
IF(SC3H«(i) ,LE. 0,0) CALL PUTERR (26,NUME^» ,HLOERR)
....CRITICAL SHEAR STRESS FOR DEPOSITION («G/M**2)
CQL. l-lo.,..oSH«n).., .CRITICAL SHEAR STRESS FOH SAHO
11-20... ,osMR(2) . ...CRITICAL SHEAR STRESS FOR SILT
21-30. ,..OSHK(3)... .CRITICAL SHEAR STRESS FOR CL»Y
REAH(1»1) (DSHB(I),I»1»3)
IF(ECHQ)
. WRITE(6,6)
. *HITE(6,12) (DSHR(I), I»l,3)
...FI^
IF(OSHR(t) ,LP. 0,0) CALL PUTEH«(29,NUMERR,HL«E«R)
IF(n3HR(2) .LE. 0.0) CALL PUTEW
-------
2H.jAM.rt?
PAGE noOD3
oot to
OOltl
00112
00113
00114
00115
00116
00117
oolis
00119
001'CH'fc(2)....ERur>AHILITY UF SILT
21-30.,..EkiODt (3) ....EHijnAHILITY Of CLAY
IF(ECMO)
. wRITF_(6,7)
. WRI TE(<>, 12)
...FIN
IF(EHOOEd) .LE.
IF(EROi>E(2) ,LE.
IF(EKOOE(3) .Lfc.
0)
0)
0.0)
CALL
CALL
CALL
....VERTICAL DIFFUSION COEFFICIENTS (M**2/SF:C)
COL.
l-10....0F2(n..,.COEFFICIENT
11-20....0FZ(2)....COEFFICIENT
2 1-30.... OFZ( 3)..,, COEFFICIENT
3 1 -40.,.. PFZC 4)..., COEFFICIENT
FOH
FUR
FOR
FOH
SAND
SILT
CLAY
DISSOLVED
CONTAMINANT
READd»l) (DFZ
IF(F.CHO)
. w«ITE(6,8)
. WHITEC6,12)
...FIN
... .ADSOSBTION VALUES C2 CARDS)
IDFZCI>,I«1,3)
CARD *1
COL. 1-10.
11-20.
21-30,
31-40.
41-50,
51-60.
61-70.
71-80.
CARD «2
COL. 1-10..
.SORBKd).
.SORBK(2).
.S'JHHM 3) .
,30RBK(4) ,
.SORUM5) .
,SORBK(6) .
.SnKHK(7),
,SflRHK(8) .
,.SORBK(9),.
.KD VALUE
.HO VALUE
,KD VALUE
.SUSPENDED
.SUSPENDED
.SUSPENDED
.BED SAND
.BED SILT
,.BEO CLAY
WITH SAND (M**3/KG)
WITH SILT (M**3/KG)
WITH CLAY (M**3/KG)
SAND MASS TRANSFER RATE
SILT MASS TRANSFER RATE
(1/S
u/s
CLAY MASS TRANSFER RATE (1/S
MASS TRANSFER RATE (I/SEC)
MASS TRANSFER RATE (I/SEC)
MASS TRANSFER RATE tl/SEC)
READd,!) (SORHKI
r9)
....DESORPTION VALUES (2 CARDS)
CARD *
COL.
CAPO *<
COL.
1
1-10.
11-20.
21-30.
31-40.
41-50.
51-60.
61-70.
71-80,
>
1-10..
.OSORfld),
.nSOR8(2).
.DS'lPH(3).
.DSOWB(4).
.OSORH(5).
,OSOR8(6) ,
,OSOPB(7) .
.D30RHC8).
,.DSnRH(9),.,
,KO VALUE WITH SAND CM**3/Km
,KD VALUE wITH SILT CM**3/*G)
,KD VALUE «ITH CLAY (M**3/KS)
.SUSPENDED SAND MASS TRANSFER RATE (1/S
.SUSPENDED SILT MASS THANSFE" RATE (1/S
.SUSPEMOEO CLAY MASS TRANSFER RATE (1/S
,HED SAND MASS TRANSFER RATE (I/SEC)
,8ED SILT MASS TRANSFER «AT£ (I/SEC)
,8ED CLAY MASS TRANSFER HATE (1/SEO-
139
-------
(FLECS VERSION 28,06)
iQ|2Y EXCEPT
FOR RAOIOMJCLIOE DECAY, IF GIVEN, THE
REMAINING PARAMETERS ARE NOT TO BE
SUPPLIED, PESTICIDE ONLY ci/'sea
21-30.. ,.OF.CAY(6)..,.VOLATIZATION DEGRADATION RATE,
PESTICIDE ONLY (I/SEC)
31-40. ...PH... .DEGREE OF ACIDITY OF ALKALINITY
41-50.... AM. ........ SECOND OROER ACID RATE CONSTANT
FOR HYDROLYSIS
51-60....AKU SECOND OROER BASE R»TE CONSTANT
FOR HYDROLYSIS
61-70, ...AKN.,. SECOND OROER HATE CONSTANT OF NEUTROL
REACTION *ITH WATER
RADICAL OXYGEN FOR OXIOAUUW
CARD #2
COL 1-10,.,.R»2.... CONCENTRATION OP FREE RAPIC'L OXYGEN
U-20,.,.AKBIO,,,,,,,SECONO OHOER RATE CONSTANT
8100EGRAOAUON
21-30. ...8IOMAS..,,,.8IOHASS PEH U*IT VOLUME
RE AD (1,1) DEC4Y(l),nECAY(2),DECAY(6),PH,AKA,AKB,AKN,AKOK,R02,
1 AKBt(*,BIOMA9
DO (I«3,5) DEC»Y(I) * 0,0
wH£N COECAY(2) .NE. 0.0)
, IFCECHO) wr
-------
22.46)
28-,IAN.«2 10(24134 PAGF 00005
00222
00223
00224
00225
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236
00237
00239
00239
00240
00241
00242
00243
00244
00245
00246
00247
0024*
00219
00250
00251
00252
00253
00254
00255
00256
00257
00258
00259
00260
00261
00262
00263
00264
00265
00266
00267
00260
C
1
2
3
4
5
6
7
e
9
10
1
2
11
1
2
3
4
5
6
7
fl
9
1
2
3
4
5
6
7
a
9
12
1
2
13
14
15
1
2
3
4
5
6
16
C
FORMATOF10
FORMATC 1HO|
FORMAT( 1HO,
FORMATdHo,
FORMATC 1HO,
FOMMATUHO,
FOHMAT( 1HO,
FORMATUHO,
FORMAT
14X,
l«X,
(1HO,
1PE12
1PE12
.0)
13X,
13x,
I3v,
13*.'
13x,
13X,
13x,
13X,
•
•
5
,'
5, '
'PARTICLE SETTLING VELOCITY (rt/SEC)')
'DENSITY (KG/M**3)')
'DIAMETER (METEKS)')
"CRITICAL SHEAR STRESS FOR SCOMR (KG/M**2)')
'CRITICAL SHEA* STRESS FOR DEPOSITION (KG/M*»2)')
'ERPOARIHTr (KG/M**2/SEO ' )
'VENTICAL DIFFUSION COEFFICIENTS (M«*2/8EC)')
'ADSORPTION KO VALUES')
'DECAY PARAMETERS'/
•
•
9
•
9
•
HADIONIJCLIOE DECAY (I/SEC)'/
TOTAL »EC*Y (EXCEPT RAuIONUCLDE) - (I/SEC)')
FORMA T(l HO, 13X, 'DECAY PARAMETERS'/
14X,
14X,
14X,
14X,
14X,
14X,
14X,
14X,
14X.
1 4X,
14X,
14X,
14X,
1PE12
1PE12
1PE12
1PE12
1PE12
1PE12
IPE12
1PE12
1PE12
1PE12
1PE12
IPE12
1PE12
•
•
•
•
t
•
•
t
,
.
*
*
*
5
5
5
5
5
5
5
5
5
5
5
5
5
,'
,'
,'
,'
,'
.'
,'
,'
, '
, '
* '
> '
,'
•
•
•
«
•
•
«
•
»
.
0
.
•
FORMAT{ 14X, 1PE12
14X,
14X,
FORMAT
FORMAT
F.ORMAT
14X,
14*.
14X,
14X,
14X,
14X,
FORMAT
END
1PE12
1PE12
( 1 IX,
(14X,
(1HO,
1PE12
1PE12
IPE 12
1PE12
1PE12
1PE12
(1HO,
•
•
1
5
5
P
,'
,'
El
•
*
2
1PE12
13X,
*
•
;
*
•
t
i
5
5
5
5
5
5
r
»
,
r
r
f
3*,
»
•
•
•
•
»
«
•
*
•
,
,
•
t
*
•
*
•
9
•
9
9
9
9
9
•
9
t
,
,
9
KAOIONUCLIOE DECAY (I/SEC)'/
PH . DEGREE OF ACIDITY OR ALKALINITY (PH)'/
SECOND ORDER ACID HATE CONSTANT FOR HYDROLYSIS"
' (AKA)'/
3FCOMD O^DEh QASE RATE CONSTANT FOR HYDROLYSIS'
' (AK8)'/
SECOND ORDER RATE CONSTANT OF NEUTROC REACTION'
I WITH WATER (AKN)'/
CHEMICAL DEGRADATION DUE TO HYDROLYSIS'/
SECOND OWOEW RATE CONSTANT OF FREE RADICAL'
' OXYGEN (AKOX)i/
CONCENTRATION OF FREE RADICAL OXYGEN (B02)1/
CHEMICAL DEGKAIM TION DUE fO OXIWATION'/
SECOND ORDER "ATE CONSTANT FOR 4IODEGPADATION '
' (AKBIO)'/
HI01ASS PER UNIT VOLUME (^IDHAS)'/
^lODEGRAOATION' /
VOLATILIZATION' )
C t O A XI fl 1 J
J t 999^ '•U •
•
9
5
5
SILT'/
CLAY')
.'...MEDIAN RED SEDIMENT DIAMETER')
.'...DISSOLVED CONTAMINANT')
•MASS TRANSFER HATES (I/SEC)1/
•
•
•
•
•
•
•
.
.
.
.
.
SUSPENDED SAND'/
SUSPENDED SILT'/
SUSPENDED CLAY'/
SAND ATTACHED TO THE BED'/
SILT ATTACHED TO THE BED'/
CLAY ATTACHED TO THE ^EO')
•OEStJHPTIOM KO VALUES1)
(FLECS VERSION
141
-------
(FLECS VERSION 22.16)
13-MAR.ai
PACE ooooi
oonoi
00002
ooooi
o o o o H
00005
OOOD6
00007
00006
0000"'
00010
0001 1
00012
0001 3
0001U
00015
00016
00017
00018
00019
00020
00021
00022
00025
00021
00025
00026
00027
0002B
00029
00030
OOPil
OOOJ2
OOOJJ
ooo3«
000 J5
OOnib
00037
00039
OOOJ9
000*40
OOOt 1
OOC«2
OOOUJ
000««
00015
OOOIh
OOOfT
OOOIfl
00019
00050
00051
00052
0005J
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
SUBROUTINE SEDIMEUB4R, ALEN, CCIN, DELTO. If J» MELE'S
1 QHIM, QHOUT, OVf SO, SH, *LFA, BETA, V£L 1 ,
2 VEL2, V3ET, OELZ, BWlO, A^IO, DEPTH,
i 8ETAI, BtTA2)
THIS ROUTINE CALCULATES COEFFICIENTS OF CONVECTION, DIFFUSION,
DECAY AND SOURCE TERMS IN THE SEDIMENT TRANSPORT CONVECTION
-DIFFUSION EQUATION
INPUT PARAMETERSl
A8AR • AVERAGE AREA
CCIN • CONCENTRATION OF INFLOrt
OEITD • TIME STEP IN DAYS
DEPTH - DEPTH ()F RlVEP SEGMENT
I - ELEMENT INDEX
J • PARAMETER INDEX
NELEM • NUMBER OF ELEMtNTS
OHIN - INFLOW DISCHARGE
QHOUT • OUTFLOW DISCHARGE
0V • VERTICAL DISCHARGE
SO - SEDIMENT DEPOSITION RATE, ( KG/M* *3/D A Y )
SR - SEDIMENT fRQSIUN RATE, ( K G/M* *3/DA Y )
OUTPUT PAHAMETERSl
ALFA - DECAY TERM, (I/DAY)
BETA - SOURCE OR SINK TERM, (KG/M* *J/OAy )
PtTAl • INFLUENT SOURCE TERM FOR THE I-TH NODE, ( KG/M* * J/DA Y )
PETA2 - INFLUENT SOURCE TERM FOR THE Itl-TH NODE* (KG/M**J/D
VEL1 • FIRST CONVECTIVE T£RM, (M/OAY)
VEL2 - SECOND CONVECTIVE TERM, (M/DAY)
CALLED BY TRANSP.
INCLUDE 'SYiELMsU.PRMi
DIMENSION ABAR(MXELEM) , AREA (MXELEM ) , CC IN (MXELEM, MAXCON) ,
| OHIN(MXELEM), QHOUT(MAELEM), QV ( MXF.LEM) , SO ( MXELE* , 6 ) ,
2 SR(6), VSET(3), 8^ ID ( MX£L£M ) , A»UO( MXELEM)
CONVECTIVE IERM WITH CORRECTION FOR A CONTINUOUS SETTLING FLUX
AQ = QVU }
VELla(A(J-VSET(J)*0WIO(I)*ALEN)/ABARtI)
A U s 0 V ( 1 t 1)
VEL2»(AO-VSET(J)»flWiO(Itl)*ALENj/AHAR(I)
DECAY TERM
ALFA s UHoUT(I) / (ABAR(I) * DELZ)
SOURCE OR SINK TEHM
BETA=SR( J) - SD(I, J)
8£TAlsQHINU)/(AeAR(I)*OELZ}»(CCTN(I,J)/3.*CClN(I+l,J)/6.)
142
-------
CFUECS VERSION 2J.U6) 1J-M4««61 13U9JS8 PiGE 00002
00059
00055 R6TURN
00056 END
(TLECS VERSION 22,
-------
(FLECS VERSION 22.Ob)
2».JA'4«fl2 10J16I11 PAGE 00001
QQOOt
00002
00003
00004
00005
00006
00007
ooona
00009
oooio
00011
0001?
00013
00014
00015
00016
00017
oooia
00019
00020
00021
00022
00023
00024
00025
00026
00027
0002*
00029
00030
00031
OD032
00033
0003'!
nooi5
00036
00037
0003S
00039
00010
ooo«l
OOOU2
00043
r> 0 0 1 '»
00045
000'J6
000«7
00018
000«9
00050
00051
00052
00053
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 1
C
C
C
C
C**
C 2
C
C
C
C**
C 3
C
C
C**
C a
C
C
C
C
C**
C 5
C
C
C
C**
C 6
C
C
tm,n$E«ATRA.FLX
VE«SION| DRIED CHANNEL OPTION
WITH SELECTION FOR IMPLICIT SETTLING VELOCITY SCHEME
o» SUPFHPOSITION OF SETTLING VELOCITY SCHEME
r*******t****«********************************+***»***»**************
THI3 COMPUTER PWQGPAM, SERATRA, IS AN UNSTEADY, TWO-
DIMENSIONAL (LO»G1TUDINAL ANO VERTICAL) M(lt>EL TO SIMULATE
SEDIMENT-CONTAMINANT TRANSPORT IN RIVF.HS AMD RIVER-RU"
RfcSF.WVOIRS.
THE MODEL MAS GENERAL CONVECTION. DIFFUSION EUUATION3
«ITH DECAY ANO SINK/SOURCE TERMS *ITH AHPROPRIATE 80IJNOARY
cnNr>nifj'4S.
SERATWA UTILIZES THE FINITE ELEMENT COMPUTATION METHOD KITH
THE GALER*!^ WEIGHTED RFSIOUAI. TECHNMUE.
THE FOLLOWING WPPflHTS OESCHI8F. 3E«ATHA MODEL FORMULATION, IJSEH'3
MANUAL AMD giME MUOEL HESULTSi
"STUDIES OF COL"HHI* RIVE» WATER OUALITf-OEVELOPMENT OF
MATHEMATICAL MOnfLS FOH SEDIMENT ANO HAuIONUCLIOE TRANSPORT
ANALYSIS." BNwL-3-1'52. HATTELLEr PACIFIC NOHTH^EST
UAHOHATOHIES, P1CHL»"">. 'tA.
OMtSHI,r. 1977, "FINITE ELEMENT MODELS FOB SEDIMENT A^O CHNfAHlNANT
T«AM3POPT IN SURFACE WATE^S-'TRANSPORT OF SEDIMENT ANO RAOJ ONUCLI.IES
IN THE CLINCH WIVEW," HN»L-2227. HATTELLE, PACIFIC NORT"xFST
LAtJOHATURlES, PICtUAHO, *A.
ONISHI,Y. 1977, "MATttEMATICAL SIMDLATinr) UF SEDIMENT \NO RAOIH.
NUCLIOE TRANSPORT JM THE COLU^UIA RIVFH." 8NWL-2229, HATTFLLF,
PACIFIC NOKFHwEST L»HOPATnRIE3, HIGHLAND, WA,
ONI3HT,Y.. n.L. SCHkeiHEK AND R.8. CUOEIL. 19f9. -M » THEM » T I C *L
SIMULATION UF SEDtMEMT ANiJ RAO IINHCLIOE TRANSPORT IN THE CLINCH
RIVE^, TENNESSEE," ^HOCEEDTNGS OF ACS/CSJ CHEMICAL CONGRESS,
HONOLULU, H»K»H» A"«IL 1-6. 1979, "CONTAMINANTS *N° •>En IMENT 3" ,
k,4. MAKER (F.D.), ANN ARIOrf SCIENCE PllHLl5HE»S, piC,, k^N ARBOR, MI,
ONISHI,Y,, 3,*, BW01N, A.H, OLSEN» M.A, PARKHURST, S,E, '''ISEi ANO
H.H. WALTERS, 1979. "HETHUDDLOGY FOR OVERLAND AND INSTREA* ••IISRATIOM
ANU RISK »SSESS-«ef
-------
(FLECS VERSION 22.16)
19114(11 PAGE 00002
0005"
00055
00056
00057
0005f
00059
00060
00061
00062
OOCo3
00061
00065
00066
00067
00060
000b9
00070
00071
00072
00073
00071
00075
00076
00077
00079
OOP79
00080
00081
000«2
00083
OOOU4
000«5
00006
000*7
00096
00089
00090
00091
00092
00093
00094
00095
00096
P0097
00098
00099
00100
nomi
00102
00103
00101
00105
OOJ06
00117
ooioe
OOinq
c
c**
C 7
C
C
C
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
PA
ON
SF,
NO
CA
1
1
2
1
1
2
3
a
5
6
i
1
1
1
1
1
I
i
1
2
3
1
PACIFIC NORTHWEST LABORATORIES, RICHLANO, HA,
C 7 ONISHI,Y. ANO 3.E. wISE. 1979, "USER'S MANUAL FOR THE I*
SEDIMENT-CONTAMINANT TRANSPORT MODEL, SERATRA," BATTELLE, PACIFIC
WORTHIEST LA80R*T'lHTtS, RICHLANU, rt»,
CALLSl HEDOAT, COLAPS, OIMOAT, FCOOE, HYOOAT, MYDFLO, 1CFLO, 1NIPAT,
ROSFLO, aOO./
DATA SECYR /31536000/
OATA tCHO/.FALSF,/
DATA FERK'IR/, F»L8E./
DATA ZE*0/t.O*-tO/
DATA SAVECH/, FALSE, » .FALSE ., .FAl SE ., .FALSE ,. .FALSE ,,. FALSE ., .FALSR
1... FALSE.,. FALSE.,. FALSE./
CALL SrPTUPCBASE, OEV, ECHO, FNAME, FTYPE, GUIC, JNFLO, ISTRT,
MFUST, OUTFLO, SMETH, UUIC, SAVEChf J3EG,
NSTF,PS,"SEG, ITP»T, ANALY3,n£LTH, AN»L^T,DEP1I
PELEV)
*** BEAD INPUT DATA FDR BEGINNING SEGMENT ***
NUMFRR • 0
THIS IS THE Rt'STAWT OPTION.
145
-------
(FLECS VERSION 23.46)
00005
00110 >
00111 C
"0112 e
00113
00114 1
00115
00116
00117
00118
00119
00120
OQU1 1
00122 2
00123
00121 1
00125
00126
00127 C
00128 C
00129 C
00130
00131 C
00132 C
00133 C
00131
00135
00136 1
00137
00138 1
00139 C
001 '10
001 '41
001U2
001«3
OOlttiJ C
001 '15
00146 C
00147
001«8
001UAT(ANALMT, ANALYS, OELTH, ECHO, HLOER«, ITPRT,
•'SF.G, NSTEPS, NUMERR, SIMLEN, OEPMIN)
PELTO»OELTH/SEC»AY
IF(ECHO)
. CALL TiME(NOfc)
. *RITE(6,3) ISTRT, NO*
. . «F I ^
CALL OIMDAT( ALE"-*, AREA, ROIV, dEO, OLZSAV, ECHO, ELEV, HLDERR,
ISTMT, N«Ei>, NELEM, NUMERR, PELEV, FOR, RIVER,
XYSO, EL)
CALL SEOOAT(UEC*Y, DENS, OF?, 01AM, OSHR, 050, ECHO, EHOOE,
HLOEP.R, NUMERR, SCSHK, SORCK, VSET, oso^e)
CALL PH01NPCECHO, JULIAN,KAY1,KAY2,PCOEF)
CALL BtooATCti, ECHO, NBEP>, AREA, »OIV,OENS,POR, XYSO, TBEO)
READS INITIAL CONDITIONS FOH SEDIMENT AND CONTAMINANT
CALL »****************************
SEGMENT LOOP *
t*** *************************
)0 ( I SEG»ISTH1 , >tSEG)
CALL OIAG(F.CM02, ECHP3, iCHO'J, ECH05, ECHOfe, ECHO?, F.CHOH,
ECHH9, ECH010, ISEG, JSEG, SAVECH)
FOBMATCMSIOO*, 'SEGMENT NO.',I3,l«,flAt)
HESET-l'ATA-1 I M6 -CONTROLS
CAVGMxaO.O
JULSEC « JULIAN * stcn*y
UNLESS (IStG .E'J. 1)
, IF JECHO)
. . CALL TIME(MO-)
. -IHITF (6,3) ISEG, NO*
146
-------
(FLECS VERSION 22,<46)
2H.JAN.-J? 10J16I11
0000«
00166
00167
00168
001
r
f
f
]
...FIN
FFRHQP. o .FAtSE.
NUMtRH * 0
CALL OI^OATULEN, AREA, Bi'iv, BED, "U/SAy, ECHO, ELEv, HIDEPft,
ISEG, NHEO, NfLE*, MUhEHR, PPLEV, POR, «IVE9,
XVSO, Ft)
C»LL SEDDAKDF.CAY, OENS, OFZ, 01AM, DSHR, 050, ECHO, ERODE,
HLDERR, MUMER^, SCSHR, SOHBK, VSET, OSOR«)
CALL flenOaTO, £CHO, N8EU, ARE*. HOI V, DENS, P0», XYSO, THED)
REAUS INITIAL CONDITIONS FOR SEOIMgMT AND CONTAMINANT
CALL wT^DATfC, ECHO, NELEM)
CALL TRBOAT(ECHO, MLOERR, NELEM, NTRI9S, NUMEPR,
SIMLE*, TBBOPn
CALL HYOOAT
DO (K«1,MAXCON)
*** NELEM VALUES OF CL*ST MUST RE CALUCLATED FROH NELEH+l
VALUES OF C **»
IF (K .LE. 3 ,OR. K .60, 7)
PROVIDES A CONSISTENT INITIAL CONDITION FOR THE DELZ DEFINED
IN HynpAT - 80TH ^OOAL VALUES (CNOOEsC), A^n ELEMENT
AV£HA«ES (CLAST»CULO)
CALL COLLAP(ALEN,AREA,C,DELZ,EL,IELM,K,N(iLEM,CLA8T»TMASS,AWID,
V3GT,OFZ,CNOOE, XSA«EA)
..^1N
..FIN
DIAGNOSTIC «PI TES
IF (ECM05)
wRiref6, isou)
WRITE (6, 1510)
"RITE (6, 1520) (1, AREA(I),ELCI), AHAR( 1), A*ID(I), 8*10(1),
XSAREA(/),TMASS(J), 1*1, NELEM*!)
f»RlTE(6, 1530)
•«niH(6,is20) (i, (Cd.J), J»i, MA XCON ),i»i, MXELEM)
147
-------
(FLEC3 VERSION 22.46)
23-JAN.82 10116111 PAGE 00005
00222
00223
00224
00225
00226 1500
00227 1510
00226 1
00229 1520
00230 1530
00231 1540
00232 1550
00233
00234
00235
00236
00237
00239
00239
00240
002"!
00242
00243
00244
WRITE (6, 1540)
*RITE(6»1S2U) lit (CNODE (I» J) ,J«l,MAXCON),I»l,NF:LEilOX,|l-HA3S')
FORMAT(2X,I5,3X,1P7E15.4)
FORMATC'ONQDAL CONCENTRATIONS PRIOR TO COLLAP')
FORMATCONOOAL CONCENTRATIONS KOLLOKING COLLAP*)
FORHATCOELEHENT AVERAGE CONCENTRATION FOLLOWING COLLAP')
..FIN
0 (I«l,NELEM)
00 (J»1,MAXCON)
. COLO(I,J)»CLAST(I,J)
. C(I,J)BCNOOE(I,J)
...FIN
..FIN
DO (jBl,MAXCON) C(NELEM+l, j)»CNODE(NELEMti , J)
00 (jBl,MAXCON) COLO(NELEM»1»J)BO.O
IF(NELEM*1.LT,MXELEM)
. DO (I*NELE*+2,MXELEM)
. . DO (J»1,MAXCON)
• . , COLD(I.J) • 0.0
• . . C(I,J) • 0.0
FIN
00246 FIN
00249 . . ...FIN
00250 . ...FIN
00251 . CALL FCODE(FNAME, BASE, ISEG, FTYPE, OEV, GUIC, UUIC)
00252 . OPE
•<( UN IT»5, NAME »F NAME, TYPE" 'NE»',FOrt MB 'UNFORMATTEU'!!
00253 . OPE N ( UNI T»9, MAMF» 'R5TRT.F I L',TYPE»'NE*', FORM* i UNFORMATTED' )
00254 . WRI
00255 C
00256 C . ***
fE(5) ISfcG
CONVERT INPUT VALUES TO THOSE UNITS USED HY MODEL ***
00257 , 00 (J»l,3) VSETtJ)ovSET(J) * SECDAY
00258 . 00 (J«l,4) CFZCJ)BOFZ(J) * SECOAY
00259 . DO
00260 . DO
00261 . DO
00262 C***********
00263 C .
(I»l,6) DECAYCI) ' OECAY(I) * SECDAY
(184,9) SORHKCI) » SORBK(I) * SECOAY
(I»4,9) OSORH(I) • OSORB(I) * SECOAY
**********************************
THE STEP LOOP *
00264 C*********************************************
00265 C
00266 . NXE8 • NFRST
00267 . IF
(f'XEQ.E'J. 1 )
00266 , . PTOPTH » <>EPTH
00269 . . PTOELZ » I'ELZ
00270 . ,
^EL^PT B NELF.1
00271 . ...FIN
00272 . ETI
00273 . ETI
IE B NX£U - 1
^E » ETi^t * DELTH
00274 . UNTIL (NXE'J .GT, NSTEPS)
00275 . . IF (IS1 ,E'J. 1) xRITEUl.'J) ISEG.^XEO
00276 4 . . FORMATC SEGMENT »',I3.' TI*E STtP',110)
00277 . .
ETIME * ET1MF. + OELTH
148
-------
(FLECS VERSION 22,46)
2R.JAN.82 10116111 PAGE 00006
00278 C
00279
00280 1
00261 2
00282 3
00293 C
00284
00285 C
00286
002*7
00296 1
002«9
00290
00291 1100
00292
00293 1110
00294 1
00295 2
00296 3
00297
00298
00299
00300
00301 C
00302
00303 1
00304 2
00305 3
00306
00307 9999
00308
00309
00310
00311 1000
0031?
00313 1010
00314 1
00315 2
00316 3
00317 4
00313
H0319
00320
00321 1020
00322
00323
00324
00325
00326
00327
00328
00329
00330 1560
00331 1
00332 2
00333 1570
*** UPDATE THE FLOW AND CONCENTRATION ARRAYS ***
CALL HYrFLOfALEN, ARE», AK-II), OELZ, DEPTH, 050, EI.EV,
ENDHYO, ETIrtE, FERROR, HR*D, Nf.LEM, K5"*«I, PELEv»
IJHIN, aHouT, 'iv, RIVER, SLOPE» STRF.SS, TE*PR»
VEL, VOL, DEPMIN, XSAREA, 8*10, AIM", <)lr CROSEc)
IF (FERHOH) HEPORT-FATAL-ER90R-ANn.STOP
IF (NTR18S ,r,T, 0)
CALL TRBFLO(CTR8, CTHIH, ENOTfll, ETIM£, FERROR, DEPTH, NF.LEM,
NE»i'JI» NE*THH, QHIN, TWBOPT, DEP*IN)
IFfECH03)
WHIT£(6, 1 100)
FnRMATCOAFTER TRHFLO, MASS FLUX VALUES OF CTRB')
NKITE«>, 1110)
F(i«*»AT( IOELE^ENT NO, ' , 2X, 3(3X, 'CO^C. OF • , 4X) , 3( 1 X, 'CONC, •
'A3SOC.',3X),2X,'CONC*MINANf '/1 3X,' SUSPENDED SANOMx*
•SUSPENDED SILT', u, 'SUSPENDED CLAY', 3x,'*iTH SAND», &x,
'WITH SILT', 6X,'WITH CLAY', ttX, '01 SSOLVEU CONC,1)
wwtTE(6,1020)(I,(CTRH(I,K),Kat,MAXCON),lal,NELEH)
..FIN
IF (FEftROR) REPORT-FATAL-E3ROR- AND- STOP
..FIN
CALL ICFICKCCIN, DEPTH, DELZ, 050, ENOIC, ETIME, FERROW, INFtO,
ISEt'jNELC'^jNE^ICj'lHlN, 01, PEPSIN, ALEN,
UEL, UWIO, XSAREA, A^EA, ASID, OFZ, VSET,
ELfELEV.PELEVjRlVEWji'lE^OIj^EwTHB)
IF (ECHO?) XRITE (6,9999)
FO»HAT( '********»** IN SE»ATRA TIME LOOP **»******•)
IF(fcCHO,») WffUE (6,4) ISEG, NXEU
IF (ECH03)
. WRITE (6,1000)
. FORMATCOAFTER ICFLO, CCIN')
. "RJTF (6,1010)
. FQPMJT('0',2X, ''JOOE ND,',2X,3(3X,'CONC. OF ' , 4X) ,3(,1 X, 'CONC, '
. ' AS30C, ',2X),?X, 'CONTAMINANT'/lX, 'FKOrt BOTTHM'.IX,
. 'SUSPENDED SAND', IX, 'SUSPENDED SILT ', 1 X ,' SUSPENDED CLAY',
. 3X,'"tirH SAND', 6*, 'WITH SILT' >6K, 'WITH CLAY ' , 4X, ' OlSSOLVFO ' ,
. ' CO.'C1)
. NKlTF(6,1020)(J,(CCIN(J,K),Ksl,MAXCON),J»l,NF.LEH*l)
. IF (FE^wnH) REPORT-FATAL-ERPOK-ANO.STOP
...FIN
FORMAT(2X,15,2X,1P7E15,5)
IF (EC*06I
, ''»l">ELZ
, ^PJTE (6, 1570)
. «")ITE(6,101(J)
. "IBITIE (6, 1020) ( J, (C ( J,K) ,Ksl ,fi*XCO.^) , J» 1 ,NELMPT+1 )
. "WJTE(6, 15«0)
. "iRi TE (6, 1 1 1 o )
. "*RITE(6, 1020) ( J, (COLD(J,K) ,KS1,-UXCON) , J»l ,NELMPT)
. FOHHMC '01MMEOIATELY P"IiJR TO ROSFLO, FOLLO^IMG ICFLO'X
•OPTOPTH a',E12.4,» 'JEL^PT «'»IS,' PTPELZ *',E12,4/
' OFPTH «',E12.4,' N6LEM "',15,' OELZ 3',E12.<«)
, FOR! AT( IONOUAL CONCEMTMAT10NS PRIOR TO R'JSFL"')
149
-------
(FIECS VERSION 22.46)
00334 15PO . .
00335 . .
00356 . .
00337 , .
00339 . .
00359 . ,
003'tO 1. .
10116111 PAGE 00007
CONCEMTHATIONS PRIOR TO H03FL.O')
FORMATt
IF (DEPTH ,GT.
. T «
. IF r-XEJ ,GT. 1 .AND. T .ST. ZERO)
. . CAUL
.... 00(J«1,
003«2 .....
003a3
oosua FIN
003u5 UN
00346 .... IF (ECH06)
003<»7 ..... "PITE(6,
00348 WRITEC6.1010)
00349 ..... *RITE(6,1020) ( J, (C ( J, K ) , K*l , MAXCON) , .
00350
00352 II..! >»aire(6.1020)
00353 1590 ..... FOR*AT('ONDOAL CONCE'XTHATIONS FOLLOWING
0035U 1600 ..... FOWMAT('OELEM£NT AVF.M4GE CONCENTRATIONS FUCL'JuIUG S05FLU1)
00355 FIN
00356 , . . ...FIN
00357 . . . COMPIiTE-BEP-ANO.WATEM-SlJRPACf-EUEVATIONS
00358 FIH
00359 C . .
00360 C . . *** AVERAGE THE INFtO* CONCENTRATI0't8 I'^TU THE SEGMENT HY
00361 C . . INTO ACCOUNT THE T»IHUTA»Y INPUT,
00362 C . ,
00363 . . IF (DEPTH ,GT, HEPMIN)
00364 . . . «MEN (NTRIBS ,f3T. 0)
00365 .... IF (NE*TR3 .OH. NEwlC .OR. '4F.s
00366 ..... DO (K«1,1AXCON)
00367 DO (Jsl,NELEMtn
003r>S , COlJfiiY(J) a CCIM(JfK)
00369 ...... ..,'
00370 DO
00371 .......
00372
00373 C
0037U C ....... NQTEJ C^ASS 13 IN (KG/n*»3)
00375 C
00376 ....... "HEN (J.F.U.I)
00377 . WHEN (K.E-1.7)
00378 CCINU,K)aCMASS
ooseo . I ...... I..FIN
00361 EL"5E
003«2 ......... KKSK
003B3 ... IF(K.£T.3) KKaKK-3
003«4 ......... COEFaO,
003«5 . "SaVSET(KK)*AKEA(1)/(AwIO( 1 )*AL£N)
003«6 .....
00387
003»* .....
150
-------
(FLECS VERSION 22.
10116111 PAGE OOOOR
0039t)
00391
00392
00393
00394
00395
00396
00397
00 J9B
00399
00400
00401
00402
00403
04404
00405
00406
00007
oo4oa
oo
-------
IFUECS VERSION 22.46)
28-JAN.62 10116111 PMJE 00009
00446
00447
00448
00449
00450
00451
00452
00453
00454
00455
00456
00457
0045B
00459
00460
00461
00462
O0'<63
00464
00465
00466
00467
00468
00469
00470
00471
00472
00473
00474
00475
00476
00477
00478
00479
00480
00481
00482
00483
004A4
004H5
00486
00487
00468
00489
00490
00491
00492
00493
00494
00495
00496
00497
0049«
00499
00500
00501
,
1.
1700 .
1710 .
1.
1720 .
,
,
1.
C*******
C
C
*
2000 .
C
c*******
C
*
C
.
C
*
C
C
C
,
i.
2.
C
C*******
C
C
,
1.
2,
C
c*******
C
C
,
C
,
C
r.
2.
3.
.
C
.
s
•
,
,
%
,
. . WHITE(6.1720) (l,ARF«(t),ABAR(I)rA«CD(I),H(ilII)(I)f :JV(t),«HIN
. . OHOUT(I),1«1,NELEH*1)
, . FORMATCOGEOMETRY AND DISCHARGE INFfJ PRIOR TO TRANSP')
, . FOH*AT( ' OELEMENT/SOOE1 »6X» ' AREA' ,9X, ' ASAH' , 1 \X, ' A^IO* , I IX, '
, . 12X| '0V ' , 12X, 'QHIN* | 1 IX, 'UHOUT ' )
. . FC'««AT(2X,15,3X,1P7E15.4)
. . . ,F IN
. . CALL TRAN9P(FERROR,PCOEFF,BVIID,AWIP,AHAR, DEPTH, OLDC,CCHO«,
. . CROSEC,ECHQ7, ECH08, ECH09)
**********************************************************
. . CHANGE 4/2/81
, ,
. . IFIECH05}*HITE(6,2000)UC(I,J),J«1,MAXCON),I»1,NELE1M)
. . FQPMAU'OAFTER IRANSP'/50(U,IP7E15.5/))
, ,
***************************************************************
. .
. . IF (FERROR) REPORT-FATAL-ERROH-AND-STOP
, ,
. . COMPtJTE-BEO-AMO-WATER-SURFACE-EUEVATIONS
• f
. ...FIN
,
. *** SAVE THE RESULTS OF THIS TIME STEP, IT WILL BECOME INPUT TO
. THE NEXT SEGMENT ***
. WRITE(OUTFLO) OEPTH,OELZ,NELEH» (OHOUT (K) , XSARE A(K) , AWIO(K) ,
, K»l,NELEM),((C(L,K),K«l,MAXC'JN),L»l ,NEUE*+1) ,
, ((QLOC(l.»K),K»l,rtAXCQN)fL«l,'sIELEH+l)
,
A***-*****************************************************
. CHANGE 7/27/81
,
, WRITEC9) OEHTH,nELZ,N£l.EM,(':>HO(.IT(K),XSAHEA(K),AWIO(K),
. K«l,NELE^),((C(L,K),K«l,MAXCON),L»1»NELEM+t),
. ((OLr)C(L;K),K»lrMAXCON),L"l»N|ELEM*l)
.
A*******************************************************
,
,
. IF (ANALYS) SAVE-THE-RESULTS-FOR-TIME-SERTES-ANALYSIS
,
. IF (MOP(NXE'5, I1PRT) .EU. 0)
, *** SAVE THE RESULTS FOR PRINTING AND QTHER POST PROCESSING **
CALL SAVEIT(«, HOIV, BED, ELF.V, C. OCLZ, NPEO, NELE'i,
. NXEQ, RESELN, STRESS, XYSO, OLDC,
. ALE^j'^HINjUHUOT jCClN.UVjAXIDjBwlD.VSETjOENSjOELTDfDFZ,
. POH.TBEO)
..FIN
.
, NELMPT"*'ELEM
, PTOPTH«"E^TH
. PTHELZ«')ELZ
, 00 (iBlfHELEM)
. , PxSAR( I )«XSAREA( I)
. , Ptolt) ( I )aAi*Il)( I )
, ...FIN
152
-------
(FLECS VEKSION 22,46)
2H-.JAN-H2 10116111 PAljF. 00010
00502
0050J
00504
00505
00506
0050T
00508
00509
00510
00511
00512
00513
00514
00515
00516
00517
00519
00519
00520
00521
00522
00523
00524
00525
00526
00527
0052B
00529
*** CHECK SENSE SWITCH #2 Tf) SEE IF TwE HUN IS TO C
IF (IS2.EW.1)
CLf»SF.-THE-OPEN-FIL£S
STQPPEO ***
FORMATC//21X,'***** SE»ATH* *****'/
5X,'TERMINATED BY OPE«AT£H AFTER TIME PUANE tf',110/
5X,'IN SEGMENT «',I5)
STOP
..FIN
END OF TIME 3TEP LOOP
i NXEO i 1
END OF SEGMENT LOOP
PELEV « ELEV
NFVST « i
LUNTMP«INFLO
INFLOSOUTFLO
OUTFLOSLUN1MP
OUTFLO
CLOSE(UNIT«5)
CLOSE(UNIT»9)
..FIN
CLOSE" THE-OPF.N.FIUE3
STOP
00530
00551
00532
00533
TO COMPUTE-BEO«AN|J.HATE«-SU»FACE-ELEVAT10NS
. BELEV • ELEV * HED
. RESELN • DEPTH + 8ELEV
00534
00535
00536
00537
00538
00539
oosao
OOSal
00542
TO CLOSE-THE-dPEN-FILES
CLOSE(UNIT»U
00544
CLOSE(UNIT»3)
CLOSE(UMiTn4)
CLOSE(UNIT«5)
CLOSE(UNIT»9J
CL"5E(UN!T»6)
CLOSE(UN1T«7)
^0545
00546
00547
0054S
10 P-EPOHT-FATAL-ERMOR-AND-STOP
. cuiSE-rnE-oPEN. FILES
. OPEN fUNlT«l,NA«Ea'TTI')
153
-------
(PLECS VERSION 22.46) SB-JAN. «2 10U6U1 PAGE 00011
005^9 1 . FOR*AT(//10X, ***** SEHATRA — FATAL ERROR *****'/
00550 1, ' PRINT "SEU.LST" FOR DETAILS')
005S1 C . *** THE IF STATEMENT BELOW 13 A CONCESSION TO THE
00552 . IF(FE«ROR) STOP
00553 ...FIN
00550 TO RESET-OATA. TIME-CONTROLS
00555 . FNDIC « 0
00556 . ENOHYD a o
00557 . ENDTHB a o
00556 ...FI^
00559 TO 3AVE-THE-RESULTS-FOR-TIME-3ERIES-ANALY3IS
00560 C
00561 C
00562 C
00563 C
00564
00565
00566
00567
0056S
00569 C
00570 C
00571
00572
00573
OOS7'»
00575 C
00576 C
00577
00578
00579
005PO 1
00561
005*2
005*3 C
005*'* C
00585
005B6
00587
005B8 1
005*9
00590
00591 C
00592 C
00595
0059*4 C
00595
00596 C
00597
0059B
COMPILER
*** COMPUTE THE VOLUME OF EACH ELEMENT AND TH£ TOTAL VOLUME
OF THE KATEK COLUMN ***
AVOL«0.
00 (Iml,NELE-"<)
, EL^IVlJUn) a OELZ*ABAK( I)
. AVOL»AVOL+EL«VOL(I)
• • • • * **
*** AVERAGE UI3SOLVEO (Kli/H**J) ***
AVGUIS • 0.0
DO(I*lf NELF.W) AVGDIS«AVGOTSt(C(I, 7)+C(I + l,7))*ELMV
AVGOIS • AVfiDIS / AVOL
IF (AVGOIS ,C,T. ANAL^T)
*** AVERAGE SEOIMgNT (KG/M»*3) **»
AVG5ED » 0.0
00( I»l , NfcLE*1)
. AV6SF.O*AVGSEO+(C(I,1)+C(I*1»1)+C(I,Z)*C(I + 1
C(I.i)+C(I+l, 3)1/2,
...FIN
AVG3EO « AVGSEO / AVOL
OL(I)/2.
,21 +
**** AVERAGE (PARTICULATE (PC/KG)*SEOI«ENT (KG/*1**3n
PA«PCM * O.o
00 (I«1,NELE*)
. PARPCM « PARPCM *-ELMVOL(I)*(C(I,«)+C(I,5)+C
+CtI*l»«)+Ca»lf5)+C(I»lffc))/2.
...FIN
PAHPC^ a PAHJJCM / AVOL
*** AVERAGE PART1CULATE (PC/KG) ***
PARPCK • PAHPCM / AVGSEO
TOTKG « ( PAWPCM + AVGOIS 1*AVOL
CAVliMX a '1AX(CAVU,-1X, AVCOIS)
tFLU* a 0.0
(1,6)
154
-------
(FLEC8 VERSION 22.46) 2fl-JAN.d2 10H6|11 PAGE 00012
0059"?
00600
00601
00602
00603
00604 C
00605
. 00 Usl,NFL£M) TFLOw a TFLOW +
IFL'JW 8 TFtU" / SECDAY
I..F1N
PBOCEUUHE CROSS-«EFE«ENCE TABLE
00550 «E3ET"OATA-TI^E"CONTflOL3
00158
00559 S A VE-THE-HF.SULTS-FOK-T I ME-SEHIES- ANALYSIS
904B5
0053fl CLOSE-THE-OPEN.FILES
0050S 00527 00516
00530 COHPUTE-8EO-ANO-rtArER-SURFACE-ELEVAT IONS
00357 00165
00515 REPORT-FATAL-EHHOH-ANO-3TOP
0014? 00186 00284 00299 00319 00'463
(FLECS VE«SION
155
-------
CFLECS VERSION 22,
13129143 PAGE, ooooi
oonot
oooo?
00005
ooooo
00005
00006
00007
00008
0000<>
oooio
00011
00012
0001)
00010
00015
00016
0001 7
00018
00019
00020
00021
00022
000*5
00020
00025
00026
00027
00029
00029
OOOJO
OOOJI
00052
00055
OOOJ5
00056
OOOJ7
00058
OOOJ9
00010
00041
00042
00043
0004(1
000
S(fXELEM,3),
P£L(I
PEL(1
PEH2
PEL(2
SEL(I
SEL(1
StLt2
SEL(2
,1)
,2)
.1)
,2)
,1)»
,2)«
,U»
,?)»
00
00 (J«l,2)
, PfcL(I,J)
, 8tL(l,J)
1./3.
l./b.
I,/*.
l,/3.
UD(1) t 00(3) * OD(7)/3.
-00(1) » D0(«) t DO(7)/6,
-00(1) - 00(5) » 00(7)/6.
00(1) + 00(6) + 00(7)/5.
) SEL(l,I)sSEL(l,l)tveL
9EL(2,2)*SEL(2,2)+VEL
fJ) * WIDTH
SCL(I,J) * WIDTH
oo (i«p,io,i) com < ood) * WIDTH
002=OU(8)/2,0
DO (Jsl,2)
, NR*1 1 t J • 1
, 00 (*sli2)
. . MC«2»(I1 » K - 1) « NR
, , P(NR,MC)» P(NH,MC) + PEL(JrK)
^C)» S(NR,MC) t SEL(J,K)
002
...
IF
999
1000
II 00
1200
wP)U(b,999)
i ****«*«**«**«*«*««« «*****IN SETUPt)
, WRITE (*>, 1 100)
, HRlUU»f1200)U,(PU,J),J»l,J),tS(I»J),Jsl,3).R
-------
(FLECS VERSION Z?,U6) 1J.M4H-81 1JI29M3 PAGE 00002
0005« END
(FLECS VERSION ^^,
-------
(FLECS VERSION
13-MAR.8I
13129148 PAGE OOOOl
00001
00002
0000)
oooou
00005
00006
00007
00008
00009
OOnlo
OOOJt
00012
0001 J
OOOld
00015
00016
00017
00018
00019
00020
00021
0002?
00021
00020
00025
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
SUBROUTINE SHEARR(DEPTH, 050, STRESS, USTAR, V£L )
THIS SUBROUTINE CALCULATES BED SHEAR STRESS AND SHEAR VELOCITY FOR
A StfHHENt LAOEN FLOW, METHOD IS APPLICABLE FOR RESERVOIRS.
REF, HYDRAULICS OF SEDIMENT TRANSPORT BY W,H. GRAF, EO 8,«'?
FORMAL PAHAMgTERSl
06P1H • FLOW DEPTH (MgTEHS)
D5o • MEDIAN BED SEDIMENT DIAMETER (METERS)
STRESS • BED SHEAR STRESS (KG/M**2)
USTAR - SHEAR VELOCITY (M/SEC)
VtL - AVERAGE VELOCITY (M/SEC)
CALLED BYI HYDFLO, ICFLO
RHO . WATER DENSITY (KG(FORCE)/M»»J)
DATA RHO /1000,/
AKAPPA - KAHMAN CONSTANT
DATA AKAPPA /0,«/
USTAR»VEL/(l7,b«>*(ALOG10(OEPTH/(96.5*050)))*2,J/AKAPPA)
STRES3eHHO»USTAH**2,0/9,8
RETURN
END
(FLECS VERSION 22.U6)
158
-------
(FLECS VERSION 22,46)
19.MAK.flt 1<||27|15 PAGE 00001
00001
00002
OOOOJ
ooooi
00005
00006
00007
OOOOfl
00009
00010
0001 I
00012
000)3
0001U
00015
00016
0001 7
00018
ooop
00020
000*1
00022
0002J
00020
00025
00026
SUBROUTINE. SH£ARS(*LEN,ELEV, HRAO, PEl.EV,SLOPE,STRESS,USTAR]
THIS METHOD OF COMPUTING BED SHEAR STRESS AND SHEAR VELOCITY
IS APPLICABLE TO RIVERS AND STREAMS.
FORMAL PARAMETERS!
ALEN
ELEV
HRAO
PtLEV
SLOPE
STRESS
USTAR
SEGMENT LENGTH
ELEVATION OF THE CURRENT SEGMENT
HYDRAULIC RADIUS OF THE SEGMENT
ELEVAtlON OF THE PREVIOUS SEGMENT
BED SLOPE
660 SHEAR STRESS
SHEAR VELOCITY
CALLED BYJHYOFLO
G > RKAVlTr (M/S»*2)
RHO - nfcNSITY OF WATER (KG(FORCE)/M**3)
DATA RHO/1000,/
DATA G/9,801/
SLOPE « (PELEV - ELEV) / AL£N
STRESS » SLOPE * RHO * HRAD
USTAR • 3URT(G » SLOPE * HRAD)
RETURN
ENO
(FLECS VERSION 22,«6)
159
-------
(FLECS VERSION 22.16)
13-MAR-8I 13129:56 PAGE 00001
00001
00002
0000]
00004
00005
00006
OOOOT
00009
00009
00010
0001 I
00012
0001 3
000)1
00015
00016
oooi;
00019
00019
00020
00021
00022
00023
00020
00025
00026
00027
00028
00029
00030
00031
000)2
00033
00030
00035
00036
00037
00038
00039
o o () a o
oooai
OOOU2
OOOU3
00014
000«5
000«6
OOOH7
OOOU8
00049
00050
00051
00052
00053
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
SUBROUTINE S1LCLA(ARAR, B, BDIVi CC1N» OEUD, DELZ, DEPTH,
1 DENS, D8HR, ERODE, HRAD, J\, N8EO, COLD,
2 NELEM, POR, QH1N, OHOUT, 3C3HR, STRESS, V3ET,
3 XYSO, OEPO, ILAYR, 3D, SR, XNT,
n CROSEC, BWIO, ALEN, ECHO?, SCOUR)
THIS SUbROUTINE CoMPUTES THE RATE AND SOURCE TERMS FoF! THE
TRANSPORT OF 3RT (Jla2) AND CLAY (Jl«3)
INPUT PARAMETERS!
ABAR - AVERAGE AREA
8 • BED CONOIUONS
80IV • ST*N0»»l-> BED LAYER THICKNESS
c • WATER CONDITIONS
CCIN • CONCENTKA rioN or INFLOW
OELTO • TlME STEP (DAYS)
OELZ " STANDARD ELEMENT THICKNESS
DENS • DENSITY
DEPTH . DEPTH OF FLOW
DSHR • CRITICAL SHEAR STRFSS FOR DEPOSITION
EROOE - ERODAB1LITY, (KQ/M**2/SEC )
HBAD - HYDRAULIC RADIUS
Jt - »M SILT »2> CL*Y
NBEO • NIJMbf.R Of BED LAYERS
NELEM - NUMBER OF ELEMENTS
POR - POROSITY
OH1N • INFLOW DISCHARGE
gHoilT • OUTFLOW DISCHARGE
SCSHR - CRITICAL 3HEAK STRESS FUR SCOUR
STRESS - BED SMEAR STRESS
VSET • PARTICLE SETTLING VELOCITY
XYSO • THICKNESS OF TUP BED LAYER
OUTPUT PARAMETERSj
ILAYR • NO, OF HF.O LAYERS AFFECTED BY DEPOSITION AMD EROSION
SO - DEPOSITION RATE, (KG (PC J /M» * 3/DAY )
SR • EROSloN RATE, (KG[pC ]/M**3/OAY)
XNT . WEIGHT OF TOP BED SEDIMENT LA^ER, (KG/M**J)
OEPO • HEO DEPOSITION RATE ( KG f PC ) /M2/DA Y )
SCOUR - 8ED SCOUR RATE (KG IPC1 /M2/DA Y )
CALLED BYI JRANSP
CALLS: OEPCAL
INCLUDE 'SYIELMSIZ.PRM-
REAL K«FUNC,K£NS(3), DSHR(3<, EROOE(S), lLAYR(3)p
? UHIN(MXELtM), OHOUT(MxfLEH), SCSHRC3),
3 SD(MXELEM,6), SR(6), VS£T(3), XNT(3).
a CCIN(MX£LtM,MAXCON) , BWIO(MXELEM), OEPOC6J, SCOIlRl*)
160
-------
(FUECS VERSION 2?,«6)
uj2|Jt56 PAGE 00002
OOOS1
00055
00056
00057
00058
00059
00060
00061
00062
00061
00060
00065
00066
00067
00068
00069
00070
00071
00072
0007i
00074
00075
00076
00077
00078
00079
00080
00081
000«2
0006J
0008«
00085
00086
00087
OOOBA
00089
00090
00091
00092
0009J
00090
00095
00096
00097
00098
00099
001 00
OOtOI
00102
00103
00101
OOJ05
00106
OOJ07
001 08
00109
C
C
C
C
C
C
C
C
C
C
DATA SECDAr/66100./
XNTU1 )*0.
IF (NBED.GT.O)
t-,
8(NBED, J)/DEN5(}))
XNT(Jt) a XDIOP * B(NBEO,J1) * XY30
DEPO(Jt) a 0,0
OtPO(J2) s 0.0
SR(JI) * 0.0
SHfJ2) * 0.0
SCOUR(Jl) * 0,0
SCOURCJ2) » 0,0
RS s 0.0
CS * 0,0
VOLUME a CROSEC * ALEN
DO (IX « l.NELEM)
, SDdx, Jl ) « 0.0
, S0(IX,J2J * 0.0
...FIN
ILAYR(.Jl) a 0
IF(fcCH07)
, IF (STRESS ,LT. OSHR(Jl))
SEDIMENT DEPOSITION
ILAYR(Jl) x -1
AVGC * 0.0
ruiu * o.o
00 (IX » 1,N£LEM)
IT
AT
IS
OR
IMPLICITLY
NEAR UNIT*
ASSUMED THAT A
IS EMPLOYED IN
DOWNSTREAM COURANT NUMBER
THIS ANALYSIS
TOTQ * TOTO * (OHIN(IX)*UHOUTCIX})/2,
Xt Jl )+CCIN(IX*l,Jl))/U,
* ( 1 . 0- ( S THE SS/DSHR ( J 1 ) ) )
ALEN / VOLUME
AVGC s AVGC / TortJ
OEPO(Jl) s VSET(J)) • AVGC
RATE » DEPO(Jl) * SWIOt1) <
DO (K a 1
SD(K,JJ )
VOLK a AUAR(K) * OELi
SEDaQHlN(K)*(CCIN(K,Jl)*CClN(K+l,Jl))/1,tf5HOUT(K)*COLD(KfJl)/2.
CONTaQHIN(h)*(CCIN(K>J2)»CCIN(K»l»J2))/«.+QHOUT(K)*COLO(K,J2)/2,
HATEK a RATE * VULK
SO(K,J2) * OATEK * CONT / SEO / VOLK
OEPO(J2) » DEPO(J2) » 5D(K,J2) * VOLK / 8WlO(l)
..FIN
,.FJ1
(STRESS ,GT, SCSHR(J1).ANO. NpED. Gf , 0)
161
-------
(FLEC3 VERSION 22.46)
13.MAK.8I 13129156 PAGE 0000)
OOJIO
001 11
00112 C
00113 C
00114 C
0011? C
001 16
00117
001 IB
00)19
00120
00121
00122
00123
00124
00125
00126 140
00127
00)28
00129
00130
00131
00132
OOI3J
00] 34
00135
00136
00137
00138 C
00139
00110
00141
00142
00)43 C
OOJ41 155
00)15
00)46
00)17
00148
00149
00150
00151
00152
00153
00154
00155 C
00156
00157
00)5«
t
•
,
«
t
t
•
.
•
i
.
.
i
.
*
.
.
.
.
*
•
•
•
.
.
.
.
.
.
.
.
»
t
.
t
.
.
i
.
»
.
»
.
,
«
«
• .
, RS * ERODE(Jl) * SECDAY * (STRESS / 3CSHRCJ1) - 1.0)
, ILATR(Jl) B (i
*
. COMPUTE AVAILABILITY OF COHESIVE SEDIMENT IN BED LAYERS,
, MAXIMUM NUMBER OF LAYERS SCOURED IS RESTRICTED BY SAND SCOURING,
•
, RS » R3 * OELTO
. *H£N (,NOT,(R3 ,GT, XNT(Jl) .AND, ILAvR(l) ,CT. 0) )
, , KS « AMIN1 (RS,XNT(Jt ) )
. , CS « RS * B(NBED,J2)
, . XNTUI ) » XNT(J1)-RS
. ...FIN
, ELSE
, , FACs8(NB£D,J2)
. . HSUSP * RS
. . KS = 0.
. . ILAYR(JI) * ILAYRCJD + 1
. . N8 « NBEO»IL*YR{ Jl )
, , RSUSP o PSUSP-XNUJ1)
, . HS » RS + XNHJl)
. , CS » CS f FAC*XNT(J1 )
XNHJt ) « 0,0
. . FAC * 0.0
, , IF (NB.NE.O)
. . , XNOs(t,0-POR)/(B(NR,l)/OENS(l ; *B (NB, 2) /DENS ( ,2 } +B (NB , 3 } /DENS ( 3 ) )
. . . XNT(JI) * BDIV * B(NB,J1) * XNO
. , . FAC * R(NB,J2)
FIN
• •
. . IF (ILAYR(J) ),EO.ILAYR(t J) GO TO 155
, IF (RSUSP. GE, XNT(J1 ), AND. ILAYR( 1) .GT.ILAYR(Jl) )
. . , GO TO 140
• « t • f * 1 "
t (
, , CONTINUE
, . DEL * AMIN1 (RSUSP, XNT Ut ))
. , KS a RS + DEL
. . CS » CS t DEL * FAC
, , XNT(Jl) * XNT(Jl) - DEL
. ...FIN
...FIN
SCUUK(Jt) « HS / OELTO
SCOUMJ2) « CS / DELTD
SH(Jl) * RS / OELIO » BKJO(l) » ALEN / VOLUME
S«(J2) * CS / OELTO * BWlD(l) * ALEN / VOLUME
,F1N
RETURN
END
(FLECS VERSION 22.46)
162
-------
(FLECS VERSION 22.46)
13,MAR»8I I3I30H2
00001
00001
01)00?
00003
00004
00005
00006
00007
00008
00009
00010
0001 1
00012
0001)
00014
00015
00016
0001 7
00018
00019
00080
00021
0002?
00025
00024
00025
00026
00027
00028
00029
000 JO
00031
0003?
0003J
oooja
00035
00036
OOOJ7
0003H
oooiR> /
DATAYFS/'Y'/
DATA wRIStG/, FALSE./
WHITE(fl, 1 )
163
-------
(FLECS VERSION 22.4fe)
m3oii2 PAGE 00002
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
0006R
00069
00070
00071
00072
00073
00074
00075
00076
00077
00076
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00091
00094
00095
00096
00097
00098
00099
00100
OOlOt
00102
0010J
00104
00105
00106
00107
00108
00109
C
C
C
C
C
C
C
C
C
C
C
RE AD (8, 2) NCHR, (lNPFIL(I)»Ial,NCHR)
INPF IL(NCHR + 1 ) « 0
OPEN (UNI Taj , NAME" t NPF I L, TYPE" 'OLD1, READONLY)
WRITEC8»3)
REAO(8,4)AN3WER
IF (ANSWER ,EO, YES) ECHO « .TRUE,
WHITE(B,6)
READ(8,«)(FNAME(I),-J«1,E9)
CALC FOCOUE(FNAME,BASE,JSEG,FTYPE,DEV,GUIC»UUIC)
WRITE(8,7)
REAO(8,4) SMETH
WHITE ( 9» 8)
REAt>(8,4) ANSWER
IF (ANSWER. EO. YES) S»VECHU) s .TRUE.
WRITE(«,<»)
REAO(8,M) ANSWER
IF (AN3«EH. EO, YES) SAVECHt?) a .TRUE.
WHITt(8, 10)
READ(6,4) ANSWER
IF (ANgwER.EQ.YES) 8AVECH(3) a .TRUE.
WHITE(8,11)
R£AD(8,4) ANSWER
IF (ANS"£H.EQ,YES) SAVECH(4) t .TRUE.
WRITE (8,12)
REAO(8,«) ANSWER
IF (ANSWER. EO. YES) S*VECH(3) « .TRUE.
WRITE (6,13)
REAO(8,4) ANSWER
IF (ANSWEK.EO.YES) SAVECH<6) a .TRUE.
WRITE (8,14)
REAOtB,") ANgwEH
IF (ANswER.EO, rfcs) gAvECH(7) * .TRUE,
WHI TE (6,15)
REAO(6,«) »NSWER
IF (AN8WEH.EO. YtS) WRTSEG a .TRUE,
WHEN (WRT5EG) JSEG(l) B 0
ELSE
. WRITE (8,16)
. R£AO(8,17) (JSEG(J),J«1,5)
« t , F I N
CLOSE(UNIT«8)
NFRST * 1
ISTNTel
INFLOa?
OUlFLO«i
164
-------
(FLECS VERSION 22.«6)
I3-M4K-8) 13lJO|l2 PAGE OOOOJ
00110
001 1 1
00112
OOI1J
00) 19
001 IS
001 16
00117
00118
001 19
00120
00121
00142
00123
00124
00125
00126
00)27
00128
OOJ29
001 JO
00111
00132
00)33
0013«
00135
00136
00137
00138
00139
OOlUO
00141
C
C
1
2
3
U
6
7
8
9
10
1 1
12
13
IU
15
16
17
C
OPEN ( UN 1T»2, NAME t'OUMMY.DTl •, TYPE e< NEW' ,FORM" ' UNFORMATTED ' )
OPEN(UNlT»J,NAMt«iDUMMYtDT2<,TYPE«»NEw»,FORM»«lJNFO»MATTEO' >
OPEN (UN IT«fl,NAME»'HY'OKOLOGY,TMP', TYPE* 'NEW », FORM* »UNFGMA T TED ' )
OPEN(UNlls7,NAME««TRIBUTAKY.TMPf , TVPE« ' NEW ' ,FORM» "UNFORMATTED » )
OPEN(UNIT»6,NAME«i8ED,LST' , TYPE* 'NEW')
RETURN
»* FOHMATS **
FORMAH 'JtNTER NAME OF INPU1 FILE>')
FORHAI (U.30A1)
FORMAK'JUO YOU WANT THE INPUT FILE ECHOED (Y OK N)>")
FURMAJ (29A1 )
FORMAT( 'HNTER BASE FILE NAM£>')
FORMAir" WHICH SAND CAPACITY METHOD 13 TO BE U3ED7'/
1 'JENTER T (TOFFALETTI) OK C (COL8Y)>i)
FORMAlClOO YOU WANT SEfATRA HEADINGS ECHOED (Y OR N)?')
FORMAtt'lOO YOU WANT INFLUENT CONCENTR A T IONS ECHOED (Y 0" N)?l)
FORM4TCJOO YOU "ANf ELEMtNT MAfHICES ECHOED (Y OR N)7')
FORMATf'IDo YOU WANT GEOMETRY AND CONCENTRATIONS ECHOED (Y OR N)t,
1 'T')
FORMAM'SUO YOU WANT CONCENTRATION ECHOED BEFORE AND AFTER ROSFLO'
1 , ' (Y OR N)?l )
FORMATMJOO YOU WANT jCOUR/DEpOSl T I ON TO OCCUR? (Y OR N>«)
FURMAH'JDO YOU WANT COMPLETE SCOU«/DEPOS I T ION INFORMATION1,
1 « RECORDED? (Y OH Nj ' )
FORMATf'IDo YOU WANT COMPLETE ECHO* INFORMATION FoR ALL'»
1 « SEGMENTS? (Y OR N) ' )
FORMAT(i*FOR WHICH SEGMENTS DO YOU WANT COMPLETE ECHO**,
1' INFORMATION? (MAXIMUM OF 5)'l
FURMAH5I5)
END
(FLECS VERSION 22,06)
165
-------
(FLEC3 VERSION 22.«6)
13-MAR-si
PAGE ooooi
00001
00002
00003
00004
00005
00006
00007
00008
0000"'
00010
00011
00012
00013
oooi a
00015
00016
000 1 7
000 1 8
0001"'
00020
00021
00022
00023
0002«
00025
00026
00027
00028
00029
00030
00031
00032
00033
0003«
00035
00036
00037
00039
00039
00040
00n« |
00042
00043
0001U
oooas
00046
OQOu7
OOOU8
00049
00050
00051
00052
00053
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
SUBROUTINE TOFFALULEN, 050, G9If HRAO. rjTOT, SLOPE? TEMPR, VOL,
1 VSET, GSU, GSM, GSL» GSB, YU, YM, YU
THIS SUBROUTINE USES TOF FALE T T I • S METHOD TO CALCULATE THE CAPACITY
OF THE FLOW TO TRANSPORT SAND. A SUMMARY OF THIS METHOD CAN BE
FOUND IN IH£ ASCE 1975 EDITION OF "SEDIMENTATION ENGINEERING"
PAGES ?U9 - 21J.
FORMAL PAKAMETEKSi
ALEN SEGMENT U.NGTH
Dio MEDIAN BED SEDIMENT DIAMETER (METERS)
GSI TOTAL CAPACITY OF THE SEGMENT (KG/OAY/MJ
HRAD HYDRAULIC RADIUS
QTOT TOTAL FLOW WITHIN THE SEGMENT
SLOPE ENEKG* OR RIVER BED SLOPE
TEMPR WATER TEMPERATURE
VOL VOLUME
VSE1 SETTLING VELOCITY
CALLED 6YJ 3ANO
REAL K«FUNC,K7 . 0,667 » TMPR
TT=1,10 * (0,051 » 0,00009 . fMPR)
ZIsVSETU) * CONST! * V 1 (CZ * FMHAD * SLOPE)
IF(ZI ,LT ,CNv> ZI»1 ,5 * CNV
THE MANNING-STRICKLER EQUATION is USED HERE TO
DETERMINE THE HYDRAULIC RADIUS COMPONENT DUE TO
GRAIN ROUGHNESS <«'). TAKEN FROM THE 1975 ASCE
"SEDIMENTATION ENGIMEERING'iPG, 128,
166
-------
(FLECS VERSION 22.«6)
ii.MAH.8t I3i30i3<> PAGE 00002
0oo5"
00055
00056
00057
00058
0005'
00060
00061
00062
00063
OOflb'J
00065
000b6
00067
00068
00069
00070
00071
00072
0007J
00070
00075
00076
00077
00076
00079
00080
00081
00082
00083
0008U
00085
00086
00087
00089
00089
00090
00091
00092
0009J
00094
00095
00096
00097
00098
00099
00100
00101
00102
0010)
ooiou
00105
00106
00107
00108
00109
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
SUBSTITUTIONS »KE MADE IN EQUATION 2,I'M FOH SHEAR VELOCITY
AND MSU9JS, THE FORMER IS REPLACED BY EQUATION 2,1«2. AND
THE LAHEK BY 0(SUH)65.
RPRIM£«((V»*l,5) » (065 *»0,25) / (SLOPE ** 0,75))
USTARr(HPHlME • SLOPE * 32,2) ** 0,5
AFUNC»
...FIN
OCZUS1 .0 *
/ AC
CNV
OCZMsJ.O f CNV •
1 .5 *
ZI
ZI
OCZL'1.0 t CNV - 0,756 * ZI
ZINV«CNV - 0.756 * ZI
ZOs. 0.756 * ZI
ZP*0,2"0 * ZI
ZQ*0,5 * ZI
CLI HAS 0EEN MULTIPLIED 8Y 1.0E30 TO KEEP IT F»0«
EXCEEDING THE COMFUTEH OVERFLOW LIMIT
CLI»CONST2 * OCZL * (V ** 2,333) / FHRAQ ««
1 UTT • AC • K« * FDIAM) ** 1,667) / (J.O f CNV) /
i ({FHRAO / 11.2U} ** (/N) . (2,0 » FfjIAM) »* OCZL)
Pls(2.0«FOIAM/FHRAD)**Zo
C2D»CLI » PI / l.OEtJO
CHECH TO SEE IF THE CALCULATED VALUE IS REASONABLE
(< 100. 0), AND ADJUST IT IF IT IS NQT.
IF(C20. Of. 100.0) CLIs l,OEt2H / PI
CMI HAS BtEN MULTIPLIED RY 1.0E30 TO KEEP IT
EXCEEDING ' THE COMPUTE* OVERFLOW LIMIT
167
-------
(FUCS VERSION 22, Ufc)
13l30|36 PAGE 00003
00110
00111
001 12
0011J
OOtH
001 15
00116
001 17
001 18
001 19
00120
00121
00123
00|23
00121
00125
00126
00127
00126
001J9
00130
00111
00132
00133
00131
00135
00136
00137
00139
0013?
001UO
001 Ul
00102
00103
001 11
00115
OOH6
001 «7
001 18
00119
001SO
00151
00152
00153
00151
00155
1)0156
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
P(SU8)l APPEARING IN EQUATIONS 2.236H, J, K, L, M, AND
N IS IH£ WEIGHT FRACTION OF TOTAl SAND THAT THE I-TH
size FRACTION CONTAINS. SINCE WE ARE MODELING ALL
SAND AS * SINGLE SIZE FRACTION -- PCSUBH » 1,0. AND
HENCE DOES NOT APPEAR IN THE MODEL EQUATIONS,
CMIs(l3.g * CL1 * (1,0 » CNV) * V * (FHRAD ** (ZM))
CALCULATE TRANSPORT CAPACITY OF THE UPPER LAYER
FDllsFHRAU X 11.21
FO?5oFHRAU / ?,5
GSUe(CMl * (FOll ** UP)) » (FD25 ** (ZQ)J *
1 {FHHAO tt (UCZU) - (FD25 •» (oCZU)))) ' (UCZU * 1,0 E+30)
CALCULATE THE CAPACITY OF THE MIDDLE LAYER,
GSMe(C*l * (FDtl *» (ZP)) * UF025 **(OCZM)) «
1 (FOll »» (OCZM)))) / (OCZM * t.OE+30)
CALCULATE THE CAPACITY OF THE LOWER LAYER
GSLs(CMI * ((FD11 »* (ZN)) - (C2.0 * FDIAM) .» (OCZL»)))
1 / (OCZL * 1,OE»30)
CALCULATE THE CAPACITY OF THE BED LAYER
G38=(CMJ * 1(2,0 * FDIAM) ** (ZN)))/1.0E+30
TOTAL CAPACITY OF THE SEGMENT (G3l HAS UNITS OF TONj/DiY/FT)
GSI»GSU + GSM + GSL *
CONVERTING TO KG/DAY/H
GSu « G5U * CONST3
GSM * GSM * CON3TJ
GSL « G5L * CONST3
GS8 * fiSB • CONST3
YU * HRAD / 2,5
YM e HRAO / 11,21
YL = 2.0 • D50
GSltGSl * CONST3
RETURN
END
(FLECS VERSION 22,16)
168
-------
(fLECS VERSION 22.a&)
3.APR-8I I7i26t56 PAGE ooooi
oooot
oooo?
00003
ooooa
00005
00006
00007
00(108
00000
00010
OOOH
00012
00013
00014
00015
00016
00017
0001A
00017
00020
00021
00022
0002)
ooo2a
00025
00026
00027
00024
00029
00030
00031
00032
00033
00030
00015
00036
0003?
00036
00039
00010
00041
00042
00043
00044
00045
00046
00047
00040
00049
00050
00051
00052
OOOS3
C
c
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
SUBROUTINE TRANSP(TERROR,PCOEF,8*10,AW ID,ABAR,DEPTH,OUDC,ECH04,
1 CRUSEC,ECH07, ECH08, ECHO?)
THIS ROUTINE SOLVES THE MASS TRANPORT EQUATIONS 8Y AN IMPLICIT
FINITE-ELEMENT METHOD. A CRANK-NICHOLSON METHOD IS USED TO
APPROXIMAtE THE SOLUTION THROUGH TIME,
VARIABLE
ALEN
AREA
8
BOIV
BED
C
CCIN
COLD
CH03EC
DECAY
DELTH
Dt'LZ
DFZ
D1AM
OSHR
050
ERODE
FERHOR
HRAO
KAY1
KAYg
NBEO
NELE*
PCOEF
POR
QHIN
OHOUT
(3V
SCSHR
SLOPE
SMgTH
SR
STRESS
TEMPR
VOL
VSEt
XYSO
DEFtNITIONSl
SEGMENT LENGTH
ELEMENT AREAS
dEO CONDITIONS
STANDARD BED LAYER THICKNESS
BED THICKNESS
WATER CONDITIONS
CONCENTRATION OF INFLOW
CELL-CENTERED CONCENTRATION
TOTAL CROSS-SECTIONAL AREA, M>*2
FIRST ORDER DECAY
TIME STEP (SECONDS)
STANDARD ELEMENT THICKNESS
DENSlTY
DIFFUSION COEFFICIENT
DIAMETER
CRITICAL SHEAR STRESS FOR DEPOSITION
MEDIAN BED SEDIMENT DIAMETER (M)
ERODABILITY
FATAL ERROR FLAG (L*l)
HYDRAULIC RADIUS
LIGHT EXTINCTION COEFFICIENT OF WATER
LIGHT EXTINCTION OF SUSPENDED SEDIMENT IN HATER
NUMBER OF BED LAYERS
NUMBER OF ELEMENTS
FIRST TERM OF THE PHuTOLYSlS RATE OF CHANGE EQUATION
POROSITY
INFLOW DISCHARGE
OUTFLOW DISCHARGE
VERTICAL DISCHARGE
CRITICAL SHFAR 5IRESS FOR SCoUR
ENERGY OR RIVER BED SLOPE
CONTROL VARIABLE TO SELECT THE METHOD TO BE USED
WHEN COMPUTING THE SAND CARRYING CAPACITY,
SEDIMENT EROSION RATE
BED SHEAR STRESS
WATER IEMPERATURE
VOLUME
PARTICLE SETTLING VELOCITY
THICKNESS OF TOP BED LAYER
CALLED BYI SERATRA
CALLSj BEODK, BEDHIS, COLAPS, COMB, DISOLV, FALL, PARTIC. SAND.
SEOIME, SETUP, SILCAL, TRISOL
INCLUDE
169
-------
CFLECS VERSION 22,«6)
3.APR-81 17l26|56 PAGE 00002
00051
00055
00056
00057
00058
00059
00060
00061
00062
00065
00061
00065
00066
00067
00068
00069
00070
00071
00072
00073
00071
00075
00076
00077
00078
00079
00080
00081
00082
00085
00081
00085
00086
00087
00088
00089
00090
00091
00092
0009*
00091
00095
00096
00097
00099
00099
00100
00101
00102
00103
00101
00105
00106
00107
ooioa
00109
I 0 10
1015
1000
LOGICAL*! FERHOH, ECMOU, ECH07, ECHQB, ECH09
DIMENSION A(MXELEM), ABAR(MxELEH) , 00(1(1), ILAYR(3), P(HX£LEM,3),
1 R(MXELEM), 3(HXELEM,5), SD(MXELEM,6),
2 SR(6), XNf(3), Z(MXELEM), BWlD(MXELEM), AWID(MXELE*)
3 ,OLUC(MXELEM,MAXCON), DEPOC6), SCOUR(6), 8EDSO(3)
INCLUDE 'TRANS,COM»
DATA EPbl/l.OE-30/
MP1 « HELEM » 1
Mil * NELEM . 1
*** PERf-OPM CALCULATIONS OVER THE TIME STEP DELTO •**
DELTD » OtLTH / 86100.
DO (J«I,HAXCON)
DO (L * l,hPl)
. R(L) » 0,0
. 00 IN s 1,3)
, . S(L,N) a 0,0
, , P(L,N) » 0.0
DO (I«1,NELEM)
, RCDZ i 0,0
. IF (I ,EQ. 1)
. . SELECT(J)
... (1)
CALL S«NO (ABAR, ALEM, AREA, 8, 8DIV, CCIN, OELTD, DELZ.
1, . , , , DENS, 050, HRAD, NBEO, NELEM, POU, QHJN, OHQUT,
2, , . , , 3CSHR, SLOPE, SMETH, STRESS, TEMF'R, VSET.VOL, XYSO,
3, , , , , OEPO,1LAYR, SO, SR, XNT, COLO, C,
«,,,,, CRUSEC,B»<10,ECH07, SCOUR)
. . . . , 1F(ECH08)
, WHIIE(6,1000)
...,,. «M£N(jLArRt i) ,LT .0)
, , WRITEC6,1010) J,ILAYR(1),DEPO<1),DEPO(«),SD(1,1),XNT(1),
1, ,.,.,. (50(11,«),II«l»NELEM)
...... ...FIN
ELSE
,..,,,. w»ITE(6, 1015) J,ILAYRO),SCOURC1),SR(»),XNT(1),SCOUR(a),SR(fl
FIN
,,.,,, FORHAT(' Ja',l2,i ILAYRs',12,' DEPOs',E15,7, ' DEPO( + 3) = ',
1, , , . , , El5.7,i SD»',E15.7,' XNy«•,E15.7,5X,'SO(I,*3)s ' /
2, ..... (5X,8E13.7))
,.,.,, FOPHATC J»',J2,' ILAYRa',12,' SCOUH='„E15,7,' SRa',E15.7,
1, ,.,,,' XMT» ' ,E15,7, i SCOURC +J)s",E15,7, ' SR( + 3) = ',Ei5,7)
..... ...FIN
IOIN TPANSP FOLLOWING SANO>)
170
-------
(FLECS VERSION 22,«6)
J-APR-81 17126(56 PAGE 00003
oono FIN
00111 .... (2)
00112 CALL SILCLA(ABAR, 8, BOlV, CClN, DELTO, DELZ, DEPTH, DENS,
00113 1 DSHR, ERODE* HRAD, Z, NBED,COLO,NELEM, POR, QHJN,
ooiia a. .... QHuur, SCSHR, STRESS, VSET, XYSO, DEPO, ILAYR,
00115 3, .... SO, 3H, XNT, CROSEC, BWID, ALEN, ECHQ7, SCOUR)
00116 IF(ECH08)
00117 ,,,,,. "RITE(6,1030)
00118 ...... WHEN(ILAYR(2),LT.O)
OOllt , WRITE(6,1010) J,ILAYR(2)»DEPO(2),OEPU(5),SDn»2),XNT(2),
00120 I. ...... (80(11,5),II«1,NELEM)
OOm FIN
00122 .,..,, ELSE
00123 ....... WRITE(6»1015) J,ILAYR(2),SCOUR(2),SR(2),XNT(2),SCOUR(5),SB(5)
0012« FIN
00125 ..... ...FIN
00126 1020 « . . . . FOHMAT(»oIN TRANSP FOLLOWING SILT')
00127 .... ...FIN
00128 .... (3)
OOU9 CALL 3ILCLA(A9AR, 8, 80lV, CCtN, DELTD, DELZ, DEPTH, DENS,
00130 1 OSHR, EROOE> HRAp, 3, N8EO,COLD,NELEM, P0», QHJN,
00131 2, .... QHQUT, SCSHR, STRESS, VSET, XYSO, OEPO, RAYR,
00132 3, .... SO, SR, XNT, CROSEC, BWID, *LEN, ECHQ7, SCOUR)
00133 IF(ECHOS)
0013« ...... «RITE(6,1030)
00135 ...... WHEN(IL*YRU) ,LT,0)
00136 ....... WRITE(6, 1010) J,ILAYRm,OEPO(3),DEPO(6),3D(l»3),XNT(3),
00137 1. ...... (80(11,6),1I«1,NELEM)
00138 ...... ...FIN
00139 ELSE
00140 ,.,,... WRITE(6,1015) J,ILAYR(3),SCOUR(3)»SR(3),XNT{3),SCOUR(6),SR(6)
ooiui ...... ,.,FIN
001«2 .... * ...pIN
001«3 1030 ..... FOHMATCoIN TRANSP FOLLOWING (SILT) CLAY')
001«a FJN
001U5 FIN
00106 . ...FIN
00107 . CONDITIONAL
OOU8 . . (J ,LE. 3)
001U9 . , . CALL SEDIME(ABAft, ALEN, CClN, DELTD, I, J, NELEM, OHIN,
00150 1 ... OMOUT, 0V, SO, SR, ALFA, BETA, VEL1, VEL2,V3ET,
00151 2 ... OELZ, BwlD, Awlo, DEPTH, BETA1, BLTA2)
00152 , . , IF(tCHOB) wfllTE(6,)500)
00|53 (500 , . , FQRHAU'OIN TWAN3P FOLLOWING SEDIMg')
OOtS« FIN
00155 . . (J ,GE. U .AND, J ,LE, 6)
00156 . . , C*LL PAHIIC(ABAH, *LEN, B, C, CClN, COLD, DECAY, DELTD, DFZ,
00157 i ... I, J, NBED. NELEH, OHOUT, OHIN, gv, SORBK, 30, SR, ALFA, BETA,
00158 2 ... VEL1, VEL2.VSET, DELZ, DEPTH,B«ID,AwID,
00159 3 ... BtTAl, 8ETA2)
00160 , , . IF(tCHOB) WRITE(6,1510)
OOjbl 1510 , , , FORMAM'OJN TRANSP FOLLOWING PARTIC')
00162 FIN
00163 . , (J ,EQ. 7)
0016U , . , CALL OISOLV(AHAR, B, BDIV, C, CClN, COLO, DECAY, DELZ,
00165 t ... DELID, DENS, DIAM, I, MYt, KAY2, NELfiM,N8ED,
171
-------
(FLECS VERSION 22.U6)
3,APH-ai m26i56 PAGE ooooi
00166
00167
00168
00169
00170
00171
00172
00173
07
00198
00199
00200
00201
00202
0020J
oo?na
00?05
OU206
00207
0 0 ? 0 ft
00209
00?! 0
00211
2. .
3. .
• •
1 520 . .
* •
* •
ff •
(530 , .
C , .
c . .
• •
• *
• *
1 •
• •
t •
• *
• •
• *
• ft
f •
• *
t •
• *
2500 . .
1. .
ff •
2000 , .
ff 9
C . ,
ff •
1 *
t *
1 * •
c
. . PCOEF, POR, OHIN, QHOUT, ov, SOHBK, ALFA, BETA, VELI,
, , VEL2, BETAl, HETA2, OEPO, SCOUR, BEOSO, XYSO, AREA)
, IF(ECHOfl) WRITE(6, 1520) (BEDSD ( 1 1 ) , ll*\ , 3 )
. , FORMA|-('OIN TRANSP FOLLOWING DISOLV I3EDSDU-3) S|,3F18,7)
. ...FIN
, , ,F IN
IMECH08) WHITE(6,1530) VEL(,VEL2
FORM»T(' VELt «',E15.7,10X, 'VEL2 »',E15,7)
*** CONSTRUCT THE FINITE ELEMENT MATRICES FOR EACH LAYER ***
"MEN (J.LE.3) DIFUSE » l)FZ(J)
ELSE DIFUSE » DFZ(J-3)
00(1) OIFUSE/(OELZ«OELZ)
DU(2) 6.0 * OELZ
00(1) (VEL2 • «, * VELI) / 00(2)
00(«) (VELI + 2, * VEL2) / 00(2)
D0(5) (2. * VELI * VEL2) / 00(2)
OU(fc) (a, * VEL2 • VELD / 00(2)
D0(7) ALFA
00(6) BETA
D0(9) 8ETA1
00( 10) a BETA2
IF (ECHO«)
, WRITE(6,2500)
. FORMAT ( "0 II,lX,lJl,6X,'DDK*l,10)
, FORMAT (U, 13, tX,Il,U, 10(1 PE 12, «))
...FIN
VEL « VEH/DELZ
IM1,EQ,NEL£M) VEL«"VEL2/OELZ
WIDTH • AMlD(I)
CALL SETUPd, 00, P, S, R, VEL,NELEM, ECHOU , W IOTH)
.FIN
f HALPO « OELTD / 2.0
. 00
* •
f t
ff •
* ff
f •
* •
• *
• •
(L > 1.MP1)
H(L) « R(L) » OELTO
Z(L) = C(L,J)
00 (N » 1,3)
. PI 2 S(L,N)*HiLFD
. P9AR 3 P(L,N) + PI
. S(L,N) = P(L,H) - Pi
P(L,N) » P8AR
...FIN
00212
00213
0 0 ? 1 «
00215
OOP16
00?17
00218
00219
...FIN
IF
5760
E COMB***
COMB**•
TE(6,5761
00221
1st ,MPl)
172
-------
(FLECS VERSION 22.a6)
J-APR-81
17126156 PAGE 00005
00222 5762
00223
,0022"
00225
00226 6001
00227 C
00228 C
00229 C
00230
00231 C
00232
00233
00?34
00235
00236
00237 C
00239
00239 (
00240 1
00241 2
002u2
002«3 1
00244
002<»5 C
002"6 C
00247 C
00246 C
00?49
00250
00251
00252 6200
00253
0025«
00255
00?56
00257
00258
00259 6300
00260
00261 f
. FORMATC1X,I2, 1P8£16.5)
...FIN
CALL COHBfMPj, S, Z, R)
IF JECH04) WM1TE(6,6000)
-------
(FLECS VERSION 22.16)
13.MAH.81 1JI32IOO PACE 00001
00001
00002
0000*
00001
00005
00006
00007
00009
00009
oooio
OOQ1 1
00012
0001J
OOOH
00015
00016
OOOlT
00018
00019
00020
00021
00022
0002J
00021
00025
00026
00027
00029
00029
OOOJO
00031
00032
OOOSJ
00031
00035
00036
00037
0003S
00039
OOO^O
oooal
0 0 0 " 2
000" 3
00010
Q 0 ft ^5
0 0 0 ** b
OOOU7
00019
00050
00051
00052
00053
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
SUBROUTINE TRBDAT (ECHO, HLDERR, N£LEM» NTRlBS, NUMERfl,
1 SIMLEN, TR80PT)
THIS ROUTINE IS RESPONSIBLE FOR READING AND PROCESSING THg
TRIBUTARY INFLOW MASS FLUX DATA, THE DATA IS READ FROM
THE INPUT STREAM (LUN 1) AND WRITTEN TO " THlBUT ARY . TMP" (LU
FOR USE DURING THE SIMULATION,
FORMAL PAHAMETERSi
ECHO • LINE PRINTER ECHO OPTION CONTROL VARIABLE (L*l
HLDERR - HOLDING ARRAY FOR ERROR NUMBERS (BYTE)
NELE* • NUMBER OF VERTICAL ELEMENTS
NOTEl NELEM IS LATER REDEFINED IN HYOOAT
NTRJBS « NUMBER OF TRIBUTARIES (0 OR 1)
NUMERH • NUMBER OF INPUT ERRORS
3IMLEN . SIMULATION LENGTH (3F.CQNOS - 1*1)
TK8UPT - TRIBUTARY INPUT CONTROL VARIABLE
CALLED BY| SERATRA
CALLSl PUTERR
INCLUDE 'SYlELMSIZ.PRMi
BYTE HLDERRUOO)
INTEG£W»2 TRBOPT
INTEGER*1 ENOTIM,PRETIM,SIMLEN
LOGICAL*! ECHO
DIMENSION CTRB(MXELEM,MAXCON)
... .TRIBUTARY INFLOW MASS FLUX
FIRST RECORD,..,.
CoL. 1- 5, ..NTRIBS,. ..NUMBER OF TRIBUTARIES to OR l
6-10. ,.TR«OPT... .TRIBUTARY INPUT OPTION
»o> THE USER WANTS THE MODEL
DISTRIBUTE THE MASS FLUX
THE ELEMENTS,
•l; THE USER WILL SUPPLY T^E
MASS FLUX VALUES FOR EACH
REWIND 7
^
ADO, *) NTRIBS, THHOPT
IF (ECHO) WHITE «»,*) NTPIBS.TRPOPT
IF (NTFMBS ,GT. 0)
. IF (ECHO) WRITE (fe, a)
N 7)
)
)
TO
THRU
ELEMENT
174
-------
(FLECS VERSION 22,86)
13.MAR.8i i3i32ioo PAGE 00002
0005«
00055
00056
00057
00058
00059
00060
00061
00062
00063
00061
00065
00066
00067
00068
00069
00070
00071
00072
0007S
0007"
00075
00076
00077
00078
00079
00080
00081
00082
00083
00080
00085
00086
00087
oooee
oooe9
00090
00091
00092
00093
0009U
00095
00096
00097
00098
00099
ootoo
00101
00102
00103
0010«
00103
00106
00107
00108
00109
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
t
2
3
a
REPEAT UNTIL (ENOTIM ,EQ. -9999)
CARD 11-8
COL. 1-10...ENDTIM,,,.ENDING TIME FOR THE DATA THAT FOLLOWS,
A VALUE OF -9999 TERMINATES THE DATA
(SECONDS).
READ(1,6) ENOTIM
UNLESS (ENDT1H .EG. -9999)
RECORD Two..,.,,.TRIBUTARY MASS FLUX AND DEPTH
********** CAUTION **********
THE MASS FLUX UNITS ARE DIFFERENT FROM THOSE OF
INITIAL WATER AND UPSTREAM WATER CONCENTRATIONS,
PADIONUCUDE IS PC/SEC, PESTICIDE IS KG/SEC,
SEDIMENT IS KG/5EC
********** CAUTION ****.******
COL, 1-IO...CTR8(1,1),,.MAS5 FLUX OF SAND (KG/M**3 ) *(M**3/SEC)
11-20, ..CTRBft,2)...MASS FLUX
21-30, ,,CTR8(1,3),,.MASS FLUX
31-«0,,,CTRB(J,«)..,MASS FLUX
WITH SAND
al-50. ,,CTRB(J,5),,.MASS FLUX
WITH SILT
51-60,..CTRB(J,6),,.MASS FLUX
WITH CLAY
«»!-70,,,CTRB(J,7) , . .MASS FLUX
OF SILT
OF CL*V
OF CONTAMINANT ASSOCIATED
(PC/KG)*(XG/M**3)*(M**VSEC)
OF CONTAMINANT ASSOCIATED
OF CONTAMINANT ASSOCIATED
OF DISSOLVED CONTAMINANT
(PC/M**3)*(M*«3/3£C)
WHEN (TRHQPT ,EO. 0) N s 1
ELSE N B NELEM
DO
-------
(FLECS VERSION 32,H6) 13-MAR-61 I3l32l00 PACE 00003
OOltO 3 'WIIH 3ILP,6X,'HITH CLAY',«X,'DISSOLVED CONC'.ZX,
OOJ11 « '(M..3/5EC)')
00112 5 FORM»T(2X,110,2X,7UPE12,5,SX),1PE12,5)
001IJ 6 FORMAT(IJO)
00114 C
00115 END
(FLECS VERSION 22.06)
176
-------
(FLECS VERSION 22.16)
1}.MAR-81 13132123 PAGE 00001
oonoi
oooo?
00003
0000«
oooos
00006
00007
00008
00009
00010
0001 1
0001?
00011
OOOlo
00015
00016
0001 7
0001 8
00019
00020
00021
00022
00023
00021
00025
00026
00027
00029
00029
OOOJO
00031
00032
00033
00031
00035
00036
OOOJ7
00038
00039
00040
oooai
00012
00013
OOoau
OOOU?
00006
00017
OOOflB
0004Q
00050
00051
00052
00053
C
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
SUBROUTINE TRBFLOtCTRB, CTRIB, ENDTR8, ETIMf., FERR0R, DEPTH,
i NELEM, NCWUI, NE*TRB, QHIN, TRBOPT, DEPMIN)
WHEN THERE is A TRIBUTARY TO THE SEGMENT THIS SUBROUTINE is
CALLED tACH TIME STEP TO READ THE DATA FRO* LUN 7
WHICH WAS WRITTEN BY SUBROUTINE TRBOAT,
FORMAL PARAMETERS!
CT«8 REDISTRIBUTED CONCENTRATIONS
CTRlb ORIGINAL TRIBUTARY MA$S FLUX
ENDTRB ENDING TIME OF THE CURRENT TRIBUTARY DATA (1*0)
ETiMt ELAPSED TIME OF THE SIMULATION u*
-------
(FCECS VERSION 22.
-------
(FLECS VERSION 22.
-MAR.8| 1JIJ2I36 PARE 00001
0000)
00002
OOOOJ
OOOOU
00005
00006
00007
00008
00009
00010
00011
00012
00015
00011
00015
00016
00017
oooje
00019
00030
00021
00032
00023
00030
00025
00026
00027
00028
C
C
C
C
C
C
C
C
C
C
C
SUBROUTINE TRISOL(*P1 ,01 ,D,02,R)
CALLED BY TRANSP.
USED IN ORIGINAL VERSION OF SERATRA
INCLUDE 'UMSIZ.PRM*
DIMENSION D(MXELEM), Ol(MXELCM), D2(*XELEM), R(MX£L£M)
NBMPt
Ml *N«1
FORWARD ELIMINATION
D0( Ia| ,N1 )
, 01 0«0 1 ( I ) /0( I )
. O(iti)=0(i*n-o2(n*oio
, R(Ifl)9"(I>l)*R(I)*010
...FIN
BACKWARD SUBSTITUTION
a (N) *a (N)/0(N)
00 (I«I,N1 )
, K»N-I
, P(K )«(R(K)«D2(K)*ff (K»l ) )/0(K)
...FIN
RETURN
END
(FLEC3 VERSION 22,a6)
179
-------
(FlECS VERSION 22,<<6)
13-MAH-81
I3I32J«B PAGE 00001
00001
00002
OOOOJ
ooooa
0000?
00006
00007
OOOOB
00009
00010
00011
00012
00013
oooiu
00015
00016
0001 7
00019
00019
00020
00021
00022
0002J
0002U
00025
00026
00027
00028
00029
000 JO
OOOJt
00032
00053
OOOJ«
OOOJ5
OOOJ6
00017
OOOjp
OOOi'
000«0
o o o m
0 0 0 a 2
000LEN
HTtAREA(J»l)*9LEN
NT«J
IELM(I)*J
GO TO 10
..FIN
..FIN
CONTINUE
LINFAHLY INTERPOLATE *10\» AT ELFMENT'S T"P NODE
wToP*«a » (ELTQP-f B)*(WT-«8) /(ET-EB)
ASSUME TRAPZOIDAL SH&PES TO FIND C«OSS-SFCT IUN^I. IREAS
ARATS(WTvWTuP)*{ET-ELTuP>'2,
AHBTa{l|lTOPt'*6)*(ELInP-£B)/^,
DElEWMINE IF NE« ELEMENT SURFACES HAyR «EEN FQuND TO
(A) LIE WITHIN A SINGLE DATA SET
IB) LIE IN SEQUENTIAL DATA SETS, QP
(C) BE 3EPARAUO BY QN£ OR MO&E DATA 3ETS
FINALLY, FORM THE CROSS-SECTIONAL ARE*
180
-------
(FLECS VERSION 22.«6)
IJ.MAR-61 13l32l«9 PAGE 00002
00050
00055
00056
00057
00058
00059
00060
0006]
00062
0006J
00060 C
00065 C
00066 C
00067
00068
00069
00070
00071 C
00072 C
0007J C
00070
00075
00076
00077
00078
00079 f
00080 E
INDIC»NT«NB
IFflNDIC ,EQ, 0) XSAREAU )»AR8T"AR88
IF (INUIC ,GE, t) XSAREA(I)3AR8TtARAB
IF (INOIC ,GE. ?)
XSAKao,
00(lI«N8tl«NT-| )
, xSAP«xSAR»(AREA(I)»XSAREA(l )/OELZ
A8AR(I)«AW10(I)*ALEN
BW I0( I + 1 )*WTGP
VOL«VOL»X8AREA(I)»ALEN
OVERWRITE INITIAL INFORMATION FOR NEXT ELEMENT
£LTOP«ELTOP*DELZ
ARABcAKAT
ARRH»*H8T
NB«NT
..FIN
!£TURN
NO
(FLECS VERSION 22.«6)
181
-------
(FLECS VERSION 22,
10-APR-S1
I3i56iu9 PAGE ooooi
OOOOi
00002
00003
oooou
00005
00006
00007
00008
00009
00010
00011
0001?
00011
0001«
00015
QOQ1!)
0001 7
00019
00019
00020
00021
00022
00023
000?"
0002?
00026
00027
00028
00029
00030
00031
OOOJ2
00033
00035
00036
00037
0003«
00039
000«0
o o o 'j i
000'42
000«3
0 0 0 '4 U
000«6
00017
000«8
00n<49
00050
00051
00052
0005J
SUBROUTINE UPSOAT(ECHO, HLOERR, NUMERSi S 1MLEN, OW JO , UEL )
THIS ROUTINE IS RESPONSIBLE FOR READING THE UPSTREAM INFLOW
CONDITIONS TO SEGMENT I. TH£ DATA IS READ FROM THE INPUT STREAM
(LUN 1) AND THEN WRITTEN TO "DUMMr. DTP (LUN 2) FOR SUBSEQUENT
USE DURING THE SIMULATION,
FORMAL PAHAMETEHSi
ECHO - LINE PRINTER ECHO OPTION CONTROL VARIABLE ( L » 1 )
HLOERR • HOLDING ARRAY FOR ERROR NUMBERS (BYTE)
NUMEPR . NUMBER OF INPUT ERRORS
3IMLEN • SIMULATION LENGTH (SECONDS • I***)
CALLED BV| SERATRA
CALLSl PUIERR
INCLUDE ''
BYTE
LOGICAL*! ECHO
DIMENSION CCIN(MXELE*,MAXCON), UWI D( MXELEM) , UEL(MXELEM)
REWIND 2
IF (ECHO) WRITE(6,6)
., .... UPSTREAM INFLOW CONDITIONS TO SEGMENT 1
REPEAT UNTIL (ENOTIM .£0, -9909)
RECOKD ONE .....
COL, 1-IO...ENOTIM, f ,, ENDING TIM£ FOR DATA THAT FOLLOWS,
A VALUE OF -9999 TERMINATES THE DATA,
(SECONDS)
11-15. ,,NH, ,.,.., .NUMBER OF ELEMENTS
lb-25,..UOEPTH,, .ELEVATION OF FpEE SURFACE ABOVE BED
HEADM.S) ENDTIM, NM,
,NE. -9
-------
(FLECS VERSION 22,06)
13156109 PARE 00002
00050
00055
00056
00057
00058
00059
00060
00061
00062
00063
00060
00065
00066
00067
00068
00069
00070
00071
00072
00073
00070
00075
00076
00077
0007B
00079
00080
00081
00082
00083
00090
00085
00086
00087
00088
00069
00090
00091
00092
00093
00090
00095
00096
00097
00"99
00099
ooioo
00101
00102
00103
00100
00105
00106
00107
00(08
00)09
C 3. ..CONCENTRATION OF CLAy (KG/M*«3)
£
C
c
£
c
c
c
c
c
«. ..CONCENTS* flON OF CONTAMINANT ASSOCIATED WITH
SAND
5, ..CONCENTRATION OF CONTAMINANT ASSOCIATED WITH
SILT
6, ..CONCENTRATION OF CONTAMINANT ASSOCIATED WtTH
CLAY
7, ..CONCENTRATION 0^ DISSOLVED CONTAMINANT
WHEN(NM ,GT. 0)
. DO U»1,NM)
, . READ (1,1) (CCiN(I,K), K»l, MAXCON)
, . , .FIN
...FIN
ELSE
, NMsMXELEM
. HEAO (l,|) (CCIN(l.K), K»l, MAxCoN)
. DO (I»2,NM)
, DO (Ksl, MAXCUN) CCIN(I,K)sCCIN(l,K)
, , . .FIN
, , , F I N
IF (ECHO)
. WRjTE(b,0) ENDTIM
, WRI TE (6, 3 )
. 00 (1*1, NM)
. WRITE(6,2) I, (CCIN(I,K) ,Kst, MAXCON)
. . . .p J^
...FIN
W«ITE (2) ENDTIM ,NM,UDEPTH, ( (CC IN( 1 , K ) , K« t , MAXCDN) , I a 1 ,NM )
PKETIM « ENIJTIM
, ,F JN
..FIN
IF (PREUM tuT, S1MLEN) CALL PUT£RR(20, NUM£RR, HLDERR)
C
REwiNt' 2
C RECORD 1HHEE..,.. CHANNEL CROSS-SECTION DATA
c uwio -• wfinH Or SEGMENT AT NODES
c ,...,utL •• ELEVATION OF NODE ABOVE BOTTOM
C
C***************************************************************
C CAIJTJUNJ
C IF CONCENTRATIONS «RF INPUT AT VARIOUS DEPTHS,
C THE CHANNEL CROSS-SECTIONAL DATA MUST BE INPUT
C AT THQSE SAME DEPTHS.
C
C**************************************************************
READU,7) NWJO, UWJDTH, OEu
WHEN (NMIO.EQ.O)
e««.
DO ( I * 1 1 M^ELfcM )
Ut»lD( J)»UW IOFH
! UtL(I)aE
, F*E+OEL
...FIN
..FIN
183
-------
(FLEC3 VERSION 32,06)
t3i56i«9 PAGE ooooe
0005«
00055
00056
00057
00058
00059
00060
00061
00062
00061
0006«
00065
00066
00067
00068
00069
00070
00071
00072
00071
00071
00075
00076
00077
0007B
00079
oooeo
00091
OOQ82
00083
0008U
00085
00096
00087
00088
00089
00090
00091
00092
00n?j
0009U
00095
00096
00097
00098
00099
001 00
ooiot
00102
0010}
OOtOM
00105
00106
00107
00108
00109
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C*
C
C
C
C
C
c »
, . J,., CONCENTRATION OF CLAy (KG/M«*J)
, . «... CONCENTRATION OF CONTAMINANT ASSOCIATED WITH
. . SAND
. . 5, ..CONCENTRATION OF CONTAMINANT ASSOCIATED WITH
. . SILT t
. . 6, ..CONCENTRATION OF CONTAMINANT ASSOCIATED WITH
. . CLAY
. . 7, ..CONCENTRATION OF DISSOLVED CONTAMINANT
. .
, , ^H£N(NM ,GT. o)
... 00 (1»J,NM)
.... READ (1,1) (CCIN(I,K), K»l> MAXCON)
FIN
. . ...FIN
. . ELSE
. . . NMsMXELEM
. . , HEAD (1,1) (CCIN(1,K), K.I, MAxCQN)
... DO U'2,NM)
.... 00 (Ksl, MAXCON) CCIN(T.»K)«CCIN(1,K)
FIN
• • t f t " * ™
. . IF (ECHO)
. . . w«lTE(b,'O ENDTIM
. . . WRIT£(6,J>
... 00 (Isl,NM)
.... WRITE{6,2) I, (CC1N(I,K),K*1, MAXCON)
FIN
fin
t •
, . W*ITE<2) ENDTIM, NM, UDEPTH, ( (CCIN(1 ,K) ,K«l ,M*XCON) , Ia| ,NM)
, , PHETIM a ENIJTJM
. ...FIN
...FIN
IF (PRE'IM tLT, SIMLEN) CALL PUTERR(2«> NUMgRR, HLDERR)
RtWINU 2
RECORD IHHEE.....CHANNEL CROSS-SECTION DATA
,,..,uKio — wfnTH OF SEGMF.NT AT NQOES
. ....utL •• ELEVATION OF NODE ABOVE BOTTOM
'
•»•**•«»»»*»<•****•«»*»•»»»*•»»«»«*«**•«»»*»•***»•***«»*******
CA1IT1UNI
IF CONCENTRATIONS »PE INPUT AT VARIOUS DEATHS,
THE CM»NNtL CROSS-SECTIONAL DATA MUST HE. INPUT
AT THQSE SAME DtPTHS.
**»»**»«»***<»**»**»»»t*******»»*»»»*«»**»**.»»***»*»***»*i>**
READ(1,7) NWID, UWfOTH, DEL
WHEN (NWJO.EQ.O)
. ES<>,
, 00 ( I - I , M
-------
(FUCS VfHSION
PAGE oooos
00110
00111
001 12
00111
09114
001 15
001 16
OOU7
00116
001 19
00120
00121
00122
0012J
00120
00125
00126
00127
00128
00129
00130
001 Jl
00112
oom
001 3«
001i<5
00116
00137
OOUfl
00119
00140
ooiai
ELSE
, READ(l,l) (UWID( I) ,Ial «NWID)
, HEAO(1,1)
, if (wwio.tr.
, . DO (I » NW1D » 1,MX£IEM)
, . . UWIO(I) * UWIDTH
. . . UEL(I) « UEUI-I) + OEL
..... FIN
. ...FIN
...FIN
IF(ECHO)
, WRITE(6,8) (I,UWID(I),I»l,MxELEM)
, WRIT£(6,9) (I
RETURN
FORMAT (8FIO.O)
*SSOC,',2X),2X,ICONTAMINANT'X
2 9x, 'SUSPENDED SAND t, | x ,' SUSPENDED SIM ', 1 X, » SUSPENDED CL*Y»,
i JX.'WITH SANo',6x,»t*rrH SILT', 6x, 'WITH cLAr»,«x, »oissoi. VED '
« 'CONC,')
FO»H»? { 1H0.8X, IIO» ' ,,,0»TA SET E^olNO TlM£')
FORM*! (110, 15, F10, 0)
FORMAT (1HO,««X, "UPSTREAM INFLOW WATER CONO I T IUHS ' )
FORMAT (I5,2FJO,0)
FORMAT (IHO,56X, IX. SECTION M IDT H3 I /It ( J7X , 5 ( I 3, 1 P E I 2 . 5) /) )
FORMAT (IHO,58X, IX- SECTION ELEVAT IONS ' /O (27X , 5( I 3 , 1 PC 12 ,5) /))
END
(FUECS VERSION 22,
185
-------
(FLECS VERSION 22,06)
29.MAK-8I 07l22»56 PAGE OQOOl
oooot
00002
OOOOJ
00004
00005
00006
00007
OOOOP
oooo1?
00010
0001 1
00012
OOOtJ
000 11
00015
00016
0001 7
00016
00019
00020
00021
00022
00023
00020
00025
00026
00027
00028
00029
00030
00031
OOOS2
0003J
0003"
00035
000j6
00037
00038
00039
00000
00001
0 0 0
-------
(FLEC3 VERSION 2?.«6)
29.MAH-81 0/123156 PAGE 00003
00051
00055
00056
00057
0005«
00059
OOOjO
00061
00062
00065
00060
00065
00066
00067
00068
00069
00070
00071
00072
00075
0007U
00075
00076
00077
If (ECHO)
. WRirt(6,3>
, WR1TM6.5)
, 00 (J«l,NELEMf1)
,(C(J,K),Ks!,MAXCON)
...FIN
C
RETURN
t
I FQRMAI(8F10.0)
2 FORMATUS.FIO.OJ
3 FORMAT OH«iS2X» ' INITI*L «4TER CONDITIONS')
5 FORMAT(1HO,'ELEMENT',IX,
1 3(3x,'CONC, OF),ux),3(lX,'CONC. AS30C,',2x),2X,'CONTAMJN&NT ' -
2 9X,'SUSPENDED SANOMX, '3USPENPEO SILT ', IX, 'SUSPENDED CL*Y',
3 3x,'«irH s*Noi,6x, »WITH siLr»,6x, IWITH CLAV^X,'DISSOLVED '
-------