We consider the functional form of the boundary condition for the Boltzmann equation (kinetic boundary condition) at an interface between a vapor and its condensed phase on the basis the of nonequilibrium molecular dynamics simulation. The molecular dynamics system treated here consists of two parallel condensed phases with different tempreratures and the vapor between them. We examine the velocity distribution function of molecules outgoing from the condensed phase in the above system, and thereby we establish a physically appropriate form of kinetic boundary condition at the interface.