source: lmdz_wrf/WRFV3/external/fftpack/fftpack5/cfft2i.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: cfft2i.f,v 1.2 2004/06/15 21:08:32 rodney Exp $               
8!                                                                       
9!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
10                                                                       
11      SUBROUTINE CFFT2I (L, M, WSAVE, LENSAV, IER)
12      INTEGER L, M, IER
13      REAL WSAVE(LENSAV)
14!                                                                       
15! Initialize error return                                               
16!                                                                       
17      IER = 0
18!                                                                       
19      IF (LENSAV .LT. 2*L + INT(LOG(REAL(L))) +                         &
20     &                    2*M + INT(LOG(REAL(M))) +8) THEN             
21        IER = 2
22        CALL XERFFT ('CFFT2I', 4)
23        GO TO 100
24      ENDIF
25!                                                                       
26      CALL CFFTMI (L, WSAVE(1), 2*L + INT(LOG(REAL(L))) + 4, IER1)
27      IF (IER1 .NE. 0) THEN
28        IER = 20
29        CALL XERFFT ('CFFT2I',-5)
30        GO TO 100
31      ENDIF
32      CALL CFFTMI (M, WSAVE(2*L+INT(LOG(REAL(L))*LOG(2.)) + 3),         &
33     &            2*M + INT(LOG(REAL(M))) + 4, IER1)                   
34      IF (IER1 .NE. 0) THEN
35        IER = 20
36        CALL XERFFT ('CFFT2I',-5)
37      ENDIF
38!                                                                       
39  100 CONTINUE
40      RETURN
41      END                                           
Note: See TracBrowser for help on using the repository browser.