Calculate a large spectrum

Calculate a full-range spectrum on several partial spectral ranges, and combine them.

plot merged large spectrum

Out:

parse_local_quanta not implemented for molecules of HITRAN group 6
parse_global_quanta not implemented for molecules of HITRAN class 3
Created ~/radis.json in /home/docs
Added HITEMP-OH database in /home/docs/radis.json
Spectrum Name:  OHspectrum140135626572688//OHspectrum140135626572808//OHspectrum140135626572928
Spectral Quantities
----------------------------------------
   abscoeff     (99,502 points)
   emisscoeff   (99,502 points)
   absorbance   (99,502 points)
   transmittance_noslit         (99,502 points)
   emissivity_noslit    (99,502 points)
   radiance_noslit      (99,502 points)
Populations Stored
----------------------------------------
Physical Conditions
----------------------------------------
   Tgas                 1000 K
   Trot                 1000 K
   Tvib                 1000 K
   isotope              N/A
   mole_fraction        0.1
   molecule             {'OH'}
   path_length          1 cm
   pressure_mbar        10000.0 mbar
   rot_distribution     boltzmann
   self_absorption      True
   state                X
   thermal_equilibrium  True
   vib_distribution     boltzmann
   wavenum_max          N/A cm-1
   wavenum_min          N/A cm-1
Computation Parameters
----------------------------------------
   Tref                 296 K
   add_at_used          numpy
   broadening_max_width  10 cm-1
   broadening_method    fft
   cutoff               1e-27 cm-1/(#.cm-2)
   dbformat             hitemp-radisdb
   dbpath               /home/docs/.radisdb/OH-13_HITEMP2020.h5
   dlm_log_pG           0.1375350788016573
   dlm_log_pL           0.20180288881201608
   folding_thresh       1e-06
   include_neighbouring_lines  True
   optimization         min-RMS
   parfuncfmt           hapi
   pseudo_continuum_threshold  0
   radis_version        0.9.29
   wavenum_max_calc     N/A cm-1
   wavenum_min_calc     N/A cm-1
   waveunit             cm-1
   wstep                0.1 cm-1
   zero_padding         N/A
Information
----------------------------------------
   calculation_time     0.8400000000000001 s
   db_use_cached        True
   export_lines         False
   lines_calculated     5,710
   lines_cutoff         21,940
   lines_in_continuum   0
   load_energies        False
   lvl_use_cached       True
   total_lines          27650
   warning_broadening_threshold  0.01
   warning_linestrength_cutoff  0.01
----------------------------------------


<matplotlib.lines.Line2D object at 0x7f73de39dac8>

from radis import MergeSlabs, calc_spectrum

spectra = []
for (wmin, wmax) in [(50, 3000), (3000, 7000), (7000, 10000)]:

    spectra.append(
        calc_spectrum(
            wmin,
            wmax,
            Tgas=1000,
            pressure=10,  # bar
            molecule="OH",
            path_length=1,
            mole_fraction=0.1,
            wstep=0.1,
            databank="hitemp",
            verbose=False,
        )
    )

s = MergeSlabs(*spectra, resample="full", out="transparent")
print(s)
s.plot("transmittance_noslit", wunit="nm")

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

Gallery generated by Sphinx-Gallery