Anomalous scaling behaviors such as significantly large growth exponent (<TEX>${\beta}$</TEX>) and small reciprocal of dynamic exponent (1/z) values for many molecular crystalline thin films have been reported. In this study, the variation of scaling exponent values and consequent growth behaviors of molecular thin films were more quantitatively analysed using a (1+1)-dimensional surface lateral diffusion model. From these simulations, influence of step edge barriers and grain boundaries of molecular thin films on the various scaling exponent values were elucidated. The simulation results for the scaling exponents were also well consistent with the experimental data for previously reported molecular thin film systems.