source: lmdz_wrf/WRFV3/external/fftpack/fftpack5/sinq1b.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: sinq1b.f,v 1.2 2004/06/15 21:29:20 rodney Exp $               
8!                                                                       
9!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
10                                                                       
11      SUBROUTINE SINQ1B ( N, INC, X, LENX, WSAVE, LENSAV,               &
12     &                   WORK, LENWRK, IER)                             
13      INTEGER    N, INC, LENX, LENSAV, LENWRK, IER
14      REAL       X(INC,*), WSAVE(LENSAV), WORK(LENWRK)
15!                                                                       
16      IER = 0
17!                                                                       
18      IF (LENX .LT. INC*(N-1) + 1) THEN
19        IER = 1
20        CALL XERFFT ('SINQ1B', 6)
21      ELSEIF (LENSAV .LT. 2*N + INT(LOG(REAL(N))) +4) THEN
22        IER = 2
23        CALL XERFFT ('SINQ1B', 8)
24      ELSEIF (LENWRK .LT. N) THEN
25        IER = 3
26        CALL XERFFT ('SINQ1B', 10)
27      ENDIF
28!                                                                       
29      IF (N .GT. 1) GO TO 101
30      X(1,1) = 4.*X(1,1)
31      RETURN
32  101 NS2 = N/2
33      DO 102 K=2,N,2
34         X(1,K) = -X(1,K)
35  102 END DO
36      CALL COSQ1B (N,INC,X,LENX,WSAVE,LENSAV,WORK,LENWRK,IER1)
37      IF (IER1 .NE. 0) THEN
38        IER = 20
39        CALL XERFFT ('SINQ1B',-5)
40        GO TO 300
41      ENDIF
42      DO 103 K=1,NS2
43         KC = N-K
44         XHOLD = X(1,K)
45         X(1,K) = X(1,KC+1)
46         X(1,KC+1) = XHOLD
47  103 END DO
48  300 RETURN
49      END                                           
Note: See TracBrowser for help on using the repository browser.