Changes between Version 32 and Version 33 of ImplementAnalysisseik
- Timestamp:
- Sep 15, 2011, 2:28:28 PM (13 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
ImplementAnalysisseik
v32 v33 62 62 63 63 64 === `U_collect_state` (collect_state .F90) ===64 === `U_collect_state` (collect_state_pdaf.F90) === 65 65 66 66 This routine is independent of the filter algorithm used. … … 68 68 69 69 70 === `U_init_dim_obs` (init_dim_obs .F90) ===70 === `U_init_dim_obs` (init_dim_obs_pdaf.F90) === 71 71 72 72 This routine is used by all global filter algorithms (SEEK, SEIK, EnKF, ETKF). … … 86 86 87 87 88 === `U_obs_op` (obs_op .F90) ===88 === `U_obs_op` (obs_op_pdaf.F90) === 89 89 90 90 This routine is used by all global filter algorithms (SEEK, SEIK, EnKF, ETKF). … … 109 109 110 110 111 === `U_init_obs` (init_obs .F90) ===111 === `U_init_obs` (init_obs_pdaf.F90) === 112 112 113 113 This routine is used by all global filter algorithms (SEEK, SEIK, EnKF, ETKF). … … 128 128 129 129 130 === `U_prepoststep` (prepoststep_seik .F90) ===130 === `U_prepoststep` (prepoststep_seik_pdaf.F90) === 131 131 132 132 The routine has already been described on the [ModifyModelforEnsembleIntegration#U_prepoststepprepoststep_seik.F90 page on modifying the model code for the ensemble integration]. For completeness, the description is repeated: … … 163 163 164 164 165 === `U_prodRinvA` (prodrinva .F90) ===165 === `U_prodRinvA` (prodrinva_pdaf.F90) === 166 166 167 167 This routine is used by all filter algorithms that use the inverse of the observation error covariance matrix (SEEK, SEIK, and ETKF). … … 188 188 * The interface has a difference for SEIK and ETKF: For ETKF the third argument is the ensemble size (`dim_ens`), while for SEIK it is the rank (`rank`) of the covariance matrix (usually ensemble size minus one). In addition, the second dimension of `A_p` and `C_p` has size `dim_ens` for ETKF, while it is `rank` for the SEIK filter. (Practically, one can usually ignore this difference as the fourth argument of the interface can be named arbitrarily in the routine.) 189 189 190 === `U_init_obsvar` (init_obsvar .F90) ===190 === `U_init_obsvar` (init_obsvar_pdaf.F90) === 191 191 192 192 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).