source: trunk/WRF.COMMON/WRFV3/external/fftpack/fftpack5/cosq1i.F @ 2759

Last change on this file since 2759 was 2759, checked in by aslmd, 2 years ago

adding unmodified code from WRFV3.0.1.1, expurged from useless data +1M size

File size: 1.7 KB
Line 
1!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
2!                                                                       
3!   FFTPACK 5.0                                                         
4!   Copyright (C) 1995-2004, Scientific Computing Division,             
5!   University Corporation for Atmospheric Research                     
6!   Licensed under the GNU General Public License (GPL)                 
7!                                                                       
8!   Authors:  Paul N. Swarztrauber and Richard A. Valent               
9!                                                                       
10!   $Id: cosq1i.f,v 1.2 2004/06/15 21:08:32 rodney Exp $               
11!                                                                       
12!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
13                                                                       
14      SUBROUTINE COSQ1I (N, WSAVE, LENSAV, IER)
15      INTEGER    N, LENSAV, IER
16      REAL       WSAVE(LENSAV)
17!                                                                       
18      IER = 0
19      IF (LENSAV .LT. 2*N + INT(LOG(REAL(N))) +4) THEN
20        IER = 2
21        CALL XERFFT ('COSQ1I', 3)
22        GO TO 300
23      ENDIF
24!                                                                       
25      PIH = 2.*ATAN(1.)
26      DT = PIH/FLOAT(N)
27      FK = 0.
28      DO 101 K=1,N
29         FK = FK+1.
30         WSAVE(K) = COS(FK*DT)
31  101 END DO
32      LNSV = N + INT(LOG(REAL(N))) +4
33      CALL RFFT1I (N, WSAVE(N+1), LNSV, IER1)
34      IF (IER1 .NE. 0) THEN
35        IER = 20
36        CALL XERFFT ('COSQ1I',-5)
37      ENDIF
38  300 CONTINUE
39      RETURN
40      END                                           
Note: See TracBrowser for help on using the repository browser.