Add more accurate heliocentric angle calculation function, move get_solar_radius to calculations module