From 421524d88cc682efa33acee2bb3dbc0c45957934 Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Thu, 7 Jul 2016 19:27:03 +0000 Subject: [PATCH] Fix the broken S/N calculation for QRA64. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6888 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 --- lib/qra64a.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/lib/qra64a.f90 b/lib/qra64a.f90 index 6812deaf8..7b31df907 100644 --- a/lib/qra64a.f90 +++ b/lib/qra64a.f90 @@ -45,16 +45,16 @@ subroutine qra64a(dd,nf1,nf2,nfqso,ntol,mycall_12,sync,nsnr,dtx,nfreq, & ss(1:NZ,j)=s savg=savg+s enddo - savg=savg/nsteps - call pctile(savg,NZ,45,base) - savg=savg/base - 1.0 - ss=ss/base fa=max(nf1,nfqso-ntol) fb=min(nf2,nfqso+ntol) ia=nint(fa/df) ib=nint(fb/df) + call pctile(savg(ia),ib-ia+1,45,base) + savg=savg/base - 1.0 + ss=ss/base + red0=0. fac=1.0/sqrt(21.0) sync=0.