Add julian_date() based on NOVAS function
Simple 3D plot of the inner solar system for index