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

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

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)