for (int d = 0; d < dim; ++d)                               \

for (int d = 0; d < LC_DIM; ++d)                            \

  for (int d = 0; d < dim; ++d) {                                               \

  const int offset[LC_DIM] = {1, 1, 1};                                         \
  for (int d = 0; d < LC_DIM; ++d) {                                            \

    imax[d] = grid.tmax[d] + (grid.tmax[d] >= grid.lsh[d] ? grid.offset[d] : 0);\
    if(imax[d] > grid.lsh[d] + grid.offset[d]) {                                \
      imax[d] = grid.lsh[d] + grid.offset[d];                                   \

    imax[d] = grid.tmax[d] + (grid.tmax[d] >= grid.lsh[d] ? offset[d] : 0);     \
    if(imax[d] > grid.lsh[d] + offset[d]) {                                     \
      imax[d] = grid.lsh[d] + offset[d];                                        \

    inormal[d] = 0;                                                             \

  int inormal[3] = {0, 0, 0};                                                   \

}

} while(0)