Understanding and control of impurity transport is one of the important issues to reduce the impurity in fusion plasmas. Being based on the Binary Collision Monte-Carlo Model (BCM), a numerical model for classical/neo-classical cross field transport of impurity ions in magnetic fusion devices is being developed. The purpose of the present study is to examine, step by step, whether our proposed model correctly reproduces 1) classical and 2) neo-classical transport processes of impurity ions. The numerical results agree well with theoretical values by classical theory. Not only self-diffusion, but also impurity flow in the direction along the background density gradient has been reproduced. In addition, good agreement of diffusion coefficient with neoclassical theory has been obtained in the wide range of collisionality parameter in a simple tokamak magnetic configuration.