Inputs and outputs with astropy units

[mchalela]

Inputs should accept astropy Quantity with units. The outputs should have units only if inputs have units.