pyPDAF.PDAFomi.diag_rmsd

pyPDAF.PDAFomi.diag_rmsd(nobs: int, verbose: int) → Tuple[int, np.ndarray]

Compute RMSD between observations and observed ensemble means.

For each active OMI observation type, PDAF-OMI compares the observation vector y with the observed ensemble mean Hx and returns the root mean squared deviation sqrt(mean((y - Hx)**2)) over all observations of that type. If only the observed ensemble H(X_i) was stored for diagnostics, PDAF-OMI first computes ensemble mean.

This function is only useful after observation operators are performed.

Parameters:

• nobs (int) – Input/output number of OMI observation types. The input value can be arbitrary; PDAF-OMI overwrites it with the number of active observation types for which diagnostics are returned.

• verbose (int) – Verbosity flag. If greater than zero, PDAF-OMI prints the RMSD table.

Returns:

• nobs (int) – Number of observation-type entries returned by PDAF-OMI.

• rmsd (ndarray[np.float64, ndim=1]) – PDAF-owned RMSD values. rmsd[i] is the RMSD for observation type i + 1.