[Mesa-users] Change abundances of C12 and O16 after He-burning stage

[Athanasiadou, S. (Sofia)]

Hello everybody,

I am trying to change the central abundances of C12 and O16 after central He-burning stage, customizing the reaction rates "c12ag" and "3a" from star_job defaults as it follows:

    num_special_rate_factors = 1
    reaction_for_special_factor(1) = 'c12ag'
    special_rate_factor(1) = 0.1

The problem is that when I plot these abundances in pyplot, nothing changes really (the path to history_data is correct I double checked it!). What should I do?

Extra info:

The initial mass of the star is 10Msolar, without convection overshooting and mass loss. Stopping condition:

    xa_central_lower_limit_species(1) = 'he4'
    xa_central_lower_limit(1) = 1d-4

Thank you very much,

Sofia.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20171214/eadd1da4/attachment.html>