Changes between Version 6 and Version 7 of LnDevel
- Timestamp:
- Sep 3, 2010, 3:33:33 PM (14 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
LnDevel
v6 v7 58 58 59 59 60 == `U_init_dim_obs` (init_dim_obs.F90) == 60 == User-supplied routines == 61 62 Here all user-supplied routines are described that are required in the call to `PDAF_put_state_seik`. For some of the generic routines, we link to the page on [ModifyModelforEnsembleIntegration modifying the model code for the ensemble integration]. 63 64 To indicate user-supplied routines we use the prefix `U_`. In the template directory `templates/` these routines are provided in files with the routines name without this prefix. In the example implementation in `testsuite/src/dummymodel_1D` the routines exist without the prefix, but with the extension `_dummy_D.F90`. In the section titles below we provide the name of the template file in parentheses. 65 66 67 === `U_collect_state` (collect_state.F90) === 68 69 This routine is independent from the filter algorithm used. 70 See [ModifyModelforEnsembleIntegration#U_collect_statecollect_state.F90 here] for the description of this routine. 71 72 73 === `U_init_dim_obs` (init_dim_obs.F90) === 61 74 62 75 This routine is used by all global filter algorithms (SEEK, SEIK, EnKF, ETKF). … … 74 87 Some hints: 75 88 * It can be useful to not only determine the size of the observation vector is determined at this point. One can also already gather information about the locations of the observations, which will be used later, e.g. to implement the observation operator. An array for the locations can be defined in a module like `mod_assimilation`. 89 76 90 77 91 == `U_obs_op` (obs_op.F90) == … … 97 111 * If the observation operator involves a global operation, e.g. some global integration, while using domain-decompostion one has to gather the information from the other model domains using MPI communication. 98 112 113 99 114 == `U_init_obs` (init_obs.F90) == 100 115 … … 114 129 115 130 For a model using domain decomposition, the vector of observations that exist on the model sub-domain for the calling process has to be initialized. 131 132 133 === `U_prepoststep` (prepoststep_seik.F90) === 134 135 See [ModifyModelforEnsembleIntegration#U_prepoststepprepoststep_seik.F90 here] for the description of this routine. 116 136 117 137 … … 141 161 142 162 143 144 163 == `U_init_obsvar` (init_obsvar.F90) == 145 164