Correct a problem with estimating rms of ccf for Q65-300x submodes.

2020-12-19 10:42:27 -05:00 · 2020-12-19 10:42:27 -05:00 · 7096a548ec
parent 7005ac3c50
commit 7096a548ec
1 changed files with 2 additions and 1 deletions
--- a/lib/q65_sync.f90
+++ b/lib/q65_sync.f90
@ -226,9 +226,10 @@ subroutine q65_sync(nutc,iwave,nmax,mode_q65,codewords,ncw,nsps,nfqso,ntol, &
  xdt=jpk*dtstep
  sq=0.
  nsq=0
+  jd=(lag2-lag1)/4
  do i=-ia,ia
     do j=lag1,lag2
-        if(abs(j-jpk).gt.8 .and. abs(i-ipk).gt.ia/2) then
+        if(abs(j-jpk).gt.jd .and. abs(i-ipk).gt.ia/2) then
           sq=sq + ccf(i,j)**2
           nsq=nsq+1
        endif