A calculation for spontaneous plasma rotation based on a full ion distribution function (full-f) approach is carried out for H-mode plasmas. It is found that both equilibrium pressure gradient and fluctuation-induced Reynolds stress contribute to the spontaneous toroidal rotation. While the equilibrium diamagnetic drift overwhelmingly dominates the rotation in the edge plasma regions, the turbulent Reynolds stress generated rotations distribute mostly in the intermediate and core plasma regions. Numerical results based on this approach accord reasonably well with the experimentally measured profiles.