The restitution of the trajectory of near-Earth orbiting satellites with accuracies to some centimeters by evaluating precise tracking data, is essential in satellite geodesy for applications like station positioning, earth orientation parameter determination, altimetry, SAR interferometry, and gravity field recovery from observed satellite orbit perturbations. Besides the stationary gravitational geopotential and non-gravitational forces, tidal induced temporal variations of the potential have to be considered in the numerical integration of the satellite's equation of motion.