Refine the equilibrium using the ESC equilibrium code in p and q mode
integer, intent(in) :: mpol     number of Fourier harmonics
integer, intent(out) :: ier