Fixed wrong finite difference in div(phi*vi*v)