Changes between Version 36 and Version 37 of InitPdaf
- Timestamp:
- Apr 12, 2012, 11:44:33 AM (13 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
InitPdaf
v36 v37 50 50 * 4: ETKF 51 51 * 5: LETKF 52 * 6: ESTKF 53 * 7: LESTKF 52 54 * `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. 53 55 * `step_null`: An integer defining the initial time step. For some cases it can use useful to set `step_null` larger to 0. … … 111 113 === Initialization for ensemble-based filters === 112 114 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.115 Generally, 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. 114 116 115 117 The 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. … … 169 171 }}} 170 172 171 The correctness of the ensemble initialization in `U_ ensemble_init` should be checked by the user.173 The correctness of the ensemble initialization in `U_init_ens` should be checked by the user.