2012-10-11 14:33:50 -04:00
|
|
|
subroutine pctile(x,npts,npct,xmedian)
|
|
|
|
|
|
|
|
|
|
real x(npts)
|
|
|
|
|
integer hist(0:1000)
|
|
|
|
|
|
2012-10-31 14:33:56 -04:00
|
|
|
if(npts.le.0) then
|
|
|
|
|
xmedian=1.0
|
|
|
|
|
go to 900
|
|
|
|
|
endif
|
|
|
|
|
|
2012-10-11 14:33:50 -04:00
|
|
|
ave=sum(x)/npts
|
2012-10-31 14:33:56 -04:00
|
|
|
s=0.
|
|
|
|
|
ns=0
|
|
|
|
|
do i=1,npts
|
|
|
|
|
if(x(i).lt.3.0*ave) then
|
|
|
|
|
s=s+x(i)
|
|
|
|
|
ns=ns+1
|
|
|
|
|
endif
|
|
|
|
|
enddo
|
|
|
|
|
ave=s/ns
|
|
|
|
|
|
2012-10-11 14:33:50 -04:00
|
|
|
hist=0
|
|
|
|
|
do i=1,npts
|
|
|
|
|
j=nint(100.0*x(i)/ave)
|
|
|
|
|
if(j.lt.0) j=0
|
|
|
|
|
if(j.gt.1000) j=1000
|
|
|
|
|
hist(j)=hist(j)+1
|
|
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
nsum=0
|
|
|
|
|
ntest=nint(npts*float(npct)/100.0)
|
|
|
|
|
do j=0,1000
|
|
|
|
|
nsum=nsum+hist(j)
|
|
|
|
|
if(nsum.ge.ntest) exit
|
|
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
xmedian=j*ave/100.0
|
2012-10-31 14:33:56 -04:00
|
|
|
if(j.lt.10) xmedian=1.0
|
2012-10-11 14:33:50 -04:00
|
|
|
|
2012-10-31 14:33:56 -04:00
|
|
|
900 continue
|
2012-10-11 14:33:50 -04:00
|
|
|
return
|
|
|
|
|
end subroutine pctile
|
