RUNT VE«SION FEB 74 B 17«12 04/23/76 



SUBROUTINE COEFF 

 C THIS SUBROUTINE CALCULATES THE PAPTS 0^ I(T,J) AND K(I,J) 

 C WHICH ARE INDEPtNOENT OF FkEQUENCY MiJMeER K, 

 2 CnnMON RI12(25.25),RKr6(25,25), P0T(»';,'5)f HOW< 25, 25) . FE ( ?5t 6) , 



lFIt25.6)t RI(25.25), RJ(25»25). RK(?5,4), '>L(25,<.), 

 2R«U(3»3»13)» RLA»1(3.3,lO), FB(3,10), nELFBOtlO), Hwn( 25»6. 10 ) . 

 3 DELW(25.6.10) . xaL(25,10) 

 2 C0MM0N/QNE/X(25), Y(25),X8(25),YB(?5),ANG('5),0EL(25) .VV(25) 



1»FEIN(25). FIIN(25). RNQRM<25,3), JCCS) 

 2 COMMON /Two/ N,NNW, NWAVEL* ISYM, TSKIP. NC. P IE. GAM?i A. M, TK, TP 



2 C0MM0N/ONE2/CC3(25),SS3(25) 



2 N2 • N/2 



6 00 i I = l.M 



10 IF(I .GT. N) GO TO 7 



13 IFdSYM .EO. 1 .AND. I .GT. N2) 60 TO 7 



26 Xll • X(I) - xe(i) 

 34 Y11-Y(I)-YB(1) 



41 X21 • Xll ♦ XR(1» 



45 Y21-Y(I)+Y)(1) 



53 PP1-AL06(X11**2+Y11**2) 



66 P0l»AL0G(XH**2 + Y21**2) 



101 TPl = ATAN2(yH,XH) 



105 T01»ATAN2(Y21tXll) 



111 DO 1 J » 1»N 



112 X12'X(I)-XB(J+1) 

 120 Y12«Y(I)-YB(J+1) 

 124 Y22»Y(I)+Y8(J+1) 



131 PP2»AL0G(X12**2tY12**2) 



144 PQ2«AL0G(Xi2**2+Y22**2) 



157 TP2=ATAN2(Y12»X12) 



163 TQ2«ATAN2t Y22»X12) 



C CORRECTION FOR DISCONTINUITY IN ATAN2 ^T PTF 

 167 IF(X11 .GT. 0. .OR. X12 .GT. 0.) GD TO ^ 



201 IF(TP2 .GT. 



214 1F(TP2 .LT. 



227 o C3 * CC3(J) 

 232 S3'SS3(J) 

 235 Al-PIE 



237 1F(I-J)2»3.2 

 241 2 A1»TP1-TP2 

 243 3 A2«T02-TC1 

 245 A5«C3*(-X8( J + 1)+XB( J)-X12*0.5*PP2 + X11*0.5*I>'>1 



1+Y12*TP2-Y11*TP1)+S3*(YB( J)-Y8(J+1)-X12*TP2 



1-Y12*0.5*PP2+X11*TP1+Y11*0.5*PP1) 

 307 A6«C3*(-X8( J*i)+X«( J)-X12*0.5*P02 + Xlt*0.5«'»Ol 



1+Y22*T02-Y21*T01)-S3*(-YB( J)+Y8(J+1)-X1?*T02 



1-Y22*0.5*P02+X11*T01+Y21*C,5*P01) 

 351 4 >11«X12 



353 Y11«Y12 



354 Y21-Y22 



356 PP1-PP2 



357 P01-PQ2 



361 TPI-TP2 



362 T01-Ta2 

 364 RI12(I.J) » Al - A2 



Tah,le D-2. Continued 



0. 



.AND. 



TPl 



.LT. 0.1 



1 Toi . TPl * TP 



0. 



.AND. 



TPl 



.GT. 0.1 



1 TPl = TPl - TP 



120 



