PDAF_diag_CRPS Subroutine

subroutine PDAF_diag_CRPS(dim, dim_ens, element, oens, obs, CRPS, reli, resol, uncert, status)

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

Arguments

Type IntentOptional 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