CarpetRegrid: allow refined regions to move

darcs-hash:20041008214349-891bb-c7bfec9052c47d14c70ed7c9c620e97be1b959f6.gz

1 files changed, 10 insertions, 1 deletions

diff --git a/Carpet/CarpetRegrid/src/manualcoordinatelist.cc b/Carpet/CarpetRegrid/src/manualcoordinatelist.cc
index 361c979a3..1a2a4b871 100644
--- a/Carpet/CarpetRegrid/src/manualcoordinatelist.cc
+++ b/Carpet/CarpetRegrid/src/manualcoordinatelist.cc
@@ -162,9 +162,18 @@ namespace CarpetRegrid {
         // rvect const spacing = base_spacing * pow(CCTK_REAL(mgfact), basemglevel) / ipow(reffact, rl);
         rvect const spacing = base_spacing / ipow(reffact, rl);
         assert (all(abs(ext.stride() - spacing) < spacing * 1.0e-10));
+        
+        rvect offset = rvect(0);
+        if (c < num_offsets) {
+          if (rl >= offset_firstlevel) {
+            assert (dim==3);
+            offset = rvect(offsetx[c], offsety[c], offsetz[c]);
+          }
+        }
+        
         ManualCoordinates_OneLevel
           (cctkGH, hh, rl, refinement_levels,
-           ext.lower(), ext.upper(), ob, bbs, obs);
+           ext.lower() + offset, ext.upper() + offset, ob, bbs, obs);
       }
       
       // make multiprocessor aware