source: lmdz_wrf/WRFV3/external/fftpack/fftpack5/cfftmb.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.0 KB
Line 
1!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
2!                                                                       
3!   FFTPACK 5.0                                                         
4!                                                                       
5!   Authors:  Paul N. Swarztrauber and Richard A. Valent               
6!                                                                       
7!   $Id: cfftmb.f,v 1.2 2004/06/15 21:08:32 rodney Exp $               
8!                                                                       
9!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
10                                                                       
11      SUBROUTINE CFFTMB (LOT, JUMP, N, INC, C, LENC, WSAVE, LENSAV,     &
12     &                  WORK, LENWRK, IER)                             
13      INTEGER  LOT, JUMP, N, INC, LENC, LENSAV, LENWRK, IER
14      COMPLEX       C(LENC)
15      REAL     WSAVE(LENSAV)     ,WORK(LENWRK)
16      LOGICAL XERCON
17!                                                                       
18      IER = 0
19!                                                                       
20      IF (LENC .LT. (LOT-1)*JUMP + INC*(N-1) + 1) THEN
21        IER = 1
22        CALL XERFFT ('CFFTMB ', 6)
23      ELSEIF (LENSAV .LT. 2*N + INT(LOG(REAL(N))) + 4) THEN
24        IER = 2
25        CALL XERFFT ('CFFTMB ', 8)
26      ELSEIF (LENWRK .LT. 2*LOT*N) THEN
27        IER = 3
28        CALL XERFFT ('CFFTMB ', 10)
29      ELSEIF (.NOT. XERCON(INC,JUMP,N,LOT)) THEN
30        IER = 4
31        CALL XERFFT ('CFFTMB ', -1)
32      ENDIF
33!                                                                       
34      IF (N .EQ. 1) RETURN
35!                                                                       
36      IW1 = N+N+1
37      CALL CMFM1B (LOT,JUMP,N,INC,C,WORK,WSAVE,WSAVE(IW1),              &
38     &                           WSAVE(IW1+1))                         
39      RETURN
40      END                                           
Note: See TracBrowser for help on using the repository browser.