Abstract : The paper indicates how a second order accurate two step method used in the numerical computation of hydrodynamic flows may be derived directly from the integral conservation laws. The result of Anderson, Preiser and Rubin makes it immediately applicable to arbitrary orthogonal curvilinear coordinate systems. A necessary and sufficient condition for stability is derived for the three-dimensional Cartesian coordinate case. (Author)