|
Last change
on this file was
2759,
checked in by aslmd, 2 years ago
|
|
adding unmodified code from WRFV3.0.1.1, expurged from useless data +1M size
|
|
File size:
1.7 KB
|
| Line | |
|---|
| 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: cosq1i.f,v 1.2 2004/06/15 21:08:32 rodney Exp $ |
|---|
| 11 | ! |
|---|
| 12 | !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
|---|
| 13 | |
|---|
| 14 | SUBROUTINE COSQ1I (N, WSAVE, LENSAV, IER) |
|---|
| 15 | INTEGER N, LENSAV, IER |
|---|
| 16 | REAL WSAVE(LENSAV) |
|---|
| 17 | ! |
|---|
| 18 | IER = 0 |
|---|
| 19 | IF (LENSAV .LT. 2*N + INT(LOG(REAL(N))) +4) THEN |
|---|
| 20 | IER = 2 |
|---|
| 21 | CALL XERFFT ('COSQ1I', 3) |
|---|
| 22 | GO TO 300 |
|---|
| 23 | ENDIF |
|---|
| 24 | ! |
|---|
| 25 | PIH = 2.*ATAN(1.) |
|---|
| 26 | DT = PIH/FLOAT(N) |
|---|
| 27 | FK = 0. |
|---|
| 28 | DO 101 K=1,N |
|---|
| 29 | FK = FK+1. |
|---|
| 30 | WSAVE(K) = COS(FK*DT) |
|---|
| 31 | 101 END DO |
|---|
| 32 | LNSV = N + INT(LOG(REAL(N))) +4 |
|---|
| 33 | CALL RFFT1I (N, WSAVE(N+1), LNSV, IER1) |
|---|
| 34 | IF (IER1 .NE. 0) THEN |
|---|
| 35 | IER = 20 |
|---|
| 36 | CALL XERFFT ('COSQ1I',-5) |
|---|
| 37 | ENDIF |
|---|
| 38 | 300 CONTINUE |
|---|
| 39 | RETURN |
|---|
| 40 | END |
|---|
Note: See
TracBrowser
for help on using the repository browser.
