Changes between Version 36 and Version 37 of InitPdaf


Ignore:
Timestamp:
Apr 12, 2012, 11:44:33 AM (12 years ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • InitPdaf

    v36 v37  
    5050  * 4: ETKF
    5151  * 5: LETKF
     52  * 6: ESTKF
     53  * 7: LESTKF
    5254 * `subtype`: An integer defining the sub-type of the filter algorithm (see the example code in `testsuite/src/dummymodel_1D` for choices). If `PDAF_init` is called with `subtype=-1` the available options are shown for the selected filter algorithm.
    5355 * `step_null`: An integer defining the initial time step. For some cases it can use useful to set `step_null` larger to 0.
     
    111113=== Initialization for ensemble-based filters ===
    112114
    113 Generally, the filters SEIK, LSEIK, EnKF, ETKF, and LETKF are ensemble-based filters.For these filters only the array `ens_p` needs to be initialized by the ensemble of model states. If a model with domain decomposition is used, the full ensemble for the local sub-domain of the MPI process has to be initialized.
     115Generally, the filters SEIK, LSEIK, EnKF, ETKF, LETKF, ESTKF, and LESTKF are ensemble-based filters. For these filters only the array `ens_p` needs to be initialized by the ensemble of model states. If a model with domain decomposition is used, the full ensemble for the local sub-domain of the MPI process has to be initialized.
    114116
    115117The arrays `state_p` and `Uinv` are allocated to their correct sizes because they are used during the assimilation cycles. They are not yet initialized and it is allowed to use these arrays in the initialization. An exception from this is EnKF for which `Uinv` is allocated only with size (`1`,`1`), because `Uinv` is not using for EnKF.
     
    169171}}}
    170172
    171 The correctness of the ensemble initialization in `U_ensemble_init` should be checked by the user.
     173The correctness of the ensemble initialization in `U_init_ens` should be checked by the user.