pyrt.rayleigh_co2#

pyrt.rayleigh_co2(column_density, wavelength)[source]#

Compute the Rayleigh CO ₂ Column.

Parameters:

column_density (ArrayLike) – 1-dimensional array of the column density [\(\frac{particles}{\text{m^2}}\)] in each layer.
wavelength (ArrayLike) – 1-dimensional array of the wavelengths [microns] to compute the optical depth at.

Returns:

A Rayleigh column.

Return type:

Column

Notes

This algorithm computes the optical depth in each layer i by using

The molecular cross-section is given by laboratory measurements from Sneep and Ubachs, 2005, JQSRT.

Examples

Compute the Rayleigh scattering optical depth from (very simplified) column densities and wavelengths.

>>> import numpy as np
>>> import pyrt
>>> column_density = np.linspace(1, 10, num=15) * 10**26
>>> wavs = np.linspace(0.2, 1, num=5)
>>> rayleigh_column = pyrt.rayleigh_co2(column_density, wavs)
>>> rayleigh_column.optical_depth.shape
(15, 5)

Simply sum over the layers to get the column integrated optical depth due to Rayleigh scattering.

>>> np.sum(rayleigh_column.optical_depth, axis=0)
array([0.78456184, 0.03590366, 0.00672394, 0.00208873, 0.00084833])