diff --git a/nes/nc_projections/default_nes.py b/nes/nc_projections/default_nes.py
index 5e6269656a36f3bd1119ccc0c1fd98b5d43d0441..06c8ce76ce1709ab1e444adfe85703a83ebbeef2 100644
--- a/nes/nc_projections/default_nes.py
+++ b/nes/nc_projections/default_nes.py
@@ -976,8 +976,10 @@ class Nes(object):
         """
 
         if self.balanced:
+            print("balance")
             return self.get_read_axis_limits_balanced()
         else:
+            print("unbalance")
             return self.get_read_axis_limits_unbalanced()
 
     def get_read_axis_limits_unbalanced(self):
@@ -997,6 +999,7 @@ class Nes(object):
                        't_min': None, 't_max': None}
         
         idx = self.get_idx_intervals()
+        print("IDX", idx)
         
         if self.parallel_method == 'Y':
             y_len = idx['idx_y_max'] - idx['idx_y_min']
@@ -1005,7 +1008,7 @@ class Nes(object):
                     self.size, y_len))
             axis_limits['y_min'] = ((y_len // self.size) * self.rank) + idx['idx_y_min']
             if self.rank + 1 < self.size:
-                axis_limits['y_max'] = ((y_len // self.size) * (self.rank + 1)) + idx['idx_y_max']
+                axis_limits['y_max'] = ((y_len // self.size) * (self.rank + 1)) + idx['idx_y_min']
             else:
                 axis_limits['y_max'] = idx['idx_y_max']
 
@@ -1023,7 +1026,7 @@ class Nes(object):
                     self.size, x_len))
             axis_limits['x_min'] = ((x_len // self.size) * self.rank) + idx['idx_x_min']
             if self.rank + 1 < self.size:
-                axis_limits['x_max'] = ((x_len // self.size) * (self.rank + 1)) + idx['idx_x_max']
+                axis_limits['x_max'] = ((x_len // self.size) * (self.rank + 1)) + idx['idx_x_min']
             else:
                 axis_limits['x_max'] = idx['idx_x_max']
 
@@ -1041,7 +1044,7 @@ class Nes(object):
                     self.size, t_len))
             axis_limits['t_min'] = ((t_len // self.size) * self.rank) + idx['idx_t_min']
             if self.rank + 1 < self.size:
-                axis_limits['t_max'] = ((t_len // self.size) * (self.rank + 1)) + idx['idx_t_max']
+                axis_limits['t_max'] = ((t_len // self.size) * (self.rank + 1)) + idx['idx_t_min']
 
             # Non parallel filters
             axis_limits['y_min'] = idx['idx_y_min']
@@ -1063,7 +1066,7 @@ class Nes(object):
         else:
             self.last_level += 1
         axis_limits['z_max'] = self.last_level
-
+        print("Rank {0:03d}: Axis limits: {1}".format(self.rank, axis_limits))
         return axis_limits
 
     def get_read_axis_limits_balanced(self):