             dg::blas2::gemv( LaplacianM_perp, x[i], temp);
             dg::blas2::gemv( LaplacianM_perp, temp, y[i]);
             dg::blas1::axpby( -p.nu_perp, y[i], 0., y[i]); // - nu_perp lapl_RZ (lapl_RZ (lnN,U)) //factor MISSING!?!
+            
+            dg::blas1::pointwiseDot( lapiris_, temp, temp); //N_i U_i
+            dg::blas1::axpby( -0.1, temp, 1., y[i]); // - nu_perp lapl_RZ (lapl_RZ (lnN,U)) //factor MISSING!?!
+
 //             dg::blas2::gemv( LaplacianM_para, x[i], temp);
 //             dg::blas1::axpby(  -p.nu_parallel, temp, 1., y[i]);
         }

        w3d_( dg::create::w3d(g)), v3d_(dg::create::v3d(g)),
        temp( g.size()),
//         pupil_( dg::evaluate( solovev::Pupil( gp), g))
        pupil_( dg::evaluate( solovev::TanhDampingProf(gp ), g))
    {
        LaplacianM_perp = dg::create::laplacianM_perp( g, dg::normed, dg::symmetric);
//         LaplacianM_para = dg::create::laplacianM_parallel( g, dg::PER);

#include "dg/backend/timer.cuh"
#endif //DG_BENCHMARK

#define APAR
namespace eule
{
//diffusive terms (add mu_hat?)
