Molecular reaction dynamics of Cl + H2 (D2) has been studied on the latest analytical potential energy surface called BW3 using the Monte Carlo quasi-classical trajectory method. Excitation functions, differential cross sections and angular distributions of HCl and DCl products have been calculated. The excitation functions of the Cl (2 P 3/2) + n-H2 and Cl(2 P 3/2) + n-D2 reactions are also studied. The results are compared with those of quasi-classical trajectory [M. Alagia et al.: Phys. Chem. Chem. Phys. 2 (2000); F. J. Aoiz et al.: J. Phys. Chem. 100 (1996)], quantum mechanical (QM) calculations [F. J. Aoiz et al.:J. Chem. Phys. 115 (2001)] and experimental data [S. H. Lee et al.: J. Chem. Phys. 110 (1999); F. Dong et al.: J. Chem. Phys. 115 (2001)]. Discussions are given to some new results.