Ionic abundances using Pyneb

[Celso Junior]

Hi,

I'm starting to work with PyNeb. Thank you for this wonderful code. 

I know that we can easily derive oxygen ionic abundances in relation to hydrogen, using something like 

O3 = pn.Atom('O',3)

O3a = O3.getIonAbundance(int_ratio=6.95, tem=13879, den=1153, wave=5007, Hbeta=1.0)

However, I need to estimate the carbon abundance (CIV 1550 and CIII 1909) in relation to He. Is it possible? 

Best regards, 

CELSO B.O. JUNIOR

[Christophe Morisset]

Hi Celso,

The getIonAbundance method only works with hydrogen line normalisation. 

But you can always obtain ionic abundance ratios from observed line ratios and line emissivities ratios, using the relation between I(line) and n_e, n_ion and Emissivity(line).

Saludos,

Ch.