model_sundials.f90

Go to the documentation of this file.

!********************************beginning*************************************
subroutine  fcvfun(t, y, ydot, ipar, rpar, ier)
    use iso_c_binding
    use constants, only:dp
    use integration, only:neq
    use model, only: odesystem
    implicit none
    integer :: ier
    integer(c_long) :: ipar(1)
    real(dp) :: t
    real(dp) :: y(neq)
    real(dp) :: ydot(neq)
    real(dp) :: rpar(1)
    call odesystem(neq, t, y, ydot)
end subroutine fcvfun
!***********************************end****************************************