This work directory is a prototype of using the other_build_initial_model hook to make a WD.

The code’s not great, but mostly seems to work. It was written against MESA r12115 and later updated for r12778. (Look at the git tags to see which commits should work for these versions).

Revised and improved code targeting the MESA development version is under construction elsewhere.

This sets the initial central temperature of the WD.

If true, this reads the composition from relax_composition_filename in star_job.

If false, it the subroutine get_xa_from_code in src/model_builder.f90 is called. This routine must provide the composition as a function of q.

If false, this assumes the coordinate in relax_composition_filename is xq (= 1 - q). That is the format typically assumed by relax_initial_composition.

If true, this assumes the coordinate in relax_composition_filename is q. The included file composition_test.dat uses q and mesa_49.net.