anchor.diff:      integer :: i1, i2, i3
       type(vector3d_dp) :: r
       real(dp) :: min_width
  +
  +    type(grid_cell) :: cell
   
       ! Read geometry file
       call mp_read_keyword(group, '.', "geometry", geo%id)
  
positive.code:        cell = find_position_in_grid(r, ilevel_new, igrid_new)
      end if
    end function find_position_in_grid
  
  
    subroutine grid_geometry_debug(debug_flag)
      implicit none
  
negative.code:     integer(hid_t) :: g_level
      character(len=100) :: level_name
  
      integer :: igrid, igrid1, igrid2, icell, iv
  
      integer :: start_id
      type(level_desc), pointer :: level, level1, level2
  