From 56dfb99a5d2e0689a25c62ab47ac92bad4502af3 Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Sat, 18 Feb 2023 08:46:13 -0500
Subject: [PATCH] Correct the logic that determines when the Q65 red sync curve
is plotted.
---
lib/qra/q65/q65.f90 | 12 ++++--------
1 file changed, 4 insertions(+), 8 deletions(-)
diff --git a/lib/qra/q65/q65.f90 b/lib/qra/q65/q65.f90
index fab3516c0..12dc4f08e 100644
--- a/lib/qra/q65/q65.f90
+++ b/lib/qra/q65/q65.f90
@@ -125,10 +125,7 @@ subroutine q65_dec0(iavg,iwave,ntrperiod,nfqso,ntol,lclearave, &
lclearave=.false.
endif
ccf1=0.
- if(iavg.eq.0) then
- ccf2=0.
- ccf2_avg=0.
- endif
+ if(iavg.eq.0) ccf2=0.
dtstep=nsps/(NSTEP*12000.0) !Step size in seconds
lag1=-1.0/dtstep
lag2=1.0/dtstep + 0.9999
@@ -694,16 +691,15 @@ subroutine q65_write_red(iz,xdt,ccf2_avg,ccf2)
call q65_sync_curve(ccf2_avg,1,iz,rms1)
call q65_sync_curve(ccf2,1,iz,rms2)
- rewind 17
- write(17,1000) xdt,minval(ccf2_avg),maxval(ccf2_avg)
i1=max(1,nint(nfa/df))
i2=min(iz,int(nfb/df))
y0=minval(ccf2(i1:i2))
y0_avg=minval(ccf2_avg(i1:i2))
g=0.4
g_avg=0.
- if(maxval(ccf2_avg(i1:i2)).ne.y0_avg) g_avg=g
-
+ if(navg(iseq).ge.2) g_avg=g
+ rewind 17
+ write(17,1000) xdt,g_avg*minval(ccf2_avg),g_avg*maxval(ccf2_avg)
do i=i1,i2
freq=i*df
y1=g_avg*(ccf2_avg(i)-y0_avg)