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. |