source: lmdz_wrf/WRFV3/external/fftpack/fftpack5/cmfm1f.F @ 1

Last change on this file since 1 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

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: cmfm1f.f,v 1.2 2004/06/15 21:08:32 rodney Exp $               
8!                                                                       
9!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
10                                                                       
11      SUBROUTINE CMFM1F (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 CMF2KF (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
29         GO TO 120
30   62    CALL CMF2KF (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
31         GO TO 120
32   53    CALL CMF3KF (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
33         GO TO 120
34   63    CALL CMF3KF (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
35         GO TO 120
36   54    CALL CMF4KF (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
37         GO TO 120
38   64    CALL CMF4KF (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
39         GO TO 120
40   55    CALL CMF5KF (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
41         GO TO 120
42   65    CALL CMF5KF (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
43         GO TO 120
44   56    CALL CMFGKF (LOT,IDO,IP,L1,LID,NA,C,C,JUMP,INC,CH,CH,          &
45     &     1,LOT,WA(IW))                                               
46         GO TO 120
47   66    CALL CMFGKF (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.