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

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

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
>
>