Context Navigation

r1f2kb.F @ 1

Last change on this file since 1 was 1, checked in by lfita, 10 years ago

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 2.0 KB

Line
1	!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
2	!
3	! FFTPACK 5.0
4	!
5	! Authors: Paul N. Swarztrauber and Richard A. Valent
6	!
7	! $Id: r1f2kb.f,v 1.2 2004/06/15 21:29:20 rodney Exp $
8	!
9	!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
10
11	SUBROUTINE R1F2KB (IDO,L1,CC,IN1,CH,IN2,WA1)
12	REAL CC(IN1,IDO,2,L1), CH(IN2,IDO,L1,2), WA1(IDO)
13	!
14	DO 101 K=1,L1
15	CH(1,1,K,1) = CC(1,1,1,K)+CC(1,IDO,2,K)
16	CH(1,1,K,2) = CC(1,1,1,K)-CC(1,IDO,2,K)
17	101 END DO
18	IF (IDO-2) 107,105,102
19	102 IDP2 = IDO+2
20	DO 104 K=1,L1
21	DO 103 I=3,IDO,2
22	IC = IDP2-I
23
24	CH(1,I-1,K,1) = CC(1,I-1,1,K)+CC(1,IC-1,2,K)
25	CH(1,I,K,1) = CC(1,I,1,K)-CC(1,IC,2,K)
26
27	CH(1,I-1,K,2) = WA1(I-2)*(CC(1,I-1,1,K)-CC(1,IC-1,2,K)) &
28	& -WA1(I-1)*(CC(1,I,1,K)+CC(1,IC,2,K))
29	CH(1,I,K,2) = WA1(I-2)*(CC(1,I,1,K)+CC(1,IC,2,K))+WA1(I-1) &
30	& *(CC(1,I-1,1,K)-CC(1,IC-1,2,K))
31
32	103 CONTINUE
33	104 END DO
34	IF (MOD(IDO,2) .EQ. 1) RETURN
35	105 DO 106 K=1,L1
36	CH(1,IDO,K,1) = CC(1,IDO,1,K)+CC(1,IDO,1,K)
37	CH(1,IDO,K,2) = -(CC(1,1,2,K)+CC(1,1,2,K))
38	106 END DO
39	107 RETURN
40	END

Note: See TracBrowser for help on using the repository browser.