Compute mean observation error variance
This routine will only be called, if the adaptive forgetting factor feature is used. Please note that this is an experimental feature.
The routine is called in the loop over all local analysis domains during each analysis by the routine PDAF_set_forget_local that estimates a local adaptive forgetting factor. The routine has to initialize the mean observation error variance for the current local analysis domain. (See init_obsvar_f for a global variant)
The routine assumed a diagonal observation error covariance matrix.
If the observation counter is zero the computation of the mean variance is initialized. The output is always the mean variance. If the observation counter is >0 first the variance sum is computed by multiplying with the observation counter.
Revision history: * 2019-09 - Lars Nerger - Initial code from restructuring observation routines * Later revisions - see repository log
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(obs_l), | intent(inout) | :: | thisobs_l | |||
| type(obs_f), | intent(inout) | :: | thisobs | |||
| real, | intent(inout) | :: | meanvar_l | |||
| integer, | intent(inout) | :: | cnt_obs_l |