Calculate a large spectrum

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

Uses the MergeSlabs() function

plot merged large spectrum

Out:

Spectrum Name:  OHspectrum140197093785408//OHspectrum140197093650496//OHspectrum140197093652800
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)
Physical Conditions
----------------------------------------
   Tgas                 1000 K
   Trot                 1000 K
   Tvib                 1000 K
   isotope              N/A
   mole_fraction        0.1
   molecule             {'OH'}
   overpopulation       None
   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          cython
   broadening_method    voigt
   cutoff               1e-27 cm-1/(#.cm-2)
   dbformat             hitemp-radisdb
   dbpath               /home/docs/.radisdb/hitemp/OH-13_HITEMP2020.h5
   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
   profiler             N/A
   pseudo_continuum_threshold  0
   radis_version        0.10.2
   spectral_points      N/A
   truncation           50 cm-1
   wG                   N/A
   wL                   4
   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.3700845241546631 s
   chunksize            None
   db_use_cached        True
   export_lines         False
   export_populations   None
   lines_calculated     5,672
   lines_cutoff         21,775
   lines_in_continuum   0
   load_energies        False
   lvl_use_cached       True
   total_lines          27447
   warning_broadening_threshold  0.01
   warning_linestrength_cutoff  0.01
----------------------------------------


<matplotlib.lines.Line2D object at 0x7f822de77970>

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 1.993 seconds)

Gallery generated by Sphinx-Gallery