add two coarser mesh files and switch current dense mesh to coarser one.
The new mesh is not working yet. Use the old mesh for now.
The calculated xi_kin_bar values need to be taken into account.
The number of unknowns in the local problem should be the same as number of components.
print each time step result for debugging purpose.
give a bit loose criteria.
switch on the p1F and p2F calculation.