oh_einstein_temp_convert is a Python package which translates temperatures originally calculated from OH*(6-2) transiton spectra using a certain set of Einstein coefficients to a new set of Einstein coefficients with no need of the original data.
It is published on Pypi and kept here, this repo also includes a Jupyter notebook with the same code in a more prose-style with plots to explain the steps.
Rowan Dayton-Oxland University of Southampton Github
Use the package manager pip to install oh_einstein_temp_convert.
pip install oh_einstein_temp_convert
- Python >= 3.10
- numpy, pandas
- OS independent
Ain (int) - Original Einstein coefficient source index (see table)
Aout (int) - Desired output Einstein coefficient source index (see table)
Einstein A Source | Index |
---|---|
Mies et al., 1974 | 0 |
Loo and Groenenboom, 2008 | 1 |
Langhoff et al., 1986 | 2 |
Goldman et al., 1998 | 3 |
Turnbull and Lowe, 1989 | 4 |
import oh_einstein_temp_convert as oh
Temperatures = [] # List of temperatures
# Temperatures originally calculated from e.g. Mies et al., 1974 coefficients
Ain_index = 0 # Choose index of input temperatures
# Output temperatures calculated from e.g. Loo and Groenenboom, 2008 coefficients.
Aout_index = 1 # Choose index of output temperatures
# returns list of converted temperatures
Result = oh.convert_temperatures(Temperatures, Ain_index, Aout_index)
- Get the original temperature and original Einstein coefficient set
- Calculate
$ln(\frac{I}{Ain \cdot 2(2J' + 1)})$ and add a correction term for the new Einstein coefficient set$ln(\frac{Ain}{Aout})$ for each$F(J')$ value in the OH*(6-2) P-branch. - Plot the corrected
$ln(\frac{I}{Aout \cdot 2(2J' + 1)})$ term against$F(J')$ - Extract the output temperature from the gradient of the line by linear fit
This comes from the following equation; $$ ln(\frac{I}{A \cdot 2(2J' + 1)}) = \frac{-(h c)}{(k T)} * F(J') + ln(\frac{N}{QR}) $$
Pull requests are welcome. For major changes, please open an issue first to discuss what you would like to change.
Einstein and other quantum coefficients, and general inspiration from the Synthetic Hydroxyl Spectrum Generator Sigernes, F., Shumilov, N., Deehr, C.S., Nielsen, K.P., Svenøe, T., and Havnes, O., The Hydroxyl rotational temperature record from the Auroral Station in Adventdalen, Svalbard (78°N, 15°E) , Journal of Geophysical Research, Vol 108 (A9), 1342, doi 1029/2001JA009023, 2003.