xarray would allow more comprehendible calls to data for slicing versus n-dimensional arrays.
The project is already outgrowing the n-dimensional array paradigm.
Converting to xarray would be beneficial, but should be layered as an optional way of interacting so as to still allow array interaction with calls.
This means "duplicating" functions (choosing an interface) or multiple ways to call or somethign similar. The model of the python control (looks like matlab, or doesn't depending no desire) is a paradigm.