|
<br />KO=U3/N3-MO*Ul/N3
<br />KO=EXP(KO)
<br />C
<br />C CHECK GOODNESS OF FIT BETWEEN GIVEN AND FITTED X-SECTION
<br />IUB=O
<br />IU7=0
<br />245 IF(IUB.GT.l} WRITE(IOUTI9220) KO,MO
<br />9220 FORMAT(/IXJ DEPTH' ,4X,'bIVEN TW' ,4X,'FIT TW' ,4X,'K =' ,F6.2,
<br />13X 'M =' Fo.2)
<br />DO'24B Jl=2 N2
<br />T4=KO*H(I,Jl)**MO
<br />IF(IU8.GT.1) WRITE(IOUT,9221) H(I,Jl),T(I,Jl),T4
<br />9221 !~~I~{ift~~i5~~*})'Ttt~~1)~Gi~0.5) IU7=N2-1
<br />248 CON'l;INUE
<br />IF (IUB.GT.l) GO TO 250
<br />IF(IU7.LT.1) GO TO 250
<br />WRITE(IOUT,9222) D(I)
<br />9222 FORMAT(/IX, '*** WARNING *** X-SECTION FIT IS OFF AT MILE'lF6.2,
<br />1/1X,'DEPTH 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 245
<br />250 IF(NO .LT. 2 ) GO TO 175
<br />IF(I .EQ. 1) GO TO 400
<br />C
<br />C COMPUT DISTANCE WEIGHTED X-SECTION
<br />C
<br />DO 280 Jl=2,N4
<br />C
<br />C INTERP H~2,J1~ INTO ELEV TABLE TO FIND ASSOC TOPW @ X-S(I)
<br />IF(H 2,Jl .GT.(E(I Nl)-E(I 1))) JO=Nl
<br />IF(H 2,Jl .GT. (EtI,N1)-E(I,l))) GO TO 270
<br />DO 260 J2=2,Nl
<br />JO=J2
<br />IF((H(2,J1) - (E(I,J2)-E(I,I))) .LE. 0.0) GO TO 270
<br />260 CONTINUE
<br />C BRACKETING ELEVS FOUND...COMPUTE DB/DH
<br />270 J3=JO-l
<br />DO=(T(I,JO)-T(I,J3))/(E(I,JO)-E(I,J3))
<br />C
<br />C DET. ASSOC TOPW @ X-S I
<br />B(1,Jl)=T(I,J3)+(E(I,1)+H(2,J1)-E(I,J3))*DO
<br />C CHECK IF HV IS EXCEEDED
<br />IF(H(2,Jl) .GE. C(2,I)) B(1,J1) = T(I,N2)
<br />C B(3,J1) IS THE WEIGHTED TOPWIDTH
<br />B(3tJl)=B(31J1)+((B(2,Jl)+B(I,Jl))/2.)*(D(I)-D(I-1))
<br />C SAVE IN ERPOLAT~D TOPW AS B(2,Jl)
<br />B( 2tJl)=BO,Jl)
<br />280 CON lNUE
<br />C
<br />C DIVIDE WEIGHTED TOPW BY DISTANCE
<br />DO 290 Jl=2,N4
<br />Bl4tJ1)=B(3,J1)/(D(I)-D(1))
<br />290 CON lNUE
<br />C
<br />C COMPUTE DISTANCE WEIGHTED MANNING /I
<br />N=O.
<br />DO 300 J1=2 I
<br />N=N+((C(1,Jl)+C(l,Jl-1))/2)*lD(Jl)-D(Jl-l))
<br />300 CONTINUE
<br />N = N/(D(I)-D(I))
<br />C
<br />C FIT K & M PARMS TO DIST WEIGHTED X-S
<br />C
<br />C
<br />
<br />IF(N1 .GT. 2) GO TO 310
<br />CONLY 2 X-S POINTS GIVEN. SET M = 1
<br />HI=1 0
<br />K1=B(4ANl)/H(2,Nl)
<br />GO TO J30
<br />C INITIALIZE LEAST SQUARES ALGORITHM
<br />310 Ul=O
<br />U2=0.
<br />U3=0.
<br />
<br />1I-6
<br />
|