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:                  nx,ny,nz,xmn,ymn,zmn,xsiz,ysiz,zsiz, extve,  &         ! input grid parameters (no rotation allowed)
                   nxdis,nydis,nzdis,                          &
                   radius,radius1,radius2,sang1,sang2,sang3, &            ! input search parameters
                   ndmax, ndmin, noct, bhidmax, MAXSBX, MAXSBY, MAXSBZ,  &
                   nst,c0,it,cc,aa,aa1,aa2,ang1,ang2,ang3,   &            ! input variogram
                   ktype, skmean, itrend, UNEST, &                        ! input kriging parameters 
                   idrift, &                                              ! input drift parameters                         
  
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)
  
  