 #include <fstream>
 #include <sstream>
 #include "file/read_input.h"
+#include <cmath>
 
 #include "draw/host_window.h"
 int main()
 {
-    unsigned Nx=100, Ny=100,polcoeff=3;
+    unsigned Nx,Ny,n;
 //     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";
+    std::cin >> n>> Nx>>Ny;
     
     std::vector<double> v;
     try{ v = file::read_input( "geometry_params.txt"); }

    dg::blas2::gemv( equigrid, hvisual4, visual4);
    dg::blas2::gemv( equigrid, hvisual5, visual5);
    dg::blas2::gemv( equigrid, hvisual6, visual6);

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

    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
