1*-----------------------------------------------------------------------
2      SUBROUTINE RADB3 (IDO,L1,CC,CH,WA1,WA2)
3      DIMENSION       CC(IDO,3,L1)           ,CH(IDO,L1,3)           ,
4     1                WA1(*)     ,WA2(*)
5      DATA TAUR,TAUI /-.5,.866025403784439/
6      DO 101 K=1,L1
7         TR2 = CC(IDO,2,K)+CC(IDO,2,K)
8         CR2 = CC(1,1,K)+TAUR*TR2
9         CH(1,K,1) = CC(1,1,K)+TR2
10         CI3 = TAUI*(CC(1,3,K)+CC(1,3,K))
11         CH(1,K,2) = CR2-CI3
12         CH(1,K,3) = CR2+CI3
13  101 CONTINUE
14      IF (IDO .EQ. 1) RETURN
15      IDP2 = IDO+2
16      DO 103 K=1,L1
17         DO 102 I=3,IDO,2
18            IC = IDP2-I
19            TR2 = CC(I-1,3,K)+CC(IC-1,2,K)
20            CR2 = CC(I-1,1,K)+TAUR*TR2
21            CH(I-1,K,1) = CC(I-1,1,K)+TR2
22            TI2 = CC(I,3,K)-CC(IC,2,K)
23            CI2 = CC(I,1,K)+TAUR*TI2
24            CH(I,K,1) = CC(I,1,K)+TI2
25            CR3 = TAUI*(CC(I-1,3,K)-CC(IC-1,2,K))
26            CI3 = TAUI*(CC(I,3,K)+CC(IC,2,K))
27            DR2 = CR2-CI3
28            DR3 = CR2+CI3
29            DI2 = CI2+CR3
30            DI3 = CI2-CR3
31            CH(I-1,K,2) = WA1(I-2)*DR2-WA1(I-1)*DI2
32            CH(I,K,2) = WA1(I-2)*DI2+WA1(I-1)*DR2
33            CH(I-1,K,3) = WA2(I-2)*DR3-WA2(I-1)*DI3
34            CH(I,K,3) = WA2(I-2)*DI3+WA2(I-1)*DR3
35  102    CONTINUE
36  103 CONTINUE
37      RETURN
38      END
39