[mesa-users] Using neu module

[Robert Farmer]

Hi
A few points, you'll want to have a
use const_def
at the top of your file to import the dp constant

When you get the negative value what is the value of info? If not zero then
there was an error

>First, am I correctly understanding how neu_get outputs? My understanding
is that when calling neu_get, the variable "loss" becomes populated with
>both the total energy loss rate due to thermal neutrinos, and its
derivatives. Furthermore, is loss(1) the loss rate, while loss(2) through
loss(5) are its
>derivatives? Also, is the "sources" array similarly organized, such that
sources(1,1) is the pair loss rate, sources(2, 1) is the plasmon loss rate,
etc?

Yes, sources is size (num_neu_types, num_neu_rvs) hen you can do a look up
against the constants in the neu_def file to determine what each index is

Rob


On Wed, Sep 28, 2016 at 4:23 PM, Jonathan Baker <jonat92 at vt.edu> wrote:

> Hi all,
>
> I'm trying to grab the neu module and use it separately from MESA, but I'm
> having some troubles. My question is twofold.
>
> First, am I correctly understanding how neu_get outputs? My understanding
> is that when calling neu_get, the variable "loss" becomes populated with
> both the total energy loss rate due to thermal neutrinos, and its
> derivatives. Furthermore, is loss(1) the loss rate, while loss(2) through
> loss(5) are its derivatives? Also, is the "sources" array similarly
> organized, such that sources(1,1) is the pair loss rate, sources(2, 1) is
> the plasmon loss rate, etc?
>
> Second, assuming I'm correct in understanding the outputs, for some values
> of T and Rho, I'm getting negative loss rate due to bremsstrahlung
> neutrinos. For example, if I input T=23095520474 and Rho=20092329582, I get
> a brems loss rate of -79390045138215.922 ergs/g/s. This doesn't seem
> correct. Am I using the module outside normal operating conditions? Am I
> misunderstanding how the neutrinos subroutine in mod_neu works?
>
> I've attached a copy of the sample/src/sample.f I've used to produce the
> above result.
>
> Thanks in advance,
> Jonathan Baker
>
> ------------------------------------------------------------
> ------------------
>
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160928/c5a5b105/attachment.html>