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

    dg::HVec hvisual4 = dg::evaluate( curvatureR, grid);
    dg::HVec hvisual5 = dg::evaluate( curvatureZ, grid);
    dg::HVec hvisual6 = dg::evaluate( gradLnB, grid);
    //allocate mem for visual
    dg::HVec visual1( grid.size());
    dg::HVec visual2( grid.size());

        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(),  -1.0, thrust::maximum<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(), -1.0, thrust::maximum<double>()   );
        colors.scalemin() =  (float)thrust::reduce( visual5.begin(), visual5.end(), colors.scalemax() ,thrust::minimum<double>() );
//                 colors.scalemax() = -0.0058;
//         colors.scalemin() = -0.009;
        title <<"K^Z / "<<colors.scalemin()<<"  " << colors.scalemax()<<"\t";
        render.renderQuad( visual5, grid.n()*grid.Nx(), grid.n()*grid.Ny(), colors);
        
