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