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



The page document the user-supplied call-back routine U_add_obs_err.

The routine U_add_obs_err is a call-back routine that has to be provided by the user. In the simplified interface the predefined name of the routine is add_obs_error_pdaf, but in the full interface, the user can choose the name of the routine. The routine is used by the EnKF and called during the filter analysis step. The purpose of the routine is to add the observation error covariance matrix to the input/output matrix C. This matrix is the projection of the ensemble covariance matrix onto the observation space that is computed during the analysis step of the EnKF.

The interface is the following:

SUBROUTINE U_add_obs_err(step, dim_obs, C)


  • step : integer, intent(in)
    Current time step
  • dim_obs : integer, intent(in)
    Size of observation vector
  • C : real, intent(inout), dimension(dim_obs,dim_obs)
    Matrix to which the observation error covariance matrix is added

The routine is called during the analysis step.

The operation is for the global observation space. Thus, it is independent of whether the filter is executed with or without parallelization.


  • The routine does not require that the observation error covariance matrix is added as a full matrix. If the matrix is diagonal, only the diagonal elements have to be added.