issue_comments

If people think this would be useful addition to the API then we could add it - or it could just be a cookbook recipe.

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.

@jthielen you should be able to adapt my example to check for being within a tolerance rather than strict equality.