turning some functions to internal, and making sure that these functions return astropy Quantities.
made the true obliquity an astropy Angle Quantity
more edits to differential rotation
implementing astropy units in the rotation code