If requested by the checkpointing code pass the 'with_ghostzones' flag in the

parameter table to HyperslabDefineLocalMapping().
Note: you will need to update IOUtil also.


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

1 files changed, 17 insertions, 2 deletions

diff --git a/src/DumpVar.c b/src/DumpVar.c
index 7c9ea2e..2ce2995 100644
--- a/src/DumpVar.c
+++ b/src/DumpVar.c
@@ -13,6 +13,7 @@
 #include <string.h>
 
 #include "cctk.h"
+#include "util_Table.h"
 #include "cctk_Parameters.h"
 #include "CactusPUGH/PUGH/src/include/pugh.h"
 #include "CactusBase/IOUtil/src/ioGH.h"
@@ -322,7 +323,7 @@ static int WriteGS (const cGH *GH, const ioRequest *request, IOFile file)
 @@*/
 static int WriteGA (const cGH *GH, const ioRequest *request, IOFile file)
 {
-  int i, mapping, hdatasize, retval;
+  int i, mapping, hdatasize, table, retval;
   void *hdata;
   char *fullname;
   const ioGH *ioUtilGH;
@@ -338,16 +339,30 @@ static int WriteGA (const cGH *GH, const ioRequest *request, IOFile file)
 
 
   /* define the hyperslab mapping */
+  table = -1;
+  if (request->with_ghostzones)
+  {
+    table = Util_TableCreateFromString ("with_ghostzones = 1");
+    if (table < 0)
+    {
+      CCTK_WARN (1, "Failed to hyperslab parameter create table from string");
+    }
+  }
   mapping = Hyperslab_DefineLocalMappingByIndex (GH, request->vindex,
                                                  request->hdim,
                                                  request->direction,
                                                  request->origin,
                                                  request->extent,
                                                  request->downsample,
-                                                 -1, NULL,
+                                                 table, NULL,
                                                  request->hsize_chunk,
                                                  request->hsize,
                                                  request->hoffset);
+  if (table >= 0)
+  {
+    Util_TableDestroy (table);
+  }
+
   if (mapping < 0)
   {
     fullname = CCTK_FullName (request->vindex);