Drag modeling is the greatest uncertainty in the dynamics of low Earth satellite orbits where ballistic coefficient and density errors dominate drag errors. This paper examines fitted drag coefficients found as part of a precision orbit determination process for Stella, Starlette, and the GEOSAT Follow-On satellites from 2000 to 2005. The drag coefficients for the spherical Stella and Starlette satellites are assumed to be highly correlated with density model error. The results using MSIS-86, NRLMSISE-00, and NRLMSISE-00 with dynamic calibration of the atmosphere (DCA) density corrections are compared. The DCA corrections were formulated for altitudes of 200–600 km and are found to be inappropriate when applied at 800 km. The yearly mean fitted drag coefficients are calculated for each satellite for each year studied. The yearly mean drag coefficients are higher for Starlette than Stella, where Starlette is at a higher altitude. The yearly mean fitted drag coefficients for all three satellites decrease as solar activity decreases after solar maximum.