[mesa-users] Reaction Rates vs Radius

Bill Paxton paxton at kitp.ucsb.edu
Fri Jul 31 15:15:16 EDT 2015


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!

-------------- next part --------------
A non-text attachment was scrubbed...
Name: run_star_extras.f
Type: application/octet-stream
Size: 17889 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150731/9a886df7/attachment.obj>
-------------- next part --------------


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

-------------- next part --------------
A non-text attachment was scrubbed...
Name: log_rate_r_c12_ag_o16.pdf
Type: application/pdf
Size: 30466 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150731/9a886df7/attachment.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log_rate_r_he4_he4_he4_to_c12.pdf
Type: application/pdf
Size: 29854 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150731/9a886df7/attachment-0001.pdf>
-------------- next part --------------

-------------- next part --------------
A non-text attachment was scrubbed...
Name: Abundances.pdf
Type: application/pdf
Size: 50022 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150731/9a886df7/attachment-0002.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Power.pdf
Type: application/pdf
Size: 55213 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150731/9a886df7/attachment-0003.pdf>
-------------- next part --------------

-------------- next part --------------
A non-text attachment was scrubbed...
Name: logT.pdf
Type: application/pdf
Size: 30528 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150731/9a886df7/attachment-0004.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: logRho.pdf
Type: application/pdf
Size: 30641 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150731/9a886df7/attachment-0005.pdf>
-------------- next part --------------




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
> 
> ------------------------------------------------------------------------------
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users



More information about the Mesa-users mailing list