[Mesa-users] Detect zones with H fusion

[Josiah Schwab]

Hi Roque,

> eps_nuc in the profile file allows you to determine the zones where
> there is energy generation from nuclear reactions.

And if you need to get even more specific (i.e. distinguish different
types of burning) there are different "categories".

See add_reaction_categories in profile_columns.list.  These have names
like pp, cno, etc.

If you're doing this in run_star_extras, there's eps_nuc_categories.
which is a (num_categories, nz) array.  You can get the category index
constants from chem_def.

So the specific energy generation rate by hydrogen burning in cell k of
the model is

  s% eps_nuc_categories(ipp, k) + s% eps_nuc_categories(icno, k)

Hope that helps,
Josiah