EMPIRE DA
v1.9.1
Data assimilation codes using EMPIRE communication
|
Go to the source code of this file.
Functions/Subroutines | |
subroutine | update_state (state, fpsi, kgain, betan) |
Subroutine to update the state. More... | |
subroutine update_state | ( | real(kind=rk), dimension(state_dim), intent(out) | state, |
real(kind=rk), dimension(state_dim), intent(in) | fpsi, | ||
real(kind=rk), dimension(state_dim), intent(in) | kgain, | ||
real(kind=rk), dimension(state_dim), intent(inout) | betan | ||
) |
Subroutine to update the state.
This subroutine is here because, mathematically, in a particle filter \(x^{k+1} = f(x^k) + A^k + \xi^k\)
However sometimes the result needs to be bounded, some variables need to be exactly related or maybe even something else.
This can be changed for the specific model if it needs to be, in order to bound variables etc.
NOTE this the theory not mathematically correct.
[in] | fpsi | deterministic model update \(f(x^{n-1})\) |
[in] | kgain | nudging term |
[in,out] | betan | Stochastic term |
[out] | state | The updated state vector |
Definition at line 47 of file update_state.f90.