You have to use a constraint and a Lagrange multiplier,
as in Error convergence on Stokes problem
The Lagrange multiplier \lambda is a constant in space, and \mu is the associated test constant, that sets the constraint. To do it correctly you have to define a block matrix as in the FreeFem doc p.671 “Lagrange Multipliers”.