diff --git a/harmonicmaps-eoc.cc b/harmonicmaps-eoc.cc
index 269e88615d1c0d6eadbfd6f5922cda1f69a7bbc0..1db6c2c996d894291839d105beef083818a6675a 100644
--- a/harmonicmaps-eoc.cc
+++ b/harmonicmaps-eoc.cc
@@ -174,8 +174,10 @@ int main (int argc, char *argv[]) try
 
     array<unsigned int,dim> elements;
     elements.fill(numBaseElements);
-    shared_ptr<GridType> referenceGrid = StructuredGridFactory<GridType>::createSimplexGrid(FieldVector<double,dim>(0),
-                                                                                            FieldVector<double,dim>(1),
+    FieldVector<double,dim> lowerLeft  = parameterSet.get<FieldVector<double,dim> >("lowerLeft");
+    FieldVector<double,dim> upperRight = parameterSet.get<FieldVector<double,dim> >("upperRight");
+    shared_ptr<GridType> referenceGrid = StructuredGridFactory<GridType>::createSimplexGrid(lowerLeft,
+                                                                                            upperRight,
                                                                                             elements);
     referenceGrid->globalRefine(numLevels-1);
 
@@ -208,8 +210,8 @@ int main (int argc, char *argv[]) try
 
         array<unsigned int,dim> elements;
         elements.fill(numBaseElements);
-        shared_ptr<GridType> grid = StructuredGridFactory<GridType>::createSimplexGrid(FieldVector<double,dim>(0),
-                                                                                       FieldVector<double,dim>(1),
+        shared_ptr<GridType> grid = StructuredGridFactory<GridType>::createSimplexGrid(lowerLeft,
+                                                                                       upperRight,
                                                                                        elements);
 
         grid->globalRefine(i-1);
diff --git a/harmonicmaps-eoc.parset b/harmonicmaps-eoc.parset
index ce10259d114217dc5b3c176b83f63099b8d008fe..051081d76559aec8f6306438901be76bade22c13 100644
--- a/harmonicmaps-eoc.parset
+++ b/harmonicmaps-eoc.parset
@@ -29,3 +29,7 @@ baseTolerance = -1
 # Number of elements of the rod grid
 numBaseElements = 4
 
+# Bounding box
+lowerLeft = -2 -1
+upperRight = 2 1
+