Implement a switching mechanism for time integration
(for 1st-order implicit Euler and 2nd-order Crank-Nicolson)
Remove unnecessary subroutines in bcut
Merge branch 'time_integration'