Reinstated only updating domain every so many iterations. Also, deleted begin and end messages which were previously used for debugging.
Fixed all solve modes for Stokes and Navier-Stokes without the stability terms. A quick recheck of the adjoints should be made in case using grad_u_old_old is no longer necessary. Stability terms are definitely still broken and need total rehaul of all parts of assembly.