Changes between Version 9 and Version 10 of ImplementAnalysisseik


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

--

Legend:

Unmodified
Added
Removed
Modified
  • ImplementAnalysisseik

    v9 v10  
    9292The 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`.
    9393
     94For a model with domain decomposition, `A_p` contains the part of the matrix that resides on the model sub-domain of the calling process. The product has to be computed for this sub-domain, too.
    9495
    9596Hints:
     
    101102== `U_init_obsvar` (init_obsvar.F90) ==
    102103
     104This routine is used by the global filter algorithms SEIK and  ETKF as well as the local filters LSEIK and LETKF. The routine is only called if the adaptive forgetting factor is used (`type_forget=1` in the example impementation).
     105
     106The interface for this routine is:
    103107{{{
    104108SUBROUTINE init_obsvar(step, dim_obs_p, obs_p, meanvar)
     
    110114}}}
    111115
    112 ! User-supplied routine for PDAF (SEIK/LSEIK/ETKF/LETKF)
    113 ! with adaptive forgetting factor. This routine will
    114 ! only be called, if the adaptive forgetting factor
    115 ! feature is used. Please note that this is an
    116 ! experimental feature.
    117 !
    118 ! The routine is called in SEIK during the analysis or
    119 ! in LSEIK before the loop over local analysis domains
    120 ! by the routine PDAF\_set\_forget that estimates an
    121 ! adaptive forgetting factor.  The routine has to
    122 ! initialize the mean observation error variance. 
    123 ! For SEIK this should be the global mean, while for
    124 ! LSEIK it should be the mean for the PE-local
    125 ! sub-domain.  (See init\_obsvar\_local() for a
    126 ! localized variant for LSEIK.)
     116The routine is called in the global filters during the analysis or
     117by the routine that computes an adaptive forgetting factor (PDAF_set_forget).
     118The routine has to initialize the mean observation error variance. 
     119For the global filters this should be the global mean.
     120
     121Hin