WSJT-X/lib/spec64.f90
Joe Taylor e7bd8383fb QRA64 decoder is now basically functional with qra64_decode_fastfading(). Currently only submode A. More to come...
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@7299 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
2016-11-07 16:47:40 +00:00

27 lines
821 B
Fortran

subroutine spec64(c0,npts2,mode64,jpk,s3,LL,NN)
parameter (NSPS=2304) !Samples per symbol at 4000 Hz
complex c0(0:360000) !Complex spectrum of dd()
complex cs(0:NSPS-1) !Complex symbol spectrum
real s3(LL,NN) !Synchronized symbol spectra
nfft6=nsps
fsample=4000.0
df=fsample/nfft6
do j=1,63
jj=j+7 !Skip first Costas array
if(j.ge.32) jj=j+14 !Skip middle Costas array
ja=jpk + (jj-1)*nfft6
jb=ja+nfft6-1
cs(0:nfft6-1)=1.3e-8*c0(ja:jb)
call four2a(cs,nfft6,1,-1,1)
do ii=1,LL
i=ii-65
if(i.lt.0) i=i+nfft6
s3(ii,j)=real(cs(i))**2 + aimag(cs(i))**2
enddo
enddo
return
end subroutine spec64