|
<br />READ (IN , 70Z) QO
<br />C
<br />150 WRITE(IOUT,913)
<br />913 FORMAT(/lX~'INPUT NUMBER OF X-SECTIONS (STARTING AT DAM) USED TO' ,
<br />1/lX1 'DESCKIBE ROUTING REACH. (MIN=l, MAX=B)')
<br />REAIl(IN 703) NO
<br />IF(NO.Lt.l ) GO TO 150
<br />C
<br />160 WRITE(IOUT 914)
<br />914 FORMAT(/lit 'INPUT NUMBER OF PAIRS OF ELEV-TOPWIDTHS USED TO' ,
<br />l/lXL 'DESCR BE EACH X-SECTION (MIN=2, MAX=B)')
<br />IF(vERS.EQ.SH) WRITE(IOUT 9140)
<br />9140 FORMAT(lX!'MODEL ASSUMES DEPTH #1 = O. AND TOPW #1 = 0.')
<br />READ(IN 7u3) Nl
<br />IF(N1.Lt.Z ) GO TO 160
<br />C
<br />C BEGIN READING IN X-SECTION INFORMATION
<br />C
<br />D(1)=O.
<br />DO 400 I=I,NO
<br />T(I,l)=O.
<br />IF(I.GT 1) WRITE(IOUT 915) I
<br />915 FORMAT(/lX,'INPUT DIStANCE DOWNSTREAM (MILES) TO X-SECTION #',IZ)
<br />IF(I.GT.l) READ(INI70Z) D(I)
<br />IF(VERS.NE.SH) GO Tu 165
<br />E(1J.l)=O.
<br />CO=112
<br />80=0.
<br />ELEV=STG
<br />FTHSL(l)=FEET(l)
<br />FTMSL(Z)=FEET(Z)
<br />GO TO 167
<br />C
<br />165 WRITE(IOUT,916) I
<br />916 FORMAT(/lX1'INPUT LOWEST CHANNEL ELEV.(FT.MSL) AT X-SECTION #',12)
<br />IF(I .EQ. ) WRITE(IOUT,917)
<br />917 FORMAT(IX! '(ELEVATION OF BASE OF DAM (FT. MSL))')
<br />READ(IN 7uZ) E(I,l)
<br />167 IF(I .Gt. 1 ) GO TO 170
<br />C
<br />C IF SURFACE AREA NOT GIVEN, ESTIMATE AO USING VO
<br />IF(AO .LT. 0.01) AO=2.~VO/(CO-E(I,I))
<br />C IF VOLUME NOT GIVEN, ESTIMATE VO USING AO
<br />IF(VO .LT. 0.01) VO=(AO/Z.)*(CO-E(l,l))
<br />IF(BO .LE. (E(l 1)+5.)) GO TO 170
<br />C PERCHED BREACH...RECOMPUTE VO BASED ON AVAILABLE VOLUME
<br />VO=H2/2.*AO*(1.+((BO-E(1,1))/(CO-E(I,1)))**2.)
<br />C
<br />C ONCE FIRST X-S HAS BEEN READ IN & FITTED & SLOPE FOUND COMPUTE PEAK
<br />170 IF(I .NE. Z ) GO TO 1BO
<br />175 H=HO
<br />K=KO
<br />N=C(1 ~1)
<br />:a:g~3 'B .
<br />CALL DAMBRK(K,M,N~SZ,H2tH3,~4,Q3,Q,H1)
<br />C SAVE FLOW~DEPTH AND T~EAK A DAM
<br />P~l 1 =Q
<br />P l:Z =Hl
<br />P 1t3 =TO/60.
<br />Tl= 0/60.
<br />CALL TIMFLD(K~M,N,SZ,H3,H4,Q,Q3,T1,T4,T5)
<br />C SAVE TIMES OF FLuODING AND DEFLOODING AT DAM
<br />P(I,4l=T4 .
<br />P(1,5 =T5
<br />IF(NO .LT. 2) GO TO 400
<br />
<br />C
<br />
<br />C
<br />
<br />C
<br />IBO
<br />C
<br />
<br />918
<br />
<br />DO 190 J=Z ,N1
<br />WRITE(IOUT,91B) ELEVAFTMSL(l),FTMSL(Z),J,I
<br />FORMAT(/lX,'INPUT ' ,~A4,' #',12,' AND ASSOC.
<br />l' r IZ)
<br />READ(IN,70Z) ElI,J)
<br />
<br />TOPWIDTH FOR X-S '
<br />
<br />II-4
<br />
|