jacobianAddAbsSpecies

Workspace.jacobianAddAbsSpecies(self: pyarts.arts._Workspace, jacobian_quantities: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.ArrayOfRetrievalQuantity]] = self.jacobian_quantities, jacobian_agenda: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Agenda]] = self.jacobian_agenda, atmosphere_dim: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Index]] = self.atmosphere_dim, p_grid: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Vector]] = self.p_grid, lat_grid: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Vector]] = self.lat_grid, lon_grid: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Vector]] = self.lon_grid, g1: Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Vector], g2: Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Vector], g3: Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Vector], species: Union[pyarts.arts.WorkspaceVariable, pyarts.arts.String], unit: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.String]] = "vmr", for_species_tag: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Index]] = 1, verbosity: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Verbosity]] = self.verbosity) → None

Includes an absorption species in the Jacobian.

For 1D or 2D calculations the latitude and/or longitude grid of the retrieval field should set to have zero length.

These retrieval units are at hand for all gas species:

  • "vmr": Volume mixing ratio.

  • "nd": Number density.

  • "rel": Relative unit (e.g. 1.1 means 10% more of the gas).

For water vapour, also these units are at hand:

  • "rh": Relative humidity.

  • "q": Specific humidity.

Note that for_species_tag is used to indicate if species tag VMR, rather than atmospheric gas VMR is calculated. Set it to 0 and we calculate the atmospheric gas VMR, but this only works for “analytical”.

Note that the Jacobian is set to zero where volume mixing ratio equals zero.

The number of elements added to the state vector (x) is:

n_g1 * n_g2 * n_g3

where n_g1, n_g2 and n_g3 are the length of GIN g1, g2 and g3, respectively. Here empty vectors should be considered to have a length 1. The elements are sorted with pressure as innermost loop, followed by latitude and longitude as outermost loop.

Author(s): Mattias Ekstrom, Patrick Eriksson

Parameters:

  • jacobian_quantities (ArrayOfRetrievalQuantity, optional) – The retrieval quantities in the Jacobian matrix. See jacobian_quantities, defaults to self.jacobian_quantities [INOUT]

  • jacobian_agenda (Agenda, optional) – Pure numerical Jacobian calculations. See jacobian_agenda, defaults to self.jacobian_agenda [INOUT]

  • atmosphere_dim (Index, optional) – The atmospheric dimensionality (1-3). See atmosphere_dim, defaults to self.atmosphere_dim [IN]

  • p_grid (Vector, optional) – The pressure grid. See p_grid, defaults to self.p_grid [IN]

  • lat_grid (Vector, optional) – The latitude grid. See lat_grid, defaults to self.lat_grid [IN]

  • lon_grid (Vector, optional) – The longitude grid. See lon_grid, defaults to self.lon_grid [IN]

  • g1 (Vector) – Pressure retrieval grid. [IN]

  • g2 (Vector) – Latitude retrieval grid. [IN]

  • g3 (Vector) – Longitude retreival grid. [IN]

  • species (String) – The species tag of the retrieval quantity. [IN]

  • unit (String, optional) – Retrieval unit. See above. Defaults to "vmr" [IN]

  • for_species_tag (Index, optional) – Index-bool for acting on species tags or species. Defaults to 1 [IN]

  • verbosity (Verbosity) – ARTS verbosity. See verbosity, defaults to self.verbosity [IN]