[Mesa-users] White dwarf cooling for different compositions

[Na'ama Hallakoun]

Hi Rob,

Thanks for the quick reply!
Indeed I realized I forgot to set

relax_initial_composition = .true.

But now when I do it, I get
relax_initial_composition
 do_relax_initial_composition ierr          -1
 do_star_job_controls_after ierr          -1
 before_evolve_loop ierr          -1

I'm attaching the inlist.
For the composition.dat I currently tried:

2 8
0.00000 0.0 0.0 0.0 1 0.0 0 0.0 0.0
1.00000 0.0 0.0 0.0 1 0.0 0 0.0 0.0

Thanks,
Na'ama

On Mon, Jul 9, 2018 at 6:06 PM Rob Farmer <r.j.farmer at uva.nl> wrote:

> Hi
> Showing inlists is the easiest way for us to help you. But did you also
> set
>
> relax_initial_composition = .true.
>
>
> Rob
>
> On Mon, 9 Jul 2018 at 17:36, Na'ama Hallakoun <naama at wise.tau.ac.il>
> wrote:
>
>> 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
>> _______________________________________________
>> mesa-users at lists.mesastar.org
>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180709/38233eb3/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_wd_cool_0.6M
Type: application/octet-stream
Size: 6291 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180709/38233eb3/attachment.obj>