Changes between Initial Version and Version 1 of init_dim_obs_l_pdaf

Jan 21, 2016, 3:37:37 PM (5 years ago)



  • init_dim_obs_l_pdaf

    v1 v1  
     1= init_dim_obs_l_pdaf =
     3The page document the user-supplied call-back routine `init_dim_obs_l_pdaf`.
     5The routine `init_dim_obs_l_pdaf` (called `U_init_dim_obs_l` 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 initialize the size of the observation vector used for the specified local analysis domain.
     9The interface is the following:
     11SUBROUTINE init_dim_obs_l_pdaf(domain_p, step, dim_obs_f, dim_obs_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 the full observation vector
     17 * `dim_obs_l` : `integer, intent(out)`[[BR]] Local dimension of observation vector
     19Some hints:
     20 * Usually, the observations to be considered for a local analysis are those which reside within some distance from the local analysis domain. Thus, if the local analysis domain is a single model grid point and if the model grid is a regular ij-grid, then one could use some range of i/j indices to select the observations and determine the local number of them. More generally, one can compute the physical distance of an observation from the local analysis domain and decide on this basis, if the observation has to be considered.
     21 * In the loop over the local analysis domains, the routine is always called before `init_obs_l_pdaf` is executed. Thus, as `init_dim_obs_l_pdaf` has to check which observations should be used for the local analysis domain, one can already initialize an integer array that stores the index of observations to be considered. This index should be the position of the observation in the array `observation_f`. With this, the initialization of the local observation vector in `init_obs_l_pdaf` can be speeded up.
     22 * For PDAF, we could not join the routines `init_dim_obs_l_pdaf` and `init_obs_l_pdaf`, because the array for the local observations is allocated internally of PDAF after its size has been determined in `init_dim_obs_l_pdaf`.