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: sint1i.f,v 1.2 2004/06/15 21:29:20 rodney Exp $
12 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
14 SUBROUTINE SINT1I (N, WSAVE, LENSAV, IER)
15 INTEGER N, LENSAV, IER
20 IF (LENSAV .LT. N/2 + N + INT(LOG(REAL(N))) +4) THEN
22 CALL XERFFT ('SINT1I', 3)
32 WSAVE(K) = 2.*SIN(K*DT)
34 LNSV = NP1 + INT(LOG(REAL(NP1))) +4
35 CALL RFFT1I (NP1, WSAVE(NS2+1), LNSV, IER1)
38 CALL XERFFT ('SINT1I',-5)