Changes between Initial Version and Version 1 of g2l_obs_pdaf

Jan 21, 2016, 6:38:52 PM (7 years ago)



  • g2l_obs_pdaf

    v1 v1  
     1= g2l_obs_pdaf =
     3The page document the user-supplied call-back routine `g2l_obs_pdaf`.
     5The 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.
     6The routine is used with all filter algorithms using domain-localization (LSEIK, LETKF, LESTKF) and is independent of the particular algorithm.
     7The 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.
     9The interface is the following:
     11SUBROUTINE g2l_obs_pdaf(domain_p, step, dim_obs_f, dim_obs_l, mstate_f, mstate_l)
     14 * `domain_p` : `integer, intent(in)`[[BR]] Index of current local analysis domain
     15 * `step` : `integer, intent(in)`[[BR]] Current time step
     16 * `dim_obs_f` : `integer, intent(in)`[[BR]] Size of full observation vector for model sub-domain
     17 * `dim_obs_l` : `integer, intent(in)`[[BR]] Size of observation vector for local analysis domain
     18 * `mstate_f` : `integer, intent(in), dimension(dim_p)`[[BR]] Full observation vector for model sub-domain
     19 * `mstate_l` : `integer, intent(out), dimension(dim_l)`[[BR]] Observation vector for local analysis domain
     22 * 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`.
     23 * 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`)