[mesa-users] Deuterium burning in PMS stars
michal.simon at stonybrook.edu
Fri Dec 11 15:09:45 EST 2015
Hi Aaron et al,
Your suggested fixes worked well. I learned a lot over the last few
days of often frustrating MESAing from you, Josiah, and Robert Farmer.
Thanks to all!
What MESA needs is not a thick text but a "MESA for Dummies".
On Thu, Dec 10, 2015 at 10:35 PM, Aaron Dotter <aaron.dotter at gmail.com>
> Hi Mike,
> It's safest to be explicit with things like this, i.e., specify the
> initial abundances in your inlist rather than rely on whatever MESA does by
> Default behavior may change over time (though I'm not certain that it has
> in this case).
> On Fri, Dec 11, 2015 at 2:25 PM, Michal Simon <michal.simon at stonybrook.edu
> > wrote:
>> Hi Aaron, Jayashree and I will try this fix tmw, Fri. But I am very
>> puzzled. Mesa 6596 pulled values of center h1, h2, li7 out of its table of
>> abundances, namely 0.7, 2e-5, etc correctly. Does Mesa 7624 work
>> differently? Mike
>> Sent from my iPhone
>> On Dec 10, 2015, at 5:29 PM, Aaron Dotter <aaron.dotter at gmail.com> wrote:
>> Hi Mike,
>> What do you expect for center_h2?
>> In your inlist you have not specified any initial composition information
>> so you get some default values, probably not what you want. May I suggest
>> you try setting these values in the star_job section of your inlist:
>> set_uniform_initial_composition = .true.
>> !### initial_h1
>> !### initial_h2
>> !### initial_he3
>> !### initial_he4
>> ! if `set_uniform_initial_composition` is true, then
>> ! set hydrogen and helium mass fractions according to the
>> ! If no h2 in current net, then this will be added to h1.
>> ! If no he3 in current net, then this will be added to he4.
>> (your preferred values here)
>> initial_h1 = -1
>> initial_h2 = -1
>> initial_he3 = -1
>> initial_he4 = -1
>> !### initial_zfracs
>> ! if `set_uniform_initial_composition` is true, then set metal
>> ! z fractions -- select one of the options defined in
>> `chem/public/chem_def` :
>> ! + `AG89_zfracs = 1`
>> ! + `GN93_zfracs = 2`
>> ! + `GS98_zfracs = 3`
>> ! + `L03_zfracs = 4`
>> ! + `AGS05_zfracs = 5`
>> ! + `AGSS09_zfracs = 6`
>> ! + `L09_zfracs = 7`
>> ! + `A09_Prz_zfracs = 8`
>> (from your inlist it seems you want GS98 for the solar composition so set
>> initial_zfracs = 3
>> On Fri, Dec 11, 2015 at 7:00 AM, Michal Simon <
>> michal.simon at stonybrook.edu> wrote:
>>> Dear Colleagues,
>>> I am running MESA 7624 and trying to determine the time when the very
>>> brief phase of D burning operates in PMS stars. I am having the history
>>> file list the value of center_h2 as the star contracts. The values for the
>>> fraction of deuterium are nonsense.
>>> What am I doing wrong? I attach inlist_project and
>>> history.data files.
>>> Cheers and thanks,
>>> Mike Simon
>>> mesa-users mailing list
>>> mesa-users at lists.sourceforge.net
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Mesa-users