Summary: Assuming that the gravity anomaly and disturbing potential are given on a reference ellipsoid, the author’s result in [Bull. Geod. 62, 93–101 (1988)] is applied to derive the potential coefficients on the bounding sphere of the ellipsoid to order \(e^2\) (i.e. the square of the eccentricity of the ellipsoid). By adding the potential coefficients and continuing the potential downward to the reference ellipsoid, the spherical Stokes formula and its ellipsoidal correction are obtained. The correction is presented in terms of an integral over the unit sphere with the spherical approximation of geoidal height as the argument and only three well-known kernel functions, namely those of Stokes, Vening-Meinesz and the inverse Stokes, lending the correction to practical computations. Finally, the ellipsoidal correction is presented also in terms of spherical harmonic functions. The frequently applied and sometimes questioned approximation of the constant m, a convenient abbreviation in normal gravity field representations, by \(e^2/2\), as introduced by Moritz, is also discussed. It is concluded that this approximation does not significantly affect the ellipsoidal corrections to potential coefficients and Stokes’ formula. However, whether this standard approach to correct the gravity anomaly agrees with the pure ellipsoidal solution to Stokes’ formula is still an open question.