planit 0.8.2

The PLANetary Impact Toolkit provides tools for accessing, manipulating, and analysing smoothed particle hydrodynamcis (SPH) impact simulation data, and creating planets for SPH impact simulations.

planit

The Planetary Impact Toolkit provides tools for accessing, manipulating, and analysing smoothed particle hydrodynamcis (SPH) impact simulation data, and creating planets for SPH impact simulations.

Philip J. Carter ([e-mail])
Sarah T. Stewart ([e-mail])

Overview

planit currently supports the SPH codes Gadget2-planetary and SWIFT.

makeplanet provides functions for generating adiabatic 1D planet structures and converting these into (unequilibrated) SPH planets.
snaptools provides the Snapshot class for holding input and output files, calculating remnant masses, and vapor fractions.
impacttools provides the Impact class for loading and plotting snapshots produced by an impact simulation.
eos provides classes and functions for accessing and using equations of state.
utils provides additional utility functions for analysis.

Installation

Equations of State (EoS)

planit currently includes the following EoS:

• ANEOS Forsterite (Stewart et al., 2019): ANEOSForsterite
• ANEOS Iron (Stewart, 2020): ANEOSIron
• ANEOS Fe85Si15 (Stewart, 2020): ANEOSFeSiAlloy
• ANEOS Pyrolite (Stewart et al., 2022): ANEOSPyrolite
• 5-phase Water (Senft & Stewart, 2008): 5PhaseWater
• HM80 Hydrogen-Helium (Hubbard & MacFarlane, 1980; Lock & Stewart, 2024): HM80HHe

They can be accessed using the names given above (as well as some variations) or using the corresponding SWIFT/WoMa ID number.
User supplied tables in standard SESAME format can also be used.

[!NOTE] ANEOSPyrolite has been assigned the SWIFT/WoMa ID 403 in planit but does not currently have an assigned ID in SWIFT. Using a modified version of SWIFT or using the user supplied EoS functionality will be required.

Planet generation

planit.makeplanet differs from other planet creation code in that produces planets with isentropic layers, integrates outwards from the centre of the planet, and ensures the requested mass and component mass fractions are matched (within tolerance). The default behaviour prevents the unphyiscal scenario of a core cooler than the mantle at the core-mantle boundary.

[!NOTE] makeplanet currently only supports 2-layer planets