add missing check on threholds site dimension

include/evalhyd/evalp.hpp

@@ -257,6 +257,18 @@ namespace evalhyd
                     "numbers of sites"
             );
         }
+
+        if (q_thr_.size() > 0)
+        {
+            if (q_obs_.shape(0) != q_thr_.shape(0))
+            {
+                throw std::runtime_error(
+                        "observations and thresholds feature different "
+                        "numbers of sites"
+                );
+            }
+        }
+
         if (t_msk_.size() > 0)
         {
             if (q_obs_.shape(0) != t_msk_.shape(0))