Re: [mesa-users] use mass fractions in neu.f

[Giannotti, M. <MGi...@ma...>]

Hi Bill,

Sorry, I still have some problems. I need to use z2bar in neu/private/neu.f. I tried to make the changes that you recommended and they worked (no compilation error). However, when I try to use z2bar in neu.f (in the subroutine "neutrinos") I have a compilation error since, it seems, that subroutine does not know what z2bar is. Any suggestion?

Thanks again,
maurizio

Maurizio Giannotti,
Assistant Professor of Physics,
Physical Sciences, Barry University,
11300 NE 2nd Ave., Miami Shores, FL 33161
mgi...@ma...<mailto:mgi...@ma...>
305-899-4565 (office), 305-899-3479 (fax)

________________________________
From: Bill Paxton [pa...@ki...]
Sent: Monday, September 05, 2011 4:10 PM
To: Giannotti, Maurizio
Cc: mes...@li...
Subject: Re: [mesa-users] use mass fractions in neu.f

Hi Maurizio,

We need to add z2bar to the argument list of neu_get for you.
Then we need to change star to pass that new arg.
I'll change mesa to do that, but if you'd like to try it with
your current version, you can make 2 small changes:

first, in neu/public/neu_lib, add z2bar as another arg for neu_get.

old:
      subroutine neu_get(T, log10_T, Rho, log10_Rho, abar, zbar, log10_Tlim, flags,
     >         loss, sources, info)

new:
subroutine neu_get(T, log10_T, Rho, log10_Rho, abar, zbar, z2bar, log10_Tlim, flags,
     >         loss, sources, info)
.....
         double precision, intent(in) :: zbar2 ! mean charge squared

then
cd neu/test
./mk
./ck
./export

then in star/private/micro.f

old:
         call neu_get(s% T(k), log10_T, s% rho(k), log10_rho, s% abar(k), s% zbar(k), &
                        log10_Tlim, flags, loss, sources, ierr)
new (with z2bar):
         call neu_get(s% T(k), log10_T, s% rho(k), log10_rho, s% abar(k), s% zbar(k), s% z2bar(k), &
                        log10_Tlim, flags, loss, sources, ierr)

then
cd star/test
./mk
./ck
./export



Let me know how it goes.

Cheers,
Bill







On Sep 5, 2011, at 11:05 AM, Giannotti, Maurizio wrote:

Hello,

I am trying to add a few lines to neu.f (in mesa/neu/private) and I need to know the values of the hydrogen, helium, carbon and oxygen  mass fraction. More specifically, I need to calculate

y_e+\sum z_j^2 y_j

where y_e and y_j are the number fractions per baryon of electrons and nuclei. From neu.f I can access directly only y_e.

What should I do (or where should i look)?

Thank you very much,
Maurizio

Maurizio Giannotti,
Assistant Professor of Physics,
Physical Sciences, Barry University,
11300 NE 2nd Ave., Miami Shores, FL 33161
mgi...@ma...<mailto:mgi...@ma...>
305-899-4565 (office), 305-899-3479 (fax)

------------------------------------------------------------------------------
Special Offer -- Download ArcSight Logger for FREE!
Finally, a world-class log management solution at an even better
price-free! And you'll get a free "Love Thy Logs" t-shirt when you
download Logger. Secure your free ArcSight Logger TODAY!
http://p.sf.net/sfu/arcsisghtdev2dev_______________________________________________
mesa-users mailing list
mes...@li...<mailto:mes...@li...>
https://lists.sourceforge.net/lists/listinfo/mesa-users