It would be great to enforce having units for all input parameters.
We might consider using astropy units for this.