Adjoint observation operator for data at grid points
Application of adjoint observation operator for the case that model variables are observerved at model grid points.
For this case INIT_DIM_OBS_F will prepare the index array thisobs%id_obs_p containing the information which elements of the PE-local state vector contain the observed values.
The routine is called by all filter processes. It first selects the observed elements for a PE-local domain. Afterwards, the values are gathered into the full vector using PDAFomi_gather_obsstate.
The routine has to fill the part of the full observation vector OBS_F_ALL that represents the current observation type. The routine first applied the observation operator for the current observation type and the calls PDAFomi_gather_obsstate to gather the observation over all processes and fills OBS_F_ALL.
The routine has to be called by all filter processes.
Revision history: * 2021-04 - Lars Nerger - Initial code from restructuring observation routines * Later revisions - see repository log
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(obs_f), | intent(inout) | :: | thisobs | |||
real, | intent(in) | :: | obs_f_all(:) | |||
real, | intent(inout) | :: | state_p(:) |