source: lmdz_wrf/WRFV3/external/fftpack/fftpack5/rfft2i.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.1 KB
RevLine 
[1]1!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
2!                                                                       
3!   FFTPACK 5.0                                                         
4!                                                                       
5!   Authors:  Paul N. Swarztrauber and Richard A. Valent               
6!                                                                       
7!   $Id: rfft2i.f,v 1.3 2004/06/15 21:29:20 rodney Exp $               
8!                                                                       
9!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
10                                                                       
11      SUBROUTINE RFFT2I (L, M, WSAVE, LENSAV, IER)
12      INTEGER L, M, LENSAV, IER
13      INTEGER LWSAV,MWSAV
14      REAL WSAVE(LENSAV)
15!                                                                       
16! Initialize IER                                                       
17!                                                                       
18      IER = 0
19!                                                                       
20! Verify LENSAV                                                         
21!                                                                       
22      LWSAV =   L+INT(LOG(REAL(L)))+4
23      MWSAV =   2*M+INT(LOG(REAL(M)))+4
24      IF (LENSAV .LT. LWSAV+MWSAV) THEN
25        IER = 2
26        CALL XERFFT ('RFFT2I', 4)
27        GO TO 100
28      ENDIF
29!                                                                       
30      CALL RFFTMI (L, WSAVE(1), L + INT(LOG(REAL(L))) + 4, IER1)
31      IF (IER1 .NE. 0) THEN
32        IER = 20
33        CALL XERFFT ('RFFT2I',-5)
34        GO TO 100
35      ENDIF
36      CALL CFFTMI (M, WSAVE(L+INT(LOG(REAL(L)))+5),                     &
37     &            2*M+INT(LOG(REAL(M)))+4,IER1)                         
38      IF (IER1 .NE. 0) THEN
39        IER = 20
40        CALL XERFFT ('RFFT2I',-5)
41      ENDIF
42!                                                                       
43  100 CONTINUE
44      RETURN
45      END                                           
Note: See TracBrowser for help on using the repository browser.