In step 3 the second correction to the momentum is determined. The lumped
equations are solved and the solution is added to
in routine resultsv2.f (again in fields v(1,*), v(2,*) and v(3,*)). Finally,
the fluid boundary conditions in the form
are applied and
is
changed appropriately.