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

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



    karniadakis.init( feltor, rolkar, y0, p.dt);
//     feltor.energies( y0);//now energies and potential are at time 0
    /////////////////////////////set up netcdf/////////////////////////////////////
    file::NC_Error_Handle err;
    int ncid;
    err = nc_create( argv[2], NC_NETCDF4|NC_CLOBBER, &ncid);
