Turbulent flow at high Reynolds numbers is currently not accessible on the basis of direct numerical simulation (DNS) of the Navier-Stokes equations - the computational complexity is too high to allow DNS in most realistic flow conditions. Instead, Large-Eddy Simulation (LES) offers an alternative in which the focus is on capturing the larger dynamic scales of a problem. However, the fundamental closure problem in LES induced by spatial filtering of nonlinear terms, and the role of discretization errors in the numerical treatment of the LES equations, induce a principal uncertainty in any LES prediction. This uncertainty requires quantification and control. We investigate error control capabilities of the Integral Length-Scale Approximation...