converted julian.py to use Time from astropy
Simplified test_mean_anomaly using astropy.time.Time