|
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.6 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: tables.f,v 1.2 2004/06/15 21:29:20 rodney Exp $ |
|---|
| 11 | ! |
|---|
| 12 | !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
|---|
| 13 | |
|---|
| 14 | SUBROUTINE TABLES (IDO,IP,WA) |
|---|
| 15 | REAL WA(IDO,IP-1,2) |
|---|
| 16 | ! |
|---|
| 17 | TPI = 8.*ATAN(1.) |
|---|
| 18 | ARGZ = TPI/REAL(IP) |
|---|
| 19 | ARG1 = TPI/REAL(IDO*IP) |
|---|
| 20 | DO 110 J=2,IP |
|---|
| 21 | ARG2 = REAL(J-1)*ARG1 |
|---|
| 22 | DO 100 I=1,IDO |
|---|
| 23 | ARG3 = REAL(I-1)*ARG2 |
|---|
| 24 | WA(I,J-1,1) = COS(ARG3) |
|---|
| 25 | WA(I,J-1,2) = SIN(ARG3) |
|---|
| 26 | 100 CONTINUE |
|---|
| 27 | IF (IP .LE. 5) GO TO 110 |
|---|
| 28 | ARG4 = REAL(J-1)*ARGZ |
|---|
| 29 | WA(1,J-1,1) = COS(ARG4) |
|---|
| 30 | WA(1,J-1,2) = SIN(ARG4) |
|---|
| 31 | 110 END DO |
|---|
| 32 | RETURN |
|---|
| 33 | END |
|---|
Note: See
TracBrowser
for help on using the repository browser.
