[Mesa-users] Question about neutrino loss rates

[Rob Farmer]

Before, you want to make sure the flags are set before calling (and that
you pass the flags array as an argument to neu_get) as it dictates which
neutrino rates are computed.

Rob


On Tue, 1 Sep 2020 at 09:58, Jeremy Sakstein <sakstein at hawaii.edu> wrote:

> Thanks Rob.
>
> Do I want to use the flags(plas_neu_type)=.true. when I call new get or
> before or after?
>
> On Mon, Aug 31, 2020 at 21:41 Rob Farmer <r.j.farmer at uva.nl> wrote:
>
>> Hi,
>>
>> >Hi everyone, just a follow up on this. Going through
>>
>> neu/private/mod_def.f90 it looks like I actually want
>>
>> sources(plas_neu_type, ineu) for the plasma loss rate.
>>
>> Yes, you want "ineu" to get the rate, the other elements in the second
>> column are the derivatives.
>>
>> >If so, will the units be ergs/g/s?
>>
>> As per the comment in neu_def.f90:
>>
>> integer, parameter :: ineu        = 1     ! loss rate (nonnegative) in
>> units of ergs / gram / second
>>
>> >after calling neu_get appropriately.
>>
>> As long as you pass the flags array with:
>>
>> flags(plas_neu_type)=.true.
>>
>>
>> Rob
>>
>>
>> On Tue, 1 Sep 2020 at 05:38, Jeremy Sakstein <sakstein at hawaii.edu> wrote:
>>
>>> Hi everyone, just a follow up on this. Going through
>>> neu/private/mod_def.f90 it looks like I actually want
>>> sources(plas_neu_type, ineu) for the plasma loss rate.
>>>
>>> it would be great if someone could confirm this for me.
>>>
>>> Cheers,
>>>
>>> Jeremy
>>>
>>> On Mon, Aug 31, 2020 at 3:12 PM Jeremy Sakstein <sakstein at hawaii.edu>
>>> wrote:
>>>
>>>> Hi everyone,
>>>>
>>>> I'm trying to find just the loss rate into plasma neutrinos. I think I
>>>> have this figured out but I would like to double check. Reading the neu
>>>> module, it looks like I can call:
>>>>
>>>> sources(plas_neu_type,num_neu_rvs)
>>>>
>>>>
>>>> after calling neu_get appropriately.
>>>>
>>>>
>>>> Is this correct?
>>>>
>>>> If so, will the units be ergs/g/s?
>>>>
>>>> Cheers,
>>>>
>>>> Jeremy
>>>>
>>>>
>>>>
>>>
>>>
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200901/43e4c6d0/attachment.htm>