GRADK 



83 



GRADK 



84 



GRADK 



85 



G8ADK 



86 



GRADK 



87 



GRADK 



88 



GRADK 



89 



GRADK 



90 



GRADK 



91 



GRADK 



92 



GRADK 



93 



GRADK 



94 



GRADK 



95 



GRADK 



96 



GRADK 



97 



GRADK 



98 



GRADK 



99 



GRADK 



100 



GRADK 



101 



GRADK 



102 



GRABK 



103 



GRADK 



104 



GRADK 



105 



GRADK 



106 



GRADK 



107 



GRAKK 



108 



GRADK 



109 



GRADK 



110 



GRADK 



111 



84 IBIESIMBRKM(J-1)+IBRKMU) + IBRKM(J+1)+IBRK<J-1) + IBIK(J)+ 



35 IBRKU+1) 



86 U(IBTEST.GT.O) GO 10 17 



87 C 



88 CJbbWtAM* LIMIT IHE CHANGE IN THE GRADIENT OF THE mi PHASE FUNCTION 

 39 UddMtkkk FRO* THE DISPERSION REUIION UAVE NUMBER 10 50 PERCENT 



90 C 



91 RKOLD=RKA(I.J) 



92 RKADD=0.25ARKOLD 



93 RKAIG=RK0LD**2+DUM1(I,J) 



94 TF(RKARG.LE.0.0)RKAI6=(<0NE6A*2/G)AA2) 



95 RKNEy=SQRI(RKARG) 



96 RKDIEE=RKNEy-RKOLD 



97 IF(ABS(RKDIFF).LI.(0.0O25xIK0LD)iGO TO 17 



98 IF(ABS(RKDIFF).LE.RKADD!GO TO 341 



99 IT(RKDIFF.LI.O.O)RKNEy=RKOLD-RKADD 



100 IF(RKDIEE.GE.0.O)RKNEU=RKOLD+RKADD 



101 341 G1DK(I,J)=RKNS 



102 17 CONTINUE 



103 C 



104 CfcdWdbk** SET LATERAL BOUNDARY CONDITIONS FOR IHE 6RADIENI OF IHE iAVE 



105 CJbddobbW PHASE FUNCTION 



106 C 



107 DO 38 1= 1ST AIT, KM) 



108 GRDK(I,1)=GRDK(I,2) 



109 GRDK(I,N)=6RDK(I,HM1) 



110 38 CONIINUE 



111 RETURN 



112 END 



1 CkkkkttkkkkkAkk^kkf.k&kf.&t.%tekttk&kkkk&.kkkkktxk HE IGHT - 



2 C HEIGHT 3 



3 SUBROUTINE HEIGHT! ISIARI. TEND, UMAX, yi, ALPHA, HCNV) HEIGHT 4 



4 C HEIGHT 5 



5 C*JdbbbbbbbbddbdbbbbWbbbU^ HEIGHT 6 

 h C HEIGHT 7 



7 C THIS SUBROUTINE ITERAIES 10 SOLVE FOR IHE yAVE HEIGHIS HEIGHI 3 



8 C ALONG A GIVEN RG1I HEIGHT 9 



9 C HEIGHI 10 

 10 CJUydbdddbbbbba^JbbbWbdUbbdd^^ HEIGHT 11 

 U C PARAH 2 



12 PARAMETER; 10=95,19=95) PARAM 3 



13 C PARAM 4 



14 COMMON/ ANSLES/Z(IQ,;iQ),SI(IQ,JQ).CO(IQ,JQ> HEIGHI 13 



15 CQMMON/yVH/H(IQ,JQ),CCG(IQ,JQ) HEIGHI 14 



16 COMMON/DEPIHS/D(IQ,JQ) HEIGHI 15 



17 CQMM0!WAVNUM/Ri<A<IQ,JQ),GRBK(I9,JQ) HEIGHI 16 



18 COMMON/CONST/ G,PI,PI2,RAD,HC0NVR,SC0NVR,DX,DY,DX2,DY2,I,0HEG, HEIGHI 17 



19 *M,N,Ml,M2,NMl,NM2,iyET,IDRY,iyETPl,CNIRANG,HFACI,DLEVa HEIGHI 18 



20 C0MMON/CQNS2/IGRID,IREF,ITAMX,ITHMX,IDIFE HEIGHI 19 



21 C0MMON/MUSC/XMUC(IQ).XMUS!IQ),YMUC(JQ),YMUS(JQ) HEIGHI 20 



22 C0MMON/CO0I/XX(IQ),YY(JQ) HEIGHT 21 



23 COMMON /IRANF/ DECAY, STABL.IBIK(J9),IBRKM(JQ) HEIGHI 22 



24 DIMENSION DM8 (JO), SLOPE (JQ) HEIGHT 23 



25 C HEIGHT 24 



26 C AAAAAAA A SOLVE DIFFERENCED FORM OF CONSERVATION OF yAVE ACTION EQUATION HEIGHT 25 



27 C HEIGHI 26 



28 DO 500 I=IEND,ISTARI,-1 HEIGHI 27 



29 IM=M HEIGHI 28 



30 DO 450 11=1. UMAX HEIGHT 29 



7' ' 



