Maps: Parametrized LayerΒΆ

Build a model of a parametrized layer in a wholespace. If you want to build a model of a parametrized layer in a halfspace, also use Maps.InjectActiveCell.

The model is

m = [
    'background physical property value',
    'layer physical property value',
    'layer center',
    'layer thickness'


/Users/lindseyjh/git/simpeg/simpeg/examples/03-maps/ UserWarning: Matplotlib is currently using agg, which is a non-GUI backend, so cannot show the figure.

from SimPEG import Mesh, Maps
import numpy as np
import matplotlib.pyplot as plt

def run(plotIt=True):

    mesh = Mesh.TensorMesh([50, 50], x0='CC')  # 2D tensor mesh
    mapping = Maps.ParametricLayer(mesh)  # parametric layer in wholespace

    # model
    m = np.hstack(
            1., # background value
            2., # layer value
            -0.1, # layer center
            0.2 # layer thickness
    rho = mapping * m # apply the mapping

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

if __name__ == '__main__':

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

Gallery generated by Sphinx-Gallery