This is a short test which is supposed to check that the velocity is divergence free.

It is a simple squared box with initial zero velocity and pressure, no normal flow bc on all boundaries and only driven by a circled shaped force.

Unfortunately, the divergence is NOT even close to zero even after the first timestep. This might be caused by a simple mistake how this testcase is set up or a bug. 
Anyway, I think its an essential feature and should be checked in a test case (the only testcase which does a divergence check is geostrophic_interpolation_gp and there the tolerance is set to 3.0e-5, but m$

For now the tolerances are set really high to pass, but eventually they will be corrected. 

I have no clue whats going wrong here, so feedback is welcome!

Settings:
P1P1
theta=1.0
no kmk stabilisation