CarpetRegrid: Take basemglevel into account when checking the grid spacings

Take basemglevel into account when checking the user specified grid
spacings.

darcs-hash:20060619205829-dae7b-82d163027ba08f4a0c9a4f3487e4f112c33fa270.gz

1 files changed, 3 insertions, 8 deletions

diff --git a/Carpet/CarpetRegrid/src/manualcoordinatelist.cc b/Carpet/CarpetRegrid/src/manualcoordinatelist.cc
index 7fe29fda0..bd189db0a 100644
--- a/Carpet/CarpetRegrid/src/manualcoordinatelist.cc
+++ b/Carpet/CarpetRegrid/src/manualcoordinatelist.cc
@@ -169,14 +169,9 @@ namespace CarpetRegrid {
         // TODO:
         // assert (domain_from_coordbase);
         ivect const spacereffact = spacereffacts.at(rl);
-        // Do not use basemglevel here, since this spacing is not use
-        // for Carpet (internally), but only used to check the spacing
-        // specified by the user in the parameter file.  This spacing
-        // should be independent of basemglevel (i.e.,
-        // convergence_level), so that the user does not have to
-        // change the parameter file when the convergence level is
-        // changed.
-        rvect const spacing = base_spacing / rvect(spacereffact);
+        rvect const spacing =
+          base_spacing * ipow((CCTK_REAL)mgfact, basemglevel)
+          / rvect(spacereffact);
         if (! all(abs(ext.stride() - spacing)
                   < spacing * (CCTK_REAL) 1.0e-10)) {
           assert (dim==3);