[Mesa-users] Stop at RGB tip

[Bill Wolf]

Hi Iñigo,

I don't think there is a command for that built in, but you could create a custom stopping condition based on s% power_he_burn / s% power_h_burn in run_star_extras.f. For inspiration, see Josiah's guide from the last summer school here:

https://jschwab.github.io/mesa-2018/#orgdae7a6d

Regards,

Bill


From: Iñigo Muñoz Elorza via Mesa-users <mesa-users at lists.mesastar.org>
Reply: Iñigo Muñoz Elorza <imunozelorza at gmail.com>
Date: November 25, 2018 at 4:26:18 PM
To: cmankovich at ucsc.edu <cmankovich at ucsc.edu>
Cc: mesa-users <mesa-users at lists.mesastar.org>
Subject:  Re: [Mesa-users] Stop at RGB tip  

Hi Pavel, Chris, all

Related to Adrian's question about stopping the model after the RGB tip, I would like to know if it is possible to impose a control option that works with ratios instead of the total power as suggested by Pavel, like

stop when ratio between energy from triple alpha and the one from H fusion in the shell = 100.

According to the documentation all the stop conditions seem to be single conditions, but I wonder if this kind of ratios/calculations could be defined.

Many thanks in advance for the support.

Cheers,
Iñigo

El mar., 20 nov. 2018 a las 2:21, Chris Mankovich via Mesa-users (<mesa-users at lists.mesastar.org>) escribió:
Hi Adrian,

I'd second Pavel's suggestion, and add the comment that a drop in central helium abundance might not coincide with the RGB tip for all masses and compositions. For instance, a model that will undergo helium core flashes will reach the RGB tip well before helium burning reaches the core, at least in terms of execution time.

Chris

On Mon, Nov 19, 2018 at 5:08 PM Pavel Denisenkov via Mesa-users <mesa-users at lists.mesastar.org> wrote:
Hi Adrian,

you can use the following control option:

  ! stop when total power from reactions burning helium (in Lsun units) is > this.
    power_he_burn_upper_limit = 1d3

Pavel

On Nov 19, 2018, at 16:06, Adrian Ayala <aayala at ugr.es> wrote:

Dear all,

I would like to stop some models just after the RGB tip but I didn't found this in the standard stop conditions, in the code documentation.

I wonder if were it a good criterion to stop the simulation when central helium mass fraction falls a bit (the starting of He burning). What do you think?

Thanks so much in advance!
Best wishes,

Adrián

--
Adrian Ayala
---------------------
PhD in Space Science and Astrophysics
Dep of Theoretical Physics and Astrophysics
University of Granada
_______________________________________________
mesa-users at lists.mesastar.org
https://lists.mesastar.org/mailman/listinfo/mesa-users


_______________________________________________
mesa-users at lists.mesastar.org
https://lists.mesastar.org/mailman/listinfo/mesa-users



--
Christopher R. Mankovich
PhD Candidate, Astronomy & Astrophysics
University of California Santa Cruz
_______________________________________________
mesa-users at lists.mesastar.org
https://lists.mesastar.org/mailman/listinfo/mesa-users

_______________________________________________  
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/20181127/00ecf211/attachment.html>