HITRAN/HITEMP
May (2021) Hajime Kawahara
See “Computing CO cross section using HITRAN (opacity calculator = LPF)” for tutorial to compute a cross section profile using HITRAN/HITEMP.
Molecular Database
An example to use HITRAN / HITEMP from exojax is like that.
>>> from exojax.utils.grids import wavenumber_grid
>>> from exojax.spec import moldb
>>> nus,wav,res=nugrid(22880.,23000.,1000,unit="AA")
>>> mdbCO=moldb.MdbHit('.database/05_HITEMP2019.par.bz2',nus)
bunziping
HITEMP H2O and CO2
For H2O and CO2, HITEMP provides multiple par files. To use those files, provide the directory path for moldb.MdbHit
as follows.
>>> mdbH2O=moldb.MdbHit('.database/01_HITEMP2010',nus)
>>> mdbCO2=moldb.MdbHit('.database/02_HITEMP2010',nus)
extract option
extract=True
in moldb.MdbHit
extracts the opacity data in the wavenumber range of nus
with margin
. Theforefore it can reduce the use of DRAM. It may be useful for large databases such as CH4.
>>> mdbCH4=moldb.MdbHit('.database/06_HITEMP2020.par.bz2',nus,extract=True)
This creates a new directory (such as 6101.281269066504_6108.7354917532075_1.0
) and the extracted data and header files (06_HITEMP2020.header
and 06_HITEMP2020.par
) in it:
ls .database
06_HITEMP2020.par
6101.281269066504_6108.7354917532075_1.0:
06_HITEMP2020.header 06_HITEMP2020.par
Basic Quantities
These are the basic quantities of MdbHit.
quantity |
instance |
unit |
np/jnp |
line center |
nu_lines |
cm-1 |
np |
line center |
dev_nu_lines |
cm-1 |
jnp |
lower state energy |
elower |
cm-1 |
jnp |
natural broadening |
gamma_natural |
cm-1 |
jnp |
air pressure broadening |
gamma_air |
cm-1 |
jnp |
self broadning |
gamma_self |
cm-1 |
jnp |
Einstein coefficient |
A |
s-1 |
jnp |
reference line strength |
Sij0 |
cm |
np |
log_e Sij0 |
logsij0 |
jnp |
|
statistical weight |
gpp |
jnp |
|
temperature exponent |
n_air |
jnp |
Collision Induced Absorption (CIA)
Collision Induced Absorption (CIA) is one of the important continuum in exoplanet/brown dwarfs atmosphere. HITRAN provides the CIA files . To load the CIA file, use contdb.
>>> from exojax.spec import contdb
>>> cdbH2H2=contdb.CdbCIA('.database/H2-H2_2011.cia',nus)
See “CIA coefficient” for tutorial.