 #include <fstream>
 #include <sstream>
 #include "file/read_input.h"
-#include <cmath>
 
 #include "draw/host_window.h"
 int main()
 {
-    unsigned Nx,Ny,n;
+//     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";
-    std::cin >> n>> Nx>>Ny;
-    
+     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) {  

    dg::Iris iris(gp);
    dg::Pupil pupil(gp);
    dg::Damping damping(gp);
    dg::Gradient grad(gp,0.05);
    dg::ZonalFlow zonalflow(gp,7);
//     //make dggrid
    dg::Grid2d<double> grid(Rmin,Rmax,Zmin,Zmax, n,Nx,Ny,dg::PER,dg::PER);
    

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