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

    solovev::Iris iris(gp);
    solovev::Pupil pupil(gp);
    solovev::Damping damping(gp);
    solovev::ZonalFlow zonalflow(gp,8.,0.5);
    solovev::Gradient gradient(gp,0.05);
//     //make dggrid
    dg::Grid2d<double> grid(Rmin,Rmax,Zmin,Zmax, n,Nx,Ny,dg::PER,dg::PER);
    
    dg::HVec hvisual1 = dg::evaluate( psip, grid);

    dg::blas2::gemv( equigrid, hvisual10, visual10);
    dg::blas2::gemv( equigrid, hvisual11, visual11);
    dg::blas2::gemv( equigrid, hvisual12, visual12);
    dg::blas2::gemv( equigrid, hvisual13, visual13);
    dg::blas2::gemv( equigrid, hvisual14, visual14);
    dg::blas2::gemv( equigrid, hvisual15, visual15);
    dg::blas2::gemv( equigrid, hvisual16, visual16);


    //Create Window and set window title
    GLFWwindow* w = draw::glfwInitAndCreateWindow( 1200, 1200, "");
    draw::RenderHostData render( 4, 4);
  
    //create a colormap
    draw::ColorMapRedBlueExtMinMax colors(-1.0, 1.0);
