Return the Grad-Shafranov operator (Del^*) in res
integer, intent(in) :: nt1, ns
real(i2mex_r8), intent(in) :: the(nt1)
real(i2mex_r8), intent(in) :: psi(ns)
real(i2mex_r8), intent(out) :: res(nt1, ns)
integer, intent(out) :: ier
    Return the right hand side (-p'*X^2 - g g') of the
    Grad-Shafranov equation in res
integer, intent(in) :: nt1, ns
real(i2mex_r8), intent(in) :: the(nt1)
real(i2mex_r8), intent(in) :: psi(ns)
real(i2mex_r8), intent(out) :: res(nt1, ns)
integer, intent(out) :: ier
    Return the Grad-Shafranov equation error in res. The error
    is normalized to R J_phi on axis.
integer, intent(in) :: nt1, ns
real(i2mex_r8), intent(in) :: the(nt1)
real(i2mex_r8), intent(in) :: psi(ns)
real(i2mex_r8), intent(out) :: res(nt1, ns)
integer, intent(out) :: ier