Made changes for use of Astropy units.
1. In a transform function, changed zeta to D0 - d as it should be.
2. Removed the '* 1000' in favour of units manipulation.