anchor.diff:      !     radius, radius1, radius2 - double: ellipse radius
       !     sang1, sang2, sang3      - double: ellipse rotation angles
       !     noct                     - integer: maximum per octant
  -    !     bhidmax                  - integer: maximum per drillhole
  +    !     maxbhid                  - integer: maximum per drillhole
       !     MAXSBX, MAXSBY, MAXSBZ   - integer: super-block definition
       !                                    ** make sure the search ellipse    
       !                                       fits in it. 
  
positive.code:     ! search
      integer, intent(in) :: ndmax, ndmin
      real*8, intent(in) :: radius, radius1, radius2, sang1, sang2, sang3
      integer, intent(in) :: noct, bhidmax
      integer, intent(in) :: MAXSBX, MAXSBY, MAXSBZ                        ! Don't like this, remove and use kdtree
      ! variogram
      integer, intent(in)  :: nst
  
negative.code:         ! Find the nearest samples:
      
          call srchsupr(xloc,yloc,zloc,radsqd,isrot,MAXROT,rotmat,nsbtosr, &
          ixsbtosr,iysbtosr,izsbtosr,noct, maxbhid,nd,x,y,z,bhid,tmp, &
          nisb,nxsup,xmnsup,xsizsup,nysup,ymnsup,ysizsup, &
          nzsup,zmnsup,zsizsup,nclose,close,infoct)
  
  