SimPEG.utils.surface2ind_topo#

SimPEG.utils.surface2ind_topo(mesh, topo, gridLoc='CC', method='nearest', fill_value=nan)[source]#

Get indices of active cells from topography.

For a mesh and surface topography, this function returns the indices of cells lying below the discretized surface topography.

Parameters

meshdiscretize.TensorMesh or discretize.TreeMesh: Mesh on which you want to identify active cells
topo(n, 3) numpy.ndarray: Topography data as a numpyndarray with columns [x,y,z]; can use [x,z] for 2D meshes. Topography data can be unstructured.
gridLocstr {‘CC’, ‘N’}: If ‘CC’, all cells whose centers are below the topography are active cells. If ‘N’, then cells must lie entirely below the topography in order to be active cells.
methodstr {‘nearest’,’linear’}: Interpolation method for approximating topography at cell’s horizontal position. Default is ‘nearest’.
fill_valuefloat: Defines the elevation for cells outside the horizontal extent of the topography data. Default is numpy.nan.

Returns