indexing zarr with both coordinate arrays and slice

[kanglcn]

The vectorized indexing on zarr currently do not support mixing coordinate arrays with slices. e.g.,

z = zarr.array(np.arange(3*4*5).reshape(3, 4, 5))
z.vindex[[0, 2], [1, 3]] #None in the last dim

It will have an error.

Any workaround about it?

I need to slice a big 4d array with coordinates in the first two dimensions. I don't want to load the whole data into numpy due to memory limit.

I also asked this question on stackoverflow: https://stackoverflow.com/questions/77226000/indexing-zarr-with-both-coordinate-arrays-and-slice