The low-impact-energy range of the data set used by the 3D Monte Carlo impurity transport code ERO has been supplemented by sputtering data calculated by molecular dynamics. Also reflection data for normal incidence were obtained. The computed sputtering data can be reasonably well fitted using the Bohdansky formula, yielding rather similar fitting parameter values to those used in ERO. A previously modelled case for ITER divertor has been simulated with the new data. As expected from the small change in the data set, the modelling results change only little. A method was also developed for direct coupling of ERO and the molecular dynamics code HCParcas, but the directly coupled code is not applicable to the present work.