bsort.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	bsort (n, x, y, xmin, xmax, ymin, ymax, dmax, bin, list)

Function/Subroutine Documentation

bsort()

subroutine bsort	(	integer	n,
		real*8, dimension(*)	x,
		real*8, dimension(*)	y,
		real*8	xmin,
		real*8	xmax,
		real*8	ymin,
		real*8	ymax,
		real*8	dmax,
		integer, dimension(*)	bin,
		integer, dimension(*)	list
	)

!     
      implicit none
!     
      integer list(*),bin(*),n,i,j,k,p,ndiv
!     
      real*8 x(*),y(*),factx,facty,xmin,xmax,ymin,ymax,dmax
!     
      ndiv=nint(real(n)**0.25d0)
      factx=real(ndiv)/((xmax-xmin)*1.01d0/dmax)
      facty=real(ndiv)/((ymax-ymin)*1.01d0/dmax)
      do 10 k=1,n
         p=list(k)
         i=int(y(p)*facty)
         j=int(x(p)*factx)
         if(mod(i,2).eq.0)then
            bin(p)=i*ndiv+j+1
         else
            bin(p)=(i+1)*ndiv-j
         end if
 10   continue
      call qsorti(n,list,bin)