[mesa-users] how to turn off certain categories of nuclear burning
aaron.dotter at gmail.com
Tue May 30 16:01:04 EDT 2017
I have a situation where I'd like to turn off all nuclear burning that is
neither H-burning nor He-burning in run_star_extras. In MESA/chem
terminology, I'd like reactions in these categories
ipp = 1
icno = 2
i3alf = 3
to go on as usual but to zero all other burning categories. I'm using a
very recent version (9793 to be exact). I had a way of doing this back in
7503 but that no longer works.
I have an array s% category_factors(1:num_categories); there are currently
24 categories in chem_def. This array is filled with values of 1 (first
three elements) or 0 (all others).
1. I need to associate each reaction in my net with its category. The
following seems to work:
call get_net_ptr(s%net_handle, g, ierr)
do j=1,g %num_reactions
id = g% reaction_id(j)
icat = reaction_categories(id)
2. I need to associate each reaction in my net with its s% rate_factor.
This has me stumped at the moment.
Any advice would be greatly appreciated!
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Mesa-users