USES INCLUDE HEADER: loopcontrol.h

  int imin[3], imax[3];

  for (int d = 0; d < 3; ++ d)
  {
    /*
    imin[d] = 0           + (cctk_bbox[2*d  ] ? 0 : cctk_nghostzones[d]);
    imax[d] = cctk_lsh[d] - (cctk_bbox[2*d+1] ? 0 : cctk_nghostzones[d]);
    */
    imin[d] = 0;
    imax[d] = cctk_lsh[d];
  }

#pragma omp parallel for
  for (int k = imin[2]; k < imax[2]; ++k)
  {
    for (int j = imin[1]; j < imax[1]; ++j)
    {
      for (int i = imin[0]; i < imax[0]; ++i)

  CCTK_LOOP3_ALL(TwoPunctures, cctkGH, i,j,k)

      } /* for i */
    }   /* for j */
  }     /* for k */

      } CCTK_ENDLOOP3_ALL(TwoPunctures);