module sfox_mod
  
  parameter (NMAX=15*12000)       !Samples in iwave (180,000)
  integer MM,NQ,NN,KK,ND1,ND2,ND,NFZ,NSPS,NS,NSYNC,NZ,NFFT,NFFT1

contains
  subroutine sfox_init(mm0,nn0,kk0,itu,fspread,delay)

    character*2 itu
    integer isps(30)
    integer iloc(1)
    data isps/1000,1008,1024,1050,1080,1120,1152,1176,1200,1260,   &
              1280,1296,1344,1350,1372,1400,1440,1470,1500,1512,   &
              1536,1568,1600,1620,1680,1728,1920,1944,1960,2000/


    MM=mm0              !Bits per symbol
    NQ=2**MM            !Q, number of MFSK tones
    NN=nn0              !Channel symbols, before puncture
    KK=kk0              !Information symbols, before puncture
    ND1=25              !Data symbols before sync 
    ND2=NN-2-ND1        !Data symbols after sync 
    ND=ND1+ND2          !Total data symbols
    NFZ=3               !First zero

    tsync=2.0
    jsps=nint((12.6-tsync)*12000.0/ND)
    iloc=minloc(abs(isps-jsps))
    NSPS=isps(iloc(1))  !Samples per symbol
    NS=nint(tsync*12000.0/NSPS)
    if(mod(NS,2).eq.1) NS=NS+1
    NSYNC=NS*NSPS       !Samples in sync waveform
    NZ=NSPS*(ND+NS)     !Samples in full Tx waveform
    NFFT=32768          !Length of FFT for sync waveform
    NFFT1=2*NSPS        !Length of FFTs for symbol spectra

    fspread=0.0
    delay=0.0
    if(itu.eq.'LQ') then
       fspread=0.5
       delay=0.5
    else if(itu.eq.'LM') then
       fspread=1.5
       delay=2.0
    else if(itu.eq.'LD') then
       fspread=10.0
       delay=6.0
    else if(itu.eq.'MQ') then
       fspread=0.1
       delay=0.5
    else if(itu.eq.'MM') then
       fspread=0.5
       delay=1.0
    else if(itu.eq.'MD') then
       fspread=1.0
       delay=2.0
    else if(itu.eq.'HQ') then
       fspread=0.5
       delay=1.0
    else if(itu.eq.'HM') then
       fspread=10.0
       delay=3.0
    else if(itu.eq.'HD') then
       fspread=30.0
       delay=7.0
    endif

    return
  end subroutine sfox_init

end module sfox_mod