 
                 tti.toc();
                 DG_RANK0 std::cout << " Time for internal diagnostics "<<tti.diff()<<"s\n";
+                if( abort) break;
             }
             ti.toc();
             //----------------Test if ampere equation holds
-            if( p.beta != 0)
+            if( p.beta != 0 && !abort)
             {
                 feltor.compute_lapMperpA( resultD);
                 double norm  = dg::blas2::dot( resultD, feltor.vol3d(), resultD);

                catch( dg::Fail& fail){ // a specific exception
                    DG_RANK0 std::cerr << "ERROR failed to converge to "<<fail.epsilon()<<"\n";
                    DG_RANK0 std::cerr << "Does simulation respect CFL condition?"<<std::endl;
                    abort_program();
                }
                catch( std::exception& fail) { // more generic exception
                    DG_RANK0 std::cerr << "ERROR in timestepper\n";
                    DG_RANK0 std::cerr << fail.what()<<std::endl;
                    abort_program();
                }
                dg::Timer tti;
                tti.tic();

                throw std::runtime_error( "timestepper: output-mode "+output_mode+" not recognized!\n");

        }
        bool abort = false;
        for( unsigned i=1; i<=maxout; i++)
        {
            dg::Timer ti;
