The precise regional geoid modelling requires combination of terrestrial gravity data with satellite-only Earth Gravitational Models (EGMs). In determining the geoid using the Stokes-Helmert approach, the relative contribution of terrestrial and satellite data to the computed geoid can be specified by the Stokes integration cap size defined by the spherical distance ψ ₀ and the maximum degree l ₀ of the EGM-based reference spheroid. Larger values of l ₀ decrease the role of terrestrial gravity data and increase the contribution of satellite data and vice versa for larger values of ψ ₀. The determination of the optimal combination of the parameters l ₀ and ψ ₀ is numerically investigated in this paper. A numerical procedure is proposed to find the best geoid solution by comparing derived gravimetric geoidal heights with those at GNSS/levelling points. The proposed method is tested over the Auvergne geoid computation area. The results show that despite the availability of recent satellite-only EGMs with the maximum degree/order 300, the combination of l ₀ = 160 and ψ ₀ =  45 arc-min yields the best fitting geoid in terms of the standard deviation and the range of the differences between the estimated gravimetric and GNSS/levelling geoidal heights. Depending on the accuracy of available ground gravity data and reference geoidal heights at GNSS/levelling points, the optimal combination of these two parameters may be different in other regions.