Compute likelihood for an ensemble member
The routine is called during the analysis step of the NETF or a particle filter. It has to compute the likelihood of the ensemble according to the difference from the observation (residual) and the error distribution of the observations.
In general this routine is similar to the routine prodRinvA used for ensemble square root Kalman filters. As an addition to this routine, we here have to evaluate the likelihood weight according the assumed observation error statistics.
Revision history: * 2020-03 - Lars Nerger - Initial code from restructuring observation routines * Later revisions - see repository log
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(obs_f), | intent(inout) | :: | thisobs | |||
| real, | intent(in) | :: | resid(:) | |||
| real, | intent(inout) | :: | lhood |