## Numerical calculation of equilibria with poloidal-sonic flow and finite Larmor radius effects in large aspect-ratio tokamaks

A code has been developed for calculating magnetohydrodynamic equilibria with poloidal-sonic flow and finite Larmor radius effects in high-beta tokamaks using an inverse aspect-ratio expansion and a reduced two-fluid model. The Grad-Shafranov equations governing the first- and second-order poloidal fluxes can be expressed in terms of five free profiles of the first-order poloidal flux. Sample equilibria, illustrating behaviors such as the deviation of pressure contours from the flux surfaces, and the criteria for the presence of the 'poloidal-sonic singularity' are presented.