Use astropy.units and xarray

As suggested at the plasma hack week (https://hack.plasmapy.org/), it would be nice for Aurora to use astropy.units to ensure consistent use of units. 

More widespread use of xarrays would also be valuable. 

This would be a great project for anyone hoping to begin contributing to Aurora.