Unexpected chunking of the dataset being returned by resample_in_space #1124
Description
Describe the bug
Assume we have 3d dataset (e.g. data variables have dimensions (time, y, x)) with ordinary chunking of (1, 1024, 1024). If this data set is put into resample_in_space method returns a dataset with one one chunk along the time axis, meaning that the cunksize for time is equal to the axis size of time.
Expected behavior
The chunking along the non-spatial axis should remain unchanged from the input dataset.
Screenshots
