The celestial mechanics approach (CMA) has its roots in the Bernese GPS software and was extensively used for determining the orbits of high-orbiting satellites. The CMA was extended to determine the orbits of Low Earth Orbiting satellites (LEOs) equipped with GPS receivers and of constellations of LEOs equipped in addition with intersatellite links. In recent years the CMA was further developed and used for gravity field determination. The CMA was developed by the Astronomical Institute of the University of Bern (AIUB). The CMA is presented from the theoretical perspective in (Beutler et al. 2010). The key elements of the CMA are illustrated here using data from 50 days of GPS, K-Band, and accelerometer observations gathered by the Gravity Recovery And Climate Experiment (GRACE) mission in 2007. We study in particular the impact of (1) analyzing different observables [Global Positioning System (GPS) observations only, inter-satellite measurements only], (2) analyzing a combination of observations of different types on the level of the normal equation systems (NEQs), (3) using accelerometer data, (4) different orbit parametrizations (short-arc, reduced-dynamic) by imposing different constraints on the stochastic orbit parameters, and (5) using either the