2 SUBROUTINE RADF5 (IDO, L1, CC, CH, WA1, WA2, WA3, WA4)
3 C***BEGIN PROLOGUE RADF5
5 C***PURPOSE Calculate the fast Fourier transform of subvectors of
7 C***LIBRARY SLATEC (FFTPACK)
8 C***TYPE SINGLE PRECISION (RADF5-S)
9 C***AUTHOR Swarztrauber, P. N., (NCAR)
10 C***ROUTINES CALLED (NONE)
11 C***REVISION HISTORY (YYMMDD)
13 C 830401 Modified to use SLATEC library source file format.
14 C 860115 Modified by Ron Boisvert to adhere to Fortran 77 by
15 C (a) changing dummy array size declarations (1) to (*),
16 C (b) changing definition of variables PI, TI11, TI12,
17 C TR11, TR12 by using FORTRAN intrinsic functions ATAN
18 C and SIN instead of DATA statements.
19 C 881128 Modified by Dick Valent to meet prologue standards.
20 C 890831 Modified array declarations. (WRB)
21 C 891214 Prologue converted to Version 4.0 format. (BAB)
22 C 900402 Added TYPE section. (WRB)
23 C***END PROLOGUE RADF5
24 DIMENSION CC(IDO,L1,5), CH(IDO,5,*), WA1(*), WA2(*), WA3(*),
26 C***FIRST EXECUTABLE STATEMENT RADF5
33 CR2 = CC(1,K,5)+CC(1,K,2)
34 CI5 = CC(1,K,5)-CC(1,K,2)
35 CR3 = CC(1,K,4)+CC(1,K,3)
36 CI4 = CC(1,K,4)-CC(1,K,3)
37 CH(1,1,K) = CC(1,K,1)+CR2+CR3
38 CH(IDO,2,K) = CC(1,K,1)+TR11*CR2+TR12*CR3
39 CH(1,3,K) = TI11*CI5+TI12*CI4
40 CH(IDO,4,K) = CC(1,K,1)+TR12*CR2+TR11*CR3
41 CH(1,5,K) = TI12*CI5-TI11*CI4
43 IF (IDO .EQ. 1) RETURN
45 IF((IDO-1)/2.LT.L1) GO TO 104
50 DR2 = WA1(I-2)*CC(I-1,K,2)+WA1(I-1)*CC(I,K,2)
51 DI2 = WA1(I-2)*CC(I,K,2)-WA1(I-1)*CC(I-1,K,2)
52 DR3 = WA2(I-2)*CC(I-1,K,3)+WA2(I-1)*CC(I,K,3)
53 DI3 = WA2(I-2)*CC(I,K,3)-WA2(I-1)*CC(I-1,K,3)
54 DR4 = WA3(I-2)*CC(I-1,K,4)+WA3(I-1)*CC(I,K,4)
55 DI4 = WA3(I-2)*CC(I,K,4)-WA3(I-1)*CC(I-1,K,4)
56 DR5 = WA4(I-2)*CC(I-1,K,5)+WA4(I-1)*CC(I,K,5)
57 DI5 = WA4(I-2)*CC(I,K,5)-WA4(I-1)*CC(I-1,K,5)
66 CH(I-1,1,K) = CC(I-1,K,1)+CR2+CR3
67 CH(I,1,K) = CC(I,K,1)+CI2+CI3
68 TR2 = CC(I-1,K,1)+TR11*CR2+TR12*CR3
69 TI2 = CC(I,K,1)+TR11*CI2+TR12*CI3
70 TR3 = CC(I-1,K,1)+TR12*CR2+TR11*CR3
71 TI3 = CC(I,K,1)+TR12*CI2+TR11*CI3
72 TR5 = TI11*CR5+TI12*CR4
73 TI5 = TI11*CI5+TI12*CI4
74 TR4 = TI12*CR5-TI11*CR4
75 TI4 = TI12*CI5-TI11*CI4
77 CH(IC-1,2,K) = TR2-TR5
81 CH(IC-1,4,K) = TR3-TR4
91 DR2 = WA1(I-2)*CC(I-1,K,2)+WA1(I-1)*CC(I,K,2)
92 DI2 = WA1(I-2)*CC(I,K,2)-WA1(I-1)*CC(I-1,K,2)
93 DR3 = WA2(I-2)*CC(I-1,K,3)+WA2(I-1)*CC(I,K,3)
94 DI3 = WA2(I-2)*CC(I,K,3)-WA2(I-1)*CC(I-1,K,3)
95 DR4 = WA3(I-2)*CC(I-1,K,4)+WA3(I-1)*CC(I,K,4)
96 DI4 = WA3(I-2)*CC(I,K,4)-WA3(I-1)*CC(I-1,K,4)
97 DR5 = WA4(I-2)*CC(I-1,K,5)+WA4(I-1)*CC(I,K,5)
98 DI5 = WA4(I-2)*CC(I,K,5)-WA4(I-1)*CC(I-1,K,5)
107 CH(I-1,1,K) = CC(I-1,K,1)+CR2+CR3
108 CH(I,1,K) = CC(I,K,1)+CI2+CI3
109 TR2 = CC(I-1,K,1)+TR11*CR2+TR12*CR3
110 TI2 = CC(I,K,1)+TR11*CI2+TR12*CI3
111 TR3 = CC(I-1,K,1)+TR12*CR2+TR11*CR3
112 TI3 = CC(I,K,1)+TR12*CI2+TR11*CI3
113 TR5 = TI11*CR5+TI12*CR4
114 TI5 = TI11*CI5+TI12*CI4
115 TR4 = TI12*CR5-TI11*CR4
116 TI4 = TI12*CI5-TI11*CI4
117 CH(I-1,3,K) = TR2+TR5
118 CH(IC-1,2,K) = TR2-TR5
121 CH(I-1,5,K) = TR3+TR4
122 CH(IC-1,4,K) = TR3-TR4