WSJT-X/libm65/pctile.f90

subroutine pctile(x,npts,npct,xpct)

  parameter (NMAX=100000)
  real*4 x(npts)
  real*4 tmp(NMAX)

  if(npts.le.0) then
     xpct=1.0
     go to 900
  endif
  if(npts.gt.NMAX) stop

  tmp(1:npts)=x
  call shell(npts,tmp)
  j=nint(npts*0.01*npct)
  if(j.lt.1) j=1
  if(j.gt.npts) j=npts
  xpct=tmp(j)

900 continue
  return
end subroutine pctile
Bring the sort routines up-to-date with WSJT-X. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@7497 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2017-01-14 11:52:03 -05:00			`subroutine pctile(x,npts,npct,xpct)`

			`parameter (NMAX=100000)`
			`real*4 x(npts)`
			`real*4 tmp(NMAX)`

			`if(npts.le.0) then`
			`xpct=1.0`
			`go to 900`
			`endif`
			`if(npts.gt.NMAX) stop`

			`tmp(1:npts)=x`
			`call shell(npts,tmp)`
			`j=nint(npts0.01npct)`
			`if(j.lt.1) j=1`
			`if(j.gt.npts) j=npts`
			`xpct=tmp(j)`

			`900 continue`
			`return`
			`end subroutine pctile`