PDAFomi_obs_op_gridpoint Subroutine

public subroutine PDAFomi_obs_op_gridpoint(thisobs, state_p, obs_f_all)

observation operator for data at grid points

Application of 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: * 2019-06 - Lars Nerger - Initial code from restructuring observation routines * Later revisions - see repository log

Arguments

Type IntentOptional Attributes Name
type(obs_f), intent(inout) :: thisobs
real, intent(in) :: state_p(:)
real, intent(inout) :: obs_f_all(:)