diff --git a/include/evalhyd/utils.hpp b/include/evalhyd/utils.hpp
index ff4bcd274ba74b518bcf3685dc9ced40235e484c..aad268c2f05d13cd018379cf3523a6d8cf158efc 100644
--- a/include/evalhyd/utils.hpp
+++ b/include/evalhyd/utils.hpp
@@ -11,7 +11,7 @@ namespace evalhyd
 {
     namespace utils
     {
-        /// Determine whether flows are strictly greater than given threshold(s)
+        /// Determine whether flows are greater than given threshold(s)
         ///
         /// Streamflow data is always expected as a 2D array (even if only for
         /// one time series). Threshold(s) given is (are) always expected as a
@@ -19,7 +19,7 @@ namespace evalhyd
         /// returned whose first two dimensions are of the same sizes as the
         /// streamflow data and whose third dimension is of size equal to the
         /// number of thresholds given. The returned array contains ones where
-        /// the threshold is strictly exceeded, and zeros otherwise.
+        /// the threshold is exceeded, and zeros otherwise.
         ///
         /// \param [in] q:
         ///     2D array of streamflow data
@@ -72,7 +72,7 @@ namespace evalhyd
             );
 
             // return a boolean-like array
-            return q2 <= thr;
+            return q2 < thr;
         }
     }
 }