Changeset 1389 for LMDZ4/branches/LMDZ4V5.0-dev/libf/filtrez/mod_fft_fftw.F90

Timestamp:

May 18, 2010, 9:48:01 AM (14 years ago)

Author:

Laurent Fairhead

Message:

• 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)

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/filtrez/mod_fft_fftw.F90 (modified) (2 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/filtrez/mod_fft_fftw.F90

@@ -89 +89 @@
     INTEGER,INTENT(IN)     :: nb_vect
     REAL,INTENT(IN)        :: vect(vsize+inc,nb_vect)
-    COMPLEX*16,INTENT(OUT) :: TF_vect(vsize/2+1,nb_vect)
+    COMPLEX,INTENT(OUT) :: TF_vect(vsize/2+1,nb_vect)
 
     CALL dfftw_execute_dft_r2c(plan_forward(nb_vect),vect,TF_vect)
@@ -102 +102 @@
     INTEGER,INTENT(IN)     :: nb_vect
     REAL,INTENT(OUT)       :: vect(vsize+inc,nb_vect)
-    COMPLEX*16,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect)
+    COMPLEX,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect)
 
     CALL dfftw_execute_dft_c2r(plan_backward(nb_vect),TF_vect,vect)