res Subroutine

   pure subroutine res(t, y, yprime, cj, delta, ires, rpar, ipar)
   !! Residuals routine.
      real(rk), intent(in):: t
      real(rk), intent(in):: y(neq)
      real(rk), intent(in):: yprime(neq)
      real(rk), intent(in):: cj
      real(rk), intent(out):: delta(neq)
      integer, intent(inout) :: ires
      real(rk), intent(in):: rpar
      integer, intent(in) :: ipar

      call f(t, y, delta)
      delta = yprime - delta

   end subroutine res