[mesa-users] Power in erg/s for specific nuclear reactions (not categories!)

[Josiah Schwab]

Hi Héctor,

>   Do you know if it is possible to get the power in erg/s corresponding to
> the specific nuclear reactions in MESA?

This doesn't seem trivial to me.  As far as I know, this information is
not part of the star structure, and therefore not easy to access from
run_star_extras.

>   For the reaction categories (cno, burn_c, c12_c12), I can get the power
> without problem, which I am unable to achieve for any of the specific
> reactions (e.g. r_1212). Do you know how to do it? I thought this would be
> the best way to check that the reactions are actually implemented by MESA.
> Is there any guarantee that everything is OK with the reactions?

That being said, what exactly do you want to achieve?  I'm not sure I
understand what "actually implemented" means.

If you're worried about checking the details of specific rates, I think
the way to do that is /not/ to print out reaction rates from a full star
MESA calculation, but rather to call the rates (or net) module directly.

There are examples of how to do this sort of thing in the test
directories in the rates (or net) modules.  In particular, the file

  $MESA_DIR/rates/test/src/test_rates.f

may be of interest to you.

Josiah