This code follows the thermal evolution of a neutron star crust. To compile the eos and crust modules with cython (pip install cython), use
python setup.py build_ext --inplace
An IPython notebook is included that shows how to run a model of the MXB 1659 lightcurve (see Brown & Cumming 2009). Alternatively, use python crustcool.py which plots the results in crustcool.pdf.
During accretion, the temperature at the top of the grid can be held constant, or allowed to cool and a shallow heat source included at low density.
Interactive version using bokeh:
bokeh serve --show crustcool_bokeh.py
To-do list:
- add entrainment of neutrons by the nuclei in the inner crust (no entrainment assumed at the moment)