[Mesa-users] Question about neutrino loss rates

[Jeremy Sakstein]

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/20200831/3dd31530/attachment.htm>