We demonstrate a new nonuniform mesh finite difference method to obtain accurate solutions for the elliptic partial differential equations in two dimensions with nonlinear first-order partial derivative terms. The method will be based on a geometric grid network area and included among the most stable differencing scheme in which the nine-point spatial finite differences are implemented, thus arriving at a compact formulation. In general, a third order of accuracy has been achieved and a fourth-order truncation error in the solution values will follow as a particular case. The efficiency of using geometric mesh ratio parameter has been shown with the help of illustrations. The convergence of the scheme has been established using the matrix analysis, and irreducibility is proved with the help of strongly connected characteristics of the iteration matrix. The difference scheme has been applied to test convection diffusion equation, steady state Burger’s equation, ocean model and a semi-linear elliptic equation. The computational results confirm the theoretical order and accuracy of the method.