     dg::HVec visual4( grid.size());
     dg::HVec visual5( grid.size());
     dg::HVec visual6( grid.size());
+    dg::HVec visual7( grid.size());
+    dg::HVec visual8( grid.size());
+    dg::HVec visual9( grid.size());
+    dg::HVec visual10( grid.size());
+    dg::HVec visual11( grid.size());
     //make equidistant grid from dggrid
     dg::HMatrix equigrid = dg::create::backscatter(grid);
     //evaluate on valzues from devicevector on equidistant visual hvisual vector

        title <<"1/B / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual3, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);

        colors.scalemax() = (float)thrust::reduce( visual4.begin(), visual4.end(), 0., dg::AbsMax<double>()  );
        colors.scalemin() =  (float)thrust::reduce( visual4.begin(), visual4.end(), colors.scalemax() ,thrust::minimum<double>() );
        title <<"K^R / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual4, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        
        colors.scalemax() = (float)thrust::reduce( visual5.begin(), visual5.end(), 0., dg::AbsMax<double>()  );
        colors.scalemin() =  (float)thrust::reduce( visual5.begin(), visual5.end(), colors.scalemax() ,thrust::minimum<double>() );
        title <<"K^Z / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual5, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        

        colors.scalemin() =  (float)thrust::reduce( visual6.begin(), visual6.end(), colors.scalemax() ,thrust::minimum<double>() );
        title <<"gradLnB / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual6, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        
        colors.scalemax() = (float)thrust::reduce( visual7.begin(), visual7.end(), 0., dg::AbsMax<double>()  );
        colors.scalemin() =  (float)thrust::reduce( visual7.begin(), visual7.end(), colors.scalemax() ,thrust::minimum<double>() );
        title <<"iris / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual7, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        
        colors.scalemax() = (float)thrust::reduce( visual8.begin(), visual8.end(), 0., dg::AbsMax<double>()  );
        colors.scalemin() =  (float)thrust::reduce( visual8.begin(), visual8.end(), colors.scalemax() ,thrust::minimum<double>() );
        title <<"pupil / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual8, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        
        colors.scalemax() = (float)thrust::reduce( visual9.begin(), visual9.end(), 0., dg::AbsMax<double>()  );
        colors.scalemin() =  (float)thrust::reduce( visual9.begin(), visual9.end(), colors.scalemax() ,thrust::minimum<double>() );
        title <<"damping / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual9, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        
        colors.scalemax() = (float)thrust::reduce( visual10.begin(), visual10.end(), 0., dg::AbsMax<double>()  );
        colors.scalemin() =  (float)thrust::reduce( visual10.begin(), visual10.end(), colors.scalemax() ,thrust::minimum<double>() );
        title <<"gradient / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual10, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        
        colors.scalemax() = (float)thrust::reduce( visual11.begin(), visual11.end(), 0., dg::AbsMax<double>()  );
        colors.scalemin() =  (float)thrust::reduce( visual11.begin(), visual11.end(), colors.scalemax() ,thrust::minimum<double>() );
        title <<"gradient / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual11, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        colors.scalemax() = (float)thrust::reduce( visual11.begin(), visual11.end(), 0., dg::AbsMax<double>()  );
        colors.scalemin() =  (float)thrust::reduce( visual11.begin(), visual11.end(), colors.scalemax() ,thrust::minimum<double>() );
        title <<"gradient / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual11, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        
        title << std::fixed; 
        glfwSetWindowTitle(w,title.str().c_str());
