Changes between Version 38 and Version 39 of InitPdaf
- Timestamp:
- Apr 28, 2014, 9:40:59 AM (11 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
InitPdaf
v38 v39 111 111 * `flag`: Status flag for PDAF. It is 0 upon entry and can be set by in the user-supplied routine, depending on the success of the ensemble initialization. Preferably, values above 102 should be used for failures to avoid conflicts with the error codes defined within PDAF_init. 112 112 113 In the initialization routine `U_init_ens_pdaf` one has to distinguish between ensemble-based and mode-based filters as will be described below. The only mode based filter supplied with PDAF is SEEK, while all other methods are ensemble-based.114 115 113 === Defining the state vector === 116 114 … … 122 120 === Initialization for ensemble-based filters === 123 121 124 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 parallel model with domain decomposition is used, the full ensemble for the local sub-domain of the MPI process has to be initialized.122 In the initialization routine `U_init_ens_pdaf` one has to distinguish between ensemble-based and mode-based filters. The only mode based filter supplied with PDAF is SEEK, while all other methods are ensemble-based. 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 parallel model with domain decomposition is used, the full ensemble for the local sub-domain of the MPI process has to be initialized. 125 123 126 124 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 used for EnKF.