Changes between Version 9 and Version 10 of ImplementAnalysisseik
 Timestamp:
 Aug 27, 2010, 9:39:34 AM (10 years ago)
Legend:
 Unmodified
 Added
 Removed
 Modified

ImplementAnalysisseik
v9 v10 92 92 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`. 93 93 94 For a model with domain decomposition, `A_p` contains the part of the matrix that resides on the model subdomain of the calling process. The product has to be computed for this subdomain, too. 94 95 95 96 Hints: … … 101 102 == `U_init_obsvar` (init_obsvar.F90) == 102 103 104 This 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 106 The interface for this routine is: 103 107 {{{ 104 108 SUBROUTINE init_obsvar(step, dim_obs_p, obs_p, meanvar) … … 110 114 }}} 111 115 112 ! Usersupplied 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 PElocal 125 ! subdomain. (See init\_obsvar\_local() for a 126 ! localized variant for LSEIK.) 116 The routine is called in the global filters during the analysis or 117 by the routine that computes an adaptive forgetting factor (PDAF_set_forget). 118 The routine has to initialize the mean observation error variance. 119 For the global filters this should be the global mean. 120 121 Hin