CCTK_REAL dx, dy, dz;

  vect<CCTK_REAL, dim> X;
  vect<CCTK_REAL, dim> DX;

       << "dxyz:"
       << "{" << p.dx << "," << p.dy << "," << p.dz << "}, "

      const PointDesc p{i, j, k, x, y, z, idx, dj, dk, I, inormal};

      const PointDesc p{i,     j,   k,  x,  y, z,       dx[0],     dx[1],
                        dx[2], idx, dj, dk, I, inormal, {x, y, z}, dx};