|
<br />so Q4=(1.49/N)*(SZ**.S)*(A1**1.667)/(B1**.667)
<br />tF(Q .LT. Q4) GO TO 100
<br />IF(QO.GT.Q4) GO TO 100
<br />C COMPUTE TIME OF FLOOD
<br />T4=T1-((Q-Q4)/(Q-QO))*TO/60.
<br />C COMPUTE TI~ OF D~FLOOD
<br />TS=Tl+(l(Q-Q4)/<Q-QO))*((B71Z0.*VO/(Q-QO))-TO*60.))/3600.
<br />100 RETURN
<br />END
<br />
<br />C
<br />c*************************************SUBROUTINE RPARMS
<br />SUBROUTINE RPARMS(K,M,N1SZ,Q,F1,Hl,H3,F7,V7,X1)
<br />C THIS ROUTINE COMPUTES THE RuUTINE PARAMETERS XC, FC, AND V*
<br />C
<br />C INPUT PARMS ARE
<br />C K&M=CHNL FIT PARMeAVO=RES.VOL'1CO=DAM CREST ELEV, E(I,l)=DAM BASE EL.
<br />C H1=DEPTH AT DAM~ tlj=HV, F1=THE"A, SZ=SLOPE, N=MANNING N
<br />C OUTPUT PARMS AR~
<br />C F7=FC, V7=V*, X1=XC
<br />C
<br />REAL K KO Kl,M MO,Ml,N
<br />cO~lliON1AR1y/ FtlS,6),X(6),D(8),E(8,8),T(B,8),C(3,8),H(Z,B),B(4,B),
<br />1P(B,S)
<br />CO~lliON/SCALAR/ CO,BO,VO,AO,WO,TO,QO,NO,Nl,SO
<br />C D8 IS THE DAM HEIGHT
<br />DB=CO-E(l 1)
<br />X1=((VO*41s60.*(M+l.))/(K*D8**(M+1.)))*(6./(1.+4.*(.S**(M+l.))))
<br />V064=6.*VO*43S60.
<br />IF(DB.GT.H3) Xl=V064/((K*H3**M)*(3.*D8-S.*(M*H3)/(M+l.)))
<br />IF(D8.LT.H3.0R.(DB/2.).GT.H3) GO TO SO
<br />Xl=V064/((K*H3**M)*(H3/(M+1)+D8-H3)+4*(K/(M+1))*(D8/Z)**(M+l))
<br />SO IF(Fl.LE.O.Ol) F1=0.9S
<br />Dl=F1*H1/(M+1.)
<br />IF((Fl*H1) .GT. H3) Dl = Fl*Hl-H3+(H3/(M+l.))
<br />VI = (1.49/N)*(S2**.5)*(D1**0.667)
<br />F7=Vl/ ((32 .2*Dl )**. S)
<br />A2=K*((Fl*Hl)**M)*D1
<br />IF(H1 .GT. H3) AZ=K*(H3**M)*D1
<br />V7=(VO*43S60.)/(AZ*X1)
<br />Xl =X1/S2BO .
<br />RETURN
<br />END
<br />
<br />C
<br />C*************************************** SUBROUTINE ROUTE
<br />SUBROUTINE ROUTE(F7 V7 Xl,XS,Q,Q7)
<br />C THIS SUBROUTINE COMPUTtS tHE PEAK FLOW AT A GIVEN X-S
<br />C
<br />C
<br />C
<br />C
<br />C
<br />
<br />INPUT PARMS ARE:
<br />F7 = FC, V7 = V*, XS = DIST TO X-S, XI = XC,
<br />RETURNED VALUE IS Q7
<br />REAL K,KOAKl1M,MO,M1,N .
<br />INTEGER Vo V~ FB F9 FO VS
<br />CO}~IDN/ARAY/ F(1~,6},Xt6),D(B),E(B,8),T(8,8),Cl3,8),H(2,8),B(4,8),
<br />1P(B,S)
<br />CO~illON/SCALAR/ CO,BO,VO,AO,WO,TO,QO,NO,N1,SO
<br />C SET X* VALUE
<br />X2=XS/X1
<br />C
<br />C SET
<br />
<br />V* POINTER
<br />V8=IFlX(V7)
<br />V9=IFlX(V7+1j
<br />IF{V7 .LE. 1
<br />IF V7.GT. 5
<br />IF V7 .GT. S
<br />C
<br />C SET FC POINTERS
<br />F8=IFlX (4. *F 7)
<br />F9=IFIX(l4.*F7)+1.)
<br />IF(F7 .LT. 0.2S) F8=F9
<br />C ***NOTE: WHEN FC=l.S CURVES AVAILABLE,
<br />IF(F7 .GT. 0.75l FB=3
<br />IF(F7 .GT. 0.75 F9=3
<br />DO=(F7-FLOAT(F8 /4.)/0.ZS
<br />IF(F8 .EQ. F9) DO=O.
<br />
<br />VB=l
<br />v8=S
<br />V9=5
<br />
<br />CHANGE FOLLOWING TWO LINES
<br />
<br />Il-12
<br />
|