EMPIRE DA  v1.9.1
Data assimilation codes using EMPIRE communication
 All Classes Files Functions Variables Pages
operator_wrappers.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine k (y, x)
 Subroutine to apply \(K\) to a vector y in observation space where \(K := QH^T(HQH^T+R)^{-1}\). More...
 
subroutine bprime (y, x, QHtR_1y, normaln, betan)
 subroutine to calculate nudging term and correlated random errors efficiently More...
 

Function/Subroutine Documentation

subroutine bprime ( real(kind=rk), dimension(obs_dim,pf%count), intent(in)  y,
real(kind=rk), dimension(state_dim,pf%count), intent(out)  x,
real(kind=rk), dimension(state_dim,pf%count), intent(out)  QHtR_1y,
real(kind=rk), dimension(state_dim,pf%count), intent(in)  normaln,
real(kind=rk), dimension(state_dim,pf%count), intent(out)  betan 
)

subroutine to calculate nudging term and correlated random errors efficiently

Parameters
[in]y(obs_dim,pf%count) vectors of innovations \(y-H(x^{n-1})\)
[out]x(state_dim,pf%count) vectors of \(\rho H^TR^{-1}[y-H(x^{n-1})]\)
[out]QHtR_1y(state_dim,pf%count) vectors of \(\rho QH^TR^{-1}[y-H(x^{n-1})]\)
[in]normaln(state_dim,pf%count) uncorrelated random vectors such that normaln(:,i) \(\sim \mathcal{N}(0,I)\)
[out]betan(state_dim,pf%count) correlated random vectors such that betan(:,i) \(\sim \mathcal{N}(0,Q)\)

Definition at line 109 of file operator_wrappers.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine k ( real(kind=rk), dimension(obs_dim,pf%count), intent(in)  y,
real(kind=rk), dimension(state_dim,pf%count), intent(out)  x 
)

Subroutine to apply \(K\) to a vector y in observation space where \(K := QH^T(HQH^T+R)^{-1}\).

Parameters
[in]yvector in observation space
[out]xvector in state space

Definition at line 32 of file operator_wrappers.f90.

Here is the call graph for this function:

Here is the caller graph for this function: