The plasma equilibrium plays a crucial role in nuclear fusion studies, serving as the foundation for various aspects of fusion research, including plasma stability, transport, and current drive. In this paper, a new Grad-Shafranov equation solver is developed for the fixed-boundary plasma equilibria with toroidal flow. This solver utilizes the pressure profile, safety factor profile (not current profile), and any two profiles of the toroidal angular velocity, plasma temperature, and square of the Mach number as inputs. The numerical results obtained by this solver exhibit good agreement with k...