jacobianAddSinefit

Workspace.jacobianAddSinefit(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, sensor_response_pol_grid: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.ArrayOfIndex]] = self.sensor_response_pol_grid, sensor_response_dlos_grid: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Matrix]] = self.sensor_response_dlos_grid, sensor_pos: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Matrix]] = self.sensor_pos, period_lengths: Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Vector], no_pol_variation: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Index]] = 0, no_los_variation: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Index]] = 0, no_mblock_variation: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Index]] = 0, verbosity: Optional[Union[pyarts.arts.WorkspaceVariable, pyarts.arts.Verbosity]] = self.verbosity) None

Includes sinusoidal baseline fit in the Jacobian.

Works as jacobianAddPolyfit(), beside that a series of sine and cosine terms are used for the baseline fit.

For each value in period_lengths one sine and one cosine term are included (in mentioned order). By these two terms the amplitude and “phase” for each period length can be determined. The sine and cosine terms have value 0 and 1, respectively, for first frequency.

If the simulation/retrieval deals with a single spectrum, the number of elements added to the state vector (x) is 2 * nperiods, where nperiods is the length of period_lengths. The same is true if no_pol_variation, no_los_variation and no_mblock_variation all are set to 1, even if several spectra are involved. Otherwise thenumber of elements added to x depends on the number of spectra and the settings of no_pol_variation, no_los_variation and no_mblock_variation. The sine and cosine terms for each period length are treated as a separate retrieval quantities. That is, the the elements associated with the first period length are grouped and form together a retrieval quantity, etc. Inside each retrieval quantity the pairs of sine and cosine terms are kept together, in given order.

Author(s): 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]

  • sensor_response_pol_grid (ArrayOfIndex, optional) – The “polarisation grid” associated with sensor_response. See sensor_response_pol_grid, defaults to self.sensor_response_pol_grid [IN]

  • sensor_response_dlos_grid (Matrix, optional) – The zenith and azimuth angles associated with sensor_response. See sensor_response_dlos_grid, defaults to self.sensor_response_dlos_grid [IN]

  • sensor_pos (Matrix, optional) – The sensor position for each measurement block. See sensor_pos, defaults to self.sensor_pos [IN]

  • period_lengths (Vector) – Period lengths of the fit. [IN]

  • no_pol_variation (Index, optional) – Set to 1 if the baseline off-set is the same for all Stokes components. Defaults to 0 [IN]

  • no_los_variation (Index, optional) – Set to 1 if the baseline off-set is the same for all line-of-sights (inside each measurement block). Defaults to 0 [IN]

  • no_mblock_variation (Index, optional) – Set to 1 if the baseline off-set is the same for all measurement blocks. Defaults to 0 [IN]

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