pyPDAF.PDAF.diag_CRPS

pyPDAF.PDAF.diag_CRPS()

Obtain a continuous rank probability score for an ensemble.

The implementation is based on [1].

References

Parameters:

  • element (int) – ID of element to be used. If element=0, mean values over all elements are computed

  • oens (ndarray[tuple[dim, dim_ens], np.float64]) –

    State ensemble

    The1st-th dimension dim is PE-local state dimension; the2nd-th dimension dim_ens is Ensemble size.

  • obs (ndarray[tuple[dim], np.float64]) –

    State ensemble

    The array dimension dim is PE-local state dimension.

Returns:

  • CRPS (float) – CRPS

  • reli (float) – Reliability

  • resol (float) – resolution

  • uncert (float) – uncertainty

  • status (int) – Status flag (0=success)