comicsansg
/
ft8d
mirror of https://github.com/pavel-demin/ft8d.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

fix s calculations in sync8.f90

This commit is contained in:
Pavel Demin 2018-04-02 10:37:24 +02:00
parent 16fd081814
commit afd405291e
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
sync8.f90
View File

@ -3,7 +3,7 @@ subroutine sync8(dd,nfa,nfb,syncmin,nfqso,s,candidate,ncand,sbase)
include 'ft8_params.f90' include 'ft8_params.f90'
! Search over +/- 2.5s relative to 0.5s TX start time. ! Search over +/- 2.5s relative to 0.5s TX start time.
parameter (JZ=62) parameter (JZ=62)
complex cx(0:NFFT1) complex cx(NFFT1)
real s(NFFT1,NHSYM) real s(NFFT1,NHSYM)
real savg(NFFT1) real savg(NFFT1)
real sbase(NFFT1) real sbase(NFFT1)
@ -31,10 +31,11 @@ subroutine sync8(dd,nfa,nfb,syncmin,nfqso,s,candidate,ncand,sbase)
x(1:NSPS)=fac*dd(ia:ib) x(1:NSPS)=fac*dd(ia:ib)
x(NSPS+1:)=0. x(NSPS+1:)=0.
call four2a(x,NFFT1,1,-1,1) !c2c FFT call four2a(x,NFFT1,1,-1,1) !c2c FFT
do i=1,NH1 do i=2,NH1
s(i,j)=real(cx(i+NH1))**2 + aimag(cx(i+NH1))**2 s(i-1,j)=real(cx(i+NH1))**2 + aimag(cx(i+NH1))**2
s(i+NH1,j)=real(cx(i))**2 + aimag(cx(i))**2 s(i-1+NH1,j)=real(cx(i))**2 + aimag(cx(i))**2
enddo enddo
s(NH1,j)=real(cx(1))**2 + aimag(cx(1))**2
savg=savg + s(1:NFFT1,j) !Average spectrum savg=savg + s(1:NFFT1,j) !Average spectrum
enddo enddo
call baseline(savg,nfa,nfb,sbase) call baseline(savg,nfa,nfb,sbase)