Changes between Version 2 and Version 3 of PDAF_diag_crps

Timestamp:

Mar 26, 2025, 2:34:38 PM (8 days ago)

Author:

lnerger

Comment:

--

PDAF_diag_crps

@@ -10 +10 @@
 The routine can be called in the pre/poststep routine of PDAF both before and after the analysis step to compute the CRPS.
 
-The interface is the following:
+The interface is:
 {{{
-  SUBROUTINE PDAF_diag_CRPS(dim, dim_ens, element, oens, obs, &
+  SUBROUTINE PDAF_diag_CRPS(dim_p, dim_ens, element, oens, obs, &
              CRPS, reli, resol, uncert, status)
-}}}
-with the following arguments:
-{{{
-  INTEGER, INTENT(in) :: dim                ! PE-local state dimension
-  INTEGER, INTENT(in) :: dim_ens            ! Ensemble size
-  INTEGER, INTENT(in) :: element            ! ID of element to be used
-       !< If element=0, mean values over all elements are computed
-  REAL, INTENT(in)    :: oens(dim, dim_ens) ! State ensemble
-  REAL, INTENT(in)    :: obs(dim)           ! State ensemble
-  REAL, INTENT(out)   :: CRPS               ! CRPS
-  REAL, INTENT(out)   :: reli               ! Reliability
-  REAL, INTENT(out)   :: resol              ! resolution
-  REAL, INTENT(out)   :: uncert             ! uncertainty
-  INTEGER, INTENT(out) :: status            ! Status flag (0=success)
+
+  INTEGER, INTENT(in) :: dim_p                ! PE-local state dimension
+  INTEGER, INTENT(in) :: dim_ens              ! Ensemble size
+  INTEGER, INTENT(in) :: element              ! ID of element to be used
+       ! If element=0, mean values over all elements are computed
+  REAL, INTENT(in)    :: oens(dim_p, dim_ens) ! State ensemble
+  REAL, INTENT(in)    :: obs(dim_p)           ! State ensemble
+  REAL, INTENT(out)   :: CRPS                 ! CRPS
+  REAL, INTENT(out)   :: reli                 ! Reliability
+  REAL, INTENT(out)   :: resol                ! resolution
+  REAL, INTENT(out)   :: uncert               ! uncertainty
+  INTEGER, INTENT(out) :: status              ! Status flag (0=success)
 }}}
 
@@ -34 +32 @@
  * A perfectly reliable system gives `reli=0`. An informative system gives `resol << uncert`.
  * Compared to Hersbach (2000), `resol` here is equivalent to `CPRS_pot`.
- * The routine is parallelized. It uses the PDAF-internal parallelization information is has to be called on all filter processes
+ * The routine is parallelized. For `element=0` it computes the CRPS, reli, and resol for the global state vector. It uses the PDAF-internal parallelization information is has to be called on all filter processes
  * The routine uses a rather simple sorting algorithm. Accordingly, the performance is likely suboptimal for high-dimensional cases.