Get Jacobian 'jac' on the new (the, psi) mesh
integer, intent(in) :: nt1, ns
real(i2mex_r8), intent(in) :: the(nt1), psi(ns)
real(i2mex_r8), intent(out) :: jac(nt1, ns)
integer, intent(out) :: ier
    Get d Jacobian /d theta 'jt' on the new (the, psi) mesh
integer, intent(in) :: nt1, ns
real(i2mex_r8), intent(in) :: the(nt1), psi(ns)
real(i2mex_r8), intent(out) :: jt(nt1, ns)
integer, intent(out) :: ier
    Get d Jacobian /d psi 'jp' on the new (the, psi) mesh
integer, intent(in) :: nt1, ns
real(i2mex_r8), intent(in) :: the(nt1), psi(ns)
real(i2mex_r8), intent(out) :: jp(nt1, ns)
integer, intent(out) :: ier