The Navier-Stokes equations are solved numerically for two-dimensional laminar flow of incompressible Newtonian fluid across the arrays of circular or elliptic cylinders. A square mesh is used to effect a discrete representation of the physical field. This necessitates the use of special treatment at mesh points near the curved boundary and also involves a new method of determining the surface vorticity using values of the stream function at various mesh points.