Changeset 1389 for LMDZ4/branches/LMDZ4V5.0-dev/libf/filtrez/mod_fft_mathkeisan.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_mathkeisan.F90 (modified) (3 diffs)

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

@@ -15 +15 @@
     INTEGER :: nb_vect_max
     REAL    :: rtmp=1.
-    COMPLEX*16 :: ctmp
+    COMPLEX :: ctmp
     INTEGER :: itmp=1
     INTEGER :: isign=0
@@ -37 +37 @@
     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)
     REAL                :: work(4*vsize*nb_vect)
     INTEGER             :: ierr
@@ -51 +51 @@
     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)
     REAL                :: work(4*vsize*nb_vect)
     INTEGER             :: ierr