Set max_drift = 10 for Q65-60A, 40 for Q65-15C, otherwise 0. Temporary ??

This commit is contained in:
Joe Taylor 2021-05-30 09:25:42 -04:00
parent ccf5e0daa7
commit 42ab7fcd20
1 changed files with 8 additions and 3 deletions

View File

@ -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