The numerical simulation using the Lattice Boltzmann Method in the field of computational fluid dynamics becomes wider in the engineering applications because of its simplicity of update rules compared to the conventional Navier-Stokes solvers. Here, a two-dimensional D2Q9 LB model is numerically tested with a few new computational treatment on the free surface. The single relaxation time is applied under the gravitational field where applied only in the higher density fluid because of its big density difference. At the free surface, the reconstruction techniques in combination with boundary conditions is adopted in order to get some distribution function coming into the fluid site from the air one, and surface tension, early stable test for the gravitional field is considered in it. With the implementation of the gravitational profile, conserving the overall mass and grid dependency are observed during the calculations and freesurface advance track is well captured with an experiment.