Line Survey

Plot details of every single line in a spectrum.

Uses the line_survey() function.

Out:

<frozen importlib._bootstrap>:219: RuntimeWarning: numpy.ndarray size changed, may indicate binary incompatibility. Expected 80 from C header, got 88 from PyObject
Generating cache file /home/docs/.astropy/cache/astroquery/Hitran/tempfile_CO2_1_2380.00_2400.00.h5 with metadata :
{'molecule': 'CO2', 'isotope': 1, 'wmin': 2380.0, 'wmax': 2400.0}
/home/docs/checkouts/readthedocs.org/user_builds/radis/envs/latest/lib/python3.8/site-packages/pandas/core/generic.py:2703: PerformanceWarning:
your performance may suffer as PyTables will pickle object types that it cannot
map directly to c-types [inferred_type->bytes,key->block3_values] [items->Index(['globu', 'globl', 'locu', 'locl', 'lmix'], dtype='object')]

  pytables.to_hdf(
HAPI version: 1.1.2.0
To get the most up-to-date version please check http://hitran.org/hapi
ATTENTION: Python versions of partition sums from TIPS-2017 are now available in HAPI code

           It is free to use HAPI. If you use HAPI in your research or software development,
           please cite it using the following reference:
           R.V. Kochanov, I.E. Gordon, L.S. Rothman, P. Wcislo, C. Hill, J.S. Wilzewski,
           HITRAN Application Programming Interface (HAPI): A comprehensive approach
           to working with spectroscopic data, J. Quant. Spectrosc. Radiat. Transfer 177, 15-30 (2016)
           DOI: 10.1016/j.jqsrt.2016.03.005
/home/docs/checkouts/readthedocs.org/user_builds/radis/envs/latest/lib/python3.8/site-packages/radis/misc/warning.py:350: HighTemperatureWarning: HITRAN is valid for low temperatures (typically < 700 K). For higher temperatures you may need HITEMP or CDSD. See the 'databank=' parameter
  warnings.warn(WarningType(message))
Calculating Equilibrium Spectrum
Physical Conditions
----------------------------------------
   Tgas                 1500 K
   Trot                 1500 K
   Tvib                 1500 K
   isotope              1
   mole_fraction        0.0004
   molecule             CO2
   overpopulation       None
   path_length          100 cm
   pressure_mbar        1013.25 mbar
   rot_distribution     boltzmann
   self_absorption      True
   state                X
   vib_distribution     boltzmann
   wavenum_max          2400.0000 cm-1
   wavenum_min          2380.0000 cm-1
Computation Parameters
----------------------------------------
   Tref                 296 K
   add_at_used          None
   broadening_method    voigt
   cutoff               1e-27 cm-1/(#.cm-2)
   dbformat             hitran
   dbpath               fetched from hitran
   dlm_log_pG           0.1375350788016573
   dlm_log_pL           0.20180288881201608
   export_rovib_fraction  True
   folding_thresh       1e-06
   hdf5_engine          pytables
   include_neighbouring_lines  True
   levelsfmt            None
   neighbour_lines      0 cm-1
   optimization         min-RMS
   parfuncfmt           hapi
   parfuncpath          None
   parsum_mode          full summation
   pseudo_continuum_threshold  0
   truncation           50 cm-1
   wavenum_max_calc     2400.0000 cm-1
   wavenum_min_calc     2380.0000 cm-1
   waveunit             cm-1
   wstep                0.01 cm-1
   zero_padding         -1
Information
----------------------------------------
   chunksize            None
   db_use_cached        True
   export_lines         True
   export_populations   None
   load_energies        False
   lvl_use_cached       True
   total_lines          332
   warning_broadening_threshold  0.01
   warning_linestrength_cutoff  0.01
----------------------------------------
0.02s - Spectrum calculated
288 lines

Figure({
    'data': [{'name': 'linestrength',
              'text': array(['CO2 [Q32](11`1`0 2)->(03`3`1 1)<br>int intensity at 296K: 2.12e-29  [cm-1/(molecule/cm-2)]<br>A Einstein A coefficient: 3.62e-05  [s-1]<br>El lower-state energy: 2.35e+03  [cm-1]',
                             'CO2 [P52](10`0`0 2)->(10`0`1 1)<br>int intensity at 296K: 7.49e-26  [cm-1/(molecule/cm-2)]<br>A Einstein A coefficient: 0.0867  [s-1]<br>El lower-state energy: 2.36e+03  [cm-1]',
                             'CO2 [Q41](11`1`0 2)->(03`3`1 1)<br>int intensity at 296K: 1.13e-29  [cm-1/(molecule/cm-2)]<br>A Einstein A coefficient: 5.33e-05  [s-1]<br>El lower-state energy: 2.6e+03  [cm-1]',
                             ...,
                             'CO2 [P34](20`0`0 2)->(20`0`1 1)<br>int intensity at 296K: 8.86e-28  [cm-1/(molecule/cm-2)]<br>A Einstein A coefficient: 0.0691  [s-1]<br>El lower-state energy: 3.13e+03  [cm-1]',
                             'CO2 [P5](10`0`1 2)->(10`0`2 1)<br>int intensity at 296K: 3.12e-29  [cm-1/(molecule/cm-2)]<br>A Einstein A coefficient: 0.196  [s-1]<br>El lower-state energy: 3.62e+03  [cm-1]',
                             'CO2 [P32](20`0`0 3)->(20`0`1 2)<br>int intensity at 296K: 2.21e-27  [cm-1/(molecule/cm-2)]<br>A Einstein A coefficient: 0.0789  [s-1]<br>El lower-state energy: 2.96e+03  [cm-1]'],
                            dtype=object),
              'type': 'bar',
              'width': 0.01,
              'x': array([2380.01644025, 2380.0813248 , 2380.11409764, ..., 2399.77908456,
                          2399.90028117, 2399.96245817]),
              'y': array([5.52335028e-27, 2.06837048e-23, 8.10164666e-27, ..., 5.02841105e-24,
                          1.19566096e-24, 6.37840048e-24])},
             {'name': 'radiance',
              'type': 'scatter',
              'x': array([2380.  , 2380.01, 2380.02, ..., 2399.99, 2400.  , 2400.01]),
              'y': array([nan, nan, nan, ..., nan, nan, nan]),
              'yaxis': 'y2'}],
    'layout': {'hovermode': 'closest',
               'showlegend': False,
               'template': '...',
               'title': {'text': 'Line Survey (1500K, 1.013bar, Mfrac=0.000)'},
               'xaxis': {'title': {'text': 'Wavenumber (cm-1)'}},
               'yaxis': {'range': [-26, -17],
                         'tickfont': {'color': '#1f77b4'},
                         'title': {'font': {'color': '#1f77b4'}, 'text': 'Linestrength (cm-1/(molecule/cm-2))'},
                         'type': 'log'},
               'yaxis2': {'overlaying': 'y',
                          'side': 'right',
                          'tickfont': {'color': '#ff7f0e'},
                          'title': {'font': {'color': '#ff7f0e'}, 'text': 'Radiance (mW/cm2/sr/nm)'},
                          'type': 'log'}}
})

from radis import calc_spectrum

s = calc_spectrum(
    wavenum_min=2380,
    wavenum_max=2400,
    mole_fraction=400e-6,
    path_length=100,  # cm
    Tgas=1500,
    molecule="CO2",
    isotope=[1],
    databank="hitran",
    export_lines=True,
)
s.apply_slit(2, "nm")
s.line_survey(overlay="radiance", barwidth=0.01)

Total running time of the script: ( 0 minutes 1.570 seconds)

Gallery generated by Sphinx-Gallery