Important quantities in aerodynamic flow simulations are the aerodynamic force coefficients including the pressure induced and the viscous stress induced drag, lift and moment coefficients. In addition to the exact approximation of these quantities it is of increasing importance, in particular in the field of uncertainty quantification, to estimate the error in the computed quantities. In recent years a posteriori error estimation and goal-oriented refinement approaches have been developed for the accurate and efficient computation of single target quantities. The current approaches are based on computing an adjoint solution related to each of the specific target quantities under consideration. In this paper we extend this approach to the accurate and efficient computation of multiple target quantities. Instead of computing multiple adjoint solutions, one for each target functional, the new approach is based on the solution to one discrete adjoint problem and one discrete error problem...