From d815eb759204d9012c23a1abf5aa42ea7a054e5b Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Tue, 19 Jan 2021 09:24:30 -0500
Subject: [PATCH] Replace the 'nsmo' smoothing: it is indeed useful, sometimes.
 (Can we have it both ways?)

---
 lib/qra/q65/q65.f90 | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/lib/qra/q65/q65.f90 b/lib/qra/q65/q65.f90
index 4c4773301..68d550c1e 100644
--- a/lib/qra/q65/q65.f90
+++ b/lib/qra/q65/q65.f90
@@ -11,7 +11,7 @@ module q65
                                      38,46,50,55,60,62,66,69,74,76,85/)
   integer codewords(63,206)
   integer navg,ibwa,ibwb,ncw,nsps,mode_q65,nfa,nfb
-  integer istep,lag1,lag2,npasses,nused
+  integer istep,nsmo,lag1,lag2,npasses,nused
   integer i0,j0
   real,allocatable,save :: s1a(:,:)      !Cumulative symbol spectra
   real sync(85)                          !sync vector
@@ -81,6 +81,8 @@ subroutine q65_dec0(iavg,nutc,iwave,ntrperiod,nfqso,ntol,ndepth,lclearave,  &
   if(nsps.ge.6912) jz=(txt+2.0)*12000.0/istep   !For TR 60 s and higher
   ia=ntol/df
   ia2=max(ia,10*mode_q65,nint(100.0/df))
+  nsmo=int(0.7*mode_q65*mode_q65)
+  if(nsmo.lt.1) nsmo=1
   if(first) then                         !Generate the sync vector
      sync=-22.0/63.0                     !Sync tone OFF  
      do k=1,22
@@ -251,6 +253,10 @@ subroutine q65_symspec(iwave,nmax,iz,jz,s1)
      do i=1,iz
         s1(i,j)=real(c0(i))**2 + aimag(c0(i))**2
      enddo
+! For large Doppler spreads, should we smooth the spectra here?
+     do i=1,nsmo
+        call smo121(s1(1:iz,j),iz)
+     enddo
   enddo
   s1a=s1a+s1
   navg=navg+1