Maps: Parametric Block in a Layer#

Parametric description of a block confined to a layer in a wholespace. The mapping can be applied in 2D or 3D. Here we show a 2D example.

The model is given by

m = np.r_[
   'value of the background',
   'value in the layer',
   'value in the block',
   'center of the layer (depth)',
   'thickness of the layer',
   'x-center of block',
   'width of the block'
plot block in layer
import discretize
from SimPEG import maps
import numpy as np
import matplotlib.pyplot as plt

def run(plotIt=True):
    mesh = discretize.TensorMesh([50, 50], x0="CC")  # 2D Tensor Mesh
    mapping = maps.ParametricBlockInLayer(mesh)  # mapping

    m = np.hstack(
            1.0,  # value of the background
            2.0,  # value in the layer
            3.0,  # value in the block
            -0.1,  # center of the layer (depth)
            0.2,  # thickness of the layer
            0.3,  # x-center of block
            0.2,  # width of the block

    # apply the mapping to define the physical property on the mesh
    rho = mapping * m

    if plotIt is True:
        fig, ax = plt.subplots(1, 1, figsize=(4, 6))
        mesh.plot_image(rho, ax=ax)

if __name__ == "__main__":

Total running time of the script: (0 minutes 0.262 seconds)

Estimated memory usage: 8 MB

Gallery generated by Sphinx-Gallery