I am having trouble extracting the absorption coefficients for carbon dioxide using the hitran2016 database. I am interested in the wavelengths between ca. 1200 and 2600 nm. For methane I get expected results, based on what i have seen others get. For carbon dioxide, however, i am expecting a peak of about 0.7 at 2000 nm. However, the peak i get in my plot is of just above 0.1.
My code for methane is:
from hapi import *import numpy as npimport matplotlib.pyplot as pltwl_min = 1200 #[nm] wavelegnth startwl_max = 2600 #[nm] wavelenth stopwn_min = 1/(wl_max*10**(-7)) # Wavenumber start [cm^-1]wn_max = 1/(wl_min*10**(-7)) # Wavenumber stop [cm^-1]fetch('CH4',6,1,wn_min,wn_max)describeTable('CH4')nu,coef = absorptionCoefficient_Voigt(SourceTables='CH4',HITRAN_units=False, Diluent={'self':1})And for carbon dioxide i did:
fetch('CO2',2,1,wn_min,wn_max)describeTable('CO2')nu,coef = absorptionCoefficient_Voigt(SourceTables='CO2',HITRAN_units=False, Diluent={'self':1})Not sure if I have done something differently by mistake, or if there are any bugs or misunderstandings i should fix. Any help would be very appreciated.