[Mesa-users] Question about scaling carbon rates

Fri Jul 29 16:47:20 UTC 2022

shut down c12(a,g) in approx-like networks. 

of course, this depends on what one means by "C-burning processes" 
and what mesa means by "burn_c", which is in data/rates_data/reactions.list

fxt

> On Jul 29, 2022, at 9:34 AM, Jeremy Sakstein <sakstein at hawaii.edu> wrote:
> 
> Hi,
> 
> This worked as expected but I still see some C burning.
> 
> If I compare burn_c with c12_c12 and c12_o16 then burn_c shows some burning but not the latter two.
> 
> Suppose I wanted to turn off all C-burning processes in approx21, which rates should I scale?
> 
> Cheers,
> 
> Jeremy
> 
> On Thu, Jul 28, 2022 at 11:50 PM Jeremy Sakstein <sakstein at hawaii.edu> wrote:
> Thanks Rob,
> 
> Changing the initial net worked.
> 
> Cheers,
> 
> Jeremy
> 
> On Thu, Jul 28, 2022 at 11:28 PM Rob Farmer <robert.j.farmer37 at gmail.com> wrote:
> You'll need to do a restart or manipulate the s% rate_factors array directly in your run_star_extras (See set_rate_factors in star/job/run_star_support.f90)
> 
> Or just use approx21 from the beginning and not worry about this.
> 
> Rob
> 
> On Fri, 29 Jul 2022 at 11:14, Jeremy Sakstein <sakstein at hawaii.edu> wrote:
> Thanks Rob,
> 
> So if I have the rate scaled at the start of the run in my inlist then will that scaling take place once approx21 turns on automatically later or do I need to stop and restart with a new inlist to scale the rate?
> 
> Cheers,
> 
> Jeremy
> 
> On Thu, Jul 28, 2022 at 11:10 PM Rob Farmer <robert.j.farmer37 at gmail.com> wrote:
> HI,
> 
> This is where inlists become helpful :)
> 
> >Is this because approx21 is not turned on until later?
> 
> I assume then you're not expliciting changing your network then? Then yes, you'll pick up basic.net first and then approx later, and basic.net doesn't do carbon burning.
> 
> Rob
> 
> 
> On Fri, 29 Jul 2022 at 10:20, Jeremy Sakstein <sakstein at hawaii.edu> wrote:
> Thanks Rob,
> 
> I also tried 'r1212' after I sent the email and this also did not produce the message at the start of the run
> 
> Is this because approx21 is not turned on until later?
> 
> Cheers,
> 
> Jeremy
> 
> On Thu, Jul 28, 2022 at 10:18 PM Rob Farmer <robert.j.farmer37 at gmail.com> wrote:
> Hi 
> Which rate you need depends on your choice of nuclear network. Use show_net_reactions_info = .true. in star_job to see the full detailed list.
> 
> I'm guessing you are using an approx net? in which case you want r1212
> 
> > rates/private/rates_names.f90
> 
> I'm glad you're looking into the code, unfortunately many rates aren't listed in the code explicitly as they get generated at runtime so you'll need to use show_net_reactions_info. Only reactions that need some sort of special handling are listed in the code.
> 
> Rob
> 
> 
> On Fri, 29 Jul 2022 at 03:11, Jeremy Sakstein <sakstein at hawaii.edu> wrote:
> Hi everyone,
> 
> I would like to scale the C12_C12 reactions i.e., the carbon burning processes but can't seem to find the correct rate.
> 
> I tried r_c12_c12 and didn't get the notification in the terminal output. Same thing with rc12c12.
> 
> It occurred to me that I may have to scale each rate individually so I looked through rates/private/rates_names.f90 (I'm using version 15140) and found two of the reactions so now I have:
> 
>     num_special_rate_factors=2
>     reaction_for_special_factor(1)='r_c12_c12_to_h1_na23'
>     special_rate_factor(1)=1d-10
>     reaction_for_special_factor(2)='r_c12_c12_to_he4_ne20'
>     special_rate_factor(2)=1d-10
> 
> in my inlist but this still does not give the message saying that a special rate is being used.
> 
> Is someone able to tell me the correct rates to scale to turn off C burning? 
> 
> Cheers,
> 
> Jeremy
> 
> 
> 
> 
> --
> Jeremy Sakstein (He/Him/His)
> Assistant Professor of Theoretical Physics
> Department of Physics & Astronomy
> University of Hawaiʻi at Mānoa
> Watanabe Hall
> 2505 Correa Road
> Honolulu, HI, 96822, USA
> http://www.jeremysakstein.com
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 
> 
> 
> -- 
> Jeremy Sakstein (He/Him/His)
> Assistant Professor of Theoretical Physics
> Department of Physics & Astronomy
> University of Hawaiʻi at Mānoa
> Watanabe Hall
> 2505 Correa Road
> Honolulu, HI, 96822, USA
> http://www.jeremysakstein.com
> 
> 
> -- 
> Jeremy Sakstein (He/Him/His)
> Assistant Professor of Theoretical Physics
> Department of Physics & Astronomy
> University of Hawaiʻi at Mānoa
> Watanabe Hall
> 2505 Correa Road
> Honolulu, HI, 96822, USA
> http://www.jeremysakstein.com
> 
> 
> -- 
> Jeremy Sakstein (He/Him/His)
> Assistant Professor of Theoretical Physics
> Department of Physics & Astronomy
> University of Hawaiʻi at Mānoa
> Watanabe Hall
> 2505 Correa Road
> Honolulu, HI, 96822, USA
> http://www.jeremysakstein.com
> 
> 
> -- 
> Jeremy Sakstein (He/Him/His)
> Assistant Professor of Theoretical Physics
> Department of Physics & Astronomy
> University of Hawaiʻi at Mānoa
> Watanabe Hall
> 2505 Correa Road
> Honolulu, HI, 96822, USA
> http://www.jeremysakstein.com
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 