29 





DUM1=(0MEGAA2)ADD/G 



SNELL 



28 



30 





DUH2=DUM1+1. 0/(1. 0+0. 6522ADUM1+0.4622A 



SNELL 



29 



31 





A(DUM1AA2)+0.0864A(DUM1A*4>+0.0675A(DUM1AA5)) 



SNELL 



30 



32 





DUM3=TASQRT ( GADD/DUM2 ) 



SNELL 



31 



33 





RKNUM=PI2/DUM3 



SNELL 



32 



34 





RKA(I,J)=RKNUM 



SNELL 



33 



35 





GRDK(I f J)=RKNUM 



SNELL 



34 



36 



C 





SNELL 



35 



37 



CAAAAAAAA COMPUTE IHE HEIGHTS AND ANGLES USING SNELL'S LAW 



SNELL 



36 



38 



CAAAAAAAA COMPUTE SIGMA AND CCG 



SNELL 



37 



39 



C 





SNELL 



38 



40 





DUMbRKNUMADD 



SNELL 



39 



41 





DUH2=2.0ADUM1 



SNELL 



40 



42 





DUM3=TANH(DUM1> 



SNELL 



41 



43 





0UH4=SINH(DUM2) 



SNELL 



42 



44 





3 INE=S IN ( IHETAO-CNIRANG ) ADUM3 



SNELL 



43 



45 





ZANG=PI-ASIN(SINE)+CNTRANG 



SNELL 



44 



46 





Z(I,J)=ZAMG 



SNELL 



45 



47 





SI(I,J)=SIN(ZANG! 



SNELL 



46 



48 





CO(I f J)=COS(ZANG) 



SNELL 



47 



49 





DUH5=SQRT ( COS ( IHETAO-CNTRANG) /COS ( ZANG-CNTRANG ) ) 



SNELL 



48 



50 





H ( I , J ) =HHADUH5ASQRT ( . 5/ ( . 5A ( 1 . 0+DUM2/DUM4 ) 



SNELL 



49 



51 





AADUM3)) 



SNELL 



50 



52 





CCG(I,J)=O.5AU.0+<DUM2/DUM4)}A<0MEGAA2)/ 



SNELL 



51 



53 





A(RKNUMAA2) 



SNELL 



52 



54 





1 CONTINUE 



SNELL 



53 



55 



C 





SNELL 



54 



56 



CAAAAAAAA CHECK FOR WAVE BREAKING 



SNELL 



55 



57 



C 





SNELL 



56 



58 





DO 600 MWEI.M 



SNELL 



57 



59 





DO 604 J=1,N 



SNELL 



58 



60 





DUM5=H(I,J) 



SNELL 



59 



61 





HBRK=0.78AD(I,J) 



SNELL 



60 



62 





IE(DUM5.GE.HBRK) DUM5=HBRK 



SNELL 



61 



63 



C 





SNELL 



62 



64 



CAAAAAAAACONVERT HEIGHTS TO AMPLITUDE FUNCTION (HAG/2ASIGMA) 



SNELL 



63 



65 



C 





SNELL 



64 



66 





H(I,J)=DUM5/HFACI 



SNELL 



65 



67 





604 CONTINUE 



SNELL 



66 



68 





600 CONTINUE 



SNELL 



57 



69 





RETURN 



SNELL 



68 



70 





END 



SNELL 



69 



1 



CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA ANGLE 



T 



■} 



C 





ANGLE 



3 



3 





SUBROUI INE ANGLE ( ISTARI , IEND , ITMAX , WI , ALPHA , SCNV ) 



ANGLE 



4 



4 



C 





ANGLE 



5 



5 



CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA ANGLE 



6 



6 



C 





ANGLE 



7 



7 



C 



THIS SUBROUTINE ITERATES TO SOLVE FOR IHE WAVE ANGLES 



ANGLE 



8 



8 



C 



ALONG A GIVEN ROW 



ANGLE 



9 



9 



C 





ANGLE 



10 



10 



CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA ANGLE 



11 



11 



C 





PARAM 



i 



12 





PARAMEIER(IQ=95,JQ=95> 



PARAM 



3 



13 



C 





PARAM 



4 



14 





COMHON/ANGLES/Z( IQ, JQ) ,SI( IQ, JQ) .C0( IQ, JQ) 



ANGLE 



13 



15 





C0MM0N/DEPTHS/D(IQ,JQ) 



ANGLE 



14 



16 





COMMON/WAVNUM/RKA ( IQ , JQ ) , GRDK ( IQ , JQ ) 



ANGLE 



15 



17 





COMMON/CONST/ G,PI.PI2,RAD.HC0NVR,SC0NVR,DX,DY r DX2,DY2,T,0MEG, 



ANGLE 



16 



F15 



