Carpet: Correct error in storage handling for grid arrays

When enabling or disabling storage for grid arrays, update the
pointers in the cGH structure as well.

darcs-hash:20051010210525-891bb-a01898931d0c5a1af8d1ad22905bfa559eb6f1c5.gz

1 files changed, 13 insertions, 0 deletions

diff --git a/Carpet/Carpet/src/Cycle.cc b/Carpet/Carpet/src/Cycle.cc
index 87fe7fa27..45425905d 100644
--- a/Carpet/Carpet/src/Cycle.cc
+++ b/Carpet/Carpet/src/Cycle.cc
@@ -39,10 +39,23 @@ namespace Carpet {
         case CCTK_SCALAR:
         case CCTK_ARRAY:
           if (do_global_mode) {
+	    int const numtimelevels = CCTK_NumTimeLevelsI (group);
+	    int const firstvarindex = CCTK_FirstVarIndexI (group);
             for (int var=0; var<CCTK_NumVarsInGroupI(group); ++var) {
               for (int c=0; c<arrdata.at(group).at(0).hh->components(0); ++c) {
                 arrdata.at(group).at(0).data.at(var)->cycle (0, c, mglevel);
               }
+	      {
+                int const varindex = firstvarindex + var;
+		int const c = CCTK_MyProc(cgh);
+		for (int tl=0; tl<numtimelevels; ++tl) {
+		  cgh->data[varindex][tl]
+		    = (tl < groupdata.at(group).info.activetimelevels
+		       ? ((*arrdata.at(group).at(0).data.at(var))
+			  (tl, 0, c, 0)->storage())
+		       : NULL);
+		}
+	      }
             }
           }
           break;