From 742e2ef59f8fab21bebe98d37fac7aabb007f9f8 Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Sat, 5 Dec 2020 11:09:27 -0500
Subject: [PATCH] Fix pctile for larger arrays.

---
 lib/pctile.f90 | 13 ++++---------
 1 file changed, 4 insertions(+), 9 deletions(-)

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