From e7e1b3c3cc4b8252dcc2e72cfb08fe7a55662b54 Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Mon, 2 Apr 2018 18:07:56 +0000 Subject: [PATCH] Further improvements for decoding drifting JT65 signals. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@8600 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 --- lib/afc65b.f90 | 24 ++++++++++++++++++++++-- 1 file changed, 22 insertions(+), 2 deletions(-) diff --git a/lib/afc65b.f90 b/lib/afc65b.f90 index 3078fec2f..ecdb5ae20 100644 --- a/lib/afc65b.f90 +++ b/lib/afc65b.f90 @@ -5,8 +5,28 @@ subroutine afc65b(cx,npts,fsample,nflip,mode65,a,ccfbest,dtbest) complex cx(npts) real a(5),deltaa(5) - a(1)=0. - a(2)=0. + a=0. + ccfmax=0. + istep=mode65 + do i=-30,30,istep + a(1)=i + do j=-30,30,istep + a(2)=j + chisq=fchisq65(cx,npts,fsample,nflip,a,ccf,dtmax) + if(ccf.gt.ccfmax) then + a1=a(1) + a2=a(2) + ccfmax=ccf + endif +! write(81,3081) i,j,ccf,dtmax +!3081 format(2i5,2f10.3) + enddo + enddo + +! a(1)=0. +! a(2)=0. + a(1)=a1 + a(2)=a2 a(3)=0. a(4)=0. deltaa(1)=2.0*mode65