wiki:PDAFomi_obs_op_interp_lin

Version 1 (modified by lnerger, 3 years ago) (diff)

--

PDAFomi_obs_op_interp_lin

This page documents the routine PDAFomi_obs_op_interp_lin of PDAF-OMI.

The routine provides an observation operator for observations that are computed as the linear interpolation between some grid point values in the state vector. It is usually called in obs_op_OBSTYPE in an OMI observation module.

See the page describing PDAF-OMI observation operators for a general overview.

The interface is:

  SUBROUTINE PDAFomi_obs_op_interp_lin(thisobs, nrows, state_p, obs_f_all)

    TYPE(obs_f), INTENT(inout) :: thisobs  ! Data type with full observation
    INTEGER, INTENT(in) :: nrows           ! Number of values to be averaged
    REAL, INTENT(in)    :: state_p(:)      ! Process-local model state provided by PDAF
    REAL, INTENT(inout) :: obs_f_all(:)    ! Full observed state for all observation types (array provided by PDAF)

Note:

  • The necessary index information for the observation operator was initialized in init_obs_OBSTYPE of the observation module.
  • The interpolation coefficients has to be initialized in init_ob_OBSTYPE of the observation module.