Numerical simulations for the flows containing free surface remain difficult problems because the drastic differences of physical properties of water and air, The difference of densities makes the solution instable in particular. For the stabilities of the solutions, the most typical methods to simulate free surface flows, such as Volume Of Fluid(VOF) and Level-Set(LS) methods, impose transient zones where the physical prosperities are continuously distributed. The thickness of the transient zone is the source of the numerical errors. The other side, marker-density method does not use such a transient zone. In the traditional marker-density method, however, the air velocities of free surface cells are extrapolated from the water velocity, and the pressures on the free surface are extrapolated from the air pressures for the stability of the solution. In this study, the marker-density method is modified for the decrease of such numerical errors. That is, the pressure on the free surface is determined to coincide with the pressure gradient terms of the governing equations, and the velocity of free surface cells are calculated with the governing equations. Two-dimensional steady spilling breakers behind of a submersed hydrofoil and three-dimensional spilling breaker near a wedge shaped ship model are simulated using INHAWAVE-II including the modified marker-density(MMD) method. The results are compared with the results of Fluent V6.3 including VOF method and several published research results.