The current construction relies on two separate separate EarthGeometry objects. One is a subclass that keeps track of the depth and calculates the overburden, and the other is used in the density profile. It would be nice to synchronize these.
As a corollary, the geometry used should allow for correct treatment of non-zero h_obs.