Getting Cross Sections

In this example we pull out cross sections from the ozone DBM optical property.

%matplotlib inline

import sasktran as sk
import matplotlib.pyplot as plt
import numpy as np

# wavelength grid
wavel_nm = np.arange(280, 700, 0.01)

# Used for temperature/pressure
background_atmosphere = sk.MSIS90()

# High resolution ozone cross sections
hires = sk.O3DBM()

# Pull out the absorption cross section at 20 km
hires_xsec = hires.calculate_cross_sections(background_atmosphere, latitude=0,
                                            longitude=0, altitude=20000, wavelengths=wavel_nm, mjd=54372)

plt.plot(wavel_nm, hires_xsec.absorption)
plt.yscale('log')
plt.xlabel('Wavelength [nm]')
plt.ylabel('Cross Section [cm$^2$]')

plt.show()