From 2ed4ca3aaff41fad8b164d5ef686cfae9cb3c46b Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Thu, 15 Apr 2021 13:51:01 -0400
Subject: [PATCH] Guard against base=0 when normalizing an snr2 estimate in
ccf65().
---
map65/libm65/ccf65.f90 | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/map65/libm65/ccf65.f90 b/map65/libm65/ccf65.f90
index 1bebb8d73..c5d34ec38 100644
--- a/map65/libm65/ccf65.f90
+++ b/map65/libm65/ccf65.f90
@@ -113,8 +113,9 @@ subroutine ccf65(ss,nhsym,ssmax,sync1,ipol1,jpz,dt1,flipk, &
tmp1(i)=ss(ipol2,i)
enddo
call pctile(tmp1,nhsym,40,base)
- snr2=0.398107*ccfbest2/base !### empirical
- syncshort=0.5*ccfbest2/rms - 4.0 !### better normalizer than rms?
+ snr2=0.01
+ if(base.gt.0.0) snr2=0.398107*ccfbest2/base !### empirical
+ syncshort=0.5*ccfbest2/rms - 4.0 !### better normalizer than rms?
dt2=2.5 + lagpk2*(2048.0/11025.0)
return