mesa-users] Adding a new nuclear reaction rates in MESA
Aaron Dotter
aaron.dotter at gmail.com
Mon Nov 8 14:39:15 EST 2010
Hi Zahra,
Yes, you can test the rates being used by the rates module. Bill and I have
looked at this and the easiest way to check the rate is to edit the rates
module test code in mesa/rates/test/src/test_rates.f. There is a subroutine
called "test1" that prints out the rate of one reaction for one temperature,
though you could easily modify this to do a loop over temperature.
At the very bottom of test_rates.f is the line
!call test1; stop
Uncomment that line and then modify subroutine test1 to use the rate and
temperature(s) that you want. Change the rate by setting ir to the rate ID
you've added to rate_list.txt.
Then do a ./mk and ./rn in mesa/rates/test/. Simply change the entry in
rate_list.txt and then run again to see the difference.
Since mesa/star uses mesa/rates, once you're sure that it's working as you
expect in the rates test code, it should also be working as you expect in
mesa/star.
Aaron
On Mon, Nov 8, 2010 at 3:14 AM, zahra altaha motahar <zami7965 at yahoo.com>wrote:
> Hi,
>
> Recently, I have added my new nuclear reaction rates in the MESA code. I
> produce a table of temperatures and rates and save it in the
> net_data/rates/directory. Then, I add the filename and reaction name to the
> "rate_list.txt" and run the code.
> How can I make sure that the MESA code uses our added new nuclear reaction
> rates instead of the NACRE rates for its modeling?
>
> Thank you,
> Zahra
>
>
> ------------------------------------------------------------------------------
> The Next 800 Companies to Lead America's Growth: New Video Whitepaper
> David G. Thomson, author of the best-selling book "Blueprint to a
> Billion" shares his insights and actions to help propel your
> business during the next growth cycle. Listen Now!
> http://p.sf.net/sfu/SAP-dev2dev
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20101108/647facb0/attachment.html>
More information about the Mesa-users
mailing list