modified GRmercury() to include only solar interactions. Also added a
constraint starting on line-269 which outputs a specified range of a
multi-period orbit. This allows the procession of the perihelion to be
easily tracked by only outputting the last orbit of a specified kmax.