From 77b7e714249d00871bde27ce4b9eae5a8bfcbff6 Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Fri, 31 Jul 2020 11:16:03 -0400
Subject: [PATCH] Enable mode QRA66B in qra66sim.

---
 lib/qra/qra66/qra66sim.f90 | 9 +++++++--
 1 file changed, 7 insertions(+), 2 deletions(-)

diff --git a/lib/qra/qra66/qra66sim.f90 b/lib/qra/qra66/qra66sim.f90
index 6043b3aed..934239189 100644
--- a/lib/qra/qra66/qra66sim.f90
+++ b/lib/qra/qra66/qra66sim.f90
@@ -42,8 +42,12 @@ program qra66sim
   twopi=8.d0*atan(1.d0)
   npts=NMAX                          !Total samples in .wav file
   nsps=1920
-  baud=12000.d0/nsps                 !Keying rate = 6.25 baud
   nsym=85                            !Number of channel symbols
+  if(csubmode.eq.'B') then
+     nsps=960
+     nsym=169
+  endif
+  baud=12000.d0/nsps                 !Keying rate = 6.25 baud
   h=default_header(12000,npts)
   ichk=66                            !Flag sent to genqra64
 
@@ -76,8 +80,9 @@ program qra66sim
      do i=1,npts                         !Add this signal into cdat()
         isym=i/nsps + 1
         if(isym.gt.nsym) exit
+        if(csubmode.eq.'B' .and. isym.gt.84) isym=isym-84
         if(isym.ne.isym0) then
-           freq=f0 + itone(isym)*baud*mode66
+           freq=f0 + itone(isym)*baud
            dphi=twopi*freq*dt
            isym0=isym
         endif