Numerical solutions of the Navier-Stokes equations, governing 2-dimensional, time-dependent, viscous, incompressible fluid flow past a square cylinde in an infinite region, are presented for Reynolds numbers $10^2$, $10^3$and $10^4$. Finite-difference scheme, based on SOLA-VOF is adopted and a discretization of the convection term for irregular grid is newly suggested by altering the original nonconservation form into conservation one. Distribution of finer grids around the body reveals fairly reasonable consistency with the experimental variables : drag coefficient, lift coefficient, Strouhal number, fluctuating pressure coefficient, etc.