Make pctile in libm65 the same as pctile in lib.

map65/libm65/decode0.f90

@@ -16,6 +16,7 @@ subroutine decode0(dd,ss,savg,nstandalone)
   data neme0/-99/,mcall3b/1/
   save
 
+!  write(60) ss,savg
   call timer('decode0 ',0)
 
   if(newdat.ne.0) then

map65/libm65/pctile.f90

@@ -1,22 +1,16 @@
 subroutine pctile(x,npts,npct,xpct)
 
-  parameter (NMAX=100000)
-  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
   if(j.gt.npts) j=npts
   xpct=tmp(j)
 
-900 continue
   return
 end subroutine pctile