WSJT-X/lib/pctile.f90
2021-01-07 10:18:38 -06:00

17 lines
247 B
Fortran

subroutine pctile(x,npts,npct,xpct)
real x(npts)
real,allocatable :: tmp(:)
allocate(tmp(npts))
tmp=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)
return
end subroutine pctile