Changes between Version 1 and Version 2 of likelihood_pdaf
 Timestamp:
 Mar 6, 2018, 4:52:33 PM (3 years ago)
Legend:
 Unmodified
 Added
 Removed
 Modified

likelihood_pdaf
v1 v2 4 4 5 5 The routine `likelihood_pdaf` (called `U_likelihood` inside the PDAF core routines) is a callback routine that has to be provided by the user. 6 The routine is used with the nonlinear filter NETF and called during the analysis step. The purpose of the routine is to compute the likelihood of the observation for a given ensemble member. The likelihood depends on the assumed observation error distribution. For a Gaussian observation error, the likelihood is '''exp(0.5*(yHx)^T^*R^1^*(yHx))'''. The vector '''yHX =''' `resid` is prov ed as an input argument. The likelihood has to be returned in the variable `likely`.6 The routine is used with the nonlinear filter NETF and called during the analysis step. The purpose of the routine is to compute the likelihood of the observation for a given ensemble member. The likelihood depends on the assumed observation error distribution. For a Gaussian observation error, the likelihood is '''exp(0.5*(yHx)^T^*R^1^*(yHx))'''. The vector '''yHX =''' `resid` is provided as an input argument. The likelihood has to be returned in the variable `likely`. 7 7 8 8 The interface is the following: