Changes between Version 7 and Version 8 of OMI_observation_operators

Timestamp:

Nov 30, 2020, 3:15:13 PM (4 years ago)

Author:

lnerger

Comment:

--

OMI_observation_operators

@@ -24 +24 @@
 [[PageOutline(2-3,Contents of this page)]]
 
-The observation operator is called for each observation type in the routine `obs_op_pdafomi` in the file `callback_obs_pdafomi.F90`.
+An observation operator routine is called for each observation type in the routine `obs_op_pdafomi` in the file `callback_obs_pdafomi.F90`.
+
+== Observation operators ==
 
 OMI currently provides 3 observation operators:
@@ -51 +53 @@
  Determine barycentric interpolation coefficients for triangular grids based on the coordinates of grid points and the observation
 
-An example of initializing interpolation coefficients with PDAFomi_get_interp_coeff_lin and of using PDAFomi_obs_op_interp_lin is provided in `tutorial/online_2D_serialmodel_omi/obs_C_pdafomi.F90`.
+An example of initializing interpolation coefficients with PDAFomi_get_interp_coeff_lin and of using PDAFomi_obs_op_interp_lin is provided in `tutorial/online_2D_serialmodel/obs_C_pdafomi.F90`.
 
 
@@ -133 +135 @@
 The current set of observation operators provided by PDAF-OMI is rather fundamental. However, there are also observation types which are not variables of state vector, but functions of several values. Likewise one might want to use a more sophisticated interpolation than the linear one or some interpolation that treats the horizontal and vertical directions separately. For these cases you can implement you own operators.
 
-The PDAF-package provides a template to implement new observation operators in `templates/omi/obs_op_pdafomi_TEMPLATE.F90`.
+The PDAF-package provides a template to implement new observation operators in `/templates/omi/obs_op_pdafomi_TEMPLATE.F90`.
 
 Each observation operator include the following functionality:
  1. Check whether `thisobs%doassim==1`, which indicates that the observation is assimilated
  2. initialize the observation vector `ostate_p` for the observation type for which the routine is called. For a parallel model this is done for for the process-local domain (for a model without parallelization this is the global domain)
- 3. Call `PDAFomi_gather_obsstate` to gather the full observation vector `obs_f_all`. This call is mandatory even if theh model is not parallelized.
+ 3. Call `PDAFomi_gather_obsstate` to gather the full observation vector `obs_f_all`. This call is mandatory even if the model is not parallelized.
 
-Generally one can implement and function that computes an observation from the elements of the state vector that is provided to the routine as `state_p`. The coordinate information is provided in `thisobs`. One can also modify the interface to the obsration operator routine if, e.g., additional information is required.
+Generally one can implement any function that computes an observation from the elements of the state vector that is provided to the routine as `state_p`. The coordinate information is provided in `thisobs`. One can also modify the interface to the obsration operator routine if, e.g., additional information is required.
 
 Dependent on the complexity of an observation operator it might be useful to separate the functional computations from the interpolation. For this one could consider multi-step observation operators in separate routines.