A few fixes to support astropy 4.0
This PR fixes a few errors when running with astropy
4.0, mainly to do with comparing units, all fixed by using regular numpy arrays in that particular operation.
This PR fixes a few errors when running with astropy
4.0, mainly to do with comparing units, all fixed by using regular numpy arrays in that particular operation.