Skip to main content

Rosseland And Planck Opacity Converter

Project description

PyPI version Documentation Status GitHub

RAPOC: Rosseland And Planck Opacity Converter

The RAPOC code is written by Lorenzo V. Mugnai and Darius Modirrousta-Galian and is the product of a collaboration between Sapienza Università di Roma, Università degli Studi di Palermo and INAF - Osservatorio Astronomico di Palermo. It uses molecular absorption measurements (i.e. wavelength-dependent opacities) to calculate Rosseland and Planck mean opacities that are commonly used in atmospheric modelling.

RAPOC is designed to be simple, straightforward, and easily incorporated into other codes. It is completely written in Python and documented with docstrings. In addition, a Sphinx version of the documentation with a full user guide that includes examples is available in html format.


RAPOC is under development, please report any issues or inaccuracies to the developers to support the implementation.


If you use this code or its results, please cite RAPOC: the Rosseland and Planck opacity converter by Mugnai L. V. and Modirrousta-Galian D. (submitted).


Installing from Pypi

RAPOC can be installed from the Pypi repository with the following script::

pip install rapoc

Installing from git

RAPOC may also be cloned from the main git repository::

git clone

The next step is to move into the RAPOC folder::

cd /your_path/Rapoc


pip install .

To check if one has the correct setup::

python -c "import rapoc"


RAPOC is designed to be used on its own or in conjunction with other Python codes. Given an ExoMol file in the TauREx.h5 format, Rosseland and Planck mean opacities can be calculated. For example, in order to estimate the mean opacities at a temperature (T) of 1000 K with a pressure (P) of 10,000 Pa in the wavelength range of 0.3-50 micron the following script is used,

from rapoc import Rosseland, Planck

r_model = Rosseland(input_data='exomol_file.TauREx.h5')
opacity = r_model.estimate(P_input=10000 * u.Pa, T_input=1000 * u.K, band=(0.3 *, 50*

p_model = Planck(input_data='exomol_file.TauREx.h5')
opacity = p_model.estimate(P_input=10000 * u.Pa, T_input=1000 * u.K, band=(0.3 *, 50*


To run the code you need measured data. The supported file formats are:

  • ExoMol opacities (downloadable here) with the TauREx.h5 format.


The full documentation is available here

Alternatively, RAPOC accepts user-defined documentation by using sphinx. To install it run

pip install sphinx sphinx_rtd_theme

From the Rapoc/docs folder running

cd docs
make html

This will create a html version of the documentation in Rapoc/doc/build/html/index.html.

Project details

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for rapoc, version 1.0.5
Filename, size File type Python version Upload date Hashes
Filename, size rapoc-1.0.5-py3-none-any.whl (25.2 kB) File type Wheel Python version py3 Upload date Hashes View
Filename, size rapoc-1.0.5.tar.gz (19.6 kB) File type Source Python version None Upload date Hashes View

Supported by

AWS AWS Cloud computing Datadog Datadog Monitoring DigiCert DigiCert EV certificate Facebook / Instagram Facebook / Instagram PSF Sponsor Fastly Fastly CDN Google Google Object Storage and Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Salesforce Salesforce PSF Sponsor Sentry Sentry Error logging StatusPage StatusPage Status page