No OIII] 1661,1666 in PyNeb Download

[Berg, Danielle A.]

Hi PyNeb team,

Thanks for putting together such a useful python version of all of these nebular tools!

Currently I am having one problem: I don’t seem to have access to the OIII] 1661,1666 lines.

For example: 

AST-NC165017:~ berg.249$ pip install --user PyNeb

Collecting PyNeb

  Downloading https://files.pythonhosted.org/packages/27/86/520aa037dadd1e4cf35bdc2a5357b88b0969a988c2c9fe60fdb9368a14d3/PyNeb-1.1.5.tar.gz (21.3MB)

    100% |████████████████████████████████| 21.4MB 78kB/s 

Building wheels for collected packages: PyNeb

  Running setup.py bdist_wheel for PyNeb ... done

  Stored in directory: /Users/berg.249/Library/Caches/pip/wheels/63/3b/96/0aa5b647d7109cf9aeff63244c7ff98c63a49608599bc08d84

Successfully built PyNeb

Installing collected packages: PyNeb

Successfully installed PyNeb-1.1.5

AST-NC165017:~ berg.249$ python

Python 3.6.4 |Anaconda custom (64-bit)| (default, Dec 21 2017, 15:39:08) 

[GCC 4.2.1 Compatible Clang 4.0.1 (tags/RELEASE_401/final)] on darwin

Type "help", "copyright", "credits" or "license" for more information.

>>> import pyneb as pn

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

>>> O3.lineList

array([ 883322.33863781,  326522.35505316,  518003.77007944,

          4931.22565733,    4958.9097225 ,    5006.84166184,

          2314.86703288,    2320.95013678,    2331.3973693 ,

          4363.20876917])

And so I have not been able to calculate temperatures from 1666/5007 or determine UV C/O abundances.

I also have this issue on my laptop, which is running Python2, so it does not seem to be a Python version issue.

I’ve also used the pn.atomicData.setDataFile command to check try other atomic data sets, but that does not help either.

I’d greatly appreciate your help with this issue so that I can use your code for my upcoming C/O paper!

Thank you in advance for your time.

Cheers,

Danielle

Danielle A. Berg, Ph.D. 
Postdoctoral Research Associate 
4061McPherson Laboratory
140 West 18th Avenue
Columbus, OH 43210
(P) 320.267.4014
(E) berg...@osu.edu

(w) www.danielleaberg.com 

[Christophe Morisset]

Hi Danielle,

The atomic data used by default for O++ are from the Storey et al 2014 paper, where no data are given for the level 6.
If you need to use the 1660, 1666 lines coming from this level, you need to change the atomic data to for example AK99 by pn.atomicData.setDataFile('o_iii_coll_AK99.dat')
Then you will have access to the UV lines:

In [4]: pn.atomicData.setDataFile('o_iii_coll_AK99.dat')

In [5]: o3 = pn.Atom('O',3)

In [6]: o3.lineList
Out[6]:

array([ 883322.33863781,  326522.35505316,  518003.77007944,
          4931.22565733,    4958.9097225 ,    5006.84166184,
          2314.86703288,    2320.95013678,    2331.3973693 ,

          4363.20876917,    1657.69329657,    1660.80921401,
          1666.14971595,    2496.03107   ,    5833.0097649 ])

Cheers,
Christophe

[Danielle Berg]

Thank you!