The aim of the present study is to develop the method of determining the kinetic boundary condition (KBC) at a vapor-liquid interface in net evaporation/condensation. We proposed a novel method for determining the KBC by combining the numerical simulations of the mean field kinetic theory and the molecular gas dynamics. The method was evaluated on steady vapor flow between two liquid slabs at different temperatures. A uniform net mass flux in the vapor phase induced by net evaporation and condensation is obtained from the numerical simulation of the mean field kinetic theory for both vapor and liquid phases. The KBC was specified by using the uniform net mass flux, and the numerical simulation of the molecular gas dynamics was conducted for the vapor phase. Comparing the macroscopic variables in the vapor phase obtained from both numerical simulations, we can validate the KBC whether the appropriate solutions are obtained. Moreover, the evaporation and condensation coefficients were estimated uniquely. The results showed that the condensation and evaporation coefficients were identical and constant in net evaporation. On the other hand, in net condensation, the condensation coefficient increased with the collision molecular mass flux. We also presented the applicable limit of the KBC which is assumed to be the isotropic Gaussian distribution at the liquid temperature. From these results, the KBCs in net evaporation and condensation, which enable the exact macroscopic variables to be determined, were proposed.