Here at Pitt, we are trying to solve the 2D Taylor Green Vortex solution from the tutorial and study the effect of resolution on L2 norm. As it is shown in the attached file, it appears that there is no difference at all between p=5 and p=6. It also looks like that the L2-error cannot be reduced below approx. 7e-6. I assume that the L2 norm provided at the end of simulation output is the comparison to the exact solution of the TGV. Are we missing something?