Converted satvp and test to astropy.units instead of scimath.units