Package: ctplanet
ctplanetCreate crustal thickness maps of planets from gravity and topography.Descriptionctplanet provides several functions and example scripts for generating crustal thickness maps of a planet from gravity and topography data, and the calculation of hydrostatic relief along density interfaces beneath the lithosphere.MethodspyMohoCalculate relief using a constant density crust and
mantle.pyMohoRhoCalculate relief using a constant density mantle and a
variable density crust.HydrostaticShapeLithCalculate the relief of hydrostatic interfaces beneath
the lithosphere along with the predicted gravity,
taking into account rotation and/or tides using the
approach ofWieczorek et al.(2019).HydrostaticShapeCalculate the relief of hydrostatic interfaces and
predicted gravity of a rotating hydrostatic planet
using the approach ofWieczorek et al.(2019).InertiaTensor_from_shapeCalculate the inertia tensor given a radial density
profile and shape of each interface.InertiaTensor_from_CCalculate the inertia tensor given the polar moment
of inertia and the gravitational potential
coefficients.moiCalculate the mean, normalized, moment of inertia
up to index n.ReadRefModelRead the reference interior model file.Example scriptsMoon-CrustA script that demonstrates how to calculate the thickenss
of the lunar crust using either a constant or variable density
crust. The latter can be used to reproduce the results
presented inWieczorek et al.(2013).Moon-CoreCalculate the hydrostatic relief of the lunar core accounting
for the non-hydrostatic potential that comes from the
lithosphere.Mars-CrustA script that demonstrates how to calculate the thickenss
of the Martian crust using either a constant or variable
density crust. For the variable density crust, the density is
assumed to change discontinuously across the dichotomy
boundary.Mars-Crust-hydrostatic-testsCreate a crustal thickness map of Mars from
gravity and topography and compare how results
change if hydrostatic interfaces are not taken
into account.Mars-Crust-InSightCreate a crustal thickness map of Mars from gravity
and topography, using the InSight crustal thickness
constraint.Mars-Crust-InSight-dichotomyCreate a crustal thickness map of Mars from
gravity and topography, using the InSight
crustal thickness constraint and different
densities across the dichotomy boundary.Mars-fcnCompute the free core nutation period of Mars.Mars-shapeCreate images related to Mars inWieczorek et al.(2019).Mars-j2Compute the contribution to the gravitational J2 of Mars from
hydrostatic interfaces beneath the lithosphere.Earth-shapeCompute hydrostatic relief of Earth using PREM.Ceres-shapeCalculate the hydrostatic shape of Ceres.How to install and run ctplanetDownload the ctplanet repository and install using pipgitclonehttps://github.com/MarkWieczorek/ctplanet.gitpipinstall.To execute a scriptcdexamplespythonMoon-Crust.pyDepending on how your system is set up, it might be necessary to use explicitlypython3andpip3instead ofpythonandpipin the above commands.ReferenceWieczorek, M. A., G. A. Neumann, F. Nimmo, W. S. Kiefer, G. J. Taylor,
H. J. Melosh, R. J. Phillips, S. C. Solomon, J. C. Andrews-Hanna,
S. W. Asmar, A. S. Konopliv, F. G. Lemoine, D. E. Smith, M. M. Watkins,
J. G. Williams, M. T. Zuber (2013), The crust of the Moon as seen by GRAIL,Science, 339, 671-675, doi:10.1126/science.1231530.Wieczorek, M. A.,  M. Beuthe, A. Rivoldini, and T. Van Hoolst (2019),
Hydrostatic interfaces in bodies with nonhydrostatic lithospheres,Journal of Geophysical Research: Planets, 124, doi:10.1029/2018JE005909.