A numerical model capable of simulating a 2-D airfoil flying over in the vicinity of the waves is discussed. Instead of treating the problem as a heaving oscillation one above the rigid flat wall, sources are distributed on the prescribed wave profile. The wave deformation due to the airfoil is assumed to be negligible and treated as a rigid undulated wall. The source and vortex are distributed on the surface of the foil. It is found that the variation of $C_L$ with wave steepness in severe and that the lift variation due to waves decreases as the wing height above the water surface increases.