ARTS built-in documentation server

Workspace Method yCalcAppend

Description

Replaces yCalc if a measurement shall be appended to an
existing one.

The method works basically as yCalc but appends the results to
existing data, instead of creating completely new y and its
associated variables. This method is required if your measurement
consists of data from two instruments using different observation
techniques (corresponding to different iyCalc-methods). One such
example is if emission and transmission data are combined into a
joint retrieval. The method can also be used to get around the
constrain that sensor_response is required to be the same for
all data.

The new measurement is simply appended to the input y, and the
other output variables are treated correspondingly. Data are
appended "blindly" in y_aux. That is, data of different type
are appended if iy_aux_vars differs between the two measurements,
the data are appended strictly following the order. First variable
of second measurement is appended to first variable of first
measurement, and so on. The number of auxiliary variables can differ
between the measurements. Missing data are set to zero.

The set of retrieval quantities can differ between the two
calculations. If an atmospheric quantity is part of both Jacobians,
the same retrieval grids must be used in both cases.
The treatment of instrument related Jacobians (baseline fits,
pointing ...) follows the *append_instrument_wfs* argument.

A difference to yCalc is that jacobian_quantities and
jacobian_indices are both in- and output variables. The input
version shall match the measurement to be calculated, while the
version matches the output y, the combined, measurements. Copies
of jacobian_quantities and * jacobian_indices* of the first
measurement must be made and shall be provided to the method as
*jacobian_quantities_copy* and *jacobian_indices_copy*.

Authors: Patrick Eriksson

Synopsis

yCalcAppend(

y, y_f, y_pol, y_pos, y_los, y_aux, jacobian, jacobian_quantities, jacobian_indices, atmgeom_checked, atmfields_checked, atmosphere_dim, t_field, z_field, vmr_field, cloudbox_on, cloudbox_checked, sensor_checked, stokes_dim, f_grid, sensor_pos, sensor_los, transmitter_pos, mblock_za_grid, mblock_aa_grid, antenna_dim, sensor_response, sensor_response_f, sensor_response_pol, sensor_response_za, sensor_response_aa, iy_main_agenda, jacobian_agenda, jacobian_do, iy_aux_vars, jacobian_quantities_copy, jacobian_indices_copy, append_instrument_wfs )

Variables

OUT+IN	y	(Vector)	The measurement vector.
OUT+IN	y_f	(Vector)	The frequencies associated with y.
OUT+IN	y_pol	(ArrayOfIndex)	The polarisation states associated with y.
OUT+IN	y_pos	(Matrix)	The sensor positions associated with y.
OUT+IN	y_los	(Matrix)	The line-of-sights associated with y.
OUT+IN	y_aux	(ArrayOfVector)	Data auxilary to y.
OUT+IN	jacobian	(Matrix)	The Jacobian matrix.
OUT+IN	jacobian_quantities	(ArrayOfRetrievalQuantity)	The retrieval quantities in the Jacobian matrix.
OUT+IN	jacobian_indices	(ArrayOfArrayOfIndex)	First and last column index in jacobian for each retrieval quantity.
IN	atmgeom_checked	(Index)	OK-flag for the geometry of the model atmosphere.
IN	atmfields_checked	(Index)	OK-flag for atmospheric grids and (physical) fields.
IN	atmosphere_dim	(Index)	The atmospheric dimensionality (1-3).
IN	t_field	(Tensor3)	The field of atmospheric temperatures.
IN	z_field	(Tensor3)	The field of geometrical altitudes.
IN	vmr_field	(Tensor4)	VMR field.
IN	cloudbox_on	(Index)	Flag to activate the cloud box.
IN	cloudbox_checked	(Index)	OK-flag for variables associated with the cloudbox.
IN	sensor_checked	(Index)	OK-flag for sensor related variables.
IN	stokes_dim	(Index)	The dimensionality of the Stokes vector (1-4).
IN	f_grid	(Vector)	The frequency grid for monochromatic pencil beam calculations.
IN	sensor_pos	(Matrix)	The sensor position for each measurement block.
IN	sensor_los	(Matrix)	The sensor line-of-sight (LOS) for each measurement block.
IN	transmitter_pos	(Matrix)	Transmitter positions.
IN	mblock_za_grid	(Vector)	The zenith angle grid for each measurement block.
IN	mblock_aa_grid	(Vector)	The azimuthal angle grid for each measurement block.
IN	antenna_dim	(Index)	The dimensionality of the antenna pattern (1-2).
IN	sensor_response	(Sparse)	The matrix modelling the total sensor response.
IN	sensor_response_f	(Vector)	The frequencies associated with the output of sensor_response.
IN	sensor_response_pol	(ArrayOfIndex)	The polarisation states associated with the output of sensor_response.
IN	sensor_response_za	(Vector)	The relative zenith angles associated with the output of sensor_response.
IN	sensor_response_aa	(Vector)	The relative azimuth angles associated with the output of sensor_response.
IN	iy_main_agenda	(Agenda)	See agendas.cc.
IN	jacobian_agenda	(Agenda)	See agendas.cc.
IN	jacobian_do	(Index)	Flag to activate jacobian calculations.
IN	iy_aux_vars	(ArrayOfString)	Selection of quantities for iy_aux and when applicable also y_aux.
GIN	jacobian_quantities_copy	(ArrayOfRetrievalQuantity)	Copy of jacobian_quantities of first measurement.
GIN	jacobian_indices_copy	(ArrayOfArrayOfIndex)	Copy of jacobian_indices of first measurement.
GIN	append_instrument_wfs	(Index, Default: 0)	Flag controlling if instrumental weighting functions are appended or treated as different retrieval quantities.