|
<br />IP(B,5)
<br />CO~lliON/SCALAR/ C01BOIVO,AOLWO,TO,QO,NO,N1,SO
<br />C CHECK FOR INSTANTAN~OU:; FaILuRE
<br />IF(TO.LT.(HZ/1000)) GO TO lZ5
<br />C1=( 23 .4*AO )tWO
<br />C H8 IS HEAD OVER WEIR
<br />HB=lCl/llTO/60.)+C1/(HZ**.5)))**2.
<br />Q=00+3.1*WO*H8**1.5
<br />C COMPUTE gv AT Dn
<br />IF(N .LT. Z S2=SO/5ZBO.
<br />IF(NO .LT Z GO TO 50
<br />S2=(E(I,l)-E(2,1))/((D(Z)-D(1))*5Z80.)
<br />50 CALL COMPQV(K1M,N,SZ,H3,Q3)
<br />C CmWUTE DEPTH AT JJAM
<br />CALL STAGE(KIM,N1S2,Q,H3,Q3,Hl)
<br />C CHECK FOR SUBMERI>ENC~
<br />C NO CHECK IF FLOW AT DAM IS SUPERCRITICAL
<br />Dl=Hl/(M+l.)
<br />IF(Hl.GT .H3) D1=H1-H3+H3/ (M+l.)
<br />V9=(1.49/N)*SQRT(S2)*D1**(2./3.)
<br />FRD=V9/(32.Z*D1)**0.5
<br />IF(FRD.GE.l.) GO TO 1~0
<br />IF(((Hl+E(1,1)-BO)/H8) .LT. 0.67) GO TO 150
<br />C6=1.
<br />DO 100 II=l 10
<br />C7=1-27.B*(((Hl+E(1,1)-BO)/HB)-0.67)**3.
<br />IF(C7 .LT. .8) C7 = .B
<br />C7=(C6+Cn/Z.
<br />C6=C7
<br />Q6=Q
<br />O-C 7*Q
<br />!F(ABS(Q6-Q) .LT. Q6/20.) GO TO 150
<br />C RECOMPUTE DEPTH AND HEAD WITH NEW fLOW
<br />CALL STAGE(K M N SZ,Q H3 Q3 HI)
<br />H8=H8+((Q6-Q}*(60.*TO}/(435~0.*Z.*AO))
<br />100 CONTINUE
<br />GOTO 150
<br />
<br />C
<br />
<br />C
<br />C INSTANTANEOUS FAILURE CALCULATION
<br />125 M1 =M+l.
<br />Bl=WO/T(1,N1)
<br />C6=(1.0+0.148*(Bl**Z.)*(M1**Z.)-0.083*(Bl**3.)*(Ml**Z.5))**1.5
<br />C7=1.-.5467*Bl*Ml**.5+.2989*(B1**2.)*Ml**.Z5
<br />C7=C7-.1634*(Bl**3.)*M1**.125
<br />C7=C7+.0893*(B1**4.)*Ml**.0625-.04B6*(Bl**5.)*M1**.03125
<br />C7=C7**1.5
<br />Q=QO+3.1*WO*(C6*C7)*HZ**1.5
<br />SZ=SO/5280
<br />IF(NO.GT.l) SZ=(E(1,1)-E(2,1))/llDlZ)-D(1))*5280)
<br />CALL COMPQV(K,M N S2 H3 Q3)
<br />CALL STAGE(K,M,N,S2,Q,H1,Q3,Hl)
<br />150 RETURN
<br />END
<br />
<br />C
<br />C***************************************** TIME OF FLOOD ROUTINE
<br />SUBROUTINE TI~~LD(K M,N S2 H3,H4,Q,Q3 T1 T4,T5)
<br />C THIS SUBROUTINE COMPUTES THE tIME OF FLOOD AND DEFLOOD
<br />C INPUT PARMS ARE
<br />C K,M = CHNL FIT PARMS, S2=SLOPE, Q=FLOW, H4=FLOOD DEPTH, Q3=QV, H3=HV
<br />C
<br />REAL K KO Kl,M MO,Ml,N
<br />COMMON'ARAY/ Ft15,6),X(6),D(8),E(8,8),T(8,8),C(3,8),H(2,8),B(4,B),
<br />1pl8,5)
<br />CO~lliON/SCALAR/ CO,BO,VO,AO,WO,TO,QO,NO,Nl,SO
<br />04=-999.
<br />14=-999.
<br />T5=-999.
<br />IF(H4 .LT. 0.01) GO TO 100
<br />IF(H4 .GT. H3) GO TO 25
<br />Bl=K*H4**M
<br />Al=lK/(M+l.))*H4**(K+1.)
<br />GO TO 50
<br />25 B1=K*ij3**M
<br />Al=(K/lM+l.))*H3**(M+1.)+(H4-H3)*Bl II-l1
<br />
|