A multigrid method for an invariant formulation of the incompressible Navier-Stokes equations in general coordinates