1 |
|
|
MODULE SUFFT_MOD |
2 |
|
|
CONTAINS |
3 |
|
|
SUBROUTINE SUFFT |
4 |
|
|
|
5 |
|
|
USE PARKIND1 ,ONLY : JPIM ,JPRB |
6 |
|
|
|
7 |
|
|
USE TPM_DIM |
8 |
|
|
USE TPM_GEN |
9 |
|
|
USE TPM_DISTR |
10 |
|
|
USE TPM_GEOMETRY |
11 |
|
|
USE TPM_FFT |
12 |
|
|
|
13 |
|
|
IMPLICIT NONE |
14 |
|
|
|
15 |
|
|
INTEGER(KIND=JPIM) :: JGL,IGLG |
16 |
|
|
LOGICAL :: LLP1,LLP2 |
17 |
|
|
|
18 |
|
|
! ------------------------------------------------------------------ |
19 |
|
|
|
20 |
|
|
LLP1 = NPRINTLEV>0 |
21 |
|
|
LLP2 = NPRINTLEV>1 |
22 |
|
|
IF(LLP1) WRITE(NOUT,*) '=== ENTER ROUTINE SUFFT ===' |
23 |
|
|
|
24 |
|
|
ALLOCATE(T%TRIGS(R%NDLON,D%NDGL_FS)) |
25 |
|
|
IF(LLP2)WRITE(NOUT,9) 'T%TRIGS ',SIZE(T%TRIGS),SHAPE(T%TRIGS) |
26 |
|
|
ALLOCATE(T%NFAX(19,D%NDGL_FS)) |
27 |
|
|
IF(LLP2)WRITE(NOUT,9) 'T%NFAX ',SIZE(T%NFAX),SHAPE(T%NFAX) |
28 |
|
|
|
29 |
|
|
DO JGL=1,D%NDGL_FS |
30 |
|
|
IGLG = D%NPTRLS(MYSETW)+JGL-1 |
31 |
|
|
CALL SET99(T%TRIGS(1,JGL),T%NFAX(1,JGL),G%NLOEN(IGLG)) |
32 |
|
|
ENDDO |
33 |
|
|
|
34 |
|
|
! ------------------------------------------------------------------ |
35 |
|
|
|
36 |
|
|
9 FORMAT(1X,'ARRAY ',A10,' ALLOCATED ',8I8) |
37 |
|
|
|
38 |
|
|
END SUBROUTINE SUFFT |
39 |
|
|
END MODULE SUFFT_MOD |