RUNT VERSION FEB 74 8 17:12 04/23/76 



SUflRQUTINE CQMP(K) 



C THIS SUBROUTINE COMPUTES THE CneFCTCI 5»<TS O'^pENOENT ON K 



C AND CALLS Oi^ LNEQF TO SOLVE THE SIMUL TAMFOUS EQUATIONS 



C FOR THE VELOCITY POTENTIALS FE(I»JJ) AMI FKT.JJJtFns FE? 



6 COMMON Rn2(25.25),RK5t)(25,25)f PnT(?^,'f>, 4nw( 25 . 25 ) . FF ( ?5, "^ ) , 

 1FI(25»6). RI(2!j»25). RJ(25t25), "KCJ.i), "JL(25»4), 

 2RMU(3»3.10). RLAM(3»3flOI. PB(3»10). nFLFRO.lO). HWB ( 2 5, 6» 10 ) . 

 3 OtLW(25,6»10) t XQL(25,10) 



6 COMf10N/ONE/X(25)»Y(25),XB(25),YB(?5),AMG(?5),DFL(2'. ).VV(25) 



l,FtINt2i). FIIN(2S). PNORH(25,3). JCC!) 



6 CaMM0N/0NE2/CC3(2'j) .553(25) 



6 COMMON /Two/ N.NNW. NV(AVEL» ISYM, T<^KIp» MT, P I F . GAMMA, M, TK j TP 



6 DIMENSION A(50.50)»8(50>4)iERASEt5T) 



6 REAL K 



6 N2 « N/2 



12 DO 1 I.l.M 



13 DO 6 12' 1.4 



14 RK(I,I2) • 0.0 

 21 6 RL(1,I2) « 0.0 

 27 DO 4 IC • l>5 



31 4 IF ( I .EO. JC(IC)) GO TO 1 



37 IFdSYM .Nt. 1) GO TO 8 



42 1F(I .6T. N2 .ANl). I .LE. N) GU TH <» 



65 8 Xll - X(I) - X5(l) 

 63 X21 • xn +XB(1) 



66 Y2l«Y<I)*Y3(l) 



74 P01-ALOG(X11**2+Y21**£) 



107 T01-ATAN2(Y21,Xii) 



113 CALL CPVtXll,Y21,£2l.Cll,Sll»A911»AlOn,K) 



124 C21-C11 



126 S21-S11 



127 AQ21'=A911 



131 Al021»A10n 



132 DO 7 J»1»N 



135 X12»X(1)-X8( J+1) 



143 Y22«Y(I)+YB( J*l) 



147 PQ2«ALOG(X12**2+Y22**2> 



163 TQ2-ATAN2(Y22,X12) 



167 S3-SS3(J) 



172 C3«CC3(J) 



175 CALL CPV(X12,Y22»E22,C12,S12»A912.ilT12,K) 



206 DO 13 IC = 1,3 



211 13 IF(J .EQ. JCdO) GO TO 41 



217 A3=A1011-A1012 



221 A4=£21*S11-E22*S12 



22 5 A7=»S3*(0.5*(P31-PO2)+A912-A'511)+C3*(TQl-TO? + 4 10n-A1012) 



243 A8 = E21*SIN(K*X11-ANG( J) )-F22*S IN ( K*X12-ANG ( J) ) 



266 5 Rid, J) • 2. * A3 + RU2(I,J1 



300 IF (I .NE. J) GO TO 3 



303 Rid, J) = Rid, J) - TP 



314 3 RJd.J) » -TP*A4 



322 POT(I,J) ' TK+A7 ♦ RK56(I,JI 



334 HDW(!,J) • -TK*PIE*A6 



342 DO 10 L ■ i,3 



344 RKd,L) ■ RKd.L) + POTd,J) ♦ RNHCMCJ.L) 



Table D-2. Continued 



122 



