issue_comments: 610519628

Hi @lanougue , thanks for the suggestion!

If I understand correctly, you want to check that all elements are close along one dimension, and if so, then select only one index from that dimension? That seems to me to be two consecutive operations, the first of which is a reduction, and the second is just .isel: ```python da = xr.DataArray([[1.,2.],[1.,2.]], dims=('x','y'))

def reduce_if_constant_along_dim(da, dim): first = da.isel(**{dim: 0}) constant_along_dim = (da == first).all(dim)

true = xr.full_like(da, fill_value=True).isel(**{dim: 0}, drop=True)
if constant_along_dim.equals(true):
    return da.isel(**{dim: 0}, drop=True)
else:
    return da

print(reduce_if_constant_along_dim(da, dim='x')) bash <xarray.DataArray (y: 2)> array([1., 2.]) ```

or are you imagining something that applies the above function to every dim, more like: ```python def drop_constant_dims(da): for dim in da.dims: da = reduce_if_constant_along_dim(da, dim) return da

print(drop_constant_dims(da)) bash <xarray.DataArray (y: 2)> array([1., 2.]) `` There might be a slightly neater way usingreduce` somehow though.