SUMMARYA description is given of the effects of the discretization of the domain of integration on the accuracy of the results obtained with boundary-fitted finite differences. Three-dimensional unsteady heat conduction problems and two dimensional Navier-Stokes equations are considered. Comparisons with analytic solutions are given for all cases. It is shown that grid shapes influence the accuracy of the results and quantitative error evaluations are provided for some interesting cases.