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

U_likelihood_l
v1 v2 5 5 The routine `U_likelihood_l` is a callback routine that has to be provided by the user. In the simplified interface the predefined name of the routine is `likelihood_l_pdaf`, but in the full interface, the user can choose the name of the routine. 6 6 The routine is used in the localized nonlinear filter LNETF and has to compute the likelihood of the observation for a given ensemble member according to the observations used for the local analysis. 7 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))''' for the local observations. The vector '''yHX =''' `resid_l` is prov ed as an input argument. The likelihood has to be returned in the variable `likely_l`.7 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))''' for the local observations. The vector '''yHX =''' `resid_l` is provided as an input argument. The likelihood has to be returned in the variable `likely_l`. 8 8 9 9 This routine is also the place to perform observation localization. To initialize a vector of weights, the routine `PDAF_local_weight` can be called. The procedure is used in the example implementation and also demonstrated in the template routine.