Changes between Version 4 and Version 5 of prepoststep_pdaf

Timestamp:

Jan 21, 2015, 2:56:00 PM (9 years ago)

Author:

lnerger

Comment:

--

prepoststep_pdaf

@@ -26 +26 @@
  * If parallelization with domain decomposition is used, the variables `state_p` and `ens_p` will only contain the state information for the process-local domain. Accordingly, also `dim_p` and `dim_obs_p` only contain the state vector size for the local domain. 
  * For the local filters (LESTKF, LETKF, LSEIK) `dim_obs_p` is the observation dimension for the `full` observations.
+ * The array `Uinv` is only needed for SEEK filter. Also one needs its values for the ESTKF, ETKF, and SEIK filter to compute the analysis error covariance matrix when these filters used with a fixed-covariance matrix (subtype=3). In all other cases the ensemble array `ens_p` contains all required information.
 
 Hints:
  * If a user considers to perform adjustments to the estimates (e.g. for balances), this routine is the right place for it.
- * Only for the SEEK filter the state vector (`state_p`) is initialized. For all other filters, the array is allocated, but it can be used freely during the execution of `prepoststep_pdaf`.
+ * Only for the SEEK filter the state vector (`state_p`) is initialized. For all other filters, the array is allocated, but it can be used freely during the execution of `prepoststep_pdaf`. FOr example one can use the array to store the ensemble mean without allocating a new array.
  * The interface through which `prepoststep_pdaf` is called does not include the array of smoothed ensembles. In order to access the smoother ensemble array one has to set a pointer to it using a call to the routine `PDAF_get_smootherens` (see page on [AuxiliaryRoutines auxiliary routines])