subroutine pctile(x,tmp,nmax,npct,xpct)
  real x(nmax),tmp(nmax)

  do i=1,nmax
     tmp(i)=x(i)
  enddo
  call sort(nmax,tmp)
  j=nint(nmax*0.01*npct)
  if(j.lt.1) j=1
  xpct=tmp(j)

  return
end subroutine pctile