Corrected to properly handle staggered grids.

This may have to be modified when the new coordinate infrastructure is
working properly.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/IOASCII/trunk@128 94b1c47f-dcfd-45ef-a468-0854c0e9e350

1 files changed, 8 insertions, 0 deletions

diff --git a/src/Write1D.c b/src/Write1D.c
index c7876ea..46cb6b9 100644
--- a/src/Write1D.c
+++ b/src/Write1D.c
@@ -317,6 +317,12 @@ fprintf (stderr, "%d/%d\n", origin[dir], extent);
                     "IOASCII_Write1D: failed to get physical coordinate ranges "
                     "for variable '%s'", fullname);
       }
+#else
+    /* correct extent in the case of staggered grids */
+    if (CCTK_StaggerDirIndex(dir,gdata.stagtype)==1)
+    {
+      --extent;
+    }
 #endif
     }
     else    /* origin for CCTK_ARRAYS is always (0, 0, 0) */
@@ -324,6 +330,8 @@ fprintf (stderr, "%d/%d\n", origin[dir], extent);
       memset (origin, 0, gdata.dim * sizeof (CCTK_INT));
     }
 
+
+
     mapping = Hyperslab_DefineGlobalMappingByIndex (GH, vindex, 1,
                                                     direction,
                                                     origin,