[Mesa-users] Question about scaling carbon rates

Jeremy Sakstein sakstein at hawaii.edu
Fri Jul 29 09:14:00 UTC 2022


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


More information about the Mesa-users mailing list