[mesa-users] reaction rate of deuterium fusion

[Masanobu Kunitomo]

Dear Frank and Josiah,

Thank you very much for your help!


With "list_net_reactions = .true.", I found that the name of the deuterium burning is “r_h2_pg_he3”. However, this seems not to be related to the energy production rate of deuterium burning.


I found “eps_h”, which is the energy production rate by hydrogen burning defined in “mesa/star/private/report.f” as 
            eps_h(k) = s% eps_nuc_categories(i_rate,ipp,k) + &
                       s% eps_nuc_categories(i_rate,icno,k).

I guess "eps_nuc_categories(ipp)” is the energy production rate by pp chain reactions. If “eps_nuc_categories” only for deuterium burning is defined, I would like to plot it, but I could not find it.


>Josiah
If you do not mind, could you please tell me how to add an "other_net” routine for output the energy production rate of a specific reaction? Is it included in the subroutine “extras_controls”?


Best wishes,
Masanobu



> On Mar 20, 2015, at 00:41, Josiah Schwab <jwschwab at berkeley.edu> wrote:
> 
> Hi Frank,
> 
>> perhaps one way to get started is to list the reaction rate names via
>> 
>> list_net_reactions = .true.
>> 
>> http://mesa.sourceforge.net/star_job_defaults.html#output
>> 
>> in your &star_job namelist.
> 
> Can you get at a single rate from inside run_star_extras.f?
> 
> I wanted to do something like this a while ago.  I ended up patching
> MESA and adding an "other_net" routine.^ Otherwise, I didn't seem like I
> could get deep enough.
> 
> Josiah
> 
> ^ It's on my TODO list to rebase this and commit it upstream.