program ft4d

   include 'ft4_params.f90'

   character*8 arg
   character*17 cdatetime 
   character*512 data_dir
   character*37 decodes(100)
   character*16 fname
   character*12 mycall
   character*12 hiscall
   character*80 infile
   character*61 line
   
   real*8 fMHz

   integer ihdr(11)
   integer*2 iwave(NMAX)                 !Generated full-length waveform

   fs=12000.0/NDOWN                       !Sample rate
   dt=1/fs                                !Sample interval after downsample (s)
   tt=NSPS*dt                             !Duration of "itone" symbols (s)
   baud=1.0/tt                            !Keying rate for "itone" symbols (baud)
   txt=NZ*dt                              !Transmission length (s)

   nargs=iargc()
   if(nargs.lt.1) then
      print*,'Usage:   ft4d [-a <data_dir>] [-f fMHz] [-n nQSOProgress] file1 [file2 ...]'
      go to 999
   endif
   iarg=1
   data_dir="."
   call getarg(iarg,arg)
   if(arg(1:2).eq.'-a') then
      call getarg(iarg+1,data_dir)
      iarg=iarg+2
   endif
   call getarg(iarg,arg)
   if(arg(1:2).eq.'-f') then
      call getarg(iarg+1,arg)
      read(arg,*) fMHz
      iarg=iarg+2
   endif
   nQSOProgress=0
   if(arg(1:2).eq.'-n') then
      call getarg(iarg+1,arg)
      read(arg,*) nQSOProgress 
      iarg=iarg+2
   endif
   nfa=0
   nfb=4200
   ncontest=4
   ndecodes=0
   nfqso=1500
   mycall="K9AN"
   hiscall="K1JT"

   do ifile=iarg,nargs
      call getarg(ifile,infile)
      j2=index(infile,'.wav')
      open(10,file=infile,status='old',access='stream')
      read(10,end=999) ihdr,iwave
      close(10)
      cdatetime=infile(1:13)//'.000'

      call ft4_decode(cdatetime,0.0,nfa,nfb,nQSOProgress,ncontest,nfqso,iwave, &
          ndecodes,mycall,hiscall,nrx,line,data_dir)
      
      do idecode=1,ndecodes
         call get_ft4msg(idecode,nrx,line)
         write(*,'(a61)') line
      enddo
   enddo !files

   write(*,1120)
1120 format("<DecodeFinished>")

999 end program ft4d