diff --git a/lib/superfox/sfox_mod.f90 b/lib/superfox/sfox_mod.f90 new file mode 100644 index 000000000..9e5273238 --- /dev/null +++ b/lib/superfox/sfox_mod.f90 @@ -0,0 +1,71 @@ +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(23) + 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/ + + 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 (125) + 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 (24,576) + NZ=NSPS*(ND+NS) !Samples in full Tx waveform (151,552) + 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