     heavi.set_origin( p.lx*p.posX, p.ly*p.posY);
     dg::HVec heavy = dg::evaluate( heavi, grid);
     double normalize = dg::blas2::dot( heavy, w2d, input0);
-    double gamma = sqrt( p.kappa*p.n0*(1+p.tau)/p.sigma);
+    double gamma = sqrt( p.kappa*p.n0/(1+p.n0)*(1+p.tau)/p.sigma);//global gamma!!
 
     unsigned idx = (unsigned)(T/gamma/p.dt/p.itstp) + 1;
     {

        posY_max = hy*(1./2. + (double)(position/Nx));
        //init Heaviside with this
        heavi.set_origin( posX_max, posY_max);
        dg::HVec heavy = dg::evaluate( heavi, grid);
        std::cout << sqrt(p.tau)/p.sigma <<" ";
        std::cout << std::fixed<< t5file.get_time( idx)*gamma << " ";//(1)
        //std::cout << " "<<posX_max<<" "<<posY_max <<" ";
        std::cout <<std::scientific<< dg::blas2::dot( heavy, w2d, input0)/normalize << "\n";


#include "galerkin/parameters.h"

const double T = 10;
//is sigma the radius or the diameter

double X( double x, double y) {return x;}
