Return in res the difference 'delta' between the
    PEST1 straight field line poloidal angle and 'the'
    on the (the, psi) grid.
integer, intent(in) :: nt1, ns
real(i2mex_r8), intent(in) :: the(nt1), psi(ns)
real(i2mex_r8), intent(out) :: res(nt1, ns)
integer, intent(out) :: ier
    Return in res the derivative d /dthe of
    'q*delta' where delta is the difference between
    the PEST1 straight field line poloidal angle and
    'the' the poloidal angle on the (the, psi) grid.
integer, intent(in) :: nt1, ns
real(i2mex_r8), intent(in) :: the(nt1), psi(ns)
real(i2mex_r8), intent(out) :: res(nt1, ns)
integer, intent(out) :: ier
    Return in res the derivative d /dpsi of
    'q*delta' where delta is the difference between
    the PEST1 straight field line poloidal angle and
    'the' the poloidal angle on the (the, psi) grid.
integer, intent(in) :: nt1, ns
real(i2mex_r8), intent(in) :: the(nt1), psi(ns)
real(i2mex_r8), intent(out) :: res(nt1, ns)
integer, intent(out) :: ier