|
<br />.
<br />
<br />READ(IN1702) T(I,J)
<br />190 CONTINU~
<br />C
<br />WRITE(IOUT,919) I
<br />919 FOWT(/IX"IN~UT ~NNING N AT X-S #',12,' (DEFAULT - O.OSO)')
<br />READ IN 704) Cll,I)
<br />IF{C l,i) .LT. 0.0001) C(I,I) = O.OSO
<br />IF(VERS.NE.SH) GO TO 19S
<br />g~H ~:H:
<br />GO to 197
<br />C
<br />19S WRITE(IOUT 920) I
<br />920 FORMAT(/IX~ 'INPUT VALLEY WALL HEIGHT ABOVE CHANNEL BOTTOM AT "
<br />I'X-S #',121 IX, 'IF N/A INPUT 0.')
<br />READ(IN 70", C\2,I)
<br />IF(C(2,!) .LT. 0.01) C(2,I) = (E(I,Nl)-E(I,l)) + H2
<br />WRITE(IOUT 921) I
<br />FO~T(/l*, 'INPUT FLOOD DEPTH AT X-S #',12,
<br />l/lX IF N7A ENTER 0.')
<br />READ(IN,702) C(3,I)
<br />C
<br />197 IF(NO .LT. 2) WRITE(IOUT 922)
<br />922 FORMAT(/lX1'INPUT CHANNtL SLOPE (FT/MI) BELOW DAM')
<br />IF(NO .LT. L) READ(IN 702) SO
<br />IF(SO.EQ.O) SO=H2**2*(f(1,Nl)/(24.7S*VO))
<br />C
<br />C FIT X-SECTION W/ K & M PARMS
<br />IF(Nl .GT. 2) GO TO 200
<br />CONLY 2 X-S POINTS GIVEN. SET M = 1
<br />MO=1 0
<br />KO=T(I Nl)/(ElI Nl)-E(I 1))
<br />IF(I .l.T. 2) Ht2,Nl)=EtI,Nl)-E(I,1)
<br />IF(I .LT. 2) B(2,Nl)-T(I,Nl)
<br />N4=N1
<br />GO TO 2S0
<br />C INITIALIZE LEAST SQUARES ALGORITHM
<br />200 U1=0.
<br />U2-0.
<br />U3-0.
<br />U4=0.
<br />N2=Nl
<br />DO 230 J1=2,N1
<br />, H(l Jl)=E(I J1)-E(I 1)
<br />C SAVE H'S OF FIRST X-S FOR WEIGHTED X-S ROUTINE
<br />IF(I .GT. 1) GO TO 210
<br />H(2 J1)=H(l,J1)
<br />IF(H(1,J1) .GEI C(2,1)) H(2,Jl)=C(2,1)
<br />B(2,J1J-T(I JIJ
<br />C SEE IF HV IS EXCEEDED
<br />210 IF(H(l Jl-1) .GE. C(2 I)) N2=J1-1
<br />IF(H(l,jl) .GE. C(2,I)} H(11Jl)=C(2,I)
<br />IF(H(lJJ1-1) .GEI C\2,I)) Gu TO 240
<br />Ul=U1+LOG(H(1 J1J)
<br />U2=U2+(LOG(H(I,JI)))**2.
<br />IF(H(l Jl) .LT. C(2 I)) GO TO 220
<br />C HV MET OR EXCEEDED COM~UTE TW ASSOC. W/ HV
<br />DO=(T{I,J1)-T(i J1-1))/(E(I,J1)-E(I Jl-1))
<br />T(IiJ1)=T(I,Jl-I)+(H(l,J1)+E(I,1)-E(I,Jl-l))*DO
<br />IF( .EO. IJ B(21l>=T\I,JI)
<br />220 U3-U3+J:OG(tO JI )
<br />U4=U4+LOG(Hl1,jl} *LOG(T(I,J1))
<br />230 CONTINUE
<br />240 N3=N2-1
<br />IF(I .EQ. 1) N4=N2
<br />IF(N2 .CT. 2) GO TO 242
<br />CONLY 2 X-S POINTS BELOW HV. SET M = 1
<br />MO=l 0
<br />KO=T(I1N2)/H(1,N2)
<br />GO TO ",SO
<br />242 MO=U4-Ul*U3/N3
<br />MO=MO/(U2-(U1**2.)/N3)
<br />IF(MO .LE. 0.0) MO = 0.01
<br />
<br />C
<br />
<br />C
<br />921
<br />
<br />II-S
<br />
|