ctplanet.HydrostaticShapeLith

ctplanet.HydrostaticShapeLith#

ctplanet.HydrostaticShapeLith(radius, rho, ilith, potential, topo, rho_surface, r_sigma, omega, lmax, rp=None, mp=None, nmax=7)#

Calculate the shape of hydrostatic relief in a rotating planet or moon with a non-hydrostatic lithosphere, along with the total gravitation potential of the hydrostatic interfaces. For the case of a moon in synchronous rotation, optionally include the tidal potential.

Returns:

  • hlm (array of SHCoeffs class instances, size(n+1)) – Array of SHCoeffs class instances of the spherical harmonic coefficients of the hydrostatic relief at each interface.

  • clm_hydro (SHCoeffs class instance containing the gravitational potential) – resulting from all hydrostatic interfaces.

  • mass (float) – Total mass of the planet, assuming a spherical shape and the provided 1D density profile.

Parameters:

  • radius (ndarray, float, size(n+1)) – Radius of each density interface, where index 0 corresponds to the center of the planet and n corresponds to the surface.

  • density (ndarray, float, size(n+1)) – Density of layer i between radius[i] and radius[i+1]. The density at index 0 is from the center of the planet to radius[1], whereas the the density at index n should be zero.

  • ilith (int) – Index of the interface that corresponds to the base of the lithosphere.

  • potential (SHGravCoeffs class instance) – Observed gravitational potential coefficients.

  • topo (SHCoeffs class instance) – Observed shape of the planet.

  • rho_surface (float) – Effective density of the surface relief.

  • r_sigma (float) – Radius of the mass sheet source in the lithosphere.

  • omega (float) – Angular rotation rate of the planet.

  • lmax (int) – Maximum spherical harmonic degree to compute for the hydrostatic relief at each interface.

  • rp (float, optional, default = None) – If specified, include the tidal potential acting on a synchronously rotating moon, where rp is the average distance between the planet and satellite.

  • mp (float, optional, default = None) – The mass of the host planet, at a distance rp from the satellite.

  • nmax (int, optional, default = 7) – The order of the approximation when computing the gravitational potential of the surface.