  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.
