Skip to content

synth_telluric

Generate a synthetic telluric transmission spectrum.
Category Corrections
Backend numpy - implemented here on top of numpy primitives
Version 1.0.1
Reads -
Writes extras.telluric_template, spectrum

Each band is modelled as an inverted Gaussian (1 - depth * G) with hand-curated centre / depth / FWHM — a quick-look approximation, not a line-by-line radiative-transfer model (for that use molecfit, Smette et al. 2015, A&A 576, A77, or an observed telluric standard). Depths and widths are unit-airmass values; the airmass parameter scales them via T = T0**airmass, the Bouguer / Beer-Lambert relation. The output spectrum is stored in ctx.extras['telluric_template'] and, when store_as_spectrum=True (default False), also on ctx.spectrum. remove_telluric_division falls back to this key (after 'telluric_template_scaled') when its default key is absent.

Parameters

Parameter Default Required Description
wavelength_min 5800.0 - Lower bound of the synthetic axis (Å).
wavelength_max 9500.0 - Upper bound of the synthetic axis (Å).
n_points 4000 - Number of samples in the synthetic axis.
airmass 1.0 - Atmospheric path length (1.0 = zenith).
store_as_spectrum False - If true, also write the template to ctx.spectrum.

Use it

from spectro_kernel import run_algorithm

output = run_algorithm("synth_telluric", ctx, {
    "wavelength_min": 5800.0,
    "wavelength_max": 9500.0,
    "n_points": 4000,
    "airmass": 1.0,
    "store_as_spectrum": False
})

spectro run synth_telluric --input spectrum.fits \
  --param wavelength_min=5800.0 \
  --param wavelength_max=9500.0 \
  --param n_points=4000 \
  --param airmass=1.0 \
  --param store_as_spectrum=false

{
  "tool": "synth_telluric",
  "arguments": {
    "session_id": "<session_id>",
    "params": {
      "wavelength_min": 5800.0,
      "wavelength_max": 9500.0,
      "n_points": 4000,
      "airmass": 1.0,
      "store_as_spectrum": false
    }
  }
}

Every algorithm is an MCP tool of the same name; describe_algorithm returns this page's metadata as JSON.

References

  • Hinkle, Wallace & Livingston 2003, BAAS 35, 1260 — Kitt Peak atmospheric transmission atlas 0.5-5.5 µm (positions of the O2 γ/B/A bands and the H2O bands used here).
  • Hardie 1962, in Astronomical Techniques (Stars and Stellar Systems II, ed. Hiltner), Univ. Chicago Press, p. 178 — Bouguer / Beer-Lambert extinction law: transmission scales as T = T(1)**airmass.