     std::map<std::string, const dg::DVec* > v4d;
     v4d["ne-1 / "] = &y0[0][0],               v4d["ni-1 / "] = &y0[0][1];
     v4d["Ue / "]   = &feltor.fields()[1][0],  v4d["Ui / "]   = &feltor.fields()[1][1];
-    v4d["Omega / "] = &feltor.potential()[0]; v4d["Apar / "] = &feltor.induction();
+    v4d["Phi / "] = &feltor.potential()[0]; v4d["Apar / "] = &feltor.induction();
     const feltor::Quantities& q = feltor.quantities();
     double dEdt = 0, accuracy = 0, dMdt = 0, accuracyM  = 0;
     std::map<std::string, const double*> v0d{

                visual.begin(), visual.end(), 0., thrust::maximum<double>() );
            colors.scalemin() = -colors.scalemax();
            title <<pair.first << colors.scalemax()<<"\t";
            for( unsigned k=0; k<p.Nz/red;k++)
            {
                unsigned size=grid.n()*grid.n()*grid.Nx()*grid.Ny();
                dg::HVec part( visual.begin() +  k*red   *size,
                               visual.begin() + (k*red+1)*size);
                render.renderQuad( part, grid.n()*grid.Nx(),
                                         grid.n()*grid.Ny(), colors);
            }
            dg::blas1::scal(avisual,0.);
            toroidal_average(visual,avisual);
            render.renderQuad( avisual, grid.n()*grid.Nx(),
                                        grid.n()*grid.Ny(), colors);
        }
        title << std::fixed;
        title << " &&   time = "<<time;

    {
        for( auto pair : v4d)
        {
            if(pair.first == "Phi / ")
            {
                //dg::blas2::gemv( laplacianM, *pair.second, dvisual);
                //dg::assign( dvisual, hvisual);
                dg::assign( *pair.second, hvisual);
            }
            else if(pair.first == "ne-1 / " || pair.first == "ni-1 / ")
            {
