Global models of the Earth gravity field and topographic/bathymetric data can be used for the gravimetric determination of the Moho discontinuity based on the Vening Meinesz-Moritz theory. In this paper, we mathematically develop this method in such a way that the local data can be used for Moho modelling. Two integral formulae are presented, one for integrating the data and one for their inversion. The kernels of both integrals are well-behaving meaning that the contribution of far-zone quantities being integrated are not very significant in the results. Both of these methods are applied for computing the Moho model of Iran and their results are compared to the Moho model determined based on the global models. Consistency of the computed Moho models from the simulated data and the global models verifies the correctness of both approaches. The presented methods are consistent even for the case of using real data. Numerical results show that the minimum value of the Moho models derived by the simulated data and global models are about 31 km, whilst those derived from the real data are about 3 km smaller. Similarly, the mean value of Moho depths derived from real data is about 1 km smaller than that from the global models.