PDAFomi_gather_obs Subroutine

public subroutine PDAFomi_gather_obs(thisobs, dim_obs_p, obs_p, ivar_obs_p, ocoord_p, ncoord, lradius, dim_obs_f)

Gather full observational information

This routine uses different gather routines from PDAF to obtain combined full observational information.

The observation-type specific variables that are initialized here are * thisobs\%dim_obs_p - PE-local number of module-type observations * thisobs\%dim_obs_f - full number of module-type observations * thisobs\%off_obs_f - Offset of full module-type observation in overall full obs. vector * thisobs\%off_obs_g - Offset of global module-type observation in overall full obs. vector * thisobs\%obs_f - full vector of module-type observations * thisobs\%ocoord_f - coordinates of observations in OBS_MOD_F * thisobs\%ivar_obs_f - full vector of inverse obs. error variances of module-type

If the full observations are restricted to those relevant for a process domain, there are further initialized * thisobs\%dim_obs_g - Number of global observations * thisobs\%id_obs_f_lim - Ids of full observations in global observations

Revision history: * 2020-03 - Lars Nerger - Initial code from restructuring observation routines * Later revisions - see repository log

Arguments

Type IntentOptional Attributes Name
type(obs_f), intent(inout) :: thisobs
integer, intent(in) :: dim_obs_p
real, intent(in) :: obs_p(:)
real, intent(in) :: ivar_obs_p(:)
real, intent(in) :: ocoord_p(:,:)
integer, intent(in) :: ncoord
real, intent(in) :: lradius
integer, intent(out) :: dim_obs_f