26 

 27 

 28 

 29 

 30 

 31 

 32 

 33 

 34 

 35 

 36 

 37 

 38 

 39 

 40 

 41 

 42 

 43 

 44 

 45 

 46 

 47 

 48 

 49 

 50 

 51 

 52 

 53 

 54 

 55 

 56 

 57 

 58 

 59 

 60 

 61 

 62 

 63 

 64 

 65 

 66 



Ckkkkkkkk COMPUTE THE DIFFRACIWE TERMS 



C 



Ckkkkkkkk USE 4 POINT BACKYARDS DIFFERENCE FOR THE X CURVATURE OF H' 



Ckkkkkkkk USE 3 POINT BACKWARDS DIFFERENCE FOR THE X GRADIENT OF 'H.CCG' 



Ckkkkkkkk USE CENIRAL DIFFERENCES FOR IHE Y CURVATURE OF 'H' 



CAAAAAAAA USE CENIRAL DIFFERENCES FOR THE Y GRADIENT OF 'H.CCG' 



C 



DO 2 I=ISTART,IEND 



QLX2=DX2AXMUC(I) 



DLXSQ=(DLX2A0.5)AA2 



DO 2 J=2.NM1 



DLY2=DY2AYHUC(J) 



DLYSQ=<DLY2A0.5)AA2 



CCGU=CCG(I,J) 



HIJ=H(I,J) 



HJP1=H(I,J+1) 



HJM1=H(I,J-1) 



HIPl=H(I+l r J> 



HIP2=H(I+2.J) 



DUM6=(-3.0AHIJ+4.0AHIP1-HIP2)/DLX2 



DUM4=(2.0AHIJ-5.0*HIP1+4.0AHIP2-H( 1+3, J) >/DLXSQ 



BUM4=BUM4-0 . 5ADUM6A ( XMUC ( 1+1 ) -XHUC ( 1-1 ) ) / ( ( XMUC ( I) AA2) ADX ) 



DU«2=(-3.0ACCGIJ+4.0ACCG<I+1,J)-CCG(I+2,J))/DLX2 



DUM7=(HJP1-HJM1)/DLY2 



DIM5=(HJP1-2.0AHIJ+HJM1)/DLYSQ 



DUH5=DUH5-0.5ADUh7A(YHUC(J+l)-YMUC(M))/((Y«UC(J)*A2)*nY) 



DUM3= ( CCG < I . J+l ) -CCG ( I . J- 1 ) ) /DLY2 



DIM ( I, ])=( DUM4+DUM5+ ( DUM2ADUM6+DUM3ADUM7 ) /CCG IJ ) /H I J 

 2 CONTINUE 

 C 



Ckkkkkkkk CHECK FOR POINT TO POINT OSCILLATIONS IN yAVE PHASE GRADIENT 

 C 



DO 410 I=ISIART.IEND 



DUH1(I,1)=DUH1(I,2) 



DUM1(I.N)=DUMHI,NM1) 



SMU=4.0 



SBEIA=0.25 



DO 409 3=1,N 



X(J)=DUM1(I.J) 



409 CONTINUE 



DO 408 J=3.NM2 

 BUH2=ABS(X(J+1)-X(J)) 

 DU«3=ABS(X(J)-X(J-D) 

 DUM4=ABS(X(J+1)-X(]-1))A0.5 

 IFC{DUM2+DUH3).LI.(SMUADUM4))GQ TO 408 

 DUH5=XU+1>-2.0AX(J)+X(J-1! 

 DUM6=X(J+2)-2.0AX(J+l)+X(J) 

 DUM7=X(J)-2.0AX(J-l)+X(J-2) 



IF((DUM5ADUM6).GE.0.0.AND.(DUH5ADUM7>.GE.0.O>GO 10 408 

 DUM1(I.J)=X(J)+SBETAA(X(J+1)-2.0AX(J>+X<J-1>) 

 408 CONTINUE 



410 CONTINUE 

 C 



DO 17 I=ISTART.IEND 



DO 17 J=2,NH1 

 C 



Ckkkkkkkk DO NOT INCLUDE DIFFRACIWE EFFECTS INSIDE OR IMMEDIATELY 

 Ckkkkkkkk ADJACENT TO IHE BREAKER ZONE 

 C 



GRADK 



24 



GRADK 



ir 



GRADK 



26 



GRADK 



27 



GRADK 



28 



GRADK 



29 



GRADK 



30 



GRADK 



31 



GRADK 



32 



GRADK 



33 



GRADK 



34 



GRADK 



35 



GRADK 



36 



GRADK 



37 



GRADK 



38 



GRADK 



39 



GRADK 



40 



GRADK 



41 



GRADK 



42 



GRADK 



43 



GRADK 



44 



GRADK 



45 



GRADK 



46 



GRADK 



47 



GRADK 



48 



GRADK 



49 



GRADK 



50 



GRADK 



51 



GRADK 



52 



GRADK 



53 



GPADK 



54 



GRADK 



55 



GRADK 



56 



GRADK 



57 



GRADK 



58 



GRADK 



59 



GRADK 



60 



GRADK 



61 



GRADK 



62 



GRADK 



63 



GRADK 



64 



GRADK 



65 



GRADK 



66 



GRADK 



67 



GRADK 



63 



GRADK 



69 



GRADK 



70 



GRADK 



71 



GRADK 



72 



GRADK 



73 



GRADK 



"4 



GRADK 



75 



GRADK 



76 



GRADK 



77 



GRADK 



73 



GRADK 



79 



GRADK 



30 



GRADK 



81 



GRADK 



82 



F10 



