[mesa-users] net_iso?

[jingluan at caltech.edu]

Hello :-)

I am confused with the returned values of s% net_iso(s% chem_id(j))
(j=1,...,s% species);

I ask mesa to print out these information:

'chem_check,j'  j  s% chem_id(j), s% net_iso(s% chem_id(j)), s% xa(j,1)
 chem_check,j   1      2             1                        0.71
 chem_check,j   2      3             2                        2E-005
 chem_check,j   3      5             3                        2E-005
 chem_check,j   4      6             4                        0.27
 chem_check,j   5      13            5                        1.25E-008
 chem_check,j   6      18            6                        1.E-099
 chem_check,j   7      26            7                        1.E-099
 chem_check,j   8      38            8                        4.1E-003
 chem_check,j   9      51            9                        1.2E-003
 chem_check,j   10     69            10                       1.1E-002
 chem_check,j   11     114           11                       2.5E-003
 chem_check,j   12     168           12                       9.5E-004

I am confused by column 3: s% chem_id(j). From the initial abundance
setting, I can figure out chem_id(1)=2 stands for h1, chem_id(2)=3 stands
for h2, chem_id(3)=5 for he3, chem_id(4)=6 for he4. But how about
chem_id=13,18,26,38,51,69,114,168? What isotopes do they stand for please?
I know that I should refer to mesa/chem/public/chem_def.f, but I could not
figure out which part of this file explains the correspondence between
isotope name and the value of chem_id(j).

Thank you very much!!

Sincerely,
Jing