[mesa-users] how to turn off certain categories of nuclear burning

[Aaron D. <aar...@gm...>]

Hi,

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)
      enddo

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!


Thanks,
Aaron