[Mesa-users] White dwarf cooling for different compositions

[Na'ama Hallakoun]

Hi,

I'm new to MESA, and I'm trying to produce white dwarf cooling tracks for
different compositions (i.e. changing the width of the He and H layers).
I'm using the wd_cool_0.6M example from the test suite, where I added to
the inlist_wd_cool_0.6M under &star_job:

num_steps_to_relax_composition = 1 ! use this many steps to do conversion
relax_composition_filename = 'composition.dat' ! file holding composition
profile information

and in composition.dat I have, for example:
6 8
0.00000 0.0 0.0 0.0 0.5 0.0 0.5 0.0 0.0
0.98999 0.0 0.0 0.0 0.5 0.0 0.5 0.0 0.0
0.99000 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0
0.99989 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0
0.99990 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
1.00000 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0

The resulting abundance plot seems fine, but I don't see any difference in
the cooling track when I edit the composition.dat file. What am I doing
wrong?

In addition, the wd_cool_0.6M test uses a wd_0.6.mod model that starts from
an effective temperature of ~17000K, so it messes up the age. How can I
start from a hotter WD?

Thanks!
Na'ama
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180709/c90d6b51/attachment.html>