Context Navigation

cmfm1b.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.4 KB

Line
1	!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
2	!
3	! FFTPACK 5.0
4	!
5	! Authors: Paul N. Swarztrauber and Richard A. Valent
6	!
7	! $Id: cmfm1b.f,v 1.2 2004/06/15 21:08:32 rodney Exp $
8	!
9	!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
10
11	SUBROUTINE CMFM1B (LOT,JUMP,N,INC,C,CH,WA,FNF,FAC)
12	COMPLEX C(*)
13	REAL CH(), WA(), FAC(*)
14	!
15	! FFTPACK 5.0 auxiliary routine
16	!
17	NF = FNF
18	NA = 0
19	L1 = 1
20	IW = 1
21	DO 125 K1=1,NF
22	IP = FAC(K1)
23	L2 = IP*L1
24	IDO = N/L2
25	LID = L1*IDO
26	NBR = 1+NA+2*MIN(IP-2,4)
27	GO TO (52,62,53,63,54,64,55,65,56,66),NBR
28	52 CALL CMF2KB (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
29	GO TO 120
30	62 CALL CMF2KB (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
31	GO TO 120
32	53 CALL CMF3KB (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
33	GO TO 120
34	63 CALL CMF3KB (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
35	GO TO 120
36	54 CALL CMF4KB (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
37	GO TO 120
38	64 CALL CMF4KB (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
39	GO TO 120
40	55 CALL CMF5KB (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
41	GO TO 120
42	65 CALL CMF5KB (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
43	GO TO 120
44	56 CALL CMFGKB (LOT,IDO,IP,L1,LID,NA,C,C,JUMP,INC,CH,CH,1, &
45	& LOT,WA(IW))
46	GO TO 120
47	66 CALL CMFGKB (LOT,IDO,IP,L1,LID,NA,CH,CH,1,LOT,C,C, &
48	& JUMP,INC,WA(IW))
49	120 L1 = L2
50	IW = IW+(IP-1)*(IDO+IDO)
51	IF(IP .LE. 5) NA = 1-NA
52	125 END DO
53	RETURN
54	END

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