2021-06-30 10:40:38 -04:00
|
|
|
subroutine spec64(c0,npts,nsps,mode_q65,jpk,s3,LL,NN)
|
2016-07-26 13:51:13 -04:00
|
|
|
|
2020-10-22 19:07:21 -04:00
|
|
|
parameter (MAXFFT=20736)
|
2021-06-30 10:40:38 -04:00
|
|
|
complex c0(0:npts-1) !Complex spectrum of dd()
|
2020-10-08 16:48:11 -04:00
|
|
|
complex cs(0:MAXFFT-1) !Complex symbol spectrum
|
2016-11-07 11:47:40 -05:00
|
|
|
real s3(LL,NN) !Synchronized symbol spectra
|
2016-12-15 13:42:33 -05:00
|
|
|
real xbase0(LL),xbase(LL)
|
2020-10-22 14:05:20 -04:00
|
|
|
! integer ipk1(1)
|
2020-10-08 16:48:11 -04:00
|
|
|
integer isync(22) !Indices of sync symbols
|
|
|
|
|
data isync/1,9,12,13,15,22,23,26,27,33,35,38,46,50,55,60,62,66,69,74,76,85/
|
2016-07-26 13:51:13 -04:00
|
|
|
|
2016-12-15 13:42:33 -05:00
|
|
|
nfft=nsps
|
2021-01-27 11:12:57 -05:00
|
|
|
j=0
|
|
|
|
|
n=1
|
|
|
|
|
do k=1,84
|
|
|
|
|
if(k.eq.isync(n)) then
|
|
|
|
|
n=n+1
|
|
|
|
|
cycle
|
|
|
|
|
endif
|
|
|
|
|
j=j+1
|
|
|
|
|
ja=(k-1)*nsps + jpk
|
|
|
|
|
jb=ja+nsps-1
|
2021-06-30 11:28:03 -04:00
|
|
|
if(ja.lt.0) ja=0
|
|
|
|
|
if(jb.gt.npts-1) jb=npts-1
|
2021-06-30 14:49:00 -04:00
|
|
|
nz=jb-ja
|
|
|
|
|
cs(0:nz)=c0(ja:jb)
|
|
|
|
|
if(nz.lt.nfft-1) cs(nz+1:)=0.
|
2021-01-27 11:12:57 -05:00
|
|
|
call four2a(cs,nsps,1,-1,1) !c2c FFT to frequency
|
|
|
|
|
do ii=1,LL
|
|
|
|
|
i=ii-65+mode_q65 !mode_q65 = 1 2 4 8 16 for Q65 A B C D E
|
|
|
|
|
if(i.lt.0) i=i+nsps
|
|
|
|
|
s3(ii,j)=real(cs(i))**2 + aimag(cs(i))**2
|
2020-10-08 16:48:11 -04:00
|
|
|
enddo
|
2021-01-27 11:12:57 -05:00
|
|
|
enddo
|
2016-07-26 13:51:13 -04:00
|
|
|
|
2016-12-15 13:42:33 -05:00
|
|
|
df=6000.0/nfft
|
|
|
|
|
do i=1,LL
|
|
|
|
|
call pctile(s3(i,1:NN),NN,45,xbase0(i)) !Get baseline for passband shape
|
|
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
nh=25
|
|
|
|
|
xbase(1:nh-1)=sum(xbase0(1:nh-1))/(nh-1.0)
|
|
|
|
|
xbase(LL-nh+1:LL)=sum(xbase0(LL-nh+1:LL))/(nh-1.0)
|
|
|
|
|
do i=nh,LL-nh
|
|
|
|
|
xbase(i)=sum(xbase0(i-nh+1:i+nh))/(2*nh+1) !Smoothed passband shape
|
|
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
do i=1,LL
|
|
|
|
|
s3(i,1:NN)=s3(i,1:NN)/(xbase(i)+0.001) !Apply frequency equalization
|
|
|
|
|
enddo
|
|
|
|
|
|
2016-07-26 13:51:13 -04:00
|
|
|
return
|
|
|
|
|
end subroutine spec64
|
