wrfv2_fire/external/fftpack/fftpack5/cosqmi.F

   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: cosqmi.f,v 1.2 2004/06/15 21:14:57 rodney Exp $
  11 !
  12 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  13
  14       SUBROUTINE COSQMI (N, WSAVE, LENSAV, IER)
  15       INTEGER    N, LENSAV, IER
  16       REAL       WSAVE(LENSAV)
  17 !
  18       IER = 0
  19 !
  20       IF (LENSAV .LT. 2*N + INT(LOG(REAL(N))) +4) THEN
  21         IER = 2
  22         CALL XERFFT ('COSQMI', 3)
  23         GO TO 300
  24       ENDIF
  25 !
  26       PIH = 2.*ATAN(1.)
  27       DT = PIH/FLOAT(N)
  28       FK = 0.
  29       DO 101 K=1,N
  30          FK = FK+1.
  31          WSAVE(K) = COS(FK*DT)
  32   101 END DO
  33       LNSV = N + INT(LOG(REAL(N))) +4
  34       CALL RFFTMI (N, WSAVE(N+1), LNSV, IER1)
  35       IF (IER1 .NE. 0) THEN
  36         IER = 20
  37         CALL XERFFT ('COSQMI',-5)
  38       ENDIF
  39   300 CONTINUE
  40       RETURN
  41       END