Maximum recursion calling `np.divide` with `DataArray` as `where` argument

What happened?

Calling np.divide() passing a bool expression involving a DataArray results in a RecursionError starting at apply_ufunc. Passing DataArray.data works fine.

What did you expect to happen?

No response

Minimal Complete Verifiable Example

import numpy as np
import xarray as xr

N = 10
array = np.random.random((N,))
da = xr.DataArray(np.random.random((N,)), dims=("x",), coords={"x": np.arange(N)})

# This works
print(np.divide(array, da, where=(da > 0.2).data))

# This gives recursion error
print(np.divide(array, da, where=da > 0.2))

MVCE confirmation

• Minimal example — the example is as focused as reasonably possible to demonstrate the underlying issue in xarray.
• Complete example — the example is self-contained, including all data and the text of any traceback.
• Verifiable example — the example copy & pastes into an IPython prompt or Binder notebook, returning the result.
• New issue — a search of GitHub Issues suggests this is not a duplicate.

Relevant log output

Traceback (most recent call last):
  File "divide_recursion_mvce.py", line 7, in <module>
    print(np.divide(array, da, where=da > 0.2))
          ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/arithmetic.py", line 86, in __array_ufunc__
    return apply_ufunc(
           ^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/computation.py", line 1197, in apply_ufunc
    return apply_dataarray_vfunc(
           ^^^^^^^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/computation.py", line 304, in apply_dataarray_vfunc
    result_var = func(*data_vars)
                 ^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/computation.py", line 761, in apply_variable_ufunc
    result_data = func(*input_data)
                  ^^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/arithmetic.py", line 86, in __array_ufunc__
    return apply_ufunc(
           ^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/computation.py", line 1210, in apply_ufunc
    return apply_array_ufunc(func, *args, dask=dask)
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

...

  File "venv/lib64/python3.11/site-packages/xarray/core/computation.py", line 1210, in apply_ufunc
    return apply_array_ufunc(func, *args, dask=dask)
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/computation.py", line 815, in apply_array_ufunc
    if any(is_chunked_array(arg) for arg in args):
       ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/computation.py", line 815, in <genexpr>
    if any(is_chunked_array(arg) for arg in args):
           ^^^^^^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/pycompat.py", line 99, in is_chunked_array
    return is_duck_dask_array(x) or (is_duck_array(x) and hasattr(x, "chunks"))
           ^^^^^^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/pycompat.py", line 95, in is_duck_dask_array
    return is_duck_array(x) and is_dask_collection(x)
                                ^^^^^^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/pycompat.py", line 87, in is_dask_collection
    if module_available("dask"):
       ^^^^^^^^^^^^^^^^^^^^^^^^
  File "venv/lib64/python3.11/site-packages/xarray/core/utils.py", line 1160, in module_available
    return importlib.util.find_spec(module) is not None
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "<frozen importlib.util>", line 103, in find_spec
  File "<frozen importlib._bootstrap>", line 1080, in _find_spec
  File "<frozen importlib._bootstrap>", line 924, in find_spec
  File "<frozen importlib._bootstrap>", line 241, in _call_with_frames_removed
RecursionError: maximum recursion depth exceeded while calling a Python object

Anything else we need to know?

No response

Environment

Run the minimal example and compare np.divide with DataArray as where against where=(da > 0.2).data. Read xarray/core/arithmetic.py and the apply_ufunc paths in xarray/core/computation.py shown in the traceback. Done means the DataArray where form no longer recurses and has regression coverage for the reproduced case.