Observation operator for averaging grid point values
Application of observation operator for the case that the observation value is given as the average of model grid point values. For this case INIT_DIM_OBS_F will prepare the index array thisobs%id_obs_p that contains several rows holding the indices of the state vector elements which are to be averaged to represent an observation.
The routine is called by all filter processes, and the operation has to be performed by each these processes for its PE-local domain before the information from all PEs is gathered.
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
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(obs_f), | intent(inout) | :: | thisobs | |||
integer, | intent(in) | :: | nrows | |||
real, | intent(in) | :: | state_p(:) | |||
real, | intent(inout) | :: | obs_f_all(:) |