From ed957a03904fd8157e0de3f8cc03c2b749eca99c Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Sat, 8 May 2021 10:11:12 -0400 Subject: [PATCH] Protect against rms=0 in ccf65. --- map65/libm65/ccf65.f90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/map65/libm65/ccf65.f90 b/map65/libm65/ccf65.f90 index 38ba1a959..592fe95ed 100644 --- a/map65/libm65/ccf65.f90 +++ b/map65/libm65/ccf65.f90 @@ -110,7 +110,8 @@ subroutine ccf65(ss,nhsym,ssmax,sync1,ipol1,jpz,dt1,flipk, & if(abs(lag-lagpk).gt.1) sq=sq + (ccf(lag,ipol1)-base)**2 enddo rms=sqrt(sq/49.0) - sync1=ccfbest/rms - 4.0 + sync1=-4.0 + if(rms.gt.0.0) sync1=ccfbest/rms - 4.0 dt1=lagpk*(2048.0/11025.0) - 2.5 ! Find base level for normalizing snr2.