search

vancesteven / PlanetProfile

Construct 1D interior structure models based on planetary properties. Self-consistent thermodynamics are used for fluid, rock, and mineral phases. Sound speeds, attenuation, and electrical conductivities are computed as outputs.
https://vancesteven.github.io/PlanetProfile/
Other
18 stars 13 forks source link

Gravity inside silicate and core layers are set to constants #11

Open itsmoosh opened 6 years ago

itsmoosh commented 6 years ago

According to Vance et al. 2018, inside the silicate layer, gravity is set to one mean value, and in the core, gravity is set to zero. This seems to be born out by output plots of the gravitational acceleration.

g(r) is calculated within the ice and ocean layers here: https://github.com/vancesteven/PlanetProfile/blob/b751522098010750f25da8fcc4472169e1203468/PlanetProfile.m#L177

There appears to be commented lines propagating g to the deeper layers. This seems like a beneficial feature for profiles considering the deeper interior.

itsmoosh commented 4 years ago

This appears to cause up to 30% too-large values for Callisto's mass.

itsmoosh commented 2 years ago

This has been fixed in the Python version; future revamping of the Matlab version will incorporate the same algorithms.