wiki:g2l_obs_pdaf

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

--

g2l_obs_pdaf

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

The routine g2l_obs_pdaf (called U_g2l_obs inside the PDAF core routines) is a call-back routine that has to be provided by the user. The routine is used with all filter algorithms using domain-localization (LSEIK, LETKF, LESTKF) and is independent of the particular algorithm. The routine is called during the loop over the local analysis domains in the analysis step. It has to provide a local observation vector mstate_l for the observation domain that corresponds to the local analysis domain with index domain_p. Provided to the routine is the full observation vector mstate_f from which the local part has to be extracted.

The interface is the following:

SUBROUTINE g2l_obs_pdaf(domain_p, step, dim_obs_f, dim_obs_l, mstate_f, mstate_l)

with

  • domain_p : integer, intent(in)
    Index of current local analysis domain
  • step : integer, intent(in)
    Current time step
  • dim_obs_f : integer, intent(in)
    Size of full observation vector for model sub-domain
  • dim_obs_l : integer, intent(in)
    Size of observation vector for local analysis domain
  • mstate_f : integer, intent(in), dimension(dim_p)
    Full observation vector for model sub-domain
  • mstate_l : integer, intent(out), dimension(dim_l)
    Observation vector for local analysis domain

Hints:

  • The vector mstate_f that is provided to the routine is one of the observed state vectors that are produced by obs_op_f_pdaf.
  • Some operations performed here are analogous to those required to initialize a local vector of observations in init_obs_l_pdaf. If that routine reads first a full vector of observations (e.g. in init_dim_obs_f_pdaf), this vector has to be restricted to the relevant observations for the current local analysis domain. For this operation, one can for example initialize an index array when init_dim_obs_l_pdaf is executed. (Which happens before g2l_obs_pdaf)