A full end-to-end modular software suite for radiative transfer and related atmospheric analysis
42 Stars
Updated Last
9 Months Ago
Started In
July 2020


An end-to-end modular software suite for vectorized atmospheric radiative transfer calculations, based on the Matrix Operator Method.

Written in Julia.

Tests Docs License Github Commit Frequency

InstallationModules (vSmartMOM, Absorption, Scattering) • SupportLicense

This project aims to revamp and modernize key atmospheric remote sensing tools. Specifically, it will enable the fast computation of atmospheric optical properties, full-polarized radiative transfer simulations, and commonly-used inversion routines.

By taking advantage of modern software tools, such as GPU acceleration and HPC computing, the software suite significantly accelerates computationally-intensive calculations and models, while keeping the interface easy-to-use for researchers and students.


vSmartMOM can be installed using the Julia package manager. From the Julia REPL, type ] to enter the Pkg REPL mode and run

pkg> add vSmartMOM


Note: This section provides only a quick overview of the available modules in vSmartMOM.jl.

For in-depth examples, tutorials, and implementation details, please see the complete Documentation.


The vSmartMOM module allows end-to-end simulation of radiative transfer (RT) throughout Earth's atmosphere and surface. Specifically, it:

  1. Enables 1D vectorized plane-parallel RT modeling based on the Matrix Operator Method.
  2. Incorporates fast, high fidelity simulations of scattering atmospheres containing haze and clouds – including pressure- and temperature-resolved absorption profiles of gaseous species in the atmosphere.
  3. Enables GPU-accelerated computations of the resulting hyperspectral reflectances/transmittances.

Key functions:

  • parameters_from_yaml(filepath::String): Load a custom set of RT parameters from a YAML file.
  • default_parameters(): Load a default set of RT parameters.
  • model_from_parameters(parameters::vSmartMOM_Parameters): Using the parameters, calculate derived parameters that will be used in the main RT code. The derived parameters include cross-section profiles, scattering phase functions, etc.
  • rt_run(model::vSmartMOM_Model): Use the defined model to perform RT simulation.


This module enables absorption cross-section calculations of atmospheric gases at different pressures, temperatures, and broadeners (Doppler, Lorentzian, Voigt). It uses the HITRAN energy transition database for calculations. While it enables lineshape calculations from scratch, it also allows users to create and save an interpolator object at specified wavelength, pressure, and temperature grids. It can perform these computations either on CPU or GPU.

Key functions:

  • read_hitran(filepath::String): Creates a HitranTable struct from the fixed-width HITRAN file with transitions.
  • make_hitran_model(hitran::HitranTable, broadening::AbstractBroadeningFunction, ...): Create a HitranModel struct that holds all of the model parameters needed to perform a absorption cross-section (transitions, broadening type, wing_cutoff, etc.)
  • make_interpolation_model(hitran::HitranTable, broadening::AbstractBroadeningFunction, ): Similar to creating a HitranModel, but this will perform the interpolation at the given wavelength, pressure, and temperature grids and store the interpolator in InterpolationModel.
  • absorption_cross_section(model::AbstractCrossSectionModel, grid::AbstractRange{<:Real}, pressure::Real, temperature::Real, ...): Performs an absorption cross-section calculation with the given model (HitranModel or InterpolationModel), at a given wavelength grid, pressure and temperature


This module enables scattering phase-function calculation of atmospheric aerosols with different size distributions, incident wavelengths, and refractive indices. It can perform the calculation using either the Siewert NAI-2 or Domke PCW methods (Suniti Sanghavi 2014).

Key functions:

  • make_univariate_aerosol(size_distribution::ContinuousUnivariateDistribution, r_max, nquad_radius::Int, nᵣ, nᵢ: Create an aerosol object with size distribution and complex refractive index.
  • make_mie_model(computation_type::AbstractFourierDecompositionType, aerosol::AbstractAerosolType, λ::Real, polarization::AbstractPolarizationType, truncation_type::AbstractTruncationType, ...): Create a MieModel struct that holds all of the model parameters needed to perform a phase function calculation (computation type, aerosol, incident wavelength, etc. )
  • compute_aerosol_optical_properties(model::MieModel): Compute the aerosol optical properties using the specified model parameters

How to Contribute

vSmartMOM.jl is a growing package and thus feedback from users like you are highly appreciated. To report bugs or suggest new features in vSmartMOM.jl, please create GitHub Issues. To contribute to the package, please feel free to create a Pull Request.

If you have any questions about the methods used or would like to chat with us, please feel free to shoot us an email here.


This project is being developed in the Christian Frankenberg and Paul Wennberg labs at Caltech and is largely based on papers and ideas by Suniti Sanghavi from NASA/JPL, with support from the Schmidt Academy for Software Engineering (SASE).

Copyright Notice

Apache 2.0 License; Copyright 2022, by the California Institute of Technology. United States Government Sponsorship acknowledged.