50800 WRITE(6. 108) 



50900 108 F0RMAT(/,2X, 'THE ON-OFFSHORE TRANSPORTS, OY , FOLLOW) 



51000 DO 17 J=1,JMAX 



51100 17 WRITE(6.101) (OY ( I , J ) , I = 1 . IMAX ) 



51200 WRITE(6,109) 



51300 109 F0RMAT(/.2X. 'THE NEW CONTOUR VALUES, Y, FOLLOW) 



51400 DO 18 ij=1 , JMAX 



51500 18 WRITE(6,100) ( Y (I , J ) , I = 1 , IMAX ) 



51600 100 F0RMAT(2X, 13(F9.3)) 



51700 101 F0RMAT(2X, 13(F9.4)) 



51800 1 CONTINUE 



51900 RETURN 



52000 GO TO 446 



52100 445 STOP 



52200 446 CONTINUE 



52300 END 



52500 SUBROUTINE OTRAN 



52600 C*THIS SUBROUTINE CALCS THE BREAKER HEIGHT FOR EACH 



52700 C*OF THE I GRID LINES. METHOD- -FINDS Y-LOCATIONS BEFORE AND AFTER 



52800 C*BREAKING HAS OCCURRED BY 'REFRAC, THEN USES SHOALING TO GET THE 



52900 C*HBO.SNELL'S LAW IS USED FOR REFRACTION OVER THE SHORT DIST TO BREAKING 



53000 C* OX(I,d) IS THE TRANS BETWEEN( I - 1 , J ) AND (I,d) AT THE BLOCKCENT 



53100 COMMON/ A/ C( 60, 20) . RK( 60, 20) , Y( 60, 20) . DEEP(60, 20) , ALPHAS ( 60, 20) 



53200 C0MM0N/AA/YZER0(6O) 



53300 COMMON/B/ THETA ( 60, 20) , OXTOT ( 60) , OLDANG( 60, 20) , DY(60,20) 



53400 COMMON/C/ H( 60. 20) , CG( 60 , 20) , H0LD(6O, 20) . HB( 60, 20) , YB(60) 



535O0 COMMON/N USED/JUSE , T , CO , CGEN , CGGEN, ANGGEN , DX , BERM , THETA0( 10) ,MMAX 



53600 COMMON/D/SIGMA,G,ELO, JMAX, I MAX, PI .TWOPI ,PI02,HGEN, IJET( 10),S JETTY 



53700 C0MM0N/G/IBREAK(6O) ,HN0NBR(2O) 



53800 COMMON/E/RHO , RHOS , POROS , CONST . TKSI 



53900 COMMON/P/HB0(6O) ,DEEPB(60) 



54000 CAPPA=0.78 



54100 DO 1 1=2, IMAX 



54200 DO 2 dd=1 ,dMAX 



54300 d=JMAX-Jd+1 



54400 HDUM=(H( I ,d)+H(I-1 , J))*0.5 



54500 HBDUM=(HB(I ,d)+HB(I-1 ,d) )*0.5 



54600 C*CAN ONLY USE COND ON ONE SIDE OF STRUCT. CAN'T AVG HERE! 



54700 IF(SdETTY .EO.O.O) GO TO 4 



54800 DO 4 M=1 , MMAX 



54900 IF( I .NE . IdET(M)+1 ) GO TO 4 



55000 IF(THETAO(M) .GE.O.O) ISIDE= IdET(M) 



55100 IF(THETAO(M) .LT.O.O) ISIDE = IdET(M) + 1 



55200 C***B.C. AT STRUCT TIP ASSUMES OX COMP AS IF NO STRUCT IS PRESENT. 



55300 YSEA=0.5*(Y(ISIDE,J)+Y(ISIDE,d+1)) 



55400 IF(YSEA.GT.SdETTY) GO TO 3 



55500 HDUM=H(ISIDE,d) 



55600 HBDUM=HB(ISIDE .d) 



55700 GO TO 3 



55800 4 CONTINUE 



55900 3 IF(HDUM.'?T.HBDUM) GO TO 2 



56000 DEEPB(I)=((05*(H(I,d+1)+H(I-1,d+1)))*((0.5*(DEEP(I,d+1) 



56100 * +DEEP( I- 1 ,d+1 ) ) )**0. 25)/CAPPA)**0.B 



56200 HBOd )=CAPPA*DEEPB(I ) 



56300 C*HBO(I) AND DEEPB(I) WILL BE COMPUTED ACCORDING TO THE WAVE DIR. 



56400 C** AT THE STRUCTURE TIP,THETAO. 



56500 IF(SdETTY. EO.O.O) GO TO 1 



56600 DO 6 M=1 ,MMAX 



56700 IF(I .NE. IdET(M)+1 ) GO TO 6 



56800 C**THE TRANSPORTING WAVES WILL BE COMPUTED USING THE WAVE TO PROP SIDE. 



56900 IF(THETAO(M) .GE .0.0) GO TO 11 



57000 DEEPBd ) = (H( IdET(M)+1 ,d+1 )*DEEP( I dET ( M )+ 1 , d+ 1 ) **0 . 25/CAPPA ) * *0 . 8 



57100 I BREAK (I ) =IBREAK( IdET(M)+ 1 ) 



57200 GO TO 12 



57300 11 DEEPB( I )=(H(IdET(M),d+1 ) *DEEP(I dET( M) , d+ 1 ) **0 . 25/CAPPA ) ♦♦O . 8 



57400 IBREAKd )=IBREAK(IdET(M)) 



57500 12 HBOd )=DEEPB( I )*CAPPA 



57600 GO TO 1 



57700 6 CONTINUE 



57800 GO TO 1 



57900 2 CONTINUE 



58000 1 CONTINUE 



73 



