Re: [mesa-users] which file creates the reaction rates binary files?

[李二涛 <le...@sz...>]

 Dear Rob,

I got it. Thank you!




Er-Tao Li(李二涛)
Institute of Nuclear Technology，College of Physics & Energy, Shenzhen University
Nanhai Ave 3688, Shenzhen, Guangdong, P.R.China, 518060 
Tel: 86-755-26536232. Fax: 86-755-26534374 

From: Robert Farmer
Date: 2017-04-01 02:09
To: let
CC: mesa-users
Subject: Re: Re: [mesa-users] which file creates the reaction rates binary files?
Hi

So firstly the easiest way to see the exact rate mesa used is to do:


cd $MESA_DIR/rates/test

./show_rates ../../data/rates_data/cache/r_he3_ag_be7_1.bin



This outputs as plain text the rate mesa uses.


As to where mesa gets the rate from, by default MESA uses NACRE rates and falls back on JINA rates if NACRE is missing them/outside of JINAs temperature range



You need to set:
      set_rates_preference = .true.
      new_rates_preference = 2


in you star_job inlist to get only the JINA rates (or new_rates_preference = 1 to get only NACRE). So in this case i think we are using the rate from the  subroutine rate_he3he4_nacre in $MESA_DIR/rates/private/ratelib.f90


You can of course always set your own rates (see the custom_rate_c14ago18 test suite case as an example) for complete control over what MESA is doing



Rob



On Thu, Mar 30, 2017 at 7:22 PM, 李二涛 <le...@sz...> wrote:



 Dear Rob,

What I really want to know is the exact reaction rates expression which is used to create the binary file. I noticed that there are several sets of data for a reaction in mesa/data/rates_data/jina_reaclib_results_20130213default2. Take 3He(a,g)7Be reaction for example, there are four sets of data about the reaction:
4
       he4  he3  be7                       cd08n     1.58610e+00          
 1.770750e+01 0.000000e+00-1.282710e+01-3.812600e+00                      
 9.422850e-02-3.010180e-03 1.333330e+00                                   
4
       he4  he3  be7                       cd08n     1.58610e+00          
 1.560990e+01 0.000000e+00-1.282710e+01-3.082250e-02                      
-6.546850e-01 8.963310e-02-6.666670e-01  
2
       be7  he4  he3                       cd08nv   -1.58610e+00          
 4.083550e+01-1.840590e+01-1.282710e+01-3.812600e+00                      
 9.422850e-02-3.010180e-03 2.833330e+00                                   
2
       be7  he4  he3                       cd08nv   -1.58610e+00          
 3.873790e+01-1.840590e+01-1.282710e+01-3.082250e-02                      
-6.546850e-01 8.963310e-02 8.333330e-01 
I think that one of the first two sets of data is used to create the 3He(a,g)7Be reaction rates, and one of the last two sets of data is used to create the reverse reaction rates. But after I put the numbers into the formula: rate = Exp[ a0 + a1/T9 + a2/T91/3 + a3*T91/3 + a4*T9 + a5*T95/3 + a6*log(T9) ], I can't get the same results with the binary file. So I am wondering how can I get the exact reaction rates expression which is used by the MESA.

Thanks.




Er-Tao Li(李二涛)
Institute of Nuclear Technology，College of Physics & Energy, Shenzhen University
Nanhai Ave 3688, Shenzhen, Guangdong, P.R.China, 518060 
Tel: 86-755-26536232. Fax: 86-755-26534374 

From: Robert Farmer
Date: 2017-03-29 11:42
To: let
CC: mesa-users
Subject: Re: [mesa-users] which file creates the reaction rates binary files?
Hi


There is a binary file version of the jina rates stored in $MESA_DIR/data/rates_data/cache/jina_reaclib.bin which we use after the first run so once the text file is gone we use that (if you didn't add an empty file you would get an error message about the text file missing)


Deleting these files seems an unusual thing to do, if you explain what you want to do with the rates, we might be able to suggest a better way to do things.


Rob



On Tue, Mar 28, 2017 at 7:44 PM, 李二涛 <le...@sz...> wrote:



 Dear MESA users,

I know that the reaction rates are stored in data/rates_data/cache as binary files. From the star/defaults/star_job.defaults file,  I also know that I am using jina_reaclib_results_20130213default2 database. So, I suppose that the binary files are created by the  jina_reaclib_results_20130213default2 file. But, after I deleted a reaction binary file and  replaced the jina_reaclib_results_20130213default2 with an empty file(just have the same file name:jina_reaclib_results_20130213default2), to my surprise, the MESA still able to create a new binary file and run successfully.
  
Now, I am wondering which file creates the reaction rates binary files when running the MESA?

Thanks in advance for any help that you can offer!




Er-Tao Li(李二涛)
Institute of Nuclear Technology，College of Physics & Energy, Shenzhen University
Nanhai Ave 3688, Shenzhen, Guangdong, P.R.China, 518060 
Tel: 86-755-26536232. Fax: 86-755-26534374 
 
 

 


------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
mesa-users mailing list
mes...@li...
https://lists.sourceforge.net/lists/listinfo/mesa-users