anchor.diff:      real(dp), allocatable :: x(:), y(:), z(:)
       integer, allocatable :: neighbors(:,:)
       integer :: n_neighbors
  -    real(dp),allocatable :: points(:,:), bmin(:,:), bmax(:,:)
  +    real(dp),allocatable :: bmin(:,:), bmax(:,:)
       type(cell),pointer :: c
   
       ! Read geometry file
  
positive.code:   type(grid_geometry_desc),public,target :: geo
  
    integer :: n_filled = 0
  
  contains
  
  
negative.code:     integer :: i, n_neighbors
      type(vector3d_dp) :: r0, ri, rmid, n
      real(dp) :: t, tx, ty, tz
  
  
      icell => geo%cells(p%icell%ic)
  
  