Weighting Functions

In the SASKTRAN framework, weighting functions are returned in the form:

\[\mathrm{wf[i, j, k]} = \frac{\partial I(\lambda_\mathrm{i}, \mathrm{L_j})}{\partial \mathrm{S_k}}\]

where

Symbol

Meaning

\(\lambda_\mathrm{i}\)

Wavelength \(\mathrm{i}\)

\(\mathrm{L_\mathrm{j}}\)

Line of sight \(\mathrm{j}\)

\(\mathrm{S_\mathrm{k}}\)

Species at altitude \(\mathrm{k}\)

Some engines might support additional forms but by default must return the form given above. To calculate weighting functions the sasktran.Atmosphere.wf_species must be set.

Weighting Function Altitudes

This is fairly self explanatory; weighting functions are calculated at specific altitudes. These altitudes are usually set via an engine option.

Weighting Function Widths

In multiple Engine’s the concept of weighting function widths is important.

_images/wf_equiv_fd.svg

Visual of what is meant by the weighting function altitude, and width for an equivalent finite difference weighting function calculation.

Weighting function widths refer to the region (around a weighting function altitude) in which a perturbation is applied. A different way of say this is that the width refers to equivalent region a finite difference method would need to perturb the species of interest to obtain the same result. It should be noted that the equivalent perturbation would be linearly interpolated from 0 at the region boundaries to the magnitude of the perturbation at the weighting function altitude. Below is an illustration to clarify.

Note

The user should take special care when setting the weighting function widths and altitudes. Note that the altitudes specify the center points and the widths specify the boundaries. Special care must be taken to ensure that the perturbation region does not extend above or below the top of the atmosphere (TOA) or surface respectively.

See Also

Engine

Objects that can calculate radiances. Additionally some engines can calculate analytic or optimized Weighting Functions.