PDAFomi_prodRinvA Subroutine

public subroutine PDAFomi_prodRinvA(thisobs, ncols, A_p, C_p)

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

Arguments

Type IntentOptional Attributes Name
type(obs_f), intent(inout) :: thisobs
integer, intent(in) :: ncols
real, intent(in) :: A_p(:,:)
real, intent(out) :: C_p(:,:)