atan($neptune_orb_all_HEE[2]/sqrt(pow($neptune_orb_all_HEE[1],2) + pow($neptune_orb_all_HEE[0],2)))*RAD2DEG