     ///////////////////////////////////////Timeloop/////////////////////////////////
     const double mass0 = exp.mass(), mass_blob0 = mass0 - grid.lx()*grid.ly();
     double E0 = exp.energy(), E1 = 0, diff = 0;
+    unsigned failed_counter = 0;
+    double dt = 1e-6;
     dg::Timer t;
     t.tic();
     try

#endif//DG_BENCHMARK
        for( unsigned j=0; j<p.itstp; j++)
        {
            karniadakis.step( exp, imp, time, y1);
            //store accuracy details
            {
                DG_RANK0 std::cout << "(m_tot-m_0)/m_0: "<< (exp.mass()-mass0)/mass_blob0<<"\t";

    toefl::Implicit< dg::x::CartesianGrid2d, dg::x::DMatrix, dg::x::DVec > imp( grid, p.nu);
    /////////////////////create initial vector////////////////////////////////////
    dg::Gaussian g( p.posX*p.lx, p.posY*p.ly, p.sigma, p.sigma, p.amp);
    std::vector<dg::x::DVec> y0(2, dg::evaluate( g, grid)); // n_e' = gaussian
    dg::blas2::symv( exp.gamma(), y0[0], y0[1]); // n_e = \Gamma_i n_i -> n_i = ( 1+alphaDelta) n_e' + 1
    if( p.equations == "gravity_local" || p.equations == "gravity_global" || p.equations == "drift_global"){
        y0[1] = dg::evaluate( dg::zero, grid);
    }
    //////////////////initialisation of timekarniadakis and first step///////////////////
    double time = 0;
    //dg::ImExMultistep_s< std::vector<dg::x::DVec> > karniadakis( "ImEx-BDF-3-3", y0, y0[0].size(), p.eps_time);
    //karniadakis.init( exp, imp, time, y0, p.dt);
    dg::Adaptive<dg::ERKStep<std::vector<dg::x::DVec>>> stepper( "ARK-4-2-3 (explicit)", y0);
    /////////////////////////////set up netcdf/////////////////////////////////////
    dg::file::NC_Error_Handle err;
    int ncid;
