From 91c0e43ef645653d7acfbd4593c0f754f60f679e Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Fri, 8 Dec 2023 16:48:17 -0500
Subject: [PATCH] Don't re-define value of mode_q65; use mode_q65_tmp instead!

---
 qmap/libqmap/qmapa.f90 | 15 +++------------
 1 file changed, 3 insertions(+), 12 deletions(-)

diff --git a/qmap/libqmap/qmapa.f90 b/qmap/libqmap/qmapa.f90
index 11715747a..b1caa9413 100644
--- a/qmap/libqmap/qmapa.f90
+++ b/qmap/libqmap/qmapa.f90
@@ -61,24 +61,15 @@ subroutine qmapa(dd,ss,savg,newdat,nutc,fcenter,ntol,nfa,nfb,         &
      f0=cand(icand)%f
      ntrperiod=cand(icand)%ntrperiod
      iseq=cand(icand)%iseq
-!###
-!     if(icand.gt.2) exit
-!     f0=-31.847 + 117.602
-!     ntrperiod=30
-!     iseq=icand-1
-!###
-     mode_q65=3                           !###
-     if(ntrperiod.eq.30) mode_q65=2       !###
+     mode_q65_tmp=mode_q65
+     if(ntrperiod.eq.30) mode_q65_tmp=max(1,mode_q65-1)
      freq=f0+nkhz_center-48.0-1.27046
-!     write(*,5001) icand,ntrperiod,iseq,mode_q65,f0,f0+nkhz_center-48.0,  &
-!          cand(icand)%xdt,cand(icand)%snr
-!5001 format('a',4i5,2f10.3,2f8.1)
      ikhz=nint(freq)
      idec=-1
      call timer('q65b    ',0)
      call q65b(nutc,nqd,fcenter,nfcal,nfsample,ikhz,mousedf,ntol,       &
           ntrperiod,iseq,                                               &
-          mycall,hiscall,hisgrid,mode_q65,f0,fqso,nkhz_center,newdat,   &
+          mycall,hiscall,hisgrid,mode_q65_tmp,f0,fqso,nkhz_center,newdat,   &
           nagain2,max_drift,offset,ndepth,datetime,ndop00,idec)
      call timer('q65b    ',1)
      tsec=sec_midn() - tsec0