From f23129c33cd60b34a57eafbaa866bb7e090cfba6 Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Sat, 19 Dec 2020 11:25:18 -0500 Subject: [PATCH] Fix a problem with Q65 FTol range extending down below 0 Hz. --- lib/q65_sync.f90 | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/lib/q65_sync.f90 b/lib/q65_sync.f90 index f65f79838..c6db664a3 100644 --- a/lib/q65_sync.f90 +++ b/lib/q65_sync.f90 @@ -208,13 +208,15 @@ subroutine q65_sync(nutc,iwave,nmax,mode_q65,codewords,ncw,nsps,nfqso,ntol, & irc=-2 dat4=0 ia=ntol/df - if(i0-ia.lt.1 .or. i0+ia.gt.iz) go to 900 do lag=lag1,lag2 do k=1,85 n=NSTEP*(k-1) + 1 j=n+lag+j0 if(j.ge.1 .and. j.le.jz) then - ccf(-ia:ia,lag)=ccf(-ia:ia,lag) + sync(k)*s1(i0-ia:i0+ia,j) + do i=-ia,ia + if(i0+i.lt.1 .or. i0+i.gt.iz) cycle + ccf(i,lag)=ccf(i,lag) + sync(k)*s1(i0+i,j) + enddo endif enddo enddo