diff --git a/include/evalhyd/detail/probabilist/contingency.hpp b/include/evalhyd/detail/probabilist/contingency.hpp
index 55c17bc1d9c3e8f3ba2124495ed82185a56a2e2f..1ec44d2d3afd7bd4c1ef0d7845de0b51eb9d7271 100644
--- a/include/evalhyd/detail/probabilist/contingency.hpp
+++ b/include/evalhyd/detail/probabilist/contingency.hpp
@@ -448,32 +448,15 @@ namespace evalhyd
 
             // Compute the relative operating characteristic skill score (ROCSS).
             //
-            // \param pod
+            // \param POD
             //     Probabilities of detection.
             //     shape: (subsets, samples, levels, thresholds)
-            // \param pofd
+            // \param POFD
             //     Probabilities of false detection.
             //     shape: (subsets, samples, levels, thresholds)
-            // \param q_thr
-            //     Streamflow exceedance threshold(s).
-            //     shape: (thresholds,)
-            // \param t_msk
-            //     Temporal subsets of the whole record.
-            //     shape: (subsets, time)
-            // \param b_exp
-            //     Boostrap samples.
-            //     shape: (samples, time slice)
-            // \param n_thr
-            //     Number of thresholds.
-            // \param n_mbr
-            //     Number of ensemble members.
-            // \param n_msk
-            //     Number of temporal subsets.
-            // \param n_exp
-            //     Number of bootstrap samples.
             // \return
-            //     Critical success indices.
-            //     shape: (subsets, samples, levels, thresholds)
+            //     ROC skill scores.
+            //     shape: (subsets, samples, thresholds)
             inline xt::xtensor<double, 3> calc_ROCSS(
                     const xt::xtensor<double, 4>& POD,
                     const xt::xtensor<double, 4>& POFD