Compute product of inverse of R with some matrix
The routine is called during the analysis step of the global square-root filters. It has to compute the product of the inverse of the process-local observation error covariance matrix with the matrix of process-local observed ensemble perturbations.
This routine assumes a diagonal observation error covariance matrix, but allows for varying observation error variances.
The routine can be applied with either all observations of different types at once, or separately for each observation type. The operation is done with all process-local observations
Revision history: * 2019-12 - Lars Nerger - Initial code from restructuring observation routines * Later revisions - see repository log
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(obs_f), | intent(inout) | :: | thisobs | |||
integer, | intent(in) | :: | ncols | |||
real, | intent(in) | :: | A_p(:,:) | |||
real, | intent(out) | :: | C_p(:,:) |