Changes between Version 2 and Version 3 of OMI_observation_operators

Timestamp:

Nov 23, 2020, 5:44:36 PM (4 years ago)

Author:

lnerger

Comment:

--

OMI_observation_operators

@@ -107 +107 @@
  - The order of the grid points in `gcoords` has to be consistent with the order of the indices specified in `thisobs%id_obs_p`
 
-I
+
+== Implementing your own observation operator ==
+
+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`.
+
+Each observation operator include the following functionality:
+ - Check whether `thisobs%doassim==1`, which indicates that the observation is assimilated
+ - 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)
+ - Call `PDAFomi_gather_obsstate` to gather the full observation vector `obs_f_all`. This call is mandatory even if theh model is not parallelized.
+
+