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: sinq1b.f,v 1.2 2004/06/15 21:29:20 rodney Exp $
12 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
14 SUBROUTINE SINQ1B ( N, INC, X, LENX, WSAVE, LENSAV, &
16 INTEGER N, INC, LENX, LENSAV, LENWRK, IER
17 REAL X(INC,*), WSAVE(LENSAV), WORK(LENWRK)
21 IF (LENX .LT. INC*(N-1) + 1) THEN
23 CALL XERFFT ('SINQ1B', 6)
24 ELSEIF (LENSAV .LT. 2*N + INT(LOG(REAL(N))) +4) THEN
26 CALL XERFFT ('SINQ1B', 8)
27 ELSEIF (LENWRK .LT. N) THEN
29 CALL XERFFT ('SINQ1B', 10)
32 IF (N .GT. 1) GO TO 101
39 CALL COSQ1B (N,INC,X,LENX,WSAVE,LENSAV,WORK,LENWRK,IER1)
42 CALL XERFFT ('SINQ1B',-5)