1 SUBROUTINE AGENR(BOAK,SB,A) 2C 3C*** GENERATES COEFFICIENTS FOR G/DELTA CALCULATIONS BY GDELTA 4C 5 COMMON /CONSNT/ PI,DEG,UNUSED,RAD 6 DIMENSION F1(4),F2(4),F3(4),F4(4),BOCO16(4),BOAK(4) 7 DIMENSION A(16),AL(16),C1(16),C2(16),C3(4) 8 DATA C1 /-.07612,.00001 ,.21677,.54120 , -.29289,-.21677,.00001 9 1 ,.32443 , -.61732,-.54120,-.32443,.00001 , 10 2 -1.0,-.92388,-.70711,-.38267/ 11 DATA C2/1.92388,1.84776,1.63099,1.30656 , 1.70711,1.63099,1.41422, 12 1 1.08979 , 1.38268,1.30655,1.08979,.76536 , 13 2 1.0,.92388,.70711,.38267/ 14 DATA F1 /5.2262,1.0360,0.,.11208/ 15 DATA F2 /1.91433,2.8284,.91418,0.0/ 16 DATA F3 /0.0,1.1944,2.1634,1.5772/ 17 DATA F4 /.14645,0.0,.85357,2.0/ 18 DATA C3 /.92388,.70711,.38268,.00001 / 19C 20 K=1 21 KP=0 22 TANSB=TAN(DEG*SB) 23 DO 1010 I=1,16 24 IF(KP.EQ.K)GO TO 1000 25 KP=K 26 BOCH=BOAK(K) 27 BOCO16(K)=-BOCH/16. 28 RCPLBC=1.0/BOCH 29 BOC2=BOCH**2 30 BOCTN=BOCH *TANSB 31 TBOCTN=2.0*BOCTN 32 DENOM=1.+C3(K)*TBOCTN 33 1000 TMP=(1.+C1(I)*BOCTN)**2 34 AL(I)=RCPLBC/C1(I)*(SQRT(TMP+BOC2*(C1(I)**2))-1.) 35 1 -RCPLBC/C2(I)*(SQRT(TMP+BOC2*(C2(I)**2))/DENOM-1.0) 36 2 -2.*TANSB*SQRT((1.+C3(K)*BOCTN)**2+BOC2*(C3(K)**2))/DENOM 37 IF(4*(I/4).EQ.I)K=K+1 38 1010 CONTINUE 39C 40 K=1 41 L=1 42 DO 1020 I=1,4 43 FST=-2.*F1(I) 44 IF(I.EQ.1)FST=-FST 45 A(L)=FST+BOCO16(I)*(2.6131*AL(K)+2.*(-.70711*AL(K+1)-.76537* 46 1 AL(K+2)+.20711*AL(K+3))) 47 L=L+1 48 FST=-2.*F2(I) 49 IF(I.EQ.2)FST=-FST 50 A(L)=FST+BOCO16(I)*(-1.4142*AL(K)+2.*(1.8478*AL(K+1)-.50000* 51 1 AL(K+2)-.76537*AL(K+3))) 52 L=L+1 53 FST=-2.*F3(I) 54 IF(I.EQ.3)FST=-FST 55 A(L)=FST+BOCO16(I)*(1.0824*AL(K)+2.*(-1.2071*AL(K+1)+1.8478 56 1 *AL(K+2)-.70711*AL(K+3))) 57 L=L+1 58 FST=-2.*F4(I) 59 IF(I.EQ.4)FST=-FST 60 A(L)=FST+BOCO16(I)*(-.5*AL(K)+1.0824*AL(K+1)-1.4142*AL(K+2) 61 1 +2.6131*AL(K+3)) 62 L=L+1 63 1020 K=L 64 RETURN 65 END 66