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