     dg::HVec hvisual10 = dg::evaluate( zonalflow, grid);
     dg::HVec hvisual11 = dg::evaluate( gradient, grid);
     dg::HVec hvisual12 = dg::evaluate( field, grid);
+    dg::HVec hvisual13 = dg::evaluate( prof, grid);
+    dg::HVec hvisual14 = dg::evaluate( damp2, grid);
+    dg::HVec hvisual15 = dg::evaluate( dampcut, grid);
+    dg::HVec hvisual16 = dg::evaluate( source, grid);
 
     //allocate mem for visual
     dg::HVec visual1( grid.size());

    dg::HVec visual10( grid.size());
    dg::HVec visual11( grid.size());
    dg::HVec visual12( grid.size());
    //make equidistant grid from dggrid
    dg::HMatrix equigrid = dg::create::backscatter(grid);
    //evaluate on valzues from devicevector on equidistant visual hvisual vector

    solovev::Iris iris(gp);
    solovev::Pupil pupil(gp);
    solovev::Damping damping(gp);
    solovev::ZonalFlow zonalflow(gp,0.5);
    solovev::Gradient gradient(gp);
    solovev::Nprofile prof(gp);
    solovev::TanhDamping damp2(gp);
    solovev::TanhDampingProf dampcut(gp);
    solovev::TanhSource source(gp,40.);
    
//     //make dggrid
    std::cout << "kR "<<curvatureR(Rmin,Zmin) <<"\n";

    dg::Grid2d<double> grid(Rmin,Rmax,Zmin,Zmax, n,Nx,Ny,dg::PER,dg::PER);
    
    dg::HVec hvisual1 = dg::evaluate( psip, grid);
