Last change
on this file since 1 was
1,
checked in by lfita, 10 years ago
|
- -- --- Opening of the WRF+LMDZ coupling repository --- -- -
WRF: version v3.3
LMDZ: version v1818
More details in:
|
File size:
1.3 KB
|
Line | |
---|
1 | !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
---|
2 | ! |
---|
3 | ! FFTPACK 5.0 |
---|
4 | ! |
---|
5 | ! Authors: Paul N. Swarztrauber and Richard A. Valent |
---|
6 | ! |
---|
7 | ! $Id: factor.f,v 1.2 2004/06/15 21:29:19 rodney Exp $ |
---|
8 | ! |
---|
9 | !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
---|
10 | |
---|
11 | SUBROUTINE FACTOR (N,NF,FAC) |
---|
12 | REAL FAC(*) |
---|
13 | INTEGER NTRYH(4) |
---|
14 | DATA NTRYH(1),NTRYH(2),NTRYH(3),NTRYH(4)/4,2,3,5/ |
---|
15 | ! |
---|
16 | NL = N |
---|
17 | NF = 0 |
---|
18 | J = 0 |
---|
19 | 101 J = J+1 |
---|
20 | IF (J-4) 102,102,103 |
---|
21 | 102 NTRY = NTRYH(J) |
---|
22 | GO TO 104 |
---|
23 | 103 NTRY = NTRY+2 |
---|
24 | 104 NQ = NL/NTRY |
---|
25 | NR = NL-NTRY*NQ |
---|
26 | IF (NR) 101,105,101 |
---|
27 | 105 NF = NF+1 |
---|
28 | FAC(NF) = NTRY |
---|
29 | NL = NQ |
---|
30 | IF (NL .NE. 1) GO TO 104 |
---|
31 | RETURN |
---|
32 | END |
---|
Note: See
TracBrowser
for help on using the repository browser.