diff --git a/bench/dolo/dolo_inner_large.R b/bench/dolo/dolo_inner_large.R
index 73b1da839..802793581 100644
--- a/bench/dolo/dolo_inner_large.R
+++ b/bench/dolo/dolo_inner_large.R
@@ -8,7 +8,7 @@ grid_setup <- list(
     pqc_in_file = "./dol.pqi",
     pqc_db_file = "./phreeqc_kin.dat", # Path to the database file for Phreeqc
     grid_def = grid_def, # Definition of the grid, containing IDs according to the Phreeqc input script
-    grid_size = c(rows, cols) / 100, # Size of the grid in meters
+    grid_size = c(cols, rows) / 100, # Size of the grid in meters
     constant_cells = c() # IDs of cells with constant concentration
 )
 
diff --git a/bench/dolo/dolo_interp.R b/bench/dolo/dolo_interp.R
index bb88e457d..f48527d00 100644
--- a/bench/dolo/dolo_interp.R
+++ b/bench/dolo/dolo_interp.R
@@ -8,7 +8,7 @@ grid_setup <- list(
     pqc_in_file = "./dol.pqi",
     pqc_db_file = "./phreeqc_kin.dat", # Path to the database file for Phreeqc
     grid_def = grid_def, # Definition of the grid, containing IDs according to the Phreeqc input script
-    grid_size = c(5, 2.5), # Size of the grid in meters
+    grid_size = c(2.5, 5), # Size of the grid in meters
     constant_cells = c() # IDs of cells with constant concentration
 )
 
diff --git a/src/Init/GridInit.cpp b/src/Init/GridInit.cpp
index a35295769..4a3acaf65 100644
--- a/src/Init/GridInit.cpp
+++ b/src/Init/GridInit.cpp
@@ -135,8 +135,8 @@ void InitialList::prepareGrid(const Rcpp::List &grid_input) {
   this->n_rows = grid_def.nrow();
   this->n_cols = grid_def.ncol();
 
-  this->s_rows = grid_size[0];
-  this->s_cols = grid_size[1];
+  this->s_cols = grid_size[0];
+  this->s_rows = grid_size[1];
 
   this->dim = n_cols == 1 ? 1 : 2;