Merge branch 'master' into ndim

2 files changed, 18 insertions, 0 deletions

diff --git a/pint/facets/numpy/numpy_func.py b/pint/facets/numpy/numpy_func.py
index 7143143..59c2f98 100644
--- a/pint/facets/numpy/numpy_func.py
+++ b/pint/facets/numpy/numpy_func.py
@@ -552,6 +552,12 @@ def _interp(x, xp, fp, left=None, right=None, period=None):
 
 @implements("where", "function")
 def _where(condition, *args):
+    if not getattr(condition, "_is_multiplicative", True):
+        raise ValueError(
+            "Invalid units of the condition: Boolean value of Quantity with offset unit is ambiguous."
+        )
+
+    condition = getattr(condition, "magnitude", condition)
     args, output_wrap = unwrap_and_wrap_consistent_units(*args)
     return output_wrap(np.where(condition, *args))
 
diff --git a/pint/testsuite/test_numpy.py b/pint/testsuite/test_numpy.py
index f2ddaf0..77d18e3 100644
--- a/pint/testsuite/test_numpy.py
+++ b/pint/testsuite/test_numpy.py
@@ -1103,6 +1103,18 @@ class TestNumpyUnclassified(TestNumpyMethods):
                 0 * self.ureg.J,
             )
 
+        helpers.assert_quantity_equal(
+            np.where([-1, 0, 1] * self.ureg.m, [1, 2, 1] * self.ureg.s, np.nan),
+            [1, np.nan, 1] * self.ureg.s,
+        )
+        with pytest.raises(
+            ValueError,
+            match=".*Boolean value of Quantity with offset unit is ambiguous",
+        ):
+            np.where(
+                self.ureg.Quantity([-1, 0, 1], "degC"), [1, 2, 1] * self.ureg.s, np.nan
+            )
+
     @helpers.requires_array_function_protocol()
     def test_fabs(self):
         helpers.assert_quantity_equal(