I’m not sure that I understand your question. The integral in that code is basically equivalent to on() with tgv=1e30. Also, you should definitely not be imposing any boundary conditions on the pressure.
I’m not sure that I understand your question. The integral in that code is basically equivalent to on() with tgv=1e30. Also, you should definitely not be imposing any boundary conditions on the pressure.