From f7464c8df61b0ee53a591a0673adb9e99819a4cc Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Thu, 3 Jun 2021 13:34:22 -0400
Subject: [PATCH] Protect against bounds error when ntol range goes below 0.

---
 lib/qra/q65/q65.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/lib/qra/q65/q65.f90 b/lib/qra/q65/q65.f90
index 0552bc8da..c3795d125 100644
--- a/lib/qra/q65/q65.f90
+++ b/lib/qra/q65/q65.f90
@@ -439,8 +439,8 @@ subroutine q65_ccf_22(s1,iz,jz,nfqso,ntol,ndepth,ntrperiod,iavg,ipk,jpk,  &
   ib=min(nfb,4900)/df
   if(nqd.ne.1 .or. iavg.ne.0) max_drift=0
   if(max_drift.ne.0) then
-     ia=nint((nfqso-ntol)/df)
-     ib=nint((nfqso+ntol)/df)
+     ia=max(100,nint((nfqso-ntol)/df))
+     ib=min(4900,nint((nfqso+ntol)/df))
   endif
 
   do i=ia,ib