We propose a general approach for the numerical approximation ofoptimal control problems governed by a linear advection–diffusionequation, based on a stabilization method applied to theLagrangian functional, rather than stabilizing the state andadjoint equations separately. This approach yields a coherentlystabilized control problem. Besides, it allows a straightforward a posteriori error estimate in which estimates of higher order termsare needless. Our a posteriori estimates stems from splitting theerror on the cost functional into the sum of an iteration errorplus a discretization error. Once the former is reduced below agiven threshold (and therefore the computed solution is “near”the optimal solution), the adaptive strategy is operated on thediscretization error. To prove the effectiveness of the proposedmethods, we report some numerical tests, referring to problems inwhich the control term is the source term of theadvection–diffusion equation.