wrfv2_fire/external/fftpack/fftpack5/cfft1b.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: cfft1b.f,v 1.2 2004/06/15 21:08:32 rodney Exp $
  11 !
  12 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  13
  14       SUBROUTINE CFFT1B (N, INC, C, LENC, WSAVE, LENSAV,                &
  15      &                  WORK, LENWRK, IER)
  16       INTEGER  N, INC, LENC, LENSAV, LENWRK, IER
  17       COMPLEX       C(LENC)
  18       REAL     WSAVE(LENSAV)     ,WORK(LENWRK)
  19 !
  20       IER = 0
  21 !
  22       IF (LENC .LT. INC*(N-1) + 1) THEN
  23         IER = 1
  24         CALL XERFFT ('CFFT1B ', 6)
  25       ELSEIF (LENSAV .LT. 2*N + INT(LOG(REAL(N))) + 4) THEN
  26         IER = 2
  27         CALL XERFFT ('CFFT1B ', 8)
  28       ELSEIF (LENWRK .LT. 2*N) THEN
  29         IER = 3
  30         CALL XERFFT ('CFFT1B ', 10)
  31       ENDIF
  32 !
  33       IF (N .EQ. 1) RETURN
  34 !
  35       IW1 = N+N+1
  36       CALL C1FM1B (N,INC,C,WORK,WSAVE,WSAVE(IW1),                       &
  37      &                           WSAVE(IW1+1))
  38       RETURN
  39       END