[mesa-users] ine20

[Jing L. <jin...@ca...>]

Hellow,

I type 'use chem_def' in run_star_extras.f. And the isotope number 
returned by chem_def for Ne20 is ine20=114; But the size of s% chem_id 
is only 9. so when I used s% chem_id(ine20), it reports error. I guess 
the 9 isotopes included in the current nuclear network are 
H1,He3,He4,C12,N14,O16,Ne20,Mg24 and Si28; So I guess I should use s% 
xa(7,k) to give the mass fraction of Ne20 at shell k?

Actually I should be able to use s% xa(s% chem_id(ine20),k) to do the 
same thing, but it does not work. What might be wrong here please?

Thank you very much :-)

-- 
Sincerely
Jing

Ph.D candidate at physics.caltech
email: jin...@ca...
address: MC350-17,Caltech,1200 E.California Blvd
          Pasadena, CA 91125