On the use of the unstable manifold correction in a Picard iteration for the solution of the velocity field in higher-order ice-flow models