resqpy.grid_surface.generate_torn_surface_for_layer_interface

resqpy.grid_surface.generate_torn_surface_for_layer_interface(grid, k0=0, ref_k_faces='top', quad_triangles=True)[source]

Returns a Surface object generated from the grid layer interface points.

Parameters:
  • grid (grid.Grid object) – the grid object from which a layer interface is to be converted to a surface

  • k0 (int) – the layer number (zero based) to be used

  • ref_k_faces (string) – either ‘top’ (the default) or ‘base’, indicating whether the top or the base interface of the layer is to be used

  • quad_triangles (boolean, optional, default True) – if True, 4 triangles are used to represent each cell k face, which gives a unique solution with a shared node of the 4 triangles at the mean point of the 4 corners of the face; if False, only 2 triangles are used, which gives a non-unique solution

Returns:

a resqml_surface.Surface object with a single triangulated patch

Notes

The resulting surface is assigned to the same model as grid, though xml is not generated and hdf5 is not written. Strictly, the RESQML business rules for a triangulated surface require a separate patch for areas of the surface which are not joined; therefore, if fault tears cut off one area of the surface (eg. a fault running fully across the grid), then more than one patch should be generated; however, at present the code uses a single patch regardless.