 
     /// Set global attributes
     std::map<std::string, std::string> att;
-    att["title"] = "Output file of feltor/diag/feltordiag.cu";
+    att["title"] = "Output file of feltor/src/feltor/feltordiag.cu";
     att["Conventions"] = "CF-1.7";
     ///Get local time and begin file history
     auto ttt = std::time(nullptr);

    err = file::define_dimension( ncid_out, "psi", &dim_ids1d[1], g1d_out);
    //Write long description
    std::string long_name = "Time at which 2d fields are written";
    err = nc_put_att_text( ncid, tvarID, "long_name", long_name.size(),
            long_name.data());
    long_name = "Flux surface label";
    err = nc_put_att_text( ncid, dim_ids1d[1], "long_name",
        long_name.size(), long_name.data());

    std::map<std::string, int> id0d, id1d, id2d;

            if( record_name[0] == 'j')
                record_name[1] = 'v';
            //1. Read toroidal average
            int dataID =0;
            err = nc_inq_varid(ncid, (record.name+"_ta2d").data(), &dataID);
            err = nc_get_vara_double( ncid, dataID,
                start2d, count2d, transferH2d.data());
            //2. Compute fsa and output fsa
            dg::blas2::symv( grid2gridX2d, transferH2d, transferH2dX); //interpolate onto X-point grid
