source: LMDZ4/branches/LMDZ4V5.0-dev/libf/filtrez/mod_fft_wrapper.F90 @ 2032

Last change on this file since 2032 was 1389, checked in by Laurent Fairhead, 15 years ago
  • Differing COMPLEX declarations were causing problems in FFT routines

compilation. The FFTs should only be used in double precision in any case

  • the ALLOCATE command for the o_trac variable was misplaced and called

several times (causing an error for some compilators)


  • Des déclarations COMPLEX différenciées causaient des problèmes de

compilation dans les routines des filtres FFT. Celles-ci ne devraient être
utilisées qu'en double précision de toutes façons.

  • L'ALLOCATE de la variable o_trac était mal placé et appelé plusieurs fois

(ce qui causait des crash pour certains compilateurs)

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.1 KB
Line 
1MODULE mod_fft_wrapper
2
3  INTEGER,SAVE             :: vsize
4  INTEGER,PARAMETER        :: inc=1
5
6CONTAINS
7 
8  SUBROUTINE Init_fft(iim,nb)
9  IMPLICIT NONE
10    INTEGER :: iim
11    INTEGER :: nb
12   
13    STOP "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique"
14  END SUBROUTINE Init_fft
15 
16 
17  SUBROUTINE fft_forward(vect,TF_vect,nb_vect)
18    IMPLICIT NONE
19    INTEGER,INTENT(IN)  :: nb_vect
20    REAL,INTENT(IN)     :: vect(vsize+inc,nb_vect)
21    COMPLEX,INTENT(INOUT) :: TF_vect(vsize/2+1,nb_vect)
22   
23    STOP "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique"
24   
25  END SUBROUTINE fft_forward
26 
27  SUBROUTINE fft_backward(TF_vect,vect,nb_vect)
28    IMPLICIT NONE
29    INTEGER,INTENT(IN)  :: nb_vect
30    REAL,INTENT(INOUT)    :: vect(vsize+inc,nb_vect)
31    COMPLEX,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect)
32 
33    STOP "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique"
34   
35  END SUBROUTINE fft_backward
36 
37END MODULE mod_fft_wrapper
Note: See TracBrowser for help on using the repository browser.