Compute the Glasser, Greene and Johnson coefficients
    E, F, and H on the set of poloidal surface psi
integer, intent(in) :: ns
real(i2mex_r8), intent(in) :: psi(ns)
real(i2mex_r8), intent(out) :: e(ns), f(ns), h(ns)
integer, intent(out) :: ier