Last change
on this file since 2759 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.