In this article, I present our (Monte Carlo) simulation results of the angular distribution of the 3-dimensional WIMP velocity, in particular, a possible annual modulation and the diurnal modulation proposed in literature, for different underground laboratories.