Changes between Version 8 and Version 9 of ImplementAnalysisseik

Timestamp:

Aug 27, 2010, 9:35:02 AM (14 years ago)

Author:

lnerger

Comment:

--

ImplementAnalysisseik

@@ -60 +60 @@
 
 The interface for this routine is:
-
 {{{
 SUBROUTINE init_obs(step, dim_obs_p, observation_p)
@@ -77 +76 @@
 == `U_prodRinvA` (prodrinva.F90) ==
 
+This routine is used by all filters whose algorithm uses the inverse of the observation error covariance matrix (SEEK, SEIK, and ETKF).
+
+The interface for this routine is:
 {{{
 SUBROUTINE prodRinvA(step, dim_obs_p, rank, obs_p, A_p, C_p)
@@ -88 +90 @@
 }}}
 
-! The routine is called during the analysis step.
-! It has to compute the product of the inverse of 
-! the observation error covariance matrix with
-! the matrix of observed EOF modes (SEEK) or 
-! observed ensemble perturbations (SEIK/ETKF).
+The routine is called during the analysis step. In the algorithms the product of the inverse of the observation error covariance matrix with some matrix has to be computed. For the SEIK filter this matrix holds the observed part of the ensemble perturbations ('''HL'''). The matrix is provided as `A_p`. The product has to be given as `C_p`.
+
+
+Hints:
+ * the routine does not require that the product is implemented as a real matrix-matrix product. Rather, the product can be implemented in its most efficient form. For example, if the observation error covariance matrix is diagonal, only the multiplication of the diagonal with matrix `A_p` has to be implemented.
+ * The observation vector `obs_p` is provided through the interface for cases where the observation error variance is relative to the actual value of the observations.