html_url,issue_url,id,node_id,user,created_at,updated_at,author_association,body,reactions,performed_via_github_app,issue
https://github.com/pydata/xarray/issues/6560#issuecomment-1118690026,https://api.github.com/repos/pydata/xarray/issues/6560,1118690026,IC_kwDOAMm_X85Crdrq,2448579,2022-05-05T15:21:22Z,2022-05-05T15:21:40Z,MEMBER,"OK it has to do with the lazy backend. Calling `data_example.load()` fixes things. @jesieleo in the mean time, if your data is small use `load_dataset` or `load` to fix it.

Though I still can't make a minimal example. The example below works fine.

``` python
import numpy as np
import xarray as xr
ds = xr.Dataset(
    {
        ""z"": (
            (""time"", ""isoBaricInhPa"", ""latitude"", ""longitude""),
            np.ones((1, 5, 721, 1440)),
        )
    },
    coords={""latitude"": np.linspace(-90, 90, 721)},
).to_netcdf(""test.nc"")
ds = xr.open_dataset('test.nc')
ds.isel(time=[0], isoBaricInhPa=1).z[:, ::10, :][:, ::-1, :]
```","{""total_count"": 1, ""+1"": 1, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,1222609313
https://github.com/pydata/xarray/issues/6560#issuecomment-1118682989,https://api.github.com/repos/pydata/xarray/issues/6560,1118682989,IC_kwDOAMm_X85Crb9t,2448579,2022-05-05T15:16:33Z,2022-05-05T15:16:52Z,MEMBER,"OK turns out my attempts at a minimal example didn't work.

There's some bug here (on xarray main)
1.  the size of `latitude` is different in `.sizes, .shape, .data.shape` and 
2. latitude really should be 73 elements. I don't see why a subsetting happens when slicing with `::-`
<img width=""844"" alt=""image"" src=""https://user-images.githubusercontent.com/2448579/166955203-ff46834a-abdf-4ea0-bbd5-4a2e0b0aebbd.png"">


-----

@jesieleo here's how I tried to create a ""minimal example""

``` python
import xarray as xr
import numpy as np

ds = xr.Dataset(
    {
        ""z"": (
            (""time"", ""isoBaricInhPa"", ""latitude"", ""longitude""),
            np.ones((1, 5, 721, 1440)), # use numpy.ones to create array of same shape.
        )
    },
    coords={""latitude"": np.linspace(-90, 90, 721)},
)
ds.isel(isoBaricInhPa=1).z[:, ::10, :][:, ::-1, :]
```

Though this doesn't replicate it!

","{""total_count"": 1, ""+1"": 1, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,1222609313