From a6779ecda3a461eaac75a250bbe2e5ce9e3e1beb Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Mon, 8 Feb 2021 11:10:02 -0500
Subject: [PATCH] Use exponential message averaging (ntc=3) rather than simple
 sums in Q65.

---
 UnitTests.txt       | 1 +
 lib/qra/q65/q65.f90 | 4 +++-
 2 files changed, 4 insertions(+), 1 deletion(-)

diff --git a/UnitTests.txt b/UnitTests.txt
index f199243bd..262d0bd39 100644
--- a/UnitTests.txt
+++ b/UnitTests.txt
@@ -9,6 +9,7 @@ ef4787:    3 10         6 14        30 33
 ada5a6:    3  6         6 10        29 36
 10f574:    2  7  97.6   6 10 95.4   21 33 80.7
 b8ea4c:    2  7 100.5   6 10 96.9   21 33 81.0
+           1  6  95.2   6 11 90.1   21 33 67.3
 
 Mode: Q65-30A
 Data: 30A_N0AN_6m_Ionoscatter (69 files, 6m ionoscatter)
diff --git a/lib/qra/q65/q65.f90 b/lib/qra/q65/q65.f90
index 195b07d3e..854b828a2 100644
--- a/lib/qra/q65/q65.f90
+++ b/lib/qra/q65/q65.f90
@@ -234,7 +234,9 @@ subroutine q65_symspec(iwave,nmax,iz,jz,s1)
      enddo
   enddo
   if(lnewdat) then
-     s1a(:,:,iseq)=s1a(:,:,iseq) + s1
+     ntc=3
+     u=1.0/ntc
+     s1a(:,:,iseq)=u*s1 + (1.0-u)*s1a(:,:,iseq)
      navg(iseq)=navg(iseq) + 1
   endif