Changes between Version 18 and Version 19 of OMI_observation_modules_PDAF3


Ignore:
Timestamp:
May 31, 2026, 9:46:40 AM (5 hours ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • OMI_observation_modules_PDAF3

    v18 v19  
    103103
    104104The main variables of **thisobs** that are directly filled in this routine are
    105  * [wiki:#thisobsdoassim thisobs%doassim]: Specify whether this observation type is assimilated
    106  * [wiki:#thisobsdisttype thisobs%disttype]: Specify the type of distance computation
    107  * [wiki:#thisobsncoord thisobs%ncoord]: Specify the number of dimensions used to compute distances
    108  * [wiki:#thisobsid_obs_p thisobs%id_obs_p]: store the indices of state vector elements that correspond to an observation (A single value for observation at grid points, or multiple values for derived quantities or interpolation; this is only used in the OMI-provided observation operators)
     105 * [wiki:OMI_observation_modules_PDAF3#thisobsdoassim thisobs%doassim]: Specify whether this observation type is assimilated
     106 * [wiki:OMI_observation_modules_PDAF3#thisobsdisttype thisobs%disttype]: Specify the type of distance computation
     107 * [wiki:OMI_observation_modules_PDAF3#thisobsncoord thisobs%ncoord]: Specify the number of dimensions used to compute distances
     108 * [wiki:OMI_observation_modules_PDAF3#thisobsid_obs_p thisobs%id_obs_p]: store the indices of state vector elements that correspond to an observation (A single value for observation at grid points, or multiple values for derived quantities or interpolation; this is only used in the OMI-provided observation operators)
    109109
    110110The next variables are filled and provided to `PDAFomi_gather_obs`. The variable names shown here are those used in the template and tutorial codes.
    111  * [wiki:#dim_obs_p dim_obs_p]: Count the number of available observations
    112  * [wiki:#obs_p obs_p]: Fill the vector of observations
    113  * [wiki:#ocoord_p ocoord_p]: store the coordinates of the observations (only actually used in case of localization)
    114  * [wiki:#ivar_obs_p ivar_obs_p]: store the inverse error variance of each observation for the default mode of OMI, which assumes a diagonal observation error covariance matrix
     111 * [wiki:OMI_observation_modules_PDAF3#dim_obs_p dim_obs_p]: Count the number of available observations
     112 * [wiki:OMI_observation_modules_PDAF3#obs_p obs_p]: Fill the vector of observations
     113 * [wiki:OMI_observation_modules_PDAF3#ocoord_p ocoord_p]: store the coordinates of the observations (only actually used in case of localization)
     114 * [wiki:OMI_observation_modules_PDAF3#ivar_obs_p ivar_obs_p]: store the inverse error variance of each observation for the default mode of OMI, which assumes a diagonal observation error covariance matrix
    115115
    116116Further possible initializations are:
    117  * When one uses PDAF-OMI's observation operator with interpolation, one further needs to initialize an array of interpolation coefficients ([wiki:OMI_observation_modules_PDAF3#thisobsicoeff_p thisobs%icoeff_p]).
    118  * For Cartesian distance computation with periodicity one also needs to set [wiki:OMI_observation_modules#thisobsdomainsize thisobs%domainsize].
    119  * When using local ensemble filters with parallelization, and can also use `thisobs%use_global_obs=0` to improve the performance of the search for local observations. This feature requires that the sub-domain limits have been provided to PDAF-OMI with `PDAFomi_set_domain_limits` (which is usually don ein `init_pdaf`).
    120  * One can also activate the [wiki:OMI_additional_functionality_PDAF3#Omittingobservationsthatarepotentialoutliers omission of observations that are too different from the ensemble mean]. This is activated by setting `thisobs%inno_omit>0.0`
    121  * One can store the name of the observation type  (thisobs%name; a string of length up to 20 characters)
    122 
    123 When parallel model with domain decomposition is used, the variables with suffix `_p` need to describe the observation information for a particular process sub-domain.
     117 * [wiki:OMI_observation_modules_PDAF3#thisobsicoeff_p thisobs%icoeff_p]: When one uses PDAF-OMI's observation operator with interpolation, one further needs to initialize this array of interpolation coefficients.
     118 * [wiki:OMI_observation_modules#thisobsdomainsize thisobs%domainsize]: This needs to be set for Cartesian distance computation with periodicity.
     119 * [wiki:OMI_observation_modules_PDAF3#thisobsuse_global_obs thisobs%use_global_obs]: When using local ensemble filters with parallelization, and can also set `thisobs%use_global_obs=0` to improve the performance of the search for local observations. This feature requires that the sub-domain limits have been provided to PDAF-OMI with [wiki:OMI_additional_functionality_PDAF3#PDAFomi_set_domain_limits PDAFomi_set_domain_limits] (which is usually done in `init_pdaf`).
     120 * [wiki:OMI_observation_modules_PDAF3#thisobsinno_omit thisobs%inno_omit]: One can also activate the [wiki:OMI_additional_functionality_PDAF3#Omittingobservationsthatarepotentialoutliers omission of observations that are too different from the ensemble mean]. This is activated by setting `thisobs%inno_omit>0.0`.
     121 * thisobs%name: One can store the name of the observation type  (thisobs%name; a string of length up to 20 characters)
     122
     123When a parallel model with domain decomposition is used, the variables with suffix `_p` need to describe the observation information for a particular process sub-domain.
    124124
    125125** 2. Call PDAFomi_gather_obs**