Improve sorting performance
At the moment, we first stack the selected data into one grand Dask array before we sort to produce the order statistics. This is simpler in some sense, but it is probably more efficient to first sort the individual arrays per day-of-year and assemble those into the overall data that we need.