fix wrong ranks in contingency metrics

include/evalhyd/detail/probabilist/contingency.hpp

@@ -46,7 +46,7 @@ namespace evalhyd
             ///     Alerts based on forecast.
             ///     shape: (sites, lead times, levels, thresholds, time)
             inline xt::xtensor<double, 5> calc_a_k(
-                    const xt::xtensor<double, 3>& sum_f_k,
+                    const xt::xtensor<double, 4>& sum_f_k,
                     std::size_t n_mbr
             )
             {
@@ -235,10 +235,6 @@ namespace evalhyd
 
         namespace metrics
         {
-            // -----------------------------------------------------------------
-            // Accuracy
-            // -----------------------------------------------------------------
-
             namespace detail
             {
                 template <class XD2>
@@ -256,7 +252,7 @@ namespace evalhyd
                 )
                 {
                     // initialise output variable
-                    xt::xtensor<double, 4> METRIC =
+                    xt::xtensor<double, 6> METRIC =
                             xt::zeros<double>({n_sit, n_ldt, n_msk, n_exp,
                                                n_mbr + 1, n_thr});