const int np = dk * (cctk_ash[2] + 1); // one extra grid point for vertex centering

#pragma omp single

                    cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2],
                    x, y, z,

                    np,
                    coordx, coordy, coordz,