[mesa-users] Deuterium burning in PMS stars
aaron.dotter at gmail.com
Thu Dec 10 22:35:51 EST 2015
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>
> 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