Fix bad iteration over z

mikekryjak · May 20, 2020 · f95d006 · f95d006
1 parent 74bb48d
commit f95d006
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/mesh/coordinates.cxx b/src/mesh/coordinates.cxx
@@ -668,7 +668,7 @@ interpolateAndNeumann(MAYBE_UNUSED(const Coordinates::metric_field_type& f), MAY
   // Set corner guard cells
   for (int i = 0; i < localmesh->xstart; i++) {
     for (int j = 0; j < localmesh->ystart; j++) {
-         for (int z=0;z<result.getNx();++z) {
+      for (int z=0;z<result.getNz();++z) {
         result(i, j, z) = BoutNaN;
         result(i, localmesh->LocalNy - 1 - j, z) = BoutNaN;
         result(localmesh->LocalNx - 1 - i, j, z) = BoutNaN;