  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
spectral-element 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].
