diff --git a/lib/pctile.f90 b/lib/pctile.f90
index 3bf2ffcf8..c0a6cdb72 100644
--- a/lib/pctile.f90
+++ b/lib/pctile.f90
@@ -1,16 +1,11 @@
 subroutine pctile(x,npts,npct,xpct)
 
-  parameter (NMAX=256*1024)
-  real*4 x(npts)
-  real*4 tmp(NMAX)
+  real x(npts)
+  real,allocatable :: tmp(:)
 
-  if(npts.le.0) then
-     xpct=1.0
-     go to 900
-  endif
-  if(npts.gt.NMAX) stop
+  allocate(tmp(npts))
 
-  tmp(1:npts)=x
+  tmp=x
   call shell(npts,tmp)
   j=nint(npts*0.01*npct)
   if(j.lt.1) j=1