[Mesa-users] Question about scaling carbon rates

Fri Jul 29 16:54:07 UTC 2022

Thanks Frank,

This is what I was looking for!

On Fri, Jul 29, 2022 at 6:47 AM Francis Timmes <fxt44 at mac.com> wrote:

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

-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20220729/45bdcf5f/attachment.htm>