rt Subroutine

   pure subroutine rt(neq, t, y, yprime, nrt, rval, rpar, ipar)
     !! Roots routine.
      integer, intent(in) :: neq
      real(rk), intent(in) :: t
      real(rk), intent(in) :: y(neq)
      real(rk), intent(in) :: yprime(neq)
      integer, intent(in) :: nrt
      real(rk), intent(out) :: rval(nrt)
      real(rk), intent(in) :: rpar
      integer, intent(in) :: ipar

      rval(1) = yprime(1)
      rval(2) = log(y(1)) - 2.2491_rk

   end subroutine rt