Changes between Version 12 and Version 13 of ImplementGenerateObs
- Timestamp:
- Feb 22, 2023, 2:07:14 PM (21 months ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
ImplementGenerateObs
v12 v13 130 130 To perform a twin experiment using the synthetic observations generated by PDAF, one runs the data assimilation as one would with real observations. If one already initializes the vector of actual observations in the routines `init_dim_obs_TYPE` in the observation modules one only needs a small modification of this routine. Namely, the only required modification is that at the end of `init_dim_obs_TYPE` one overwrites the vector of real observations with the values from the synthetic observations. If one uses the template file `readwrite_obs.F90` for this, one can use `read_syn_obs` from this file at the end of `init_dim_obs_TYPE` to overwrite the observation vector. To allow for a flexible switching between the case using real observations and the twin experiment, one can for example introduce a flag `twin_experiment` that controls whether the real observation values are overwritten. This reading is already included, but out-commented, in the templates. 131 131 132 Example implementations using `PDAFomi_put_state_generate_obs` and `readwrite_obs.F90` are provided by the test case `testsuite/src/lorenz96_omi`. These also use the flag `twin_experiment` to activate the twin experiment (Note: These two test cases always use simulated observations. Nonetheless, they allowto see how the synthetic observations are generated with PDAF and how they are used in a twin experiment).132 Example implementations using `PDAFomi_put_state_generate_obs` and `readwrite_obs.F90` are provided by the implementation of PDAF with the Lorenz-96 modle in `models/lorenz96/`. These also use the flag `twin_experiment` to activate the twin experiment (Note: The Lorenz-96 model case always use simulated observations. Nonetheless, it allows to see how the synthetic observations are generated with PDAF and how they are used in a twin experiment). 133 133