Get the pressure p on the psi mesh (p in [mu0 * Pa]=[Tesla] units)
integer, intent(in) :: ns
real(i2mex_r8), intent(in) :: psi(ns)
real(i2mex_r8), intent(out) :: p(ns)
integer, intent(out) :: ier
    Get d p /d psi on the psi mesh. psi is the poloidal flux/(2*pi)
    Get d^2 p /d psi^2 on psi mesh
integer, intent(in) :: ns
real(i2mex_r8), intent(in) :: psi(ns)
real(i2mex_r8), intent(out) :: ppp(ns)
integer, intent(out) :: ier