             radtrans = feltor.radial_transport();
             coupling= feltor.coupling();
             err = nc_open(argv[2], NC_WRITE, &ncid);
-            ns_ncid :: ncid = ncid;
             err = nc_put_vara_double( ncid, EtimevarID, Estart, Ecount, &time);
             err = nc_put_vara_double( ncid, energyID, Estart, Ecount, &E1);
             err = nc_put_vara_double( ncid, massID,   Estart, Ecount, &mass);

#include <vector>
#include <sstream>
#include <cmath>
#include <csignal>
// #define DG_DEBUG



            std::cout << "(E_tot-E_0)/E_0: "<< (E1-energy0)/energy0<<"\t";
            std::cout <<" d E/dt = " << dEdt <<" Lambda = " << diss << " -> Accuracy: "<< accuracy << "\n";
            err = nc_close(ncid);
        }
#ifdef DG_BENCHMARK
        ti.toc();
