1 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
4 ! Copyright (C) 1995-2004, Scientific Computing Division,
5 ! University Corporation for Atmospheric Research
6 ! Licensed under the GNU General Public License (GPL)
8 ! Authors: Paul N. Swarztrauber and Richard A. Valent
10 ! $Id: cost1i.f,v 1.2 2004/06/15 21:14:57 rodney Exp $
12 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
14 SUBROUTINE COST1I (N, WSAVE, LENSAV, IER)
15 INTEGER N, LENSAV, IER
20 IF (LENSAV .LT. 2*N + INT(LOG(REAL(N))) +4) THEN
22 CALL XERFFT ('COST1I', 3)
36 WSAVE(K) = 2.*SIN(FK*DT)
37 WSAVE(KC) = 2.*COS(FK*DT)
39 LNSV = NM1 + INT(LOG(REAL(NM1))) +4
40 CALL RFFT1I (NM1, WSAVE(N+1), LNSV, IER1)
43 CALL XERFFT ('COST1I',-5)