Changes between Version 15 and Version 16 of OMI_Callback_obs_pdafomi


Ignore:
Timestamp:
Nov 24, 2020, 11:05:29 AM (3 years ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • OMI_Callback_obs_pdafomi

    v15 v16  
    1111In the descriptions below we use 'TYPE' as a generic label for an observation type. When implemeting an observation type, one replace this by an actual name. For example for sea surface temperature (sst) observations, one could replace TYPE by 'sst' and write e.g. `init_dim_obs_sst`.
    1212
     13
    1314== init_dim_obs_pdafomi ==
    1415
    15 The routine is called at the beginning of each analysis step.  For PDAF, it has to initialize the size `dim_obs_p` of the observation vector according to the current time step. Apart from this routine will initialize overall observation information.
     16The routine is called at the beginning of each analysis step.  For PDAF, it has to initialize the size `dim_obs_p` of the observation vector according to the current time step. Apart from this routine will initialize overall observation information.  In this routine one just calls `init_dim_obs_TYPE` for each observation type.
    1617
    1718The '''interface''' for this routine is:
     
    6667== obs_op_pdafomi ==
    6768
    68 In this routine one just calls `obs_op_TYPE` for each observation type.
     69The routine is called during the analysis step. It has to perform the operation of the observation operator acting on a state vector that is provided as `state_p`. The observed state has to be returned in `m_state_p`. In this routine one just calls `obs_op_TYPE` for each observation type.
     70
     71The '''interface''' for this routine is:
     72{{{
     73SUBROUTINE obs_op_pdafomi(step, dim_p, dim_obs_p, state_p, m_state_p)
     74
     75  INTEGER, INTENT(in) :: step               ! Currrent time step
     76  INTEGER, INTENT(in) :: dim_p              ! PE-local dimension of state
     77  INTEGER, INTENT(in) :: dim_obs_p          ! Dimension of observed state
     78  REAL, INTENT(in)    :: state_p(dim_p)     ! PE-local model state
     79  REAL, INTENT(out) :: m_state_p(dim_obs_p) ! PE-local observed state
     80}}}
    6981
    7082The '''implementation steps''' are: