Often we have many arrays stored on disk that we want to stack together and
think of as one large array. This is common with geospatial data in which we
might have many HDF5/NetCDF files on disk, one for every day, but we want to do
operations that span multiple days.

To solve this problem we use the functions da.stack and da.concatenate.