From 6395087e6f3989808c2e78c0b14a49020b79c5e6 Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Wed, 10 Feb 2021 11:14:46 -0500
Subject: [PATCH] Tweak the averaging scheme: now straight avg up to N=4, then
 expenential with ntc=4.

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

diff --git a/lib/qra/q65/q65.f90 b/lib/qra/q65/q65.f90
index 854b828a2..072e7d0a6 100644
--- a/lib/qra/q65/q65.f90
+++ b/lib/qra/q65/q65.f90
@@ -234,11 +234,11 @@ subroutine q65_symspec(iwave,nmax,iz,jz,s1)
      enddo
   enddo
   if(lnewdat) then
-     ntc=3
+     navg(iseq)=navg(iseq) + 1
+     ntc=min(navg(iseq),4)               !Averaging time constant in sequences
      u=1.0/ntc
      s1a(:,:,iseq)=u*s1 + (1.0-u)*s1a(:,:,iseq)
-     navg(iseq)=navg(iseq) + 1
-  endif
+   endif
 
   return
 end subroutine q65_symspec