WSJT-X/lib/pctile.f90

17 lines
247 B
Fortran
Raw Normal View History

subroutine pctile(x,npts,npct,xpct)
2020-12-05 11:09:27 -05:00
real x(npts)
real,allocatable :: tmp(:)
2020-12-05 11:09:27 -05:00
allocate(tmp(npts))
2020-12-05 11:09:27 -05:00
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