From ec2d5b6233516b2b7cb5797ce35a9dc6c609d757 Mon Sep 17 00:00:00 2001
From: Steve Franke <s.j.franke@icloud.com>
Date: Mon, 25 Feb 2019 12:57:34 -0600
Subject: [PATCH] Fix the sign of one of the Nuttal window coefficients. Change
 BT to 2.0.

---
 lib/ft8/gen_ft8wave.f90 | 2 +-
 lib/symspec.f90         | 6 +++---
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/lib/ft8/gen_ft8wave.f90 b/lib/ft8/gen_ft8wave.f90
index 51c0033df..f9665ca56 100644
--- a/lib/ft8/gen_ft8wave.f90
+++ b/lib/ft8/gen_ft8wave.f90
@@ -17,7 +17,7 @@ subroutine gen_ft8wave(itone,nsym,nsps,fsample,f0,cwave,wave,icmplx,nwave)
      twopi=8.0*atan(1.0)
      dt=1.0/fsample
      hmod=1.0
-     bt=1.5
+     bt=2.0
 ! Compute the frequency-smoothing pulse
      do i=1,3*nsps
         tt=(i-1.5*nsps)/real(nsps)
diff --git a/lib/symspec.f90 b/lib/symspec.f90
index be2bd296e..20bf4dbf1 100644
--- a/lib/symspec.f90
+++ b/lib/symspec.f90
@@ -58,11 +58,11 @@ subroutine symspec(shared_data,k,ntrperiod,nsps,ingain,nminw,pxdb,s,   &
 ! discrete Hilbert Transform, by Albert Nuttall"
 !
      a0=0.3635819
-     a1=0.4891775;
+     a1=-0.4891775;
      a2=0.1365995;
-     a3=0.0106411;
+     a3=-0.0106411;
      do i=1,nfft3
-         w3(i)=a0-a1*cos(2*pi*(i-1)/(nfft3))+ &
+         w3(i)=a0+a1*cos(2*pi*(i-1)/(nfft3))+ &
                   a2*cos(4*pi*(i-1)/(nfft3))+ &
                   a3*cos(6*pi*(i-1)/(nfft3))
      enddo