     p(p),
     evec(3),
     //probe
-//     probevec(2),
-//     Xprobe(1,p.lx*p.posX), //use blob position
-//     Yprobe(1,p.ly*p.posY),//use blob position
-//     probeinterp(dg::create::interpolation( Xprobe,  Yprobe,g, dg::NEU)),
+    probe_coord_X(dg::create::linspace(0.0, p.lx, p.lx * 1. / 8.)),
+    probe_coord_Y(8, 0.5 * p.ly),
+    probeinterp(dg::create::interpolation(probe_coord_X, probe_coord_Y, g, dg::NEU)),
 //     probevalue(1,0.0),
     //damping and sources
     lh( dg::evaluate(dg::TanhProfX(p.lx*p.solb,p.solw,-1.0,0.0,1.0),g)),

    double mass_, energy_, diff_, ediff_,gammanex_,coupling_;
    std::vector<double> evec;
    //probe
//     std::vector<container> probevec;
//     const container Xprobe,Yprobe;
//     Matrix probeinterp;
//     container probevalue;    
//     dg::Grid1d<double> gy;
//     const container w1d;

    double coupling( ) {return coupling_;}

    std::vector<double> energy_vector( ) {return evec;}
    std::vector<container>& get_probe_vector( ) {return probevec;}
    void update_probes();

    double energy_diffusion( ){ return ediff_;}
    double radial_transport( ){ return gammanex_;}
