| | 1 | = PDAFomi_observation_localization_weights = |
| | 2 | |
| | 3 | This page documents the routine `PDAFomi_observation_localization_weights` of PDAF-OMI. It is used with non-diagonal observation error covariance matrices |
| | 4 | |
| | 5 | The routine returns a vector of observation localization weights computed according to the localization specifications in `thisobs_l` and observation coordinates in `thisobs`. The routine can be called in `prodRinvA_l_pdafomi` or `likelihood_l_pdafomi` in order to obtain the localization weights for observations. For diagonal observation error covariance matrices, the weighting is computed by OMI internally, but for nondiagonal R-matrices, the use had to implement the application of the localization. |
| | 6 | |
| | 7 | The interface is the following: |
| | 8 | {{{ |
| | 9 | SUBROUTINE PDAFomi_observation_localization_weights(thisobs_l, thisobs, ncols, & |
| | 10 | A_l, weight, verbose) |
| | 11 | |
| | 12 | TYPE(obs_l), INTENT(inout) :: thisobs_l ! Data type with local observation |
| | 13 | TYPE(obs_f), INTENT(inout) :: thisobs ! Data type with full observation |
| | 14 | INTEGER, INTENT(in) :: ncols ! Rank of initial covariance matrix |
| | 15 | REAL, INTENT(in) :: A_l(:, :) ! Input matrix, size (thisobs_l%dim_obs_l, ncols) |
| | 16 | INTEGER, INTENT(in) :: verbose ! Verbosity flag |
| | 17 | REAL, INTENT(out) :: weight(thisobs_l%dim_obs_l) ! Localization weights |
| | 18 | }}} |
| | 19 | |
| | 20 | |
| | 21 | Notes: |
| | 22 | * Array `A_l` is usually the input argument to the routine [wiki:prodRinvA_l_pdaf prodRinvA_l_pdaf]. It is only used when the regulated localization is used. |
| | 23 | * The regulated localization computes an alternative localization function based on the ration of state error variances to observation error variances. the regulated localization method is described in the paper Nerger et al. (Quarterly Journal of the Royal Meteorological Society, 138 (2012) 802-812; see [wiki:PublicationsandPresentations publications]. |
| | 24 | |
