C
C *** FOURIER TRNS. CHECK
C *** WAVE READ, FOURIER TRANS., INVERSE FOURIE TRANS., AND WAVE PLOT
C
      PROGRAM MFOURCK
C
C *** TITL   : TITLE OF EACH CHANNEL
C *** WV     : WAVE AMPLITUDE DATA
C *** NWV    : WAVE DATA NUMBER
C *** MWV    : NWV IN MAIN PROGRAM DIMENSION
C            : IF MWV > 2**13(=8192), CAN NOT USE IN SUBROUTINE !!!
C *** DT     : DELTA TIME ( SEC ) 
C            : DURATION = (NWV-1)*DT IN PLOTTER, NWV*DT IN FFT
C *** FOUR(I): FOURIER COEFFICIENTS
C   ( FOURIER SPECTRUM = FOURIER COEFFICIENTS * DURATION/2.0 )
C
      PARAMETER(MWV=12800)
      CHARACTER TITL*50
      DIMENSION WV(MWV)
      COMPLEX*8 FOUR(MWV)
C
      CALL WVREAD(TITL,WV,NWV,MWV,DT)
      CALL WVFOUR(WV,NWV,MWV,FOUR,NFAST)
      CALL WVFAST(NFAST,FOUR,MWV,1)
C
      DO 100 I=1,NFAST
      WV(I)=REAL(FOUR(I))
  100 CONTINUE
C
      CALL WVPLOT(TITL,WV,NWV,MWV,DT)
C
      STOP
      END
         