From afd405291e56f32accda059e847fb5bd490e0a2c Mon Sep 17 00:00:00 2001 From: Pavel Demin Date: Mon, 2 Apr 2018 10:37:24 +0200 Subject: [PATCH] fix s calculations in sync8.f90 --- sync8.f90 | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/sync8.f90 b/sync8.f90 index 49cf92d..3ad6ad8 100644 --- a/sync8.f90 +++ b/sync8.f90 @@ -3,7 +3,7 @@ subroutine sync8(dd,nfa,nfb,syncmin,nfqso,s,candidate,ncand,sbase) include 'ft8_params.f90' ! Search over +/- 2.5s relative to 0.5s TX start time. parameter (JZ=62) - complex cx(0:NFFT1) + complex cx(NFFT1) real s(NFFT1,NHSYM) real savg(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(NSPS+1:)=0. call four2a(x,NFFT1,1,-1,1) !c2c FFT - do i=1,NH1 - s(i,j)=real(cx(i+NH1))**2 + aimag(cx(i+NH1))**2 - s(i+NH1,j)=real(cx(i))**2 + aimag(cx(i))**2 + do i=2,NH1 + s(i-1,j)=real(cx(i+NH1))**2 + aimag(cx(i+NH1))**2 + s(i-1+NH1,j)=real(cx(i))**2 + aimag(cx(i))**2 enddo + s(NH1,j)=real(cx(1))**2 + aimag(cx(1))**2 savg=savg + s(1:NFFT1,j) !Average spectrum enddo call baseline(savg,nfa,nfb,sbase)