Made velocity analytic instead of numeric

Also renamed phi to psi since it wrongly named.  Now just need to use
get dfdx to find vorticity and divergence of velocity from analytic
velocity instead of velocity from analytic psi