{"id":32174859,"url":"https://github.com/equib/pyequib","last_synced_at":"2026-02-21T12:01:29.572Z","repository":{"id":60722569,"uuid":"121817983","full_name":"equib/pyEQUIB","owner":"equib","description":"pyEQUIB - Python Package for Plasma Diagnostics and Abundance Analysis","archived":false,"fork":false,"pushed_at":"2025-01-26T15:55:54.000Z","size":4353,"stargazers_count":9,"open_issues_count":0,"forks_count":2,"subscribers_count":1,"default_branch":"master","last_synced_at":"2026-01-04T11:11:03.407Z","etag":null,"topics":["atomic-level-populations","chemical-elements","line-emissivities","pyequib","python","statistical-equilibrium"],"latest_commit_sha":null,"homepage":"https://equib.github.io/pyEQUIB","language":"Python","has_issues":true,"has_wiki":null,"has_pages":null,"mirror_url":null,"source_name":null,"license":"gpl-3.0","status":null,"scm":"git","pull_requests_enabled":true,"icon_url":"https://github.com/equib.png","metadata":{"files":{"readme":"README.rst","changelog":null,"contributing":"CONTRIBUTING.md","funding":null,"license":"LICENSE","code_of_conduct":null,"threat_model":null,"audit":null,"citation":null,"codeowners":null,"security":null,"support":null,"governance":null,"roadmap":null,"authors":"AUTHORS.rst","dei":null,"publiccode":null,"codemeta":null,"zenodo":null}},"created_at":"2018-02-17T01:23:06.000Z","updated_at":"2025-04-02T22:18:31.000Z","dependencies_parsed_at":"2023-12-22T19:32:17.271Z","dependency_job_id":"797eed1a-6c12-4581-9e67-1d5aea70733e","html_url":"https://github.com/equib/pyEQUIB","commit_stats":{"total_commits":166,"total_committers":1,"mean_commits":166.0,"dds":0.0,"last_synced_commit":"40a6cecc6c76c338812816c4b20111bf5cc5217e"},"previous_names":[],"tags_count":6,"template":false,"template_full_name":null,"purl":"pkg:github/equib/pyEQUIB","repository_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/equib%2FpyEQUIB","tags_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/equib%2FpyEQUIB/tags","releases_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/equib%2FpyEQUIB/releases","manifests_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/equib%2FpyEQUIB/manifests","owner_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners/equib","download_url":"https://codeload.github.com/equib/pyEQUIB/tar.gz/refs/heads/master","sbom_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/equib%2FpyEQUIB/sbom","scorecard":null,"host":{"name":"GitHub","url":"https://github.com","kind":"github","repositories_count":286080680,"owners_count":29680147,"icon_url":"https://github.com/github.png","version":null,"created_at":"2022-05-30T11:31:42.601Z","updated_at":"2026-02-21T11:29:27.227Z","status":"ssl_error","status_checked_at":"2026-02-21T11:29:20.292Z","response_time":107,"last_error":"SSL_connect returned=1 errno=0 peeraddr=140.82.121.5:443 state=error: unexpected eof while reading","robots_txt_status":"success","robots_txt_updated_at":"2025-07-24T06:49:26.215Z","robots_txt_url":"https://github.com/robots.txt","online":false,"can_crawl_api":true,"host_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub","repositories_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories","repository_names_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repository_names","owners_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners"}},"keywords":["atomic-level-populations","chemical-elements","line-emissivities","pyequib","python","statistical-equilibrium"],"created_at":"2025-10-21T19:12:20.952Z","updated_at":"2026-02-21T12:01:29.566Z","avatar_url":"https://github.com/equib.png","language":"Python","funding_links":[],"categories":[],"sub_categories":[],"readme":"======================\npyEQUIB Python Package\n======================\n\n.. image:: https://img.shields.io/pypi/v/pyequib.svg?style=flat\n    :target: https://pypi.python.org/pypi/pyequib/\n    :alt: PyPI Version\n    \n.. image:: https://app.travis-ci.com/equib/pyEQUIB.svg?branch=master\n    :target: https://app.travis-ci.com/github/equib/pyEQUIB\n    :alt: Build Status\n    \n.. image:: https://ci.appveyor.com/api/projects/status/b3gw6vgf8s0vu8nv?svg=true\n    :target: https://ci.appveyor.com/project/danehkar/pyequib\n    :alt: Build Status\n    \n.. image:: https://coveralls.io/repos/github/equib/pyEQUIB/badge.svg?branch=master\n    :target: https://coveralls.io/github/equib/pyEQUIB?branch=master\n    :alt: Coverage Status\n    \n.. image:: https://img.shields.io/badge/license-GPL-blue.svg\n    :target: https://github.com/equib/pyEQUIB/blob/master/LICENSE\n    :alt: GitHub license\n    \n.. image:: https://img.shields.io/conda/vn/conda-forge/pyequib.svg\n    :target: https://anaconda.org/conda-forge/pyequib\n    :alt: Anaconda Cloud\n    \n.. image:: https://readthedocs.org/projects/pyequib/badge/?version=latest\n    :target: https://pyequib.readthedocs.io/en/latest/?badge=latest\n    :alt: Documentation Status\n    \n.. image:: https://img.shields.io/badge/python-2.7%2C%203.8-blue.svg\n    :alt: Support Python versions 2.7 and 3.8\n    \n.. image:: https://img.shields.io/badge/DOI-10.5281/zenodo.4287575-blue.svg\n    :target: https://doi.org/10.5281/zenodo.4287575\n    :alt: Zenodo\n    \n.. image:: http://joss.theoj.org/papers/10.21105/joss.02798/status.svg\n    :target: https://doi.org/10.21105/joss.02798\n    :alt: JOSS\n\n.. image:: https://mybinder.org/badge_logo.svg\n :target: https://mybinder.org/v2/gh/equib/pyEQUIB/HEAD?labpath=Notebooks.ipynb\n \nDescription\n===========\n\nThe **pyEQUIB** library is a collection of `Python \u003chttps://www.python.org/\u003e`_ programs developed to perform plasma diagnostics and abundance analysis using emission line fluxes measured in ionzed nebulae. It uses the `AtomNeb Python Package \u003chttps://github.com/atomneb/AtomNeb-py\u003e`_ to read collision strengths and transition probabilities for collisionally excited lines (CEL), and recombination coefficients for recombination lines (RL). This Python package can be used to determine interstellar extinctions, electron temperatures, electron densities, and ionic abundances from the measured fluxes of emission lines. It mainly contains the follwing API functions written purely in Python: \n\n* **API functions for collisionally excited lines (CEL)** have been developed based on the algorithm of the FORTRAN program `EQUIB \u003chttps://github.com/equib/EQUIB\u003e`_ originally written in FORTRAN by `Howarth \u0026 Adams (1981) \u003chttps://ui.adsabs.harvard.edu/abs/1981ucl..rept.....H/abstract\u003e`_, extended and customized by other people (`R. Clegg, D. Ruffle, X.-W. Liu, C. Pritchet, B. Ercolano, \u0026 R. Wesson \u003chttps://ui.adsabs.harvard.edu/abs/2016ascl.soft03005H/abstract\u003e`_). The program EQUIB calculates atomic level populations and line emissivities in statistical equilibrium in multi-level atoms for different physical conditions of the stratification layers where the chemical elements are ionized. Using the Python implementation of the program EQUIB, electron temperatures, electron densities, and ionic abundances are determined from the measured fluxes of collisionally excited lines.\n\n* **API functions for recombination lines (RL)** have been developed based on the algorithm of the recombination scripts by X. W. Liu and Y. Zhang from `output_mod.f90 \u003chttps://github.com/mocassin/MOCASSIN-2.0/blob/master/source/output_mod.f90\u003e`_ included in the FORTRAN program `MOCASSIN \u003chttps://github.com/mocassin/MOCASSIN-2.0\u003e`_. These API functiosn are used to determine ionic abundances from recombination lines for some heavy element ions.\n \n* **API functions for reddening and extinctions** have been developed according to the methods of the reddening law functions from `STSDAS IRAF Package \u003chttps://github.com/equib/iraf-redcorr\u003e`_, which are used to obtain interstellar extinctions and deredden measured fluxes based on different reddening laws.\n\n\nInstallation\n============\n\nDependent Python Packages\n-------------------------\n\n This package requires the following packages:\n\n    - `NumPy \u003chttps://numpy.org/\u003e`_\n    - `SciPy \u003chttps://scipy.org/scipylib/\u003e`_\n    - `AtomNeb \u003chttps://github.com/atomneb/AtomNeb-py/\u003e`_\n    \n* To get this package with the AtomNeb FITS files, you can simply use ``git`` command as follows:\n\n.. code-block::\n\n        git clone --recursive https://github.com/equib/pyEQUIB\n\nTo install the last version, all you should need to do is\n\n.. code-block::\n\n    $ python setup.py install\n\nTo install the stable version, you can use the preferred installer program (pip):\n\n.. code-block::\n\n    $ pip install pyequib\n\nor you can install it from the cross-platform package manager *conda*:\n\n.. code-block::\n\n    $ conda install -c conda-forge pyequib\n\nHow to Use\n==========\n\nThe Documentation of the Python functions provides in detail in the *API Documentation* (`equib.github.io/pyEQUIB/doc \u003chttps://equib.github.io/pyEQUIB/doc\u003e`_). \n\nSee *Jupyter Notebooks*: `Notebooks.ipynb \u003chttps://github.com/equib/pyEQUIB/blob/master/Notebooks.ipynb\u003e`_\n\nRun *Jupyter Notebooks* on `Binder \u003chttps://mybinder.org/v2/gh/equib/pyEQUIB/HEAD?labpath=Notebooks.ipynb\u003e`_:\n\n.. image:: https://mybinder.org/badge_logo.svg\n :target: https://mybinder.org/v2/gh/equib/pyEQUIB/HEAD?labpath=Notebooks.ipynb\n \nThere are three main object units:\n\n* **Collision Unit** has the API functions for plasma diagnostics and abundance analysis of collisionally excited lines. Here are some examples of using *Collision* Unit:\n\n    - *Temperature*:\n\n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_dir = os.path.join('atomic-data', 'chianti70')\n        atom_elj_file = os.path.join(base_dir,data_dir, 'AtomElj.fits')\n        atom_omij_file = os.path.join(base_dir,data_dir, 'AtomOmij.fits')\n        atom_aij_file = os.path.join(base_dir,data_dir, 'AtomAij.fits')\n        \n        atom = 's'\n        ion = 'ii'\n        s_ii_elj = atomneb.read_elj(atom_elj_file, atom, ion, level_num=5)\n        s_ii_omij = atomneb.read_omij(atom_omij_file, atom, ion)\n        s_ii_aij = atomneb.read_aij(atom_aij_file, atom, ion)\n        \n        upper_levels='1,2,1,3/'\n        lower_levels='1,5/'\n        density = np.float64(2550)\n        line_flux_ratio=np.float64(10.753)\n        temperature = pyequib.calc_temperature(line_flux_ratio=line_flux_ratio, density=density, \n                               upper_levels=upper_levels, lower_levels=lower_levels, \n                               elj_data=s_ii_elj, omij_data=s_ii_omij, aij_data=s_ii_aij)\n        print(\"Electron Temperature:\", temperature)\n\n      which gives::\n    \n        Electron Temperature:       7920.2865\n\n    - *Density*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_dir = os.path.join('atomic-data', 'chianti70')\n        atom_elj_file = os.path.join(base_dir,data_dir, 'AtomElj.fits')\n        atom_omij_file = os.path.join(base_dir,data_dir, 'AtomOmij.fits')\n        atom_aij_file = os.path.join(base_dir,data_dir, 'AtomAij.fits')\n        \n        atom = 's'\n        ion = 'ii'\n        s_ii_elj = atomneb.read_elj(atom_elj_file, atom, ion, level_num=5)\n        s_ii_omij = atomneb.read_omij(atom_omij_file, atom, ion)\n        s_ii_aij = atomneb.read_aij(atom_aij_file, atom, ion)\n        \n        upper_levels='1,2/'\n        lower_levels='1,3/'\n        temperature=np.float64(7000.0)#\n        line_flux_ratio=np.float64(1.506)#\n        density = pyequib.calc_density(line_flux_ratio=line_flux_ratio, temperature=temperature, \n                                       upper_levels=upper_levels, lower_levels=lower_levels, \n                                       elj_data=s_ii_elj, omij_data=s_ii_omij, aij_data=s_ii_aij)\n        print(\"Electron Density:\", density)\n\n      which gives::\n      \n        Electron Density:       2312.6395\n\n    - *Ionic Abundance*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_dir = os.path.join('atomic-data', 'chianti70')\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_elj_file = os.path.join(base_dir,data_dir, 'AtomElj.fits')\n        atom_omij_file = os.path.join(base_dir,data_dir, 'AtomOmij.fits')\n        atom_aij_file = os.path.join(base_dir,data_dir, 'AtomAij.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'h'\n        ion = 'ii' # H I Rec\n        hi_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        atom = 'o'\n        ion = 'iii' # [O III]\n        o_iii_elj = atomneb.read_elj(atom_elj_file, atom, ion, level_num=5) # read Energy Levels (Ej)\n        o_iii_omij = atomneb.read_omij(atom_omij_file, atom, ion) # read Collision Strengths (Omegaij)\n        o_iii_aij = atomneb.read_aij(atom_aij_file, atom, ion) # read Transition Probabilities (Aij)\n\n        levels5007='3,4/'\n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        iobs5007=np.float64(1200.0)\n        abb5007 = pyequib.calc_abundance(temperature=temperature, density=density, \n                                         line_flux=iobs5007, atomic_levels=levels5007,\n                                         elj_data=o_iii_elj, omij_data=o_iii_omij, aij_data=o_iii_aij, \n                                         h_i_aeff_data=hi_rc_data['aeff'][0])\n        print('N(O^2+)/N(H+):', abb5007)\n\n      which gives::\n      \n        N(O^2+)/N(H+):   0.00041256231 \n        \n    - *Emissivity*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_dir = os.path.join('atomic-data', 'chianti70')\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_elj_file = os.path.join(base_dir,data_dir, 'AtomElj.fits')\n        atom_omij_file = os.path.join(base_dir,data_dir, 'AtomOmij.fits')\n        atom_aij_file = os.path.join(base_dir,data_dir, 'AtomAij.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'h'\n        ion = 'ii' # H I Rec\n        hi_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        atom = 'o'\n        ion = 'iii' # [O III]\n        o_iii_elj = atomneb.read_elj(atom_elj_file, atom, ion, level_num=5) # read Energy Levels (Ej)\n        o_iii_omij = atomneb.read_omij(atom_omij_file, atom, ion) # read Collision Strengths (Omegaij)\n        o_iii_aij = atomneb.read_aij(atom_aij_file, atom, ion) # read Transition Probabilities (Aij)\n        \n        levels5007='3,4/'\n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        iobs5007=np.float64(1200.0)\n        emis = pyequib.calc_emissivity(temperature=temperature, density=density, atomic_levels=levels5007, \n                                       elj_data=o_iii_elj, omij_data=o_iii_omij, aij_data=o_iii_aij)\n        print('Emissivity(O III 5007):', emis)\n\n      which gives::\n      \n        Emissivity(O III 5007):   3.6041012e-21\n        \n\n    - *Atomic Level Population*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_dir = os.path.join('atomic-data', 'chianti70')\n        atom_elj_file = os.path.join(base_dir,data_dir, 'AtomElj.fits')\n        atom_omij_file = os.path.join(base_dir,data_dir, 'AtomOmij.fits')\n        atom_aij_file = os.path.join(base_dir,data_dir, 'AtomAij.fits')\n        \n        atom = 's'\n        ion = 'ii'\n        s_ii_elj = atomneb.read_elj(atom_elj_file, atom, ion, level_num=5)\n        s_ii_omij = atomneb.read_omij(atom_omij_file, atom, ion)\n        s_ii_aij = atomneb.read_aij(atom_aij_file, atom, ion)\n        \n        density = np.float64(1000)\n        temperature=np.float64(10000.0)#\n        nlj = pyequib.calc_populations(temperature=temperature, density=density, \n                                       elj_data=s_ii_elj, omij_data=s_ii_omij, aij_data=s_ii_aij)\n        print('Populations:', nlj)\n\n      which prints::\n      \n        Populations: 0.96992832 0.0070036315 0.023062261 2.6593671e-06 3.1277019e-06\n\n    - *Critical Density*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_dir = os.path.join('atomic-data', 'chianti70')\n        atom_elj_file = os.path.join(base_dir,data_dir, 'AtomElj.fits')\n        atom_omij_file = os.path.join(base_dir,data_dir, 'AtomOmij.fits')\n        atom_aij_file = os.path.join(base_dir,data_dir, 'AtomAij.fits')\n        \n        atom = 's'\n        ion = 'ii'\n        s_ii_elj = atomneb.read_elj(atom_elj_file, atom, ion, level_num=5)\n        s_ii_omij = atomneb.read_omij(atom_omij_file, atom, ion)\n        s_ii_aij = atomneb.read_aij(atom_aij_file, atom, ion)\n        \n        temperature=np.float64(10000.0)\n        n_crit = pyequib.calc_crit_density(temperature=temperature, \n                                           elj_data=s_ii_elj, omij_data=s_ii_omij, aij_data=s_ii_aij)\n        print('Critical Densities:', n_crit)\n\n      which gives::\n      \n        Critical Densities: 0.0000000 5007.8396 1732.8414 1072685.0 2220758.1\n\n    - *All Ionic Level Information*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_dir = os.path.join('atomic-data', 'chianti70')\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_elj_file = os.path.join(base_dir,data_dir, 'AtomElj.fits')\n        atom_omij_file = os.path.join(base_dir,data_dir, 'AtomOmij.fits')\n        atom_aij_file = os.path.join(base_dir,data_dir, 'AtomAij.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'h'\n        ion = 'ii' # H I Rec\n        hi_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        atom = 'o'\n        ion = 'iii' # [O III]\n        o_iii_elj = atomneb.read_elj(atom_elj_file, atom, ion, level_num=5) # read Energy Levels (Ej)\n        o_iii_omij = atomneb.read_omij(atom_omij_file, atom, ion) # read Collision Strengths (Omegaij)\n        o_iii_aij = atomneb.read_aij(atom_aij_file, atom, ion) # read Transition Probabilities (Aij)\n        \n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        pyequib.print_ionic(temperature=temperature, density=density,\n                    elj_data=o_iii_elj, omij_data=o_iii_omij, aij_data=o_iii_aij,\n                    h_i_aeff_data=hi_rc_data['aeff'][0])\n\n      which gives::\n      \n        Temperature =   10000.0 K\n        Density =    1000.0 cm-3\n        \n        Level    Populations   Critical Densities \n        Level 1:   3.063E-01   0.000E+00\n        Level 2:   4.896E-01   4.908E+02\n        Level 3:   2.041E-01   3.419E+03\n        Level 4:   4.427E-05   6.853E+05\n        Level 5:   2.985E-09   2.547E+07\n          \n         2.597E-05  \n             88.34um \n             (2--\u003e1) \n         2.859E-22  \n        \n         0.000E+00   9.632E-05  \n             32.66um      51.81um \n             (3--\u003e1)     (3--\u003e2) \n         0.000E+00   7.536E-22  \n        \n         2.322E-06   6.791E-03   2.046E-02  \n           4932.60A    4960.29A    5008.24A \n            (4--\u003e1)     (4--\u003e2)     (4--\u003e3) \n         4.140E-25   1.204E-21   3.593E-21  \n        \n         0.000E+00   2.255E-01   6.998E-04   1.685E+00  \n           2315.58A    2321.67A    2332.12A    4364.45A \n            (5--\u003e1)     (5--\u003e2)     (5--\u003e3)     (5--\u003e4) \n         0.000E+00   5.759E-24   1.779E-26   2.289E-23  \n        \n        H-beta emissivity: 1.237E-25 N(H+) Ne  [erg/s]\n\n\n* **Recombination Unit** has the API functions for plasma diagnostics and abundance analysis of recombination lines. Here are some examples of using *Recombination* Unit:\n\n    - *He+ Ionic Abundance*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_he_i_file = os.path.join(base_dir,data_rc_dir, 'rc_he_ii_PFSD12.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'h'\n        ion = 'ii' # H I\n        h_i_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        atom = 'he'\n        ion = 'ii' # He I\n        he_i_rc_data = atomneb.read_aeff_he_i_pfsd12(atom_rc_he_i_file, atom, ion)\n        \n        h_i_aeff_data = h_i_rc_data['aeff'][0]\n        he_i_aeff_data = he_i_rc_data['aeff'][0]\n        \n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        he_i_4471_flux= 2.104\n        linenum=10# 4471.50\n        abund_he_i = pyequib.calc_abund_he_i_rl(temperature=temperature, density=density,\n                                        linenum=linenum, line_flux=he_i_4471_flux,\n                                        he_i_aeff_data=he_i_aeff_data, h_i_aeff_data=h_i_aeff_data)\n        print('N(He^+)/N(H^+):', abund_he_i)\n\n      which gives::\n      \n        N(He^+)/N(H^+):     0.040848393\n\n    - *He++ Ionic Abundance*:\n\n      .. code-block:: Python\n      \n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'h'\n        ion = 'ii' # H I\n        h_i_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        atom = 'he'\n        ion = 'iii' # He II\n        he_ii_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        h_i_aeff_data = h_i_rc_data['aeff'][0]\n        he_ii_aeff_data = he_ii_rc_data['aeff'][0]\n        \n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        he_ii_4686_flux = 135.833\n        abund_he_ii = pyequib.calc_abund_he_ii_rl(temperature=temperature, density=density,\n                                          line_flux=he_ii_4686_flux,\n                                          he_ii_aeff_data=he_ii_aeff_data, h_i_aeff_data=h_i_aeff_data)\n        print('N(He^2+)/N(H^+):', abund_he_ii)\n\n      which gives::\n      \n        N(He^2+)/N(H^+):      0.11228817\n\n    - *C++ Ionic Abundance*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_all_file = os.path.join(base_dir,data_rc_dir, 'rc_collection.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'c'\n        ion = 'iii' # C II\n        c_ii_rc_data = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion)\n        \n        atom = 'h'\n        ion = 'ii' # H I\n        h_i_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n\n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        wavelength=6151.43\n        c_ii_6151_flux = 0.028\n        abund_c_ii = pyequib.calc_abund_c_ii_rl(temperature=temperature, density=density,\n                                        wavelength=wavelength, line_flux=c_ii_6151_flux,\n                                        c_ii_rc_data=c_ii_rc_data, h_i_aeff_data=h_i_aeff_data)\n        print('N(C^2+)/N(H+):', abund_c_ii)\n\n      which gives::\n      \n        N(C^2+)/N(H+):   0.00063404650 \n      \n    - *C3+ Ionic Abundance*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_ppb91_file = os.path.join(base_dir,data_rc_dir, 'rc_PPB91.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'c'\n        ion = 'iv' # C III\n        c_iii_rc_data = atomneb.read_aeff_ppb91(atom_rc_ppb91_file, atom, ion)\n        \n        atom = 'h'\n        ion = 'ii' # H I\n        h_i_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        wavelength=4647.42\n        c_iii_4647_flux = 0.107\n        abund_c_iii = pyequib.calc_abund_c_iii_rl(temperature=temperature, density=density,\n                                          wavelength=wavelength,\n                                          line_flux=c_iii_4647_flux, c_iii_rc_data=c_iii_rc_data,\n                                          h_i_aeff_data=h_i_aeff_data)\n        print('N(C^3+)/N(H+):', abund_c_iii)\n\n      which gives::\n      \n        N(C^3+)/N(H+):   0.00017502840\n\n    - *N++ Ionic Abundance*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_all_file = os.path.join(base_dir,data_rc_dir, 'rc_collection.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'n'\n        ion = 'iii' # N II\n        n_ii_rc_data = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion)\n        n_ii_rc_data_br = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion, br=True)\n        \n        atom = 'h'\n        ion = 'ii' # H I\n        h_i_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        wavelength=4442.02\n        n_ii_4442_flux = 0.017\n        abund_n_ii = pyequib.calc_abund_n_ii_rl(temperature=temperature, density=density,\n                                        wavelength=wavelength, line_flux=n_ii_4442_flux,\n                                        n_ii_rc_br=n_ii_rc_data_br, n_ii_rc_data=n_ii_rc_data,\n                                        h_i_aeff_data=h_i_aeff_data)\n        print('N(N^2+)/N(H+):', abund_n_ii)\n\n      which gives::\n      \n        N(N^2+)/N(H+):   0.00069297541\n\n    - *N3+ Ionic Abundance*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_ppb91_file = os.path.join(base_dir,data_rc_dir, 'rc_PPB91.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'n'\n        ion = 'iv' # N III\n        n_iii_rc_data = atomneb.read_aeff_ppb91(atom_rc_ppb91_file, atom, ion)\n           \n        atom = 'h'\n        ion = 'ii' # H I\n        h_i_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        wavelength=4640.64\n        n_iii_4641_flux = 0.245\n        abund_n_iii = pyequib.calc_abund_n_iii_rl(temperature=temperature, density=density,\n                                          wavelength=wavelength, line_flux=n_iii_4641_flux,\n                                          n_iii_rc_data=n_iii_rc_data, h_i_aeff_data=h_i_aeff_data)\n        print('N(N^3+)/N(H+):', abund_n_iii)\n\n      which gives::\n      \n        N(N^3+)/N(H+):   6.3366175e-05\n\n    - *O++ Ionic Abundance*:\n    \n      .. code-block:: Python\n      \n        import pyequib\n        import atomneb\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_all_file = os.path.join(base_dir,data_rc_dir, 'rc_collection.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'o'\n        ion = 'iii' # O II\n        o_ii_rc_data = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion)\n        o_ii_rc_data_br = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion, br=True)\n                   \n        atom = 'h'\n        ion = 'ii' # H I\n        h_i_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        wavelength=4613.68\n        o_ii_4614_flux = 0.009\n        abund_o_ii = pyequib.calc_abund_o_ii_rl(temperature=temperature, density=density,\n                                        wavelength=wavelength, line_flux=o_ii_4614_flux,\n                                        o_ii_rc_br=o_ii_rc_data_br,\n                                        o_ii_rc_data=o_ii_rc_data,\n                                        h_i_aeff_data=h_i_aeff_data)              \n        print('N(O^2+)/N(H+):', abund_o_ii)\n        \n      which gives::\n      \n        N(O^2+)/N(H+):    0.0018886330\n\n    - *Ne++ Ionic Abundance*:\n    \n      .. code-block:: Python\n      \n        import pyequib\n        import atomneb\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_all_file = os.path.join(base_dir,data_rc_dir, 'rc_collection.fits')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n        \n        atom = 'ne'\n        ion = 'iii' # Ne II\n        ne_ii_rc_data = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion)\n                   \n        atom = 'h'\n        ion = 'ii' # H I\n        h_i_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n        \n        wavelength=3777.14\n        ne_ii_3777_flux = 0.056\n        abund_ne_ii = pyequib.calc_abund_ne_ii_rl(temperature=temperature, density=density,\n                                          wavelength=wavelength, line_flux=ne_ii_3777_flux,\n                                          ne_ii_rc_data=ne_ii_rc_data, h_i_aeff_data=h_i_aeff_data)\n        print('N(Ne^2+)/N(H+):', abund_ne_ii)\n\n      which gives::\n      \n        N(Ne^2+)/N(H+):   0.00043376850\n\n\n    - *He I Emissivity*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_he_i_file = os.path.join(base_dir,data_rc_dir, 'rc_he_ii_PFSD12.fits')\n        \n        atom = 'he'\n        ion = 'ii' # He I\n        he_i_rc_data = atomneb.read_aeff_he_i_pfsd12(atom_rc_he_i_file, atom, ion)\n        \n        he_i_aeff_data = he_i_rc_data['aeff'][0]\n        \n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        linenum=10# 4471.50\n        emiss_he_i = pyequib.calc_emiss_he_i_rl(temperature=temperature, density=density,\n                                        linenum=linenum, he_i_aeff_data=he_i_aeff_data)\n        print('He I Emissivity:', emiss_he_i)\n\n      which gives::\n      \n        He I Emissivity:   6.3822830e-26\n\n    - *He II Emissivity*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_sh95_file = os.path.join(base_dir,data_rc_dir, 'rc_SH95.fits')\n         \n        atom = 'he'\n        ion = 'iii' # He II\n        he_ii_rc_data = atomneb.read_aeff_sh95(atom_rc_sh95_file, atom, ion)\n\n        he_ii_aeff_data = he_ii_rc_data['aeff'][0]\n\n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        emiss_he_ii = pyequib.calc_emiss_he_ii_rl(temperature=temperature, density=density,\n                                          he_ii_aeff_data=he_ii_aeff_data)\n        print('He II Emissivity:', emiss_he_ii)\n\n      which gives::\n      \n        He II Emissivity:   1.4989134e-24\n\n    - *C II Emissivity*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_all_file = os.path.join(base_dir,data_rc_dir, 'rc_collection.fits')\n        \n        atom = 'c'\n        ion = 'iii' # C II\n        c_ii_rc_data = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion)\n        \n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        wavelength=6151.43\n        emiss_c_ii = pyequib.calc_emiss_c_ii_rl(temperature=temperature, density=density,\n                                        wavelength=wavelength, c_ii_rc_data=c_ii_rc_data)\n        print('C II Emissivity:', emiss_c_ii)\n\n      which gives::\n      \n        C II Emissivity:   5.4719511e-26\n      \n    - *C III Emissivity*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import numpy as np\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_ppb91_file = os.path.join(base_dir,data_rc_dir, 'rc_PPB91.fits')\n        \n        atom = 'c'\n        ion = 'iv' # C III\n        c_iii_rc_data = atomneb.read_aeff_ppb91(atom_rc_ppb91_file, atom, ion)\n        \n        temperature=np.float64(10000.0)\n        density=np.float64(5000.0)\n        wavelength=4647.42\n        emiss_c_iii = pyequib.calc_emiss_c_iii_rl(temperature=temperature, density=density,\n                                          wavelength=wavelength,\n                                          c_iii_rc_data=c_iii_rc_data)\n        print('C III Emissivity:', emiss_c_iii)\n\n      which gives::\n      \n        C III Emissivity:   7.5749632e-25\n\n    - *N II Emissivity*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_all_file = os.path.join(base_dir,data_rc_dir, 'rc_collection.fits')\n        \n        atom = 'n'\n        ion = 'iii' # N II\n        n_ii_rc_data = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion)\n        n_ii_rc_data_br = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion, br=True)\n        \n        wavelength=4442.02\n        emiss_n_ii = pyequib.calc_emiss_n_ii_rl(temperature=temperature, density=density,\n                                        wavelength=wavelength,\n                                        n_ii_rc_br=n_ii_rc_data_br, n_ii_rc_data=n_ii_rc_data)\n        print('N II Emissivity:', emiss_n_ii)\n\n      which gives::\n      \n        N II Emissivity:   3.0397397e-26\n\n    - *N III Emissivity*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        import atomneb\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_ppb91_file = os.path.join(base_dir,data_rc_dir, 'rc_PPB91.fits')\n        \n        atom = 'n'\n        ion = 'iv' # N III\n        n_iii_rc_data = atomneb.read_aeff_ppb91(atom_rc_ppb91_file, atom, ion)\n        \n        wavelength=4640.64\n        emiss_n_iii = pyequib.calc_emiss_n_iii_rl(temperature=temperature, density=density,\n                                          wavelength=wavelength, n_iii_rc_data=n_iii_rc_data)\n        print('N III Emissivity:', emiss_n_iii)\n\n      which gives::\n      \n        N III Emissivity:   4.7908644e-24\n\n    - *O II Emissivity*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_all_file = os.path.join(base_dir,data_rc_dir, 'rc_collection.fits')\n        \n        atom = 'o'\n        ion = 'iii' # O II\n        o_ii_rc_data = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion)\n        o_ii_rc_data_br = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion, br=True)\n        \n        wavelength=4613.68\n        emiss_o_ii = pyequib.calc_emiss_o_ii_rl(temperature=temperature, density=density,\n                                        wavelength=wavelength,\n                                        o_ii_rc_br=o_ii_rc_data_br, o_ii_rc_data=o_ii_rc_data)\n        print('O II Emissivity:', emiss_o_ii)\n        \n      which gives::\n      \n        O II Emissivity:   5.9047319e-27\n\n    - *Ne II Emissivity*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        import atomneb\n        import os\n        base_dir = 'externals/atomneb'\n        data_rc_dir = os.path.join('atomic-data-rc')\n        atom_rc_all_file = os.path.join(base_dir,data_rc_dir, 'rc_collection.fits')\n        \n        atom = 'ne'\n        ion = 'iii' # Ne II\n        ne_ii_rc_data = atomneb.read_aeff_collection(atom_rc_all_file, atom, ion)\n        \n        wavelength=3777.14\n        emiss_ne_ii = pyequib.calc_emiss_ne_ii_rl(temperature=temperature, density=density,\n                                          wavelength=wavelength, ne_ii_rc_data=ne_ii_rc_data)\n        print('Ne II Emissivity:', emiss_ne_ii)\n\n      which gives::\n      \n        Ne II Emissivity:   1.5996881e-25\n        \n* **Reddening Unit** has the API functions for estimating logarithmic extinctions at H-beta and dereddening observed fluxes based on reddening laws and extinctions. Here are some examples of using *Reddening* Unit:\n\n    - *Reddening Law Function*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        wavelength=6563.0\n        r_v=3.1\n        fl=pyequib.redlaw(wavelength, rv=r_v, ext_law='GAL')\n        print('fl(6563):', fl)\n\n      which gives::\n      \n        fl(6563):     -0.32013816\n\n    - *Galactic Reddening Law Function based on Seaton (1979), Howarth (1983), \u0026 CCM (1983)*:\n    \n      .. code-block:: Python\n      \n        import pyequib\n        wavelength=6563.0\n        r_v=3.1\n        fl=pyequib.redlaw_gal(wavelength, rv=r_v)\n        print('fl(6563):', fl)\n\n      which gives::\n      \n        fl(6563):     -0.32013816\n\n    - *Galactic Reddening Law Function based on Savage \u0026 Mathis (1979)*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        wavelength=6563.0\n        fl=pyequib.redlaw_gal2(wavelength)\n        print('fl(6563):', fl)\n\n      which gives::\n      \n        fl(6563):     -0.30925984\n\n    - *Reddening Law Function based on Cardelli, Clayton \u0026 Mathis (1989)*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        wavelength=6563.0\n        r_v=3.1\n        fl=pyequib.redlaw_ccm(wavelength, rv=r_v)\n        print('fl(6563):', fl)\n\n      which gives::\n      \n        fl(6563):     -0.29756615\n\n    - *Galactic Reddening Law Function based on Whitford (1958), Seaton (1977), \u0026 Kaler(1976)*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        wavelength=6563.0\n        fl=pyequib.redlaw_jbk(wavelength)\n        print('fl(6563):', fl)\n\n      which gives::\n      \n        fl(6563):     -0.33113684\n\n    - *Reddening Law Function based on Fitzpatrick \u0026 Massa (1990), Fitzpatrick (1999), Misselt (1999)*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        wavelength=6563.0\n        r_v=3.1\n        fmlaw='AVGLMC'\n        fl=pyequib.redlaw_fm(wavelength, fmlaw=fmlaw, rv=r_v)\n        print('fl(6563):', fl)\n\n      which gives::\n      \n        fl(6563):     -0.35053032\n\n    - *Reddening Law Function for the Small Magellanic Cloud*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        wavelength=6563.0\n        fl=pyequib.redlaw_smc(wavelength)\n        print('fl(6563):', fl)\n\n      which gives::\n      \n        fl(6563):     -0.22659261\n\n    - *Reddening Law Function for the Large Magellanic Cloud*:\n    \n      .. code-block:: Python\n    \n        import pyequib\n        wavelength=6563.0\n        fl=pyequib.redlaw_lmc(wavelength)\n        print('fl(6563):', fl)\n\n      which gives::\n      \n        fl(6563):     -0.30871187\n\n    - *Dereddening Relative Flux*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        wavelength=6563.0\n        m_ext=1.0\n        flux=1.0\n        ext_law='GAL'\n        r_v=3.1\n        flux_deredden=pyequib.deredden_relflux(wavelength, flux, m_ext, ext_law=ext_law, rv=r_v)\n        print('dereddened flux(6563)', flux_deredden)\n\n      which gives::\n      \n        dereddened flux(6563)       0.47847785\n\n    - *Dereddening Absolute Flux*:\n    \n      .. code-block:: Python\n\n        import pyequib\n        wavelength=6563.0\n        m_ext=1.0\n        flux=1.0\n        ext_law='GAL'\n        r_v=3.1\n        flux_deredden=pyequib.deredden_flux(wavelength, flux, m_ext, ext_law=ext_law, rv=r_v)\n        print('dereddened flux(6563)', flux_deredden)\n\n      which gives::\n      \n        dereddened flux(6563)      4.7847785\n\nDocumentation\n=============\n\nFor more information on how to use the API functions from the pyEQUIB libray, please read the `API Documentation  \u003chttps://equib.github.io/pyEQUIB/doc\u003e`_ published on `equib.github.io/pyEQUIB \u003chttps://equib.github.io/pyEQUIB\u003e`_.\n\n\nReferences\n==========\n* Danehkar, A. (2020). pyEQUIB Python Package, an addendum to proEQUIB: IDL Library for Plasma Diagnostics and Abundance Analysis. *J. Open Source Softw.*, **5**, 2798. doi: `10.21105/joss.02798 \u003chttps://doi.org/10.21105/joss.02798\u003e`_ ads: `2020JOSS....5.2798D \u003chttps://ui.adsabs.harvard.edu/abs/2020JOSS....5.2798D/abstract\u003e`_.\n\n* Danehkar, A. (2018). proEQUIB: IDL Library for Plasma Diagnostics and Abundance Analysis. *J. Open Source Softw.*, **3**, 899. doi: `10.21105/joss.00899 \u003chttps://doi.org/10.21105/joss.00899\u003e`_ ads: `2018JOSS....3..899D \u003chttps://ui.adsabs.harvard.edu/abs/2018JOSS....3..899D/abstract\u003e`_.\n\n* Danehkar, A. (2018). Bi-Abundance Ionisation Structure of the Wolf-Rayet Planetary Nebula PB 8, *PASA*, **35**, e005.  doi: `10.1017/pasa.2018.1 \u003chttps://doi.org/10.1017/pasa.2018.1\u003e`_ ads: `2018PASA...35....5D \u003chttps://ui.adsabs.harvard.edu/abs/2018PASA...35....5D/abstract\u003e`_.\n\n* Danehkar, A. (2021). Physical and Chemical Properties of Wolf-Rayet Planetary Nebulae, *ApJS*, **257**, 58.  doi: `10.3847/1538-4365/ac2310 \u003chttps://doi.org/10.3847/1538-4365/ac2310\u003e`_ ads: `2021ApJS..257...58D \u003chttps://ui.adsabs.harvard.edu/abs/2021ApJS..257...58D/abstract\u003e`_.\n\n\nCitation\n========\n\nUsing the **pyEQUIB** Python package in a scholarly publication? Please cite thess papers:\n\n.. code-block:: bibtex\n\n   @article{Danehkar2020,\n     author = {{Danehkar}, Ashkbiz},\n     title = {pyEQUIB Python Package, an addendum to proEQUIB: IDL Library \n              for Plasma Diagnostics and Abundance Analysis},\n     journal = {Journal of Open Source Software},\n     volume = {5},\n     number = {55},\n     pages = {2798},\n     year = {2020},\n     doi = {10.21105/joss.02798}\n   }\n\nand if you use the `proEQUIB \u003chttps://github.com/equib/proEQUIB\u003e`_ IDL library:\n\n.. code-block:: bibtex\n\n   @article{Danehkar2018,\n     author = {{Danehkar}, Ashkbiz},\n     title = {proEQUIB: IDL Library for Plasma Diagnostics and Abundance Analysis},\n     journal = {Journal of Open Source Software},\n     volume = {3},\n     number = {32},\n     pages = {899},\n     year = {2018},\n     doi = {10.21105/joss.00899}\n   }\n\nLearn More\n==========\n\n==================  =============================================\n**Documentation**   https://pyequib.readthedocs.io/\n**Repository**      https://github.com/equib/pyEQUIB\n**Issues \u0026 Ideas**  https://github.com/equib/pyEQUIB/issues\n**Conda-Forge**     https://anaconda.org/conda-forge/pyequib\n**PyPI**            https://pypi.org/project/pyequib/\n**DOI**             `10.21105/joss.02798 \u003chttps://doi.org/10.21105/joss.02798\u003e`_\n**Archive**         `10.5281/zenodo.4287575 \u003chttps://doi.org/10.5281/zenodo.4287575\u003e`_\n==================  =============================================\n","project_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fequib%2Fpyequib","html_url":"https://awesome.ecosyste.ms/projects/github.com%2Fequib%2Fpyequib","lists_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fequib%2Fpyequib/lists"}