Adding RecAtom Line Labels

[Nick Sterling]

Hello,

I am trying to compute abundances from optical recombination lines with PyNeb v1.1.17, and have a question about incomplete line label lists.

For example, if I define O3r as

DataFileDict = {'O3r': {'rec': 'o_iii_rec_P91.func'}}

O3r = pn.RecAtom('O', 3)
 
then print(O3r.labels) returns 

['3265+' '3762+' '4363.' '5007.']

However this does not the label '5592+', which is in the o_iii_rec_P91.func file.

Is it possible to add labels when for permitted lines in the ion_rec_ref.func file, if they are not part of the preset labels?  If so, how can I do this? I found that some line labels are defined in init.py in pyneb/utils and found label definitions for O1r, O2r, N2r, N3r, and C2r, but not for other CNO ions with permitted lines. Are those defined elsewhere?

Thank you!

Nick Sterling

[Christophe Morisset]

Hi Nick,

The default case is B for recombination lines, and for this case no data is available for the 5592 line. You should use A-case to allow it:

O3r = pn.RecAtom('O', 3, case='A')

O3r.labels

Out[6]: array(['4435+', '3265+', '3762+', '5592+', '4363.', '5007.'], dtype='<U5') 

Saludos,

Ch.