3D output doesn't care about staggering as it did in Cactus 3.x

where it did nothing if (pughGH->stagger != PUGH_VERTEXCTR)

Also use PUGH_NO_STAGGER instead of PUGH_VERTEXCTR for the highest
ownership[][][] index which should have the same meaning.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusPUGHIO/IOFlexIO/trunk@74 ebee0441-1374-4afa-a3b5-247f3ba15b9a

1 files changed, 2 insertions, 8 deletions

diff --git a/src/DumpVar.c b/src/DumpVar.c
index e642bdf..4bfa869 100644
--- a/src/DumpVar.c
+++ b/src/DumpVar.c
@@ -373,12 +373,6 @@ void IOFlexIO_DumpGF (cGH *GH, int index, int timelevel, IOFile iof,
   myproc = CCTK_MyProc (GH);
   nprocs = CCTK_nProcs (GH);
 
-  /* If I'm not vertex centered, bail (for now) */
-  if (pughGH->stagger != PUGH_VERTEXCTR) {
-    CCTK_WARN (1, "3D Output for staggered variables not available");
-    return;
-  }
-
   /* Get the pointer to the data we want to output (includes downsampling) */
   IOFlexIO_getDumpData (GH, index, timelevel, &outme, &free_outme, bnd,
                       element_size);
@@ -777,12 +771,12 @@ void IOFlexIO_getDumpData (cGH *GH, int index, int timelevel, void **outme,
       bnd [i]   = GH->cctk_lbnd[i] / downsample;
       start [i] = bnd [i] * downsample;
       if (start [i] <
-          GH->cctk_lbnd[i] + pughGH->ownership [PUGH_VERTEXCTR][i][0]) {
+          GH->cctk_lbnd[i] + pughGH->ownership [PUGH_NO_STAGGER][i][0]) {
         start [i] += downsample;
         bnd [i] ++;
       }
       end [i]   = ((GH->cctk_lbnd [i] +
-                   pughGH->ownership [PUGH_VERTEXCTR][i][1] - 1) / downsample)
+                   pughGH->ownership [PUGH_NO_STAGGER][i][1] - 1) / downsample)
                    * downsample;
       bnd [i+3] = (end [i] - start [i]) / downsample + 1;
     }