Integral constraint for a Poiseuille flow

Thank you for the suggestion. I got the following error

Can't add expressions with different shapes.
Traceback (most recent call last):
    Fc = (u - Constant(C0)) * ksi * dx + q*phi*ds(inlet)

In this topic they also added the following command but I do not know why

u.assign(Constant([C0,C0]))