The Grad-Shafranov equation is solved using spectral elements for tokamak equilibrium with toroidal rotation. The Grad-Shafranov solver builds upon and extends the NIMEQ code [Howell and Sovinec, Comput. Phys. Commun. 185 (2014) 1415] previously developed for static tokamak equilibria. Both geometric and algebraic convergence are achieved as the polynomial degree of the spectralelement basis increases. A new analytical solution to the Grad-Shafranov equation is obtained for Solov'ev equilibrium in presence of rigid toroidal rotation, in addition to a previously obtained analytical solution for a defferent set of equilibrium and rotation profiles. The numerical solutions from the extended NIMEQ are benchmarked with the analytical solutions, with good agreements. Besides, the extended NIMEQ code is benchmarked with the FLOW code [L. Guazzotto, R. Betti, et al., Phys. Plasma 11(2004)604].