another few tweaks, a few more decodes. :-)

This commit is contained in:
Joe Taylor 2022-12-21 16:02:39 -05:00
parent e398b6a33d
commit a085b5e390

View File

@ -107,20 +107,21 @@ subroutine getcand2(ss,savg0,nts_q65,cand,ncand)
savg(ja:jb)=savg(ja:jb)/(1.015*base) savg(ja:jb)=savg(ja:jb)/(1.015*base)
savg0(ja:jb)=savg0(ja:jb)/(1.015*base) savg0(ja:jb)=savg0(ja:jb)/(1.015*base)
enddo enddo
df=96000.0/NFFT df=96000.0/NFFT
bw=65*nts_q65*1.666666667 bw=65*nts_q65*1.666666667
nbw=bw/df + 1 nbw=bw/df + 1
nb0=2*nts_q65
smin=1.4 smin=1.4
nguard=5 nguard=5
j=0 j=0
sync(1:NFFT)%ccfmax=0. sync(1:NFFT)%ccfmax=0.
do i=1,NFFT-2*nbw do i=1,NFFT-nbw-nguard
if(savg(i).lt.smin) cycle if(savg(i).lt.smin) cycle
spk=maxval(savg(i:i+nbw)) spk=maxval(savg(i:i+nb0))
ipk1=maxloc(savg(i:i+nbw)) ipk1=maxloc(savg(i:i+nb0))
i0=ipk1(1) + i - 1 i0=ipk1(1) + i - 1
fpk=0.001*i0*df fpk=0.001*i0*df
! Check to see if sync tone is present. ! Check to see if sync tone is present.
@ -130,10 +131,10 @@ subroutine getcand2(ss,savg0,nts_q65,cand,ncand)
! write(73,3073) j,fpk+32.0-2.270,snr_sync,xdt ! write(73,3073) j,fpk+32.0-2.270,snr_sync,xdt
!3073 format(i3,3f10.3) !3073 format(i3,3f10.3)
cand(j)%f=fpk cand(j)%f=fpk
cand(j)%xdt=2.8 cand(j)%xdt=xdt
cand(j)%snr=spk cand(j)%snr=snr_sync
cand(j)%iflip=0 cand(j)%iflip=0
sync(i0)%ccfmax=spk sync(i0)%ccfmax=snr_sync
ia=min(i,i0-nguard) ia=min(i,i0-nguard)
ib=i0+nbw+nguard ib=i0+nbw+nguard
savg(ia:ib)=0. savg(ia:ib)=0.