|
Planet.dzClath_km = np.max(Planet.z_m[:-1][abs(Planet.phase) == Constants.phaseClath])/1e3 \ |
This calculation (and probably all others) for layer thicknesses in SummaryTables.py uses the radii of the uppermost and lowermost layers for a particular phase to calculate layer thicknesses. The bottom layer should use the outer radius of the next material, not the outer radius of the bottom layer of the present material. The index of the bottom layer must be increased by 1, but innermost layers will need to account for the extra r = 0 index.
PlanetProfile/PlanetProfile/Utilities/SummaryTables.py
Line 157 in d99caea
This calculation (and probably all others) for layer thicknesses in SummaryTables.py uses the radii of the uppermost and lowermost layers for a particular phase to calculate layer thicknesses. The bottom layer should use the outer radius of the next material, not the outer radius of the bottom layer of the present material. The index of the bottom layer must be increased by 1, but innermost layers will need to account for the extra r = 0 index.