source: LMDZ6/branches/Amaury_dev/libf/filtrez/lmdz_fft_wrapper.f90

MODULE lmdz_fft_wrapper
  IMPLICIT NONE; PRIVATE
  PUBLIC vsize, inc, init_fft, fft_forward, fft_backward

  INTEGER,SAVE             :: vsize
  INTEGER,PARAMETER        :: inc=1

CONTAINS
  
  SUBROUTINE init_fft(iim,nb)
  IMPLICIT NONE
    INTEGER :: iim
    INTEGER :: nb
    
    PRINT *, "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique"
    stop 1
  END SUBROUTINE init_fft
  
  
  SUBROUTINE fft_forward(vect,TF_vect,nb_vect)
    IMPLICIT NONE
    INTEGER,INTENT(IN)  :: nb_vect
    REAL,INTENT(IN)     :: vect(vsize+inc,nb_vect)
    COMPLEX,INTENT(INOUT) :: TF_vect(vsize/2+1,nb_vect)
    
    PRINT *, "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique"
    stop 1
    
  END SUBROUTINE fft_forward
  
  SUBROUTINE fft_backward(TF_vect,vect,nb_vect)
    IMPLICIT NONE
    INTEGER,INTENT(IN)  :: nb_vect
    REAL,INTENT(INOUT)    :: vect(vsize+inc,nb_vect)
    COMPLEX,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect)
  
    PRINT *, "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique"
    stop 1
    
  END SUBROUTINE fft_backward
  
END MODULE lmdz_fft_wrapper