Fix to issue #579 using astropy.units