Manufactured solutions for SteadyState, axisy

Both, the 2D in cylinder coordinates test and
the corresponding 3D wedge in cartesian coordinates
are updated. Verified against analytical solution
given in both project files for initial conditions.

New references are provided for numerical comparison.

Since here two mesh resolutions were already available
the grid convergence is roughly quadratic for cylinder
coordinates and cartesian coordinates cases as expected.