To analyse an incompressible Navier4tokes flow problem in a boundary-fitted curvilinear co-odinatc system is definitely not a trivial task. In the primitive variable formulation, choices between wrhg variables and their storage points have to be made judiciously. The present work engages Contravariant velocity components and scalar pressure which stagger each other in the mesh to prevent even-odd pressure oscillations from emerging. Now that smoothness of the pnssure field is attainable, the remaining task is to ensure a disuete divergence-fire velocity field for an incompmsible flow simulation. Aside h m the flux discretiZations, the indqensable metric tensors, Jacobian and Chistoffel symbols in the transformed equations should be approximated with care. The guidmg idea is to get the property of geometric identity Pertaining to these grid-sensitive discretizatons. In addition, how to maintain the revertible one-toone equivalence at the discrete level between primitive and contravariant velocities is another theme in the present staggered formulation. A semi-implicit segregated solution algorithm felicitous for a large-scale flow simulation was utilized to solve the entire set of basic equations iteratively. Also of note is that the present segregated solution algorithm has the Virtue of reqUiring no userspecified relaxation pammeters for speedmg up the satisfaction of incompmsibility in an optimal sense. "hu benchmark problems, including an analytic problem, were investigated to j u s t i e the capability of the present formulation in handing problems with complex geometry. The test cases considered and the results obtained herein make a useful contribution in solving problems subsuming cells with arbitrary shapcs in a boundary-fitted grid system.