Changes between Version 8 and Version 9 of ImplementAnalysisseik


Ignore:
Timestamp:
Aug 27, 2010, 9:35:02 AM (10 years ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • ImplementAnalysisseik

    v8 v9  
    6060
    6161The interface for this routine is:
    62 
    6362{{{
    6463SUBROUTINE init_obs(step, dim_obs_p, observation_p)
     
    7776== `U_prodRinvA` (prodrinva.F90) ==
    7877
     78This routine is used by all filters whose algorithm uses the inverse of the observation error covariance matrix (SEEK, SEIK, and ETKF).
     79
     80The interface for this routine is:
    7981{{{
    8082SUBROUTINE prodRinvA(step, dim_obs_p, rank, obs_p, A_p, C_p)
     
    8890}}}
    8991
    90 ! The routine is called during the analysis step.
    91 ! It has to compute the product of the inverse of
    92 ! the observation error covariance matrix with
    93 ! the matrix of observed EOF modes (SEEK) or
    94 ! observed ensemble perturbations (SEIK/ETKF).
     92The 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`.
     93
     94
     95Hints:
     96 * 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.
     97 * 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.
    9598
    9699