From d0db3f74d6f61be626dcbe4b85e89d3bb107b8d3 Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Mon, 2 Apr 2018 18:57:08 +0000 Subject: [PATCH] Further improvements for decoding drifting JT65 signals. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@8601 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 --- lib/afc65b.f90 | 39 +++++++++++++++++++++++++-------------- 1 file changed, 25 insertions(+), 14 deletions(-) diff --git a/lib/afc65b.f90 b/lib/afc65b.f90 index ecdb5ae20..9599b6893 100644 --- a/lib/afc65b.f90 +++ b/lib/afc65b.f90 @@ -6,23 +6,34 @@ subroutine afc65b(cx,npts,fsample,nflip,mode65,a,ccfbest,dtbest) real a(5),deltaa(5) a=0. + i2=8*mode65 + i1=-i2 + j2=8*mode65 + j1=-j2 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) + istep=2*mode65 + do iter=1,2 + do i=i1,i2,istep + a(1)=i + do j=j1,j2,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) istep,i1,i2,j1,j2,i,j,ccf,ccfmax,dtmax,a1,a2 +!3081 format(7i4,5f8.2) + enddo enddo + i1=a1-istep + i2=a1+istep + j1=a2-istep + j2=a2+istep + istep=1 enddo - + ! a(1)=0. ! a(2)=0. a(1)=a1