Our aim is to approximate a reference velocity field solving the two-dimensional Navier–Stokes equations (NSE) in the absence of its initial condition by utilizing spatially discrete measurements of that field, available at a coarse scale, and continuous in time. The approximation is obtained via numerically discretizing a downscaling data assimilation algorithm. Time discretization is based on semiimplicit and fully implicit Euler schemes, while spatial discretization (which can be done at an arbitrary scale regardless of the spatial resolution of the measurements) is based on a spectral Galerkin method. The two fully discrete algorithms are shown to be unconditionally stable, with respect to the size of the time step, the number of time steps and the number of Galerkin modes. Moreover, explicit, uniform-in-time error estimates between the approximation and the reference solution are obtained, in both the $L^2$ and $H^1$ norms. Notably, the two-dimensional NSE, subject to the no-slip Dirichlet or periodic boundary conditions, are used in this work as a paradigm. The complete analysis that is presented here can be extended to other two- and three-dimensional dissipative systems under the assumption of global existence and uniqueness.