149 

 150 

 151 

 152 

 153 

 154 

 155 

 156 

 157 

 158 

 159 

 160 

 161 

 162 

 163 

 164 

 165 

 166 

 167 

 168 

 169 

 170 

 171 

 172 

 173 

 174 

 175 

 176 

 177 

 178 

 179 



DM7=DM8(J) 

 HH2=H(IM,J)AHFACI 

 C 



Ckkkkkkkk UPDATE THE HAVE HEIGHT 

 C 



EPSZ=0.5 



DM7=EPSZ*DM7+ ( 1 . 0-EPSZ ) AHH2 

 SUH=SUH+ABS(DM7-HH2) 

 C 



Ckkkkkkkk CONVERT HEIGHT BACK TO AHPLITUDE FUNCTION 

 C 



H(IM,J)=DH7/HEACT 

 400 CONTINUE 

 C 



Ckkkkkkkk SET LATERAL BOUNDARY CONDITIONS FOR WAVE HEIGHT 

 C 



H(IN,1)=H(IM,2) 

 H(IM,N)=H(IM.NM1) 

 C 



Ckkkkkkkk IEST FOR CONVERGENCE 

 C 



IE ( SUM . LT . ( FLOAT ( NH2 ) AHCNV ) ) GO 10 500 

 450 CONTINUE 

 C 



IFdREF.EQ.DGO 10 500 

 WRITE<6,475) IM 

 475 FORMAK IX, 'RELAXATION FOR HAVE HEIGHTS FAILED ON ROW =',I4) 

 C 

 500 CONTINUE 

 RETURN 

 END 

 Ckkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

 C 



SUBROUTINE SNELL (THETAO.HH) 

 C 

 Ckkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



c 



C THIS SUBROUTINE COMPUTES THE SNELLS LAW SOLUTION OVER IHE ENURE 



C GRID 



C 



Ckkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 



C 



PARAMETER ( IQ=95 , JQ=95 ) 

 C 



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



COMM0N/UVH/H( IQ, JQ) ,CCG( IQ , JQ! 



COMMON/DEPIHS/D<IQ,JQ) 



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



COMMON/CONST/ G,PI,PI2,RAD,HC0NVR,SC0NVR.DX,DT,DX2,DY2.I,0MEG. 

 AM,N,M1,H2,NM1,NM2,IUET.IDRY,IVETP1,CNTRANG,HFACT,DLEVEL 



C0MM0N/PRINIC/IP1,IP2,INC,JP1,JP2,DWJLI,HMULI,ZMULT,RKMULI 



C0MM0N/C0NS2/ IGR ID , IREF . ITAMX , ITHMX , ID IFF 

 C 



Ckkkkkkkk COMPUTE THE HAVE NUMBERS USING THE PADE APPROXIMATION 

 Ckkkkkkkk (SEE THE COASTAL ENGINEERING NOTEBOOK) 

 C 



DO 1 I=IUET,M 



DO 1 J=l r N 



DD=D(I,J) 



HEIGHT 



148 



HEIGHT 



149 



HEIGHT 



150 



HEIGHI 



151 



HEIGHI 



152 



HEIGHT 



153 



HEIGHI 



154 



HEIGHI 



155 



HEIGHI 



156 



HEIGHT 



157 



HEIGHT 



158 



HEIGHT 



159 



HEIGHT 



160 



HEIGHT 



161 



HEIGHT 



162 



HEIGHT 



163 



HEIGHI 



164 



HEIGHT 



165 



HEIGHT 



166 



HEIGHI 



167 



HEIGHT 



168 



HEIGHT 



169 



HEIGHI 



170 



HEIGHT 



171 



HEIGHT 



172 



HEIGHT 



173 



HEIGHT 



174 



HEIGHT 



175 



HEIGHT 



176 



HEIGHT 



177 



HEIGHT 



178 



SNELL 



2 



SNELL 



3 



SNELL 



4 



SNELL 



5 



SNELL 



6 



SNELL 



n 



SNELL 



8 



SNELL 



9 



SNELL 



10 



SNELL 



11 



PARAM 



2 



PARAM 



3 



PARAM 



4 



SNELL 



13 



SNELL 



14 



SNELL 



15 



SNELL 



16 



SNELL 



17 



SNELL 



18 



SNELL 



19 



SNELL 



20 



SNELL 



21 



SNELL 



22 



SNELL 



23 



SNELL 



24 



SNELL 



25 



SNELL 



26 



SNELL 



27 



F14 



