ctplanet.HydrostaticShape

ctplanet.HydrostaticShape(radius, rho, omega, gm, r_ref, rp=None, mp=None, i_clm_hydro=None)

Calculate the shape of hydrostatic relief in a rotating planet or moon, along with the total gravitation potential. 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. If i_clm_hydro is specified, then the potential will include only interfaces beneath index i_clm_hydro.

• 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.

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

• gm (float) – GM of the planet.

• r_ref (float) – Refernce radius for output potential coefficients.

• 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.

• i_clm_hydro (int, optional, default = None) – If specified, calculate the gravitational potential clm_hydro resulting from all interfaces below and including the radius index i_clm_hydro.