PDAFomi_obs_op Module

PDAF-OMI observation operators

This module contains generic routines for several observation operators to be used after preparation with init_dim_obs_f

The observation operators are:

• PDAFomi_obs_op_gridpoint\n Observation operator for data at grid points. The routine selects values of the state vector according to an index array
• PDAFomi_obs_op_gridavg\n Observation operator for the case that the observations are the average of grid point values. The routine computes these averages according to an index array.
• PDAFomi_obs_op_interp_lin\n Observation operator for the case that the observations are linear interpolated from the grid points. The interpolation coefficients are pre-computed.
• PDAFomi_obs_op_gatheronly\n Observation operator for the case of strongly coupled assimilation to gather an observation which only exists in other compartments.

Adjoint observation operators:

• PDAFomi_obs_op_adj_gridpoint\n Adjoint observation operator for data at grid points. The routine selects values of the state vector according to an index array
• PDAFomi_obs_op_adj_gridavg\n Adjoint observation operator for the case that the observations are the average of grid point values. The routine computes these averages according to an index array.
• PDAFomi_obs_op_adj_interp_lin\n Adjoint observation operator for the case that the observations are linear interpolatied from the grid points. The interpolation coefficients are pre-computed.
• PDAFomi_obs_op_adj_gatheronly\n Adjoint observation operator for the case of strongly coupled assimilation to gather an observation which only exists in other compartments.

Helper routines for the operators: * PDAFomi_get_interp_coeff_tri \n Routine to compute interpolation coefficients for triangular interpolation from barycentric coordinates. * PDAFomi_get_interp_coeff_lin1D \n Routine to comput linear interpolation in 1D * PDAFomi_get_interp_coeff_lin \n Routine to compute interpolation coefficients for linear interpolations (linear, bi-linear, tri-linear)