PDAFomi_obs_l Module

PDAF-OMI routines for local observations

This module contains generic routines for several observation-related operations for local filters. The routines are

• PDAFomi_set_debug_flag \n Set or unset the debugging flag for PDAFomi routines
• PDAFomi_init_dim_obs_l \n Initialize dimension of local obs. vetor and arrays for local observations
• PDAFomi_cnt_dim_obs_l \n Set dimension of local obs. vector with isotropic localization
• PDAFomi_cnt_dim_obs_l_noniso \n Set dimension of local obs. vector with nonisotropic localization
• PDAFomi_init_obsarrays_l \n Initialize arrays for the index of a local observation in the full observation vector and its corresponding distance.
• PDAFomi_init_obsarrays_l_noniso \n Initialize arrays for the index of a local observation in the full observation vector and its corresponding distance with onoisotrppic localization.
• PDAFomi_g2l_obs \n Initialize local observation vector from full observation vector
• PDAFomi_init_obs_l \n Initialize the local vector of observations
• PDAFomi_prodRinvA_l \n Multiply an intermediate matrix of the local filter analysis with the inverse of the observation error covariance matrix and apply observation localization
• PDAFomi_prodRinvA_hyb_l \n Multiply an intermediate matrix of the local filter analysis with the inverse of the observation error covariance matrix and apply observation localization. In addition apply the hybrid weight
• PDAFomi_init_obsvar_l \n Compute mean observation error variance
• PDAFomi_likelihood_l \n Compute local likelihood for an ensemble member
• PDAFomi_likelihood_hyb_l \n Compute local likelihood for an ensemble member taking into account a hybrid weight for tempering
• PDAFomi_localize_covar_iso \n Apply covariance isotropic localization in LEnKF
• PDAFomi_localize_covar_noniso \n Apply non-isotropic covariance localization in LEnKF
• PDAFomi_g2l_obs_internal \n Internal routine to initialize local observation vector from full observation vector (used by PDAFomi_init_obs_l and PDAFomi_g2l_obs)
• PDAFomi_comp_dist2 \n Compute squared distance
• PDAFomi_check_dist2 \n Compute and check distance for isotropic localization
• PDAFomi_check_dist2_noniso \n Compute and check distance for non-isotropic localization
• PDAFomi_weights_l \n Compute a vector of localization weights
• PDAFomi_deallocate_obs \n Deallocate arrays in observation type
• PDAFomi_dealloc \n Deallocate arrays in all observation types
• PDAFomi_omit_by_innovation_l !! Exclude observations if innovation is too large (thisobs%inno_exclude)

Revision history: * 2019-06 - Lars Nerger - Initial code * Later revisions - see repository log

Set debugging flag

This routine sets the debug flag for PDAF-OMI. One can set the flag dependent on the local analysis domain, the MPI rank, or the OpenMP thread ID, or and combination of them.

For debugval>0 additional information is written by the OMI routine to stdout. One should activate the debugging before calling some selected routine(s) and deactivate it with debugval=0 afterwards. This allows for a targeted checking of the functionality.

Revision history: * 2019-09 - Lars Nerger - Initial code * Later revisions - see repository log