[Mesa-users] Radius of a planet

[Josiah Schwab]

Hi Claudio,

> I am currently using Mesa to study the initial stage of the growth of a
> Giant planet. I implemented correctly the subroutine "other_wind" to
> accrete solids and gas, however in codes different from Mesa the radius
> of the planet is usually defined as the Bondi radius, the Hill radius,
> or a combination of the two. In order to compare my preliminary results
> with the ones obtained with other codes I would be interested to modify
> in Mesa the computation of the radius and compare the profile of
> temperature density and pressure that I get. However diving in the code
> I have not understood how the radius of the planet is computed and how
> to modify it. The outer boundary condition that I am using is
> "fixed_Psurf_and_Tsurf"

The radius reported by MESA will be the radius of the surface cell
(i.e. the location where your outer boundary condition is applied.)

If you would like to compute and report a different radius using the
properties of the model, you can do so during the calculation using
run_star_extras.

Josiah