ctplanet.pyMohoRho

Contents

ctplanet.pyMohoRho#

ctplanet.pyMohoRho(pot, topo, density, porosity, lmax, rho_m, thickave, filter_type=0, half=None, nmax=8, delta_max=5.0, lmax_calc=None, correction=None, quiet=False)#

Calculate the relief along the crust-mantle interface assuming a constant density mantle and a laterally varying crustal density.

Returns:

moho – The radius of the crust-mantle interface.

Return type:

SHCoeffs class instance

Parameters:

  • pot (SHGravCoeffs class instance) – Gravitational potential spherical harmonic coefficients.

  • topo (SHCoeffs class instance) – Spherical harmonic coefficients of the surface relief.

  • density (SHCoeffs class instance) – Spherical harmonic coefficients of the crustal grain density.

  • porosity (float) – Crustal porosity (from 0 to 1).

  • lmax (int) – Maximum spherical harmonic degree of the function, which determines the sampling interval of the internally computed grids.

  • rho_m (float) – Mantle density in kg / m^3.

  • thickave (float) – Average thickness of the crust in meters.

  • filter_type (int, optional, default = 0) – 0 = no filtering, 1 = minimum amplitude filter, 2 = minimum curvature filter.

  • half (float, optional, default = None) – The spherical harmonic degree where the filter is equal to 0.5. This must be set when filter_type is 1 or 2.

  • nmax (int, optional, default = 8) – The maximum order used in the Taylor-series expansion when calculating the potential coefficients.

  • delta_max (float, optional, default = 5.0) – The algorithm will continue to iterate until the maximum difference in relief between solutions is less than this value (in meters).

  • lmax_calc (int) – Maximum spherical harmonic degree when evalulating the functions.

  • correction (SHGravCoeffs class instance, optional, default = None) – If present, these coefficients will be added to the Bouguer correction (subtracted from the Bouguer anomaly) before performing the inversion. This could be used to account for the pre-computed gravitational attraction of the polar caps of Mars, which have a different density than the crust.

  • quiet (boolean, optional, default = False) – If True, suppress printing output during the iterations.

Contents