Changes between Version 12 and Version 13 of ImplementAnalysisLocal

Timestamp:

May 12, 2021, 5:32:01 PM (4 years ago)

Author:

lnerger

Comment:

--

ImplementAnalysisLocal

@@ -157 +157 @@
 
 The routine is called during the loop over the local analysis domains in the analysis step. 
-It has to provide in `dim_l` the dimension of the state vector for the local analysis domain with index `domain_p`.
+It has to provide in `dim_l` the dimension of the state vector for the local analysis domain with index `domain_p`. 
 
 Hints: 
- * If a local analysis domain is a single vertical column of the model grid, the size of the state in the local analysis domain will be just the number of vertical grid points at this location.
+ * For sharing through the module 'mod_assimilation', we further initialize and array 'coords_l' containing the coordinates that describe the local domain. These coordinates have to describe one location in space and is used in the OMI observation modules to compute the distance from observations.
+ * Any form of local domain is possible as long as it can be describe as a single location. If observations are only horizontally distributes (a common situation with satellite data in the ocean), the local analysis domain can be single vertical column of the model grid. In this case, the size of the state in the local analysis domain will be just the number of vertical grid points at this location and the horizontal coordinates are used in 'coords_l'
+ * Further, we recommend to initialize an array containing the indices of the elements of the local state vector in the global (or domain-decomposed) state vector (`id_lstatae_in_pstate` in the templates). This array is also shared through 'mod_assimilation'.