Changes between Initial Version and Version 1 of PDAFlocal_set_increment_weights

Timestamp:

Sep 9, 2024, 2:54:20 PM (2 months ago)

Author:

lnerger

Comment:

--

PDAFlocal_set_increment_weights

@@ +1 @@
+= PDAFlocal_set_increment_weights =
+
+This page documents the routine `PDAFlocal_set_increment_weights` of PDAF, which is used with the [wiki:PDAFlocal] interface. This routine was introduced with PDAF V2.3. 
+
+This routine provides the optional functionality to prescribe the assimilation increment of each element of the local state vector a different weight. This can be used, e.g. to implement a vertical localization in the case that the local state vector is a full vertical column of the model grid. In this case, one can make the increment weight depending on the height (or depth) of a grid point. Another application is to implement weakly-coupled assimilation in which the local state vector contains all variables, but only a subset of them is updated. This is achieved by givening those element that should not be updated the weight 0.
+
+The routine can be called in `init_dim_l_pdaf` when the dimension of the local state vector has to be determined. In addition to this dimension, the implementation guide recommends to also initialize the coordinates of the local analysis domain and the indices of the elements ofthe local state vector in the non-localized (global or domain-decomposed) state vector `state_p`. In addition one can initialize a vector of increment weights and provide it to PDAFlocal by calling this routine.
+
+The interface is:
+{{{
+SUBROUTINE PDAFlocal_set_increment_weights(dim_l, weights)
+
+  INTEGER, INTENT(in) :: dim_l          ! Dimension of local state vector
+  REAL, INTENT(in) :: weights(dim_l)    ! Weights array
+}}}
+
+Hints:
+ * The loop to initialize `weights` can be analogous to the loop that initialize the index array for mapping between the global and local state vector (see [wiki:PDAFlocal_set_indices].
+ * It is possible to clear the weights array by calling [wiki:PDAFlocal_clear_increment_weights]. Afterwards, unit weights are used.