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

            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);

        //////////////////////////write fields////////////////////////
        start[0] = i;
        err = nc_open(argv[2], NC_WRITE, &ncid);
        for(unsigned j = 0; j < 2; j++)
        {
            dg::blas2::symv( interpolate, y0[j], transferD);
