Added more files,
datfiles/008_S3cPoTf_angles/
datfiles/009_S3cPoTf_N128/
Working quarter circle
Initial commit for testing Po, with and without L-mult.
Files for testing Po and Tf, playing around with boundary conditions.
Check what the problem is before running the code!
Different variations of Tf.
Boundaries, goes anti-clockwise from the bottom:
0 - bottom
1 - right
2 - top (aka neumann)
3 - left (aka inflow)

P2v is pinning two different velocity node, at each ends of the neumann
boundary.

P2v1p - pin velocity at the two ends of the neumann boundary, and the
first pressure node.

P1v - pin only one end of the neumann boundary.

P1v1p - pin only one end of the neumann boundary, and the first
pressure.