WSJT-X/lib/averms.f90

subroutine averms(x,n,nskip,ave,rms)
  real x(n)
  integer ipk(1)

  ns=0
  s=0.
  sq=0.
  ipk=maxloc(x)
  do i=1,n
     if((nskip.lt.0) .or. (abs(i-ipk(1)).gt.nskip)) then
        s=s + x(i)
        sq=sq + x(i)**2
        ns=ns+1
     endif
  enddo
  ave=s/ns
  rms=sqrt(sq/ns - ave*ave)
 
  return
end subroutine averms
MOving toward use of sync64(). Not yet there! git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6939 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2016-07-22 15:52:20 -04:00			`subroutine averms(x,n,nskip,ave,rms)`
			`real x(n)`
			`integer ipk(1)`

			`ns=0`
			`s=0.`
			`sq=0.`
			`ipk=maxloc(x)`
			`do i=1,n`
Tweak subroutine averms(). Add test program test_snr to exercise new algorithm for echo snr. 2022-04-22 11:02:47 -04:00			`if((nskip.lt.0) .or. (abs(i-ipk(1)).gt.nskip)) then`
MOving toward use of sync64(). Not yet there! git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6939 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2016-07-22 15:52:20 -04:00			`s=s + x(i)`
			`sq=sq + x(i)**2`
			`ns=ns+1`
			`endif`
			`enddo`
			`ave=s/ns`
			`rms=sqrt(sq/ns - ave*ave)`

			`return`
			`end subroutine averms`