| | 26 | |
| | 27 | With OMI, the functionality to handle observations is shifted to generic routines and observation-specific modules (obs_*_pdafomi in the third column in Fig. 1, denoted obs-module below). Most routines for the observation handling are generic with OMI and included in PDAF so that only three routines with observation specific operations need to be implemented. |
| | 28 | There is one obs-module per observation type with collects the three routines. For example, one can have one obs-module for the satellite sea surface temperature from one data provider and another one for sea level anomaly data. Important is that each of these obs-modules, which are further described below, is independent from the others. This allows us to switch between different combinations of observations ensuring that their implementations don’t interfere. |
| | 29 | |
| | 30 | Since the actual opernations are done in the obs-modules obs_*_pdafomi, the remaining three call-back routines (init_dim_obs, obs_op, init_dim_obs_l) are reduced to pass-through routines. Thus, each of these routines contains only calls to one observation-specific routine from each obs-module. There is, in addition, one routine for deallocating the observation-related arrays which calls a generic PDAF-OMI routine for each observation type. As the call-back routines become very compact by this change, they are collected into the file callback_obs_pdafomi.F90. This is mainly for convenience, because all these routines are now ‘in one place’. This simplifies the addition of a new data type by adding a subroutine call to each of the routines in this file. (Note, that callback_obs_pdafomi.F90 does not define a Fortran module. This is because the call-back routines of PDAF are not contained in a module as this would require to compile the module together with the PDAF core so that the PDAF core would no longer be generic.) When adding a new observation type only a new subroutine call has to be added to each of the routines in callback_obs_pdafomi.F90. |
| | 31 | |
| | 32 | Each obs-module contains three routines: |
| | 33 | 1. `init_dim_obs` initializes all variables holding the information about one observation type. The information about the observation type is stored in a data structure (Fortran type). |
| | 34 | 1. `obs_op` applies the observation operator to a state vector. For the latter routine one can select an observation operator from the list of operators provided by OMI, or one can to implement a new operator. |
| | 35 | 1. `init_dim_obs_l` calls a PDAF-OMI routine to initialize the observation information corresponding to a local analysis domain. Here for can set a different localization radius for each observation type. |
| | 36 | |
| | 37 | For each observation type, PDAF-OMI uses a data structure that was initialized in the obs-module. The set of routines in callback_obs_pdafomi.F90 provide the observation handling for all filters and smoothers provided by PDAF. Thus once, init_dim_obs_f and obs_op_f are implemented for a particular observation and the subroutine calls in callback_obs_pdafomi.F90 for this observation type are inserted, one can use each of the provided assimilation methods. (An exception is the local EnKF, which needs one additional routine to localize the covariance matrix. The actual operation is, however, also performed by a routine provided by PDAF-OMI). |
| | 38 | |
| | 39 | |
| | 40 | |
