rt Subroutine

public pure subroutine rt(neq, t, c, cdot, nrt, rval, rpar, ipar)

Roots routine.

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: neq
real(kind=rk), intent(in) :: t
real(kind=rk), intent(in) :: c(neq)
real(kind=rk), intent(in) :: cdot(neq)
integer, intent(in) :: nrt
real(kind=rk), intent(out) :: rval(nrt)
real(kind=rk), intent(in) :: rpar(*)
integer, intent(in) :: ipar(*)

Calls

proc~~rt~3~~CallsGraph proc~rt~3 rt proc~c1_average c1_average proc~rt~3->proc~c1_average

Source Code

   pure subroutine rt(neq, t, c, cdot, nrt, rval, rpar, ipar)
   !! Roots routine.

      integer, intent(in) :: neq
      real(rk), intent(in) :: t
      real(rk), intent(in) :: c(neq)
      real(rk), intent(in) :: cdot(neq)
      integer, intent(in) :: nrt
      real(rk), intent(out) :: rval(nrt)
      real(rk), intent(in) :: rpar(*)
      integer, intent(in) :: ipar(*)

      real(rk) :: c1ave

      call c1_average(c, c1ave)
      rval(1) = c1ave - 20.0_rk

   end subroutine rt