1 | MODULE lmdz_fft_wrapper |
---|
2 | IMPLICIT NONE; PRIVATE |
---|
3 | PUBLIC vsize, inc, init_fft, fft_forward, fft_backward |
---|
4 | |
---|
5 | INTEGER,SAVE :: vsize |
---|
6 | INTEGER,PARAMETER :: inc=1 |
---|
7 | |
---|
8 | CONTAINS |
---|
9 | |
---|
10 | SUBROUTINE init_fft(iim,nb) |
---|
11 | IMPLICIT NONE |
---|
12 | INTEGER :: iim |
---|
13 | INTEGER :: nb |
---|
14 | |
---|
15 | PRINT *, "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique" |
---|
16 | stop 1 |
---|
17 | END SUBROUTINE init_fft |
---|
18 | |
---|
19 | |
---|
20 | SUBROUTINE fft_forward(vect,TF_vect,nb_vect) |
---|
21 | IMPLICIT NONE |
---|
22 | INTEGER,INTENT(IN) :: nb_vect |
---|
23 | REAL,INTENT(IN) :: vect(vsize+inc,nb_vect) |
---|
24 | COMPLEX,INTENT(INOUT) :: TF_vect(vsize/2+1,nb_vect) |
---|
25 | |
---|
26 | PRINT *, "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique" |
---|
27 | stop 1 |
---|
28 | |
---|
29 | END SUBROUTINE fft_forward |
---|
30 | |
---|
31 | SUBROUTINE fft_backward(TF_vect,vect,nb_vect) |
---|
32 | IMPLICIT NONE |
---|
33 | INTEGER,INTENT(IN) :: nb_vect |
---|
34 | REAL,INTENT(INOUT) :: vect(vsize+inc,nb_vect) |
---|
35 | COMPLEX,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect) |
---|
36 | |
---|
37 | PRINT *, "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique" |
---|
38 | stop 1 |
---|
39 | |
---|
40 | END SUBROUTINE fft_backward |
---|
41 | |
---|
42 | END MODULE lmdz_fft_wrapper |
---|