Following from the conversation here.
Having constant coefficient advection terms on the same side as the time derivative causes the wrong scheme to be used from the perspective of if that term appeared on the other side of the equation.
Casting equations to a set form, nominally eq.lhs - eq.rhs ~ 0 pre discretization and discretizing from this perspective would ensure consistency.