Invariant Discretization of the Incompressible Navier-Stokes Equations in Boundary Fitted Coordinates