diff options
Diffstat (limited to 'src/Write1D.c')
| -rw-r--r-- | src/Write1D.c | 18 |
1 files changed, 3 insertions, 15 deletions
diff --git a/src/Write1D.c b/src/Write1D.c index 7779c86..4134e06 100644 --- a/src/Write1D.c +++ b/src/Write1D.c @@ -174,8 +174,8 @@ int IOASCII_Write1D (const cGH *GH, int vindex, const char *alias) do_dir[0] = out1D_x && gdata.dim >= 1; do_dir[1] = out1D_y && gdata.dim >= 2; do_dir[2] = out1D_z && gdata.dim >= 3; - /* diagonal slice is done only if variable is non-staggered and 3D */ - do_dir[3] = out1D_d && gdata.dim == 3 && gdata.stagtype == 0; + /* diagonal slice is done only if variable is 3D */ + do_dir[3] = out1D_d && gdata.dim == 3; if (out1D_d && ! do_dir[3] && myGH->out1D_last[vindex] < 0) { CCTK_VWarn (3, __LINE__, __FILE__, CCTK_THORNSTRING, @@ -327,12 +327,6 @@ int IOASCII_Write1D (const cGH *GH, int vindex, const char *alias) origin[i] = myGH->spxyz[gdata.dim-1][dir][i]; } extent -= origin[dir]; - - /* correct extent in the case of staggered grids */ - if (CCTK_StaggerDirIndex(dir,gdata.stagtype)==1) - { - --extent; - } } else /* origin for CCTK_ARRAYS is always (0, 0, 0) */ { @@ -394,13 +388,7 @@ int IOASCII_Write1D (const cGH *GH, int vindex, const char *alias) { if (dir < 3) { - /* get the staggering offset for the xyz coordinates */ - offset = 0.5 * GH->cctk_delta_space[dir] * - CCTK_StaggerDirIndex (dir, gdata.stagtype); - for (i = 0; i < hsize; i++) - { - ((CCTK_REAL *) hdata[1])[i] += offset; - } + /* do nothing */ } else { |