Initial Conditions

Tools for the creation of initial conditions are available via the WoMa and SEAGen open-source python packages, including: creating spherical or spinning planetary (or similar) profiles; placing particles to match arbitrary profiles with precise SPH densities; and setting the initial target and impactor positions and velocities, as presented in Kegerreis et al. (2019) and Ruiz-Bonilla et al. (2020).

They are available with documentation and examples at https://github.com/srbonilla/WoMa and https://github.com/jkeger/seagen, or can be installed directly with pip (https://pypi.org/project/woma/, https://pypi.org/project/seagen/).

Settling initial conditions with fixed entropies

If the particles’ equations of state include specific entropies, and the initial conditions file includes specific entropies for each particle (in PartType0/Entropies), then configuring SWIFT with --enable-planetary-fixed-entropy will override the internal energy of each particle each step such that its specific entropy remains constant.

This should be used with caution, but may be a convenient way to maintain an entropy profile while initial conditions settle to equilibrium with their slightly different SPH densities.