     std::vector<dg::x::DVec> transferD(4, dg::evaluate(dg::zero, grid_out));
     dg::x::HVec transferH(dg::evaluate(dg::zero, grid_out));
     dg::x::IDMatrix interpolate = dg::create::interpolation( grid_out, grid);
-    dg::blas2::symv( interpolate, y1[0], transferD[0]);
-    dg::blas2::symv( interpolate, y1[1], transferD[1]);
+    dg::blas2::symv( interpolate, y0[0], transferD[0]);
+    dg::blas2::symv( interpolate, y0[1], transferD[1]);
     dg::blas2::symv( interpolate, exp.potential()[0], transferD[2]);
     dg::blas2::symv( imp.laplacianM(), exp.potential()[0], transfer);
     dg::blas2::symv( interpolate, transfer, transferD[3]);

    ///////////////////////////////////////Timeloop/////////////////////////////////
    const double mass0 = exp.mass(), mass_blob0 = mass0 - grid.lx()*grid.ly();
    double E0 = exp.energy(), E1 = 0, diff = 0;
    dg::Timer t;
    t.tic();
    try

        }
        //////////////////////////write fields////////////////////////
        start = i;
        dg::blas2::symv( interpolate, y0[0], transferD[0]);
        dg::blas2::symv( interpolate, y0[1], transferD[1]);
        dg::blas2::symv( interpolate, exp.potential()[0], transferD[2]);
        dg::blas2::symv( imp.laplacianM(), exp.potential()[0], transfer);
        dg::blas2::symv( interpolate, transfer, transferD[3]);
