added possibility to use different time integration schemes for ion concentration for exponential integrator