Changeset 1249 for LMDZ4/branches/LMDZ4-dev/libf/filtrez

Timestamp:

Oct 21, 2009, 6:04:08 PM (15 years ago)

Author:

yann meurdesoif

Message:

Corrections de Bug divers - portage vers Titane (CCRT) -
YM

Location:

LMDZ4/branches/LMDZ4-dev/libf/filtrez

Files:

• filtreg_mod.F90 (modified) (2 diffs)
• mod_fft_mkl.F90 (modified) (2 diffs)
• mod_filtre_fft.F90 (modified) (3 diffs)

LMDZ4/branches/LMDZ4-dev/libf/filtrez/filtreg_mod.F90

@@ -7 +7 @@
 
   SUBROUTINE inifilr
+  USE mod_filtre_fft
     !
     !    ... H. Upadhyaya, O.Sharma   ...
@@ -519 +520 @@
     ENDDO
 
+    IF (use_filtre_fft) THEN
+       CALL Init_filtre_fft(coefilu,modfrstu,jfiltnu,jfiltsu,  &
+                           coefilv,modfrstv,jfiltnv,jfiltsv)
+    ENDIF
+
     !   ...................................................................
 

LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_fft_mkl.F90

@@ -74 +74 @@
       ierr = DftiSetValue(FFT_Handle,DFTI_PLACEMENT,DFTI_NOT_INPLACE)
       ierr = DftiSetValue(FFT_Handle, DFTI_INPUT_DISTANCE, vsize+inc)
-      ierr = DftiSetValue(FFT_Handle, DFTI_OUTPUT_DISTANCE, vsize/2+1)
+      ierr = DftiSetValue(FFT_Handle, DFTI_OUTPUT_DISTANCE, (vsize/2+1)*2)
       ierr = DftiCommitDescriptor( FFT_Handle )
 !      Forward_handle(nb_vect)%IsAllocated=.TRUE.
@@ -114 +114 @@
       ierr = DftiSetValue(FFT_Handle,DFTI_BACKWARD_SCALE,scale_factor)
       ierr = DftiSetValue(FFT_Handle,DFTI_PLACEMENT,DFTI_NOT_INPLACE)
-      ierr = DftiSetValue(FFT_Handle, DFTI_INPUT_DISTANCE,  vsize/2+1)
+      ierr = DftiSetValue(FFT_Handle, DFTI_INPUT_DISTANCE,  (vsize/2+1)*2)
       ierr = DftiSetValue(FFT_Handle, DFTI_OUTPUT_DISTANCE, vsize+inc)
       ierr = DftiCommitDescriptor( FFT_Handle )

LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_filtre_fft.F90

@@ -118 +118 @@
 
     REAL               :: vect(iim+inc,jj_end-jj_begin+1,nbniv)
-    COMPLEX            :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv)
+!    REAL               :: vect_test(iim+inc,jj_end-jj_begin+1,nbniv)
+    COMPLEX*16         :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv)
+!    COMPLEX*16         :: TF_vect_test(iim/2+1,jj_end-jj_begin+1,nbniv)
     INTEGER            :: nb_vect
     INTEGER :: i,j,l
@@ -142 +144 @@
     CALL FFT_forward(vect,TF_vect,nb_vect)
 
+!    CALL FFT_forward(vect,TF_vect_test,nb_vect)
+!      PRINT *,"XXXXXXXXXXXXX Filtre_u_FFT xxxxxxxxxxxx"
+!      DO j=1,jj_end-jj_begin+1
+!      DO i=1,iim/2+1
+!         PRINT *,"====",i,j,"----->",TF_vect_test(i,j,1)
+!       ENDDO
+!      ENDDO
+
     DO l=1,ll_nb
       DO j=1,jj_end-jj_begin+1
@@ -149 +159 @@
       ENDDO
     ENDDO
-  
+       
     CALL FFT_backward(TF_vect,vect,nb_vect)
-      
+!    CALL FFT_backward(TF_vect_test,vect_test,nb_vect)
+          
 !      PRINT *,"XXXXXXXXXXXXX Filtre_u_FFT xxxxxxxxxxxx"
 !      DO j=1,jj_end-jj_begin+1
-!          DO i=1,iim
-!          PRINT *,"====",j,"----->",(vect_tmp(i,j,1)-vect(i,j,1))/ &
-!                                    ((vect_tmp(i,j,1)+vect(i,j,1))*0.5+1e-30)
-!          ENDDO
+!         DO i=1,iim
+!           PRINT *,"====",i,j,"----->",vect_test(i,j,1)
+!         ENDDO
 !      ENDDO
-      
+
     ll_nb=0
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)