anchor.diff:      ! search
       integer, intent(in) :: ndmax, ndmin
       real*8, intent(in) :: radius, radius1, radius2, sang1, sang2, sang3
  -    integer, intent(in) :: noct, bhidmax
  +    integer, intent(in) :: noct, maxbhid
       integer, intent(in) :: MAXSBX, MAXSBY, MAXSBZ                        ! Don't like this, remove and use kdtree
       ! variogram
       integer, intent(in)  :: nst
  
positive.code: 
      !input
      integer, intent(in) :: nxsup, nysup, nzsup
      integer, intent(in) :: irot, MAXROT, nsbtosr, noct, 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
      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
  
  
negative.code:     !     radius, radius1, radius2 - double: ellipse radius
      !     sang1, sang2, sang3      - double: ellipse rotation angles
      !     noct                     - integer: maximum per octant
      !     maxbhid                  - integer: maximum per drillhole
      !     MAXSBX, MAXSBY, MAXSBZ   - integer: super-block definition
      !                                    ** make sure the search ellipse    
      !                                       fits in it. 
  