We present a formulation of the two-dimensional complete flux (CF) scheme in terms of local orthogonal coordinates adapted to the flow, i.e., one coordinate axis is aligned with the local velocity field and the other one is perpendicular to it. This approach gives rise to an advection-diffusion-reaction boundary value problem (BVP) for the flux component in the local flow direction. For the other (diffusive) flux component we use central differences. We will demonstrate the performance of the scheme for several examples.