Essential imposition of Neumann condition in Galerkin-Legendre elliptic solvers