resqpy.grid_surface.get_boundary_dict

resqpy.grid_surface.get_boundary_dict(k_faces: Optional[ndarray], j_faces: Optional[ndarray], i_faces: Optional[ndarray], grid_extent_kji: Tuple[int, int, int]) → Dict[str, int]

Cretaes a dictionary of the indices that bound the surface (where the faces are True).

Parameters

• k_faces (np.ndarray) – a boolean array of which faces represent the surface in the k dimension

• j_faces (np.ndarray) – a boolean array of which faces represent the surface in the j dimension

• i_faces (np.ndarray) – a boolean array of which faces represent the surface in the i dimension

• grid_extent_kji (Tuple[int, int, int]) – the shape of the grid

Returns

boundary (Dict[str, int]) – a dictionary of the indices that bound the surface

Note

input faces arrays are for internal grid faces (ie. extent reduced by 1 in axis of faces); a buffer slice is included where the surface does not reach the edge of the grid; max values are not increment, ie. need to be incremented to be used as an upper end of a python range