
Hi,

This is the setup as presented in Chen & Rogers (2016). Please run the run_S.py (modified after original python scripy by Phil Arras)
for the complete simulation as documented in the paper. Basically, choose the planet mass, envelope fraction, and initial entropy you want.
There are several tweaks required for difference regimes, for instance at very low mass or very low H/He mass fractions.

Please contact me at howard@earth.northwestern.edu if you have issues running the code.


Best regards,
Howard

