18 





AM,N,M1,M2,NM1,NM2,IUET,IDRY.IWETP1,CNTRANG,HFACT,DLEVEL 



ANGLE 



17 



19 





C0MM0N/C0NS2/ IGR ID , IREE , ITAHX , IIHHX , ID IFF 



ANGLE 



18 



20 





COMMON/MUSC/XMUC ( IQ ) , XMUS ( IQ ) , YMUC ( JQ ) , THUS ( JQ ) 



ANGLE 



13 



21 





DIMENSION DUM1(IQ,JQ),X(JQ> 



ANGLE 



20 



in 



C 





ANGLE 



21 



23 



CAAAAAAAA SOLVE THE DIFFERENCED FORM OF THE IRROTAIIONALIIY OF THE 



ANGLE 



22 



24 



CAAAAAAAA GRADIENT OF THE WAVE PHASE FUNCTION EQUATION 



ANGLE 



23 



25 



C 





ANGLE 



24 



26 





DO 1 I=IEND.ISTARI.-1 



ANGLE 



25 



27 





IM=I-1 



ANGLE 



26 



28 





DO 4 11=1, UMAX 



ANGLE 



27 



29 





DO 2 J=2.NM1 



ANGLE 



28 



30 





DUM6=ALPHAA ( GRDK ( I , J+l ) AS I ( I , J+ 1 ) ) 



ANGLE 



29 



31 





A +(1.0-2.0AALPHA)A(GRDK(I,J)ASI(I,J)) 



ANGLE 



30 



32 





A +ALPHAA(GRDK(I,J-1)ASI(I,J-D) 



ANGLE 



31 



33 





DUM2=WIA(GRDK(I-1.J+1)AC0(I-1,J+1)-GRDK(I-1,J-1)A 



ANGLE 



32 



34 





AC0(I-1.M))/(DY2AYMUC(J)) 



ANGLE 



33 



35 





A +(1.0-«T)A(GRDK(I,J+1)AC0(I,J+1)-GRDK(I,J-1)A 



ANGLE 



34 



36 





ACQ(I,J-1))/(DY2AYMUCU)) 



ANGLE 



35 



37 





DUH3=DUH6- ( DXAXMUS ( 1-1 ) ) ADUM2 



ANGLE 



36 



38 





DUM4=DUM3/GRDK(I-1,J> 



ANGLE 



37 



39 





DUM1(I,J)=DUM4 



ANGLE 



38 



40 





2 CONTINUE 



ANGLE 



39 



41 





SUh=0.0 



ANGLE 



40 



42 





DO 22 J=2.NM1 



ANGLE 



41 



43 





DUH4=DUH1U,J) 



ANGLE 



42 



44 





DUh5=Z(I-l,J) 



ANGLE 



43 



45 





DUM3=SIN(DUM5) 



ANGLE 



44 



46 



C 





ANGLE 



45 



47 



CAAAAAAAA UPDATE IHE SINES 



ANGLE 



46 



48 



C 





ANGLE 



47 



49 





EPSZ=0.5 



ANGLE 



48 



50 





DUM4=EPSZADUM4+< 1 . O-EPSZ ) ADUM3 



ANGLE 



49 



51 



c 





ANGLE 



50 



52 



CAAAAAAAA LIMIT OBLIQUE WAVE ANGLES 



ANGLE 



51 



53 



C 





ANGLE 



52 



54 





IE ( DUM4 . GE . . 997 ) DUM4=0 . 997 



ANGLE 



53 



55 





IE ( DUM4 . LE . -0 . 997 ) DUH4=-0 . 997 



ANGLE 



54 



56 





DUM6=PI-ASIN(DUM4> 



ANGLE 



55 



57 





SUM=SUM+ABS(DUH6-Z(IM,J>> 



ANGLE 



56 



58 





Z(IM.J)=DUM6 



ANGLE 



57 



59 





5I(Ih,J)=DUM4 



ANGLE 



58 



60 





CO<IH.J>=C0S(DUM6) 



ANGLE 



59 



61 





22 CONTINUE 



ANGLE 



60 



62 



c 





ANGLE 



61 



63 



CAAAAAAAA SET LATERAL BOUNDARY CONDITIONS FOR WAVE ANGLES, SINES. AND 



ANGLE 



62 



64 



CAAAAAAAA COSINES 



ANGLE 



63 



65 



C 





ANGLE 



64 



66 





Z(IM,1)=Z(IM,2) 



ANGLE 



65 



67 





CO(IM,1)=CO(I«,2) 



ANGLE 



66 



68 





SHIM,1)=SI(IH,2) 



ANGLE 



67 



69 





Z(Df,N)=Z(IM,N«l) 



ANGLE 



68 



70 





C0(IM,N)=CQ(IM,NM1) 



ANGLE 



69 



71 





SKIM, N)=SI(IM, NMD 



ANGLE 



70 



72 



C 





ANGLE 



71 



73 



CAAAAAAAA CHECK FOR ANGLE CONVERGENCE 



ANGLE 



72 



74 



C 





ANGLE 



73 



75 





IF(SUM.LI.(NM2ASCNV))G0 TO 1 



ANGLE 



74 



76 





4 CONTINUE 



ANGLE 



75 



F16 



