The data I read with xarray with rasterio engine is inverted along Y axis (This is SMAPL4 data). The correct dimension values are y=ds.y*(-1).
ds = xr.open_mfdataset(file_paths, engine="rasterio", chunks=chunks, combine='nested', concat_dim='time', parallel=True)
How can I fix this, without messing up the chunks?
Debug attempt 1
ds = ds.reindex(y=ds.y*(-1))
Result: The chunking over the y-axis disappears. Plus, numeric values in the variable data field somehow disappears.
Debug attempt 2:
ds= ds.isel(y=slice(None, None, -1))
Result: This code just re-selects the grids and does not change the dimension values themselves
Debug attempt 3:
ds= ds.reindex_like(ds_template)
Result: Again, numeric values within precipitation data disappear
Debug attempt 4: Change engine from rasterio to default (auto-detect).
ds = xr.open_mfdataset(file_paths, chunks=chunks, combine='nested', concat_dim='time', parallel=True)
Result: My file has some issues with the default engine; the data is read as "attributes."