anchor.diff:      !   iysbtosr         Y offsets for super blocks to search
       !   izsbtosr         Z offsets for super blocks to search
       !   noct             If >0 then data will be partitioned into octants
  +    !   maxbhid          Maximum number of samples per BHID
       !   nd               Number of data
       !   x(nd)            X coordinates of the data
       !   y(nd)            Y coordinates of the data
       !   z(nd)            Z coordinates of the data
  +    !   bhid(nd)         BHID 
       !   tmp(nd)          Temporary storage to keep track of the squared
       !                      distance associated with each data
       !   nisb()                Array with cumulative number of data in each
  
positive.code: 
  
  subroutine srchsupr(xloc,yloc,zloc,radsqd,irot,MAXROT,rotmat, &
      nsbtosr,ixsbtosr,iysbtosr,izsbtosr,noct,nd, &
      x,y,z,tmp,nisb,nxsup,xmnsup,xsizsup, &
      nysup,ymnsup,ysizsup,nzsup,zmnsup,zsizsup, &
      nclose,close,infoct)
      !-----------------------------------------------------------------------
  
negative.code: 
      !input
      integer, intent(in) :: nxsup, nysup, nzsup
      integer, intent(in) :: irot, MAXROT, nsbtosr, noct, maxbhid, nd
      real*8, intent(in), dimension (MAXROT,3,3) :: rotmat
      real*8, intent(in) ::  xloc,yloc,zloc,radsqd 
      integer, intent(in), dimension(nsbtosr) ::ixsbtosr, iysbtosr, izsbtosr
      integer, intent(in), dimension(nd) :: bhid
      real*8, intent(in), dimension(nd) :: x, y, z 
      integer, intent(in), dimension (nxsup*nysup*nzsup) :: nisb
      real*8, intent(in) :: xmnsup,xsizsup,ymnsup, ysizsup, zmnsup, zsizsup
  
  