Re: [mesa-users] how to turn off certain categories of nuclear burning

[Robert F. <rjf...@as...>]

Hi

You need to set

s%weak_rate_factor = 0

to turn off the weak rates. when i do this and put up the pgstar power win
all i see is the cno, tri_alfa and pp as expected

Rob

On Tue, May 30, 2017 at 2:46 PM, Aaron Dotter <aar...@gm...>
wrote:

> Hi Rob,
>
> Thanks!  I've tried something like this and it doesn't work.  I also tried
> your code, just to be sure, and it doesn't work either.
>
> I *think* there is some disconnect between the reaction rates in g and
> the rate factors -- or I just don't understand how the two go together.
> The two arrays are of different length, for one thing: size(s%
> rate_factors) is greater than g% num_reactions = s% num_reactions.  Not
> sure what's going on there.  Weak reactions, maybe?
>
> I'm attaching as close to a MWE as I can make.  The attached tarball can
> be unpacked into star/work.  The code is in extras_finish_step and just
> below that.  The idea in this particular case is to turn off all but H- and
> He-burning during the post-AGB.  It should take fewer than 20 steps before
> the code to turn off burning is called.  Then, if you look at the terminal
> output for LH, LHe, and LZ, you'll see that this is not what happens.  At
> least that's what I see!
>
> Aaron
>
>
>
>
>
>
>
>
>
> On Tuesday, May 30, 2017, Robert Farmer <rjf...@as...> wrote:
>
>> Hi
>>
>>  s%rate_factors goes from 1 to g%num_reactions
>>
>> so you could do:
>>
>> call get_net_ptr(s%net_handle, g, ierr)
>> s%rate_factors=0d0
>>
>>       do j=1,g %num_reactions
>>          id = g% reaction_id(j)
>>          icat = reaction_categories(id)
>>          if(icat .eq. ipp .or. ........)then
>>                 s%rate_factors(j)=1d0
>>          end if
>>      enddo
>>
>> Rob
>>
>>
>>
>> On Tue, May 30, 2017 at 1:01 PM, Aaron Dotter <aar...@gm...>
>> wrote:
>>
>>> Hi,
>>>
>>> I have a situation where I'd like to turn off all nuclear burning that
>>> is neither H-burning nor He-burning in run_star_extras.  In MESA/chem
>>> terminology, I'd like reactions in these categories
>>>
>>> ipp = 1
>>> icno = 2
>>> i3alf = 3
>>>
>>> to go on as usual but to zero all other burning categories.  I'm using a
>>> very recent version (9793 to be exact).  I had a way of doing this back in
>>> 7503 but that no longer works.
>>>
>>> I have an array s% category_factors(1:num_categories); there are
>>> currently 24 categories in chem_def.  This array is filled with values of 1
>>> (first three elements) or 0 (all others).
>>>
>>> 1. I need to associate each reaction in my net with its category.  The
>>> following seems to work:
>>>
>>>       call get_net_ptr(s%net_handle, g, ierr)
>>>       do j=1,g %num_reactions
>>>          id = g% reaction_id(j)
>>>          icat = reaction_categories(id)
>>>       enddo
>>>
>>> 2. I need to associate each reaction in my net with its s% rate_factor.
>>> This has me stumped at the moment.
>>>
>>>
>>> Any advice would be greatly appreciated!
>>>
>>>
>>> Thanks,
>>> Aaron
>>>
>>>
>>>
>>>
>>> ------------------------------------------------------------
>>> ------------------
>>> 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
>>>
>>>
>>