Roots routine.
Type | Intent | Optional | 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(*) |
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