From 73db6eb28fa6d8208e9ed8cc7de2cc5cfc61007b Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Mon, 23 Jan 2023 17:25:31 -0500
Subject: [PATCH] Compute symbol spectra at 1/4 symbol steps, rather than 1/8
 symbol.

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

diff --git a/lib/qra/q65/q65.f90 b/lib/qra/q65/q65.f90
index 6bcd6ae97..a54807e9c 100644
--- a/lib/qra/q65/q65.f90
+++ b/lib/qra/q65/q65.f90
@@ -294,7 +294,7 @@ subroutine q65_symspec(iwave,nmax,iz,jz,s1)
   allocate(c0(0:nsps-1))
   nfft=nsps
   fac=1/32767.0
-  do j=1,jz                              !Compute symbol spectra at step size
+  do j=1,jz,2                     !Compute symbol spectra at 2*step size
      i1=(j-1)*istep
      i2=i1+nsps-1
      k=-1
@@ -314,6 +314,8 @@ subroutine q65_symspec(iwave,nmax,iz,jz,s1)
      do i=1,nsmo
         call smo121(s1(1:iz,j),iz)
      enddo
+! Interpolate to fill in the skipped-over spectra.
+     if(j.ge.3) s1(1:iz,j-1)=0.5*(s1(1:iz,j-2)+s1(1:iz,j))
   enddo
   if(lnewdat) then
      navg(iseq)=navg(iseq) + 1