Turbulence in the presence of rotation is numerically investigated with MRT-lattice Bohzmann method. Influence of rotation (Rossby number and Ekman number) on turbulence is studied in detail. Numerical results include turbulence energy, dissipation rate, iso-surface of velocity, vortex structure, dissipation length and integral length and so on. It shows that rotation slows down turbulence energy decaying rate. Moreover, the initial homogeneity of turbulence is broken down in the presence of rotation. Vortices rotating in opposite direction of the system are suppressed gradually and columnar vortices appear rotating in the same direction of the system eventually. In addition, it was found that Kolmogorov length decreases as rotation increases while integral length increases.