Added some minor changes for units again.
Multiplied by 1000 wherever there are distance, x, y or z quantities.
Changed z to the heliocentric coordinate's zeta value in HCC to HGS.