                 double max_ue = dg::blas1::reduce(
                     feltor.velocity(0), 0., dg::AbsMax<double>() );
                 std::cout << "\tMaximum ue "<<max_ue<<"\n";
-                if( p.timestepper == "adaptive")
+                if( p.timestepper == "adaptive" || p.timestepper == "adaptive-imex")
                 {
                     std::cout << "\tdt "<<dt<<"\n";
                     std::cout << "\tfailed "<<var.nfailed<<"\n";

    dg::Adaptive< dg::ARKStep< std::array<std::array<dg::x::DVec,2>,2>,
        feltor::ImplicitSolver< dg::x::CylindricalGrid3d, dg::x::IDMatrix,
            dg::x::DMatrix, dg::x::DVec> > > adapt_ark;
    double rtol = 0., atol = 0., dt = 0.;
    if( p.timestepper == "multistep")
    {

                        if( time + dt > t_output)
                            dt = t_output-time;
                        if( p.timestepper == "adaptive")
                        {
                            double dt_old;
                            do{
                                dt_old = dt;
                                adapt.step( feltor, time, y0, time, y0, dt,
                                        dg::pid_control, dg::l2norm, rtol,
                                        atol);
