Changes between Version 1 and Version 2 of U_likelihood_l


Ignore:
Timestamp:
Mar 6, 2018, 4:54:11 PM (7 years ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • U_likelihood_l

    v1 v2  
    55The routine `U_likelihood_l` is a call-back 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.
    66The 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*(y-Hx)^T^*R^-1^*(y-Hx))''' for the local observations. The vector '''y-HX =''' `resid_l` is proved as an input argument. The likelihood has to be returned in the variable `likely_l`.
     7The likelihood depends on the assumed observation error distribution. For a Gaussian observation error, the likelihood is '''exp(-0.5*(y-Hx)^T^*R^-1^*(y-Hx))''' for the local observations. The vector '''y-HX =''' `resid_l` is provided as an input argument. The likelihood has to be returned in the variable `likely_l`.
    88
    99This 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.