Four-dimensional gauge theories based on symplectic Lie groups provide elegant realisations of the microscopic origin of several new physics models. Numerical studies pursued on the lattice provide quantitative information necessary for phenomenological applications. To this purpose, we implemented Sp(2N) gauge theories using Monte Carlo techniques within Grid, a performant framework designed for the numerical study of quantum field theories on the lattice. We show the first results obtained using this library, focusing on the case-study provided by the Sp(4) theory coupled to Nas = 4 Wilson-Dirac fermions transforming in the 2-index antisymmetric representation. In particular, we discuss preliminary tests of the algorithm and we test some of its main functionalities.