2 



90 





X4=0.5A(XX(IM)+XX(I))-XX(IM> 



HEIGHT 



89 



2 



91 





Y3=0.5A(YY(J)+YY(J-1))-YY(J) 



HEIGHT 



90 



<j 



92 





Y4=YY(J-1)-YY(J) 



HEIGHT 



91 



i 



93 





ANG3=AIAN2(Y3.X3) 



HEIGHI 



92 



2 



94 





ANG4=AIAN2(Y4,X4> 



HEIGHI 



93 



? 



95 



C 





HEIGHT 



94 



2 



96 





IF(IHETA.GT.ANG3) THEN 



HEIGHI 



95 



3 



97 





KEY=IBRKM(J) 



HEIGHI 



96 



3 



98 





IKEY=I 



HEIGHT 



97 



3 



99 





JKEY=J 



HEIGHI 



98 



3 



100 





ELSE IF(IHETA.GT.ANG4) THEN 



HEIGHT 



99 



3 



101 





KEY=IBRKH(J-1) 



HEIGHI 



100 



3 



102 





IKEY=I 



HEIGHI 



101 



3 



103 





JKEY=J-1 



HEIGHI 



102 



3 



104 





ELSE 



HEIGHT 



103 



3 



105 





KEY=IBRKCJ-1) 



HEIGHT 



104 



3 



106 





IKEY=IM 



HEIGHI 



105 



3 



107 





JKEY=J-1 



HEIGHT 



106 



3 



108 





END IF 



HEIGHI 



107 



i 



109 





END IF 



HEIGHI 



108 



2 



110 



C 





HEIGHI 



109 



i 



111 



CAAAAAAAA COMPUTE INCIPIENT BREAKING WAVE HEIGHT USING MEIHOD FROM 



HEIGHT 



110 



n 



112 



CAAAAAAAA THE SHORE PROTECTION MANUAL (WEGGEL'S EMPIRICAL METHOD) 



HEIGHI 



111 



2 



113 



C 





HEIGHI 



112 



1 



114 





JJ=(J+(JKEY-J-l)/2) 



HEIGHI 



113 



1 



115 





BISI^SQRK (FLOAK IKEY-IM)ADXAXMUS( IM) )AA2+ 



HEIGHT 



114 



1 



116 





(FLQAT(JKEY-J)ADYAYHUS(JJ)>AA2) 



HEIGHI 



115 



1 



117 





SLOPE! J)=(D( IKEY, JKEY)-D( IM, J) )/DISI 



HEIGHI 



116 



1 



118 





SLOPE(J)=AMAX1<0.0 ( SLOFE(J)> 



HEIGHT 



117 



1 



119 





HSIG=1.0ADM7 



HEIGHI 



118 



1 



120 





A1=43.75A(1.0-EXP(-19.0ASL0PE(J))) 



HEIGHI 



119 



1 



121 





B1=1.56/(1.0+EXP(-19.5ASL0PE(J))) 



HEIGHT 



120 



1 



122 





DEPBRK=HSIG/(B1-<A1AHSIG/(GAIAI>)> 



HEIGHT 



121 



1 



123 





IF(DEPBRK.GI.D(IM,J)) KEY=1 



HEIGHI 



122 



1 



124 



C 





HEIGHI 



123 



1 



125 



CAAAAAAAA CHECK WHETHER 10 ACIUALLY TRANSFORM WAVE 



HEIGHI 



124 



1 



126 



C 





HEIGHI 



125 



1 



127 





IF(KEY.EQ.l) THEN 



HEIGHT 



126 



2 



128 





IBRK(J)=1 



HEIGHT 



127 



■5 



129 





DMl=V.5AiD(IKEY,JKEY>+D(IH,J)) 



HEIGHT 



128 



2 



130 





DM2=CCG ( IKEY , JKEY ) AGRDK ( IKEY , JKEY ) A ( H ( IKEY , JKEY ) AA2- 



HEIGHT 



129 



2 



131 





i STABLAD ( IKEY . JKEY ) /HFACI ) AA2 ) -CCG ( IM , J ) AGRDK ( IM , J ) A 



HEIGHI 



130 



2 



132 





(STABLAD(IM r J)/HFACI)AA2 



HEIGHT 



131 



-i 



133 





DM3=DXAXMUS(IM)ADECAY/(2.0ADM1ADH4) 



HEIGHI 



132 



•) 



134 





DM5=< DM5+DM3ADM2 > / ( 1 . 0-CCG ( IM , J ) AGRDK (IM, J) ADM3 ) 



HEIGHT 



133 



•y 



135 





DM7=DM5A(HFACTAHFACI) 



HEIGHI 



134 



2 



136 





IE(DM7.LI. ( (STABLADt IM, J) )AA2) ) DM7=(SIABLAD( IM. J))AA2 



HEIGHT 



135 



T 



137 





DM7=SQRT(DM7) 



HEIGHT 



136 



2 



138 





END IF 



HEIGHT 



137 



1 



139 





END IF 



HEIGHT 



138 



1 



140 



C 





HEIGHI 



139 





141 





DM8(J)=DM7 



HEIGHI 



140 





142 



c 





HEIGHT 



141 





143 





200 CONTINUE 



HEIGHI 



142 





144 



c 





HEIGHT 



143 





145 



CAAAAAAAA CHECK FOR WAVE HEIGHT CONVERGENCE 



HEIGHI 



144 





146 



C 





HEIGHI 



145 





147 





SUM=0.0 



HEIGHT 



146 





148 





DO 400 J=2,NM1 



HEIGHI 



147 



F13 



