A meshless model, based on the Meshless Local Petrov-Galerkin (MLPG) approach, is developed and implemented in parallel for the solution of axi-symmetric poroelastic problems. The parallel code is based on a concurrent construction of the stiffness matrix by the processors and on a parallel preconditioned iterative method of Krylov type for the solution of the resulting linear system. The performance of the code is investigated on a realistic application concerning the prediction of land subsidence above a deep compacting reservoir. The overall code is shown to obtain a very high parallel efficiency (larger than 78% for the solution phase) and it is successfully applied to the solution of a poroelastic problem with a fine discretization which produces a linear system with more than 6 million equations using up to 512 processors on the HPCx supercomputer.