Changes between Version 1 and Version 2 of U_likelihood


Ignore:
Timestamp:
Mar 6, 2018, 4:53:43 PM (6 years ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • U_likelihood

    v1 v2  
    44
    55The routine `U_likelihood` 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_pdaf`, but in the full interface, the user can choose the name of the routine.
    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*(y-Hx)^T^*R^-1^*(y-Hx))'''. The vector '''y-HX =''' `resid` is proved as an input argument. The likelihood has to be returned in the variable `likely`.
     6The 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*(y-Hx)^T^*R^-1^*(y-Hx))'''. The vector '''y-HX =''' `resid` is provided as an input argument. The likelihood has to be returned in the variable `likely`.
    77
    88The interface is the following: