Changes between Version 10 and Version 11 of OMI_nondiagonal_observation_error_covariance_matrices

Timestamp:

Sep 19, 2024, 12:34:42 PM (2 months ago)

Author:

lnerger

Comment:

--

OMI_nondiagonal_observation_error_covariance_matrices

@@ -254 +254 @@
  * The matrix **R** could be constructed e.g. from decorrelation length scales or by direct implementation. 
  * The routine [wiki:PDAF_correlation_function] can be used to obtain values of a correlation function.
+ * Th eroutine [wiki:PDAFomi_observation_localization_weights] can be used to obtain the values of the localization weights.
  * The matrix **R** cannot be stored fully for larger numbers of observations. E.g. for 100 000 observations, the full matric **R** would have 10^10^ entries, which corresponds to 80 GB memory at double precision. Thus, one would need compressed ways of storing the matrix or one can implement the products with the inverse of **R** or the addition of **R** in form of oeprations.
  * For the localized operations in `prodRinvA_l_pdafomi` and `likelihood_l_pdafomi` it should be possible to handle the part of **R** that relates only to the local observations used in a local analysis domain. However, the inverse of this local matrix part is not identical to the local part of the global inverse. The different might be mitigated by applying the localization weight to **R**.