From 42ab7fcd206c36668179170c7c0bb885b838ccf9 Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Sun, 30 May 2021 09:25:42 -0400 Subject: [PATCH] Set max_drift = 10 for Q65-60A, 40 for Q65-15C, otherwise 0. Temporary ?? --- lib/qra/q65/q65.f90 | 11 ++++++++--- 1 file changed, 8 insertions(+), 3 deletions(-) diff --git a/lib/qra/q65/q65.f90 b/lib/qra/q65/q65.f90 index 2012af83f..9f89fa2e7 100644 --- a/lib/qra/q65/q65.f90 +++ b/lib/qra/q65/q65.f90 @@ -437,11 +437,16 @@ subroutine q65_ccf_22(s1,iz,jz,nfqso,ntol,ndepth,ntrperiod,iavg,ipk,jpk, & ia=max(nfa,100)/df ib=min(nfb,4900)/df - max_drift=0 - if(nqd.eq.1 .and. iavg.eq.0 .and. ntol.le.100 .and. ntrperiod.eq.60) then + max_drift=0 !Drift units: bins/TxT +! Do we need a GUI control to set max_drift ? +! For now, turn on drift compensation only for submodes 15B amd 60A. + if(nqd.eq.1 .and. iavg.eq.0 .and. ntol.le.100 .and. ntrperiod.eq.60 .and. & + mode_q65.eq.1) max_drift=10 !Q65-60A + if(nqd.eq.1 .and. iavg.eq.0 .and. ntol.le.100 .and. ntrperiod.eq.15 .and. & + mode_q65.eq.4) max_drift=40 !Q65-15C + if(max_drift.ne.0) then ia=nint((nfqso-ntol)/df) ib=nint((nfqso+ntol)/df) - max_drift=10 !Drift units: bins/TRperiod (?) endif do i=ia,ib