DataArray.groupby drops empty coordinates

What happened?

When I run groupby on an empty DataArray, empty coordinates are dropped.
This is is xarray version 2024.10.0.

import numpy as np
import xarray as xr

data = xr.DataArray(np.empty((0, 2)), dims=['x', 'y'], coords={'x': [], 'y': [1, 1]})
print(data.groupby('y').sum())


<xarray.DataArray (x: 0, y: 1)> Size: 0B
array([], shape=(0, 1), dtype=float64)
Coordinates:
  * y        (y) int64 8B 1
Dimensions without coordinates: x

As you see, x coordinate is not longer present.

What did you expect to happen?

I expect x coordinate to be preserved.
E.g., here is the optuput of the same snippet in xarray==2024.2.0:

<xarray.DataArray (x: 0, y: 1)> Size: 0B
array([], shape=(0, 1), dtype=float64)
Coordinates:
  * x        (x) float64 0B 
  * y        (y) int64 8B 1

Minimal Complete Verifiable Example

import xarray as xr
xr.show_versions()

# your reproducer code ...
import numpy as np
import xarray as xr

data = xr.DataArray(np.empty((0, 2)), dims=['x', 'y'], coords={'x': [], 'y': [1, 1]})
assert 'x' in data.groupby('y').sum().coords

Steps to reproduce

No response

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.
• Recent environment — the issue occurs with the latest version of xarray and its dependencies.

Relevant log output

Anything else we need to know?

This issue is caused by the change introduced to DataArray._replace_maybe_drop_dims function in scope of https://github.com/pydata/xarray/issues/5361.

In this case self.dims == ('x', 'y') and variable.dims == ('y', 'x') so the condition set(self.dims) == set(variable.dims) is True.
However, the next line new_sizes = dict(zip(self.dims, variable.shape, strict=True)) assumes that the dims order is the same, otherwise it does not have any sense.

Environment

Start at DataArray._replace_maybe_drop_dims, identified in the issue as the source of the dimension-order handling, and reproduce the behavior with the provided empty DataArray example. Check the interaction between self.dims, variable.dims, and the constructed sizes, then verify the example's assertion that the empty x coordinate remains present. Add or update regression coverage for this case and run the relevant tests.