             
             //  div(K_kappa) terms
             dg::blas1::pointwiseDot(u[i],divCurvKappa,omega);              // U div(K_kappa)
-            dg::blas1::axpby( -p.tau[i], omega, 1., yp[2+i]);                //dtU = dtU -tau U div(K_kappa)
+            dg::blas1::axpby( -p.tau[i], omega, 1., yp[2+i]);                //dtw = dtw -tau U div(K_kappa)
             dg::blas1::pointwiseDot(u[i],omega,omega);                     // U^2 div(K_kappa)
             dg::blas1::pointwiseDot(npe[i],omega,omega);                     // N U^2 div(K_kappa)
             dg::blas1::axpby( -p.mu[i], omega, 1., yp[i]);                //dtw = dtw -hat(mu) N U^2 div(K_kappa)

{   
    /*  y[0] := N_e - 1
        y[1] := N_i - 1
        y[2] := w_e =  U_e + beta/mu_e A_parallel
        y[3] := w_i =  U_i + beta/mu_i A_parallel
    */
    
    dg::Timer t;


    container chi, omega,lambda;//1d container

    container binv,  curvX, curvY, curvKappaX, curvKappaY,divCurvKappa, dslnB;
    container source, damping, one;
    container profne, profNi;
    container w3d, v3d;
    
    std::vector<container> phi,apar,dsblnB, u, u2, npe, logn,un,curvphi,curvkappaphi, dsphi, dslogn,dsun,dsu2; //2d container
    std::vector<container> poissonn,poissonu,poissonw,poissonun,poissonlogn,poissonphi,poissonu2; //2d container
    std::vector<container> dsy, curvy,curvkappay;  //4d container

