Implement boolean array indexing

[tomwhite]

https://data-apis.org/array-api/latest/API_specification/indexing.html#boolean-array-indexing

The output has a data-dependent output shape, and it's difficult to implement since each chunk can be an arbitrary size, yet Zarr doesn't support non-regular chunk sizes.

It may be simplest to wait until non-regular chunking is available in Zarr 3.0 before doing this. There's a sketch of how this might be implemented here: https://github.com/zarr-developers/zarr-specs/issues/49#issuecomment-611473418

Using Zarr 3.0 non-regular chunking could be an implementation detail, by ensuring that it is only used for intermediate data, and that the final output is always a Zarr 2.0 array with regular chunking - at least until Zarr 3.0 is more widely supported.

[tomwhite]

Implementing nonzero and the unique_* functions (https://github.com/tomwhite/cubed/blob/main/api_status.md) also need support for non-regular chunking.

[hammer]

It looks like Xarray does not yet have this feature either? https://github.com/pydata/xarray/issues/1887