Computation of CRPS
This routine computes the continuous ranked probability score (CRPS) and its decomposition into uncertainty and resolution: CRPS = RELI + RESOL. In addition the uncertainty is computed. A perfectly reliable system gives RELI=0. An informative system gives RESOL << UNCERT.
The computation follows H. Hersbach, Weather and Forecasting 15(2000) 599-570. Here, RESOL is equivalent to CPRS_pot.
Revision history: * 2021-05 - Lars Nerger - Initial code based on sangoma_ComputeCRPS * Later revisions - see repository log
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | dim | |||
integer, | intent(in) | :: | dim_ens | |||
integer, | intent(in) | :: | element | |||
real, | intent(in) | :: | oens(dim,dim_ens) | |||
real, | intent(in) | :: | obs(dim) | |||
real, | intent(out) | :: | CRPS | |||
real, | intent(out) | :: | reli | |||
real, | intent(out) | :: | resol | |||
real, | intent(out) | :: | uncert | |||
integer, | intent(out) | :: | status |