The authors propose a numerical solution for the discrete exterior Dirichlet boundary value problem of gravitational potential with respect to the real Earth surface. Using a spline approach for the real Earth, the approximating spline has the same domain of harmonicity as the real gravitational potential.