CarpetReduce: Handle groups that have not all time levels allocated

darcs-hash:20050303102622-891bb-8d89ed158189c376fa1235ece02a5d8fef87fcae.gz

1 files changed, 9 insertions, 0 deletions

diff --git a/Carpet/CarpetReduce/src/reduce.cc b/Carpet/CarpetReduce/src/reduce.cc
index 6eed2a0ab..eef15d3d1 100644
--- a/Carpet/CarpetReduce/src/reduce.cc
+++ b/Carpet/CarpetReduce/src/reduce.cc
@@ -1097,7 +1097,16 @@ namespace CarpetReduce {
             for (int tl=0; tl<num_tl; ++tl) {
               myinarrays.at(tl).resize (num_invars);
               for (int n=0; n<num_invars; ++n) {
+#if 0
                 myinarrays.at(tl).at(n) = CCTK_VarDataPtrI(cgh, tl, invars[n]);
+#else
+                int const vi = invars[n];
+                int const gi = CCTK_GroupIndexFromVarI (vi);
+                int const vi0 = CCTK_FirstVarIndexI (gi);
+                myinarrays.at(tl).at(n)
+                  = ((*arrdata.at(gi).at(Carpet::map).data.at(vi-vi0))
+                     (tl, reflevel, component, mglevel)->storage());
+#endif
                 assert (myinarrays.at(tl).at(n));
               }
               inarrays.at(tl) = &myinarrays.at(tl).at(0);