|
<br />U4=0.
<br />DO 320 Jls2 N4
<br />Ul=Ul+LOG(H(2 J1)~
<br />U2sU2+(LOG(H{~jJI ))**2.
<br />U3-U3+LOG(B{4 1)
<br />U4=U4+LOG(H(2:Jl) *LOG(B(4,Jl))
<br />320 CONTINUE
<br />N3sN4-1
<br />IF(N4 .GT. 2) GO TO 322
<br />CONLY 2 X-S POINTS GIVEN. SET M = 1
<br />Ml=l 0
<br />K1=B(4AN4)/H(2,N4)
<br />GO TO -,30
<br />322 M1=U4-U1*U3/N3
<br />Ml=M1/(U2-(U1**2.)/N3)
<br />IF(M1 .LE. 0.0) Ml = 0.01
<br />Kl=U3/N3-Ml*Ul/N3
<br />Kl=EXP(Kl )
<br />
<br />C
<br />C CHECK GOODNESS OF FIT OF WEIGHTED X-SECTION
<br />IU8=0
<br />IU7s0
<br />325 IF(IUB.GT.l) WRITE(IOUT 9230) K1 M1
<br />9230 FORMAT(/lX,'DEPTH' ,4X,'WEIGHTED tw' ,4X.'FIT TW' ,4X,'K =',F6.2,
<br />13X 'M =' Fo.2)
<br />DO'32B Ji=2 N4
<br />T4=Kl*H(2,Ji)**Ml
<br />IF(IU8.GT.l) WRITE(IOUT,9231) H(1,J1),B(4,Jl),T4
<br />9231 FORMAT(/IX F6.2 6X F8.2 3X F8.2)
<br />IF (ABS(lB(4,J1}-TJ.))/B(4,jl).GT.0.5) IU7=N4-1
<br />328 CONTINUE
<br />IF (IUB.GT.1) GO TO 330
<br />IF (IU7.LT.1) GO TO 330
<br />WRITE(IOUT,9232) D(I)
<br />9232 FORMAT(/lX~'*WARNING* WEIGHTED X-SECTION FIT IS OFF AT MILE'iF6.2,
<br />l/lX,'ROUTI~G ERRORS MAY OCCUR...USER SHOULD PLOT THE X-S,' / X
<br />l'DRAW A SMOOTH CURVE THRU POINTS, AND RERUN WITH SMOOTHED POINtS')
<br />IU8=2
<br />GO TO 325
<br />C
<br />C RECOMPUTE DEPTij AT DAM USING WIGHTED X-S PROPERTIES
<br />330 S2=(E(l,l)-E(I,l))/((D(I)-D(1))*52BO.)
<br />K=Kl
<br />MsM1
<br />O=p(I 1)
<br />CALL tOMPHV(I,H3)
<br />CALL COMPQV(K,M N S2,H3 Q3)
<br />CALL STAGE(K,M,iI,S2,Q,H1,Q3,Hl)
<br />HO=Hl
<br />C Fl IS THETA
<br />Fl=0.95
<br />19=0
<br />340 IB=I9
<br />O=p(I 1)
<br />CALL STAGE(K,M,N S2 Q,H3 Q3 HI)
<br />CALL RPARMS(K,M,iI,s~,Q,Fi,Hi,H3,F7,V7,X1)
<br />X5=Xl
<br />350 !F(I9 .GT. 0) X5=D(I)
<br />Q=P (I 1)
<br />CALL ROUTE(F7,V7,Xl,X5,Q,Q7)
<br />C
<br />IF(I9 .GT. 0) GO TO 360
<br />C CHECK THETA VALUE
<br />O=Q7
<br />CALL STAGE(K1M,NAS2,Q,H3,Q3.Hl)
<br />F2=(HO+H1)/(L *Hu)
<br />IF((ABS(F1-F2j/F1) .LE. 0.02) 19=1
<br />IF(I9 .NE. 18 GO TO 350
<br />C NEW THETA MUST BE COIIPUTED
<br />Fl=(Fl+F2)/2.
<br />GO TO 340
<br />
<br />C
<br />C
<br />C SAVE PEAK FLOW AT X-S (1)
<br />360 P (I ,1 ) = Q7
<br />
<br />11-7
<br />
|