From 588e862346d38f04d001800328e9de8595083ea2 Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Fri, 25 Jun 2021 16:29:42 -0400
Subject: [PATCH] Improves setting of limits for flattening s1(), when nfqso is
 too close to 0.

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

diff --git a/lib/qra/q65/q65.f90 b/lib/qra/q65/q65.f90
index 13f254669..914932705 100644
--- a/lib/qra/q65/q65.f90
+++ b/lib/qra/q65/q65.f90
@@ -139,16 +139,17 @@ subroutine q65_dec0(iavg,nutc,iwave,ntrperiod,nfqso,ntol,ndepth,lclearave,  &
   endif
 
   i0=nint(nfqso/df)                             !Target QSO frequency
-  if(i0-64.lt.1 .or. i0-65+LL.gt.iz) go to 900  !Frequency out of range
-  call pctile(s1(i0-64:i0-65+LL,1:jz),LL*jz,45,base)
+  ii1=max(1,i0-64)
+  ii2=i0-65+LL
+  call pctile(s1(ii1:ii2,1:jz),ii2-ii1+1*jz,45,base)
   s1=s1/base
   s1raw=s1
 
 ! Apply fast AGC to the symbol spectra
   s1max=20.0                                  !Empirical choice
   do j=1,jz                                   !### Maybe wrong way? ###
-     smax=maxval(s1(i0-64:i0-65+LL,j))
-     if(smax.gt.s1max) s1(i0-64:i0-65+LL,j)=s1(i0-64:i0-65+LL,j)*s1max/smax
+     smax=maxval(s1(ii1:ii2,j))
+     if(smax.gt.s1max) s1(ii1:ii2,j)=s1(ii1:ii2,j)*s1max/smax
   enddo
 
   dat4=0