WRITES: lambdaU0GF(interior) lambdaU1GF(interior) lambdaU2GF(interior)

  # I don't really write to the outer boundary but for periodic boundaries the 
  # SYNC above took care of them (and CarpetX does not yet realize this)
  WRITES: lambdaU0GF(interiorwithboundary) lambdaU1GF(interiorwithboundary) lambdaU2GF(interiorwithboundary)

  # I do not really write to the boundary, but for periodic ones used, the SYNC does
  WRITES: BaikalX::evol_variables(boundary)

});
Loop::loop_bnd<0,0,0>(cctkGH, [&](const Loop::PointDesc &p){
    const int i0 = p.i;
    const int i1 = p.j;
    const int i2 = p.k;

    lambdaU0GF_(i0,i1,i2) = 0.;
    lambdaU1GF_(i0,i1,i2) = 0.;
    lambdaU2GF_(i0,i1,i2) = 0.;
    

Loop::GF3D<CCTK_REAL, 0, 0, 0> aDD00GF_(cctkGH, aDD00GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> aDD01GF_(cctkGH, aDD01GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> aDD02GF_(cctkGH, aDD02GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> aDD11GF_(cctkGH, aDD11GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> aDD12GF_(cctkGH, aDD12GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> aDD22GF_(cctkGH, aDD22GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> alphaGF_(cctkGH, alphaGF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> betU0GF_(cctkGH, betU0GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> betU1GF_(cctkGH, betU1GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> betU2GF_(cctkGH, betU2GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> cfGF_(cctkGH, cfGF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> hDD00GF_(cctkGH, hDD00GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> hDD01GF_(cctkGH, hDD01GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> hDD02GF_(cctkGH, hDD02GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> hDD11GF_(cctkGH, hDD11GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> hDD12GF_(cctkGH, hDD12GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> hDD22GF_(cctkGH, hDD22GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> lambdaU0GF_(cctkGH, lambdaU0GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> lambdaU1GF_(cctkGH, lambdaU1GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> lambdaU2GF_(cctkGH, lambdaU2GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> trKGF_(cctkGH, trKGF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> vetU0GF_(cctkGH, vetU0GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> vetU1GF_(cctkGH, vetU1GF);
Loop::GF3D<CCTK_REAL, 0, 0, 0> vetU2GF_(cctkGH, vetU2GF);

});
Loop::loop_bnd<0,0,0>(cctkGH, [&](const Loop::PointDesc &p){
    const int i0 = p.i;
    const int i1 = p.j;
    const int i2 = p.k;
aDD00GF_(i0, i1, i2) = 0.0;
aDD01GF_(i0, i1, i2) = 0.0;
aDD02GF_(i0, i1, i2) = 0.0;
aDD11GF_(i0, i1, i2) = 0.0;
aDD12GF_(i0, i1, i2) = 0.0;
aDD22GF_(i0, i1, i2) = 0.0;
alphaGF_(i0, i1, i2) = 0.0;
betU0GF_(i0, i1, i2) = 0.0;
betU1GF_(i0, i1, i2) = 0.0;
betU2GF_(i0, i1, i2) = 0.0;
cfGF_(i0, i1, i2) = 0.0;
hDD00GF_(i0, i1, i2) = 0.0;
hDD01GF_(i0, i1, i2) = 0.0;
hDD02GF_(i0, i1, i2) = 0.0;
hDD11GF_(i0, i1, i2) = 0.0;
hDD12GF_(i0, i1, i2) = 0.0;
hDD22GF_(i0, i1, i2) = 0.0;
lambdaU0GF_(i0, i1, i2) = 0.0;
lambdaU1GF_(i0, i1, i2) = 0.0;
lambdaU2GF_(i0, i1, i2) = 0.0;
trKGF_(i0, i1, i2) = 0.0;
vetU0GF_(i0, i1, i2) = 0.0;
vetU1GF_(i0, i1, i2) = 0.0;
vetU2GF_(i0, i1, i2) = 0.0;