[Mesa-users] Question about scaling carbon rates

[Jeremy Sakstein]

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