Function that interfaces with the ntemp functions provided by the individual integrators. This function can be called before choose_integrator is called.
The number of temporary storage levels needed.
function mol_ntmp ( ) result ( ntemp )
!! Function that interfaces with the ntemp functions provided by the
!! individual integrators. This function can be called before
!! [[choose_integrator]] is called.
use parameters, only : mol_integrator
implicit none
integer(ip) :: ntemp
!! The number of temporary storage levels needed.
select case ( mol_integrator )
case ( 'rk4' )
ntemp = rk4_int%ntemp ( )
case ( 'rk5' )
ntemp = rk5_int%ntemp ( )
case ( 'abmv5' )
ntemp = abmv5_int%ntemp ( )
case default
print*,'Invalid integrator requested'
stop
end select
end function mol_ntmp