     dg::HVec visual5( grid.size());
     dg::HVec visual6( grid.size());
     dg::HVec visual7( grid.size());
-    dg::HVec visual8( grid.size());
+    dg::HVec visual8 (grid.size());
     dg::HVec visual9( 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

        
        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());
        title.str("");

#include <fstream>
#include <sstream>
#include "file/read_input.h"

#include "draw/host_window.h"
int main()
{
//     unsigned Nx=100, Ny=100,polcoeff=3;
//     double Nxh = Nx/2.,Nyh=Ny/2.;
//     double a, elongation,triangularity;
    double Rmin,Zmin,Rmax,Zmax;
//     double A,R_0,psipmin,psipmax;
//     std::vector<double> c(13);
     std::cout << "Type n, Nx, Ny\n";
    unsigned n, Nx, Ny;
    std::cin >> n>> Nx>>Ny;   
    std::vector<double> v;
    try{ v = file::read_input( "geometry_params.txt"); }
    catch (toefl::Message& m) {  
