Better integration with astropy.units

[pllim]

Would be nice if the code can understand astropy.units.Quantity directly instead of assuming given value must have certain unit, and in general, better integration with astropy.units.