[mesa-users] Reaction Rates vs Radius

[Francis Timmes]

wow, thanks bill! 

interesting to me is the access through a do1_net.

fxt



> On Jul 31, 2015, at 12:15 PM, Bill Paxton <paxton at kitp.ucsb.edu> wrote:
> 
> 
> Here is a run_star_extras with some (non-trivial) code to show selected reaction rates in extra profile columns.  The code is perhaps beyond what the typical user would want to tackle, but it isn't too horrible and it uses only routines from public interfaces -- no cheating!
> 
> <run_star_extras.f>
> 
> Hopefully it will be clear how to modify this to show multiple rates of your choice.
> 
> These are plots from the 1st model of the agb test case showing screened rates for r_c12_ag_o16 and triple alphs (from basic.net).
> 
> -B
> 
> <log_rate_r_c12_ag_o16.pdf><log_rate_r_he4_he4_he4_to_c12.pdf>
> <Abundances.pdf><Power.pdf>
> <logT.pdf><logRho.pdf>
> 
> 
> 
> On Jul 31, 2015, at 1:15 AM, Josiah Schwab wrote:
> 
>> Hi Joe,
>> 
>>> What’s the best way to get reaction rates for various reactions versus
>>> radius?
>> 
>> The structure of the net module means that it is not easy to see the
>> value of a single reaction rate throughout the entire star.  See, for
>> example, this recent thread
>> 
>> https://sourceforge.net/p/mesa/mailman/message/33611293/
>> 
>> The profile outputs are listed in
>> 
>> star/defaults/profile_columns.list.
>> 
>> and you can use run_star_extras.f to output anything in
>> 
>> star/public/star_data.inc
>> 
>> which includes things like dxdt_nuc throughout the star.
>> 
>> http://mesa.sourceforge.net/run_star_extras.html#toc-1-3
>> 
>> Best,
>> Josiah
>>