res Subroutine

   pure subroutine res(t, y, ydot, cj, delta, ires, rpar, ipar)
   !! Residuals routine.

      real(rk), intent(in):: t
      real(rk), intent(in):: y(neq)
      real(rk), intent(in):: ydot(neq)
      real(rk), intent(in):: cj
      real(rk), intent(out):: delta(neq)
      integer, intent(out) :: ires
      real(rk), intent(inout):: rpar
      integer, intent(inout) :: ipar

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

   end subroutine res