pyPDAF.PDAFomi.diag_diffstats

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

Compare observations with observed ensemble means by observation type.

This function is the same as diag_stats().

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 one row per observation type.

Returns:

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

  • stats (ndarray[np.float64, ndim=2]) – PDAF-owned statistics matrix with shape (6, nobs). The rows are:

    stats[0, :]

    Pearson correlation between observation anomalies and observed ensemble-mean anomalies.

    stats[1, :]

    Centered RMS deviation between y and Hx.

    stats[2, :]

    Bias mean(y) - mean(Hx).

    stats[3, :]

    Mean absolute deviation mean(abs(y - Hx)).

    stats[4, :]

    Standard deviation of the observations y.

    stats[5, :]

    Standard deviation of the observed ensemble mean Hx.