resqpy.grid_surface.bisector_from_faces

resqpy.grid_surface.bisector_from_faces(grid_extent_kji: Tuple[int, int, int], k_faces: ndarray, j_faces: ndarray, i_faces: ndarray, raw_bisector: bool) Tuple[ndarray, bool][source]

Returns a numpy bool array denoting the bisection of the grid by the face sets.

Parameters:
  • grid_extent_kji (triple int) – the shape of the grid

  • k_faces (numpy bool arrays) – True where an internal grid face forms part of the bisecting surface

  • j_faces (numpy bool arrays) – True where an internal grid face forms part of the bisecting surface

  • i_faces (numpy bool arrays) – True where an internal grid face forms part of the bisecting surface

Returns:

(numpy bool array of shape grid_extent_kji, bool) where the array is set True for cells on one side of the face sets deemed to be shallower (more strictly, lower K index on average); set False for cells on othe side; the bool value is True if the surface is a curtain (vertical), otherwise False

Notes

the face sets must form a single ‘sealed’ cut of the grid (eg. not waving in and out of the grid); any ‘boxed in’ parts of the grid (completely enclosed by bisecting faces) will be consistently assigned to either the True or False part