In this paper we consider the finite element approximation of the singularities of the solution of Poisson problems in a polygonal domain with reentrant corners or changing Dirichlet-Neumann boundary conditions. We use a correction algorithm with patches of elements to improve the a priori error estimates and to obtain the same order as the optimal estimate when everything is regular. We give an application of the correction method to the problem of glacier modeling.