|
<br />IF(P(I-l,l).LT.P(I,I)) P(I,I)=.99*P(I-l,l)
<br />C COMPUTE TIME TO PEAK
<br />X5-D(I)/Z
<br />CA~L ROUTE{F7~V7,Xl,X5,Q,Q7)
<br />Q=l.3+.1*K)*QI
<br />CALL StAGE(~,M,N,SZ,Q,H3,Q3,Hl)
<br />D9-HlIlK+l.)
<br />IF(Hl .GT. H3) D9=HI-H3+H3/(M+l.)
<br />V6-{1.49/N)*(SZ**0.5)*D9**0 667
<br />C9=0.68Z*V6*(1.667-( 667*(M]lM+1.))))
<br />P(I 3)=(TO/60.)+D{I)iC9
<br />IF(P(I,3).LT.P(I-l,3)) P(I,3)-1.01*P(I-l,3)
<br />C
<br />C COMPUTE AND SAVE PEAK DEPTH AT X-S (I)
<br />K=KO
<br />M-M~
<br />SZ= E(I-l,1)-E(I,1))/((D(I)-D(I-l))*5Z80.)
<br />N-C 1 I)
<br />H3=c(i1I)
<br />Q=pO )
<br />CALL ~OMPQV(K,M,N S2 H3 Q3)
<br />CALL STAGE(K,M,N,SZ,Q,Hd,Q3,Hl)
<br />PO 2)=Hl
<br />C DET TI~ FLOODING AND DEFLOODING
<br />Tl=P 1,3)
<br />H4=C 3 I)
<br />CALL TfMFLD(K,M,N,SZ,H3,H4,Q,Q3,Tl,T4,T5)
<br />P(I,4)=T4
<br />PO, 5) =T5
<br />
<br />C
<br />C
<br />400 CONTINUE
<br />IF(NO .GT. 1) GO TO 455
<br />705 FORMAT(A4)
<br />C SET UP FORECAST POINT PROFILE DISTANCES
<br />WRITE(IOUT,9Z60)
<br />9260 FORM4T(/ZXA'ENTER NUMBER OF FORECAST POINTS BELOW DAM (MAX = 7)')
<br />READlIN, 70-,) NO
<br />NO=NO+l
<br />IF(NO .LT. 2) GO TO 455
<br />HO = P(1,2)
<br />C Fl IS THETA
<br />Fl=O .95
<br />19=0
<br />420 18=I9
<br />Q=P(1,1)
<br />CALL STAGE(K,M,N,S2 Q,H3 Q3 HI)
<br />CALL RPARMS(K,M,N,si,Q,FI,Hl,H3,F7,V7,Xl)
<br />X5=X1
<br />CALL ROUTE(F7,V7,Xl,X5,Q,Q7)
<br />C CHECK THETA VALUE
<br />Q=Q7
<br />CALL STAGE(K1M,NASZ,Q,H3,Q3,H1)
<br />FZ=(HO+Hl)/(L.*Hu)
<br />IF((ABS(FI-FZ)/F1) .LE. 0.02) 19=1
<br />IF(I9 .GT. 0) GO TO 430
<br />C NEW THETA MUST BE COMPUTED
<br />Fl=(Fl+FZ)/Z.
<br />GO TO 420
<br />
<br />C
<br />430 DO 450 I = Z,NO
<br />IM=I-l
<br />WRITE(IOUT,926) IM
<br />926 FORMAT(/lX1'INPUT DIST DOWNSTREAM (MILES) TO FORECAST PT. NO.' ,12)
<br />READ lIN,7uZ) D(I)
<br />X5 = D(IJ
<br />Q=P(1,1)
<br />CALL ROUTE(F7,V7,Xl,X5,Q,Q7)
<br />C
<br />C SAVE PEAK FLOW AT X-S (I)
<br />pO,1)= Q7
<br />C COMPUTE TIm; TO PEAK
<br />X5=D(I)/2
<br />CALL ROUTE(F7AV7,Xl,X5,Q,Q7)
<br />Q=( .3+.1*K)*QI II-8
<br />
<br />.
<br />
|