diff --git a/python/triqs_dft_tools/sumk_dft.py b/python/triqs_dft_tools/sumk_dft.py
index b1b9e036..d5fe200a 100644
--- a/python/triqs_dft_tools/sumk_dft.py
+++ b/python/triqs_dft_tools/sumk_dft.py
@@ -117,7 +117,7 @@ class SumkDFT(object):
                                                len(self.mesh))
             elif isinstance(mesh, MeshReFreq):
                 self.mesh = mesh
-                self.mesh_values = np.linspace(self.mesh.omega_min, self.mesh.omega_max, len(self.mesh))
+                self.mesh_values = np.linspace(self.mesh.w_min, self.mesh.w_max, len(self.mesh))
             else:
                 raise ValueError('mesh must be a triqs mesh of type MeshImFreq or MeshReFreq')
 
@@ -534,7 +534,7 @@ class SumkDFT(object):
             if mesh is None:
                 broadening = 0.01
             else:  # broadening = 2 * \Delta omega, where \Delta omega is the spacing of omega points
-                broadening = 2.0 * ((mesh.omega_max - mesh.omega_min) / (len(mesh) - 1))
+                broadening = 2.0 * ((mesh.w_max - mesh.w_min) / (len(mesh) - 1))
 
         # Check if G_latt is present
         set_up_G_latt = False                       # Assume not
@@ -564,7 +564,7 @@ class SumkDFT(object):
                 if isinstance(mesh, MeshImFreq):
                     mesh_values = np.linspace(mesh(mesh.first_index()), mesh(mesh.last_index()), len(mesh))
                 else:
-                    mesh_values = np.linspace(mesh.omega_min, mesh.omega_max, len(mesh))
+                    mesh_values = np.linspace(mesh.w_min, mesh.w_max, len(mesh))
             else:
                 mesh = self.mesh
                 mesh_values = self.mesh_values
@@ -574,7 +574,7 @@ class SumkDFT(object):
             if isinstance(mesh, MeshImFreq):
                 mesh_values = np.linspace(mesh(mesh.first_index()), mesh(mesh.last_index()), len(mesh))
             else:
-                mesh_values = np.linspace(mesh.omega_min, mesh.omega_max, len(mesh))
+                mesh_values = np.linspace(mesh.w_min, mesh.w_max, len(mesh))
         else:
             mesh = self.mesh
             mesh_values = self.mesh_values
diff --git a/python/triqs_dft_tools/sumk_dft_tools.py b/python/triqs_dft_tools/sumk_dft_tools.py
index 3360a7e7..c27fb275 100644
--- a/python/triqs_dft_tools/sumk_dft_tools.py
+++ b/python/triqs_dft_tools/sumk_dft_tools.py
@@ -116,7 +116,7 @@ class SumkDFTTools(SumkDFT):
                 broadening=0.001
         else:
             assert 0, "ReFreqMesh input required for calculations without real frequency self-energy"
-        mesh_val = numpy.linspace(mesh.omega_min,mesh.omega_max,len(mesh))
+        mesh_val = numpy.linspace(mesh.w_min,mesh.w_max,len(mesh))
         n_om = len(mesh)
         om_minplot = mesh_val[0] - 0.001
         om_maxplot = mesh_val[-1] + 0.001
@@ -494,7 +494,7 @@ class SumkDFTTools(SumkDFT):
                 broadening=0.001
         else:
             assert 0, "ReFreqMesh input required for calculations without real frequency self-energy"
-        mesh_val = numpy.linspace(mesh.omega_min,mesh.omega_max,len(mesh))
+        mesh_val = numpy.linspace(mesh.w_min,mesh.w_max,len(mesh))
         n_om = len(mesh)
         om_minplot = mesh_val[0] - 0.001
         om_maxplot = mesh_val[-1] + 0.001
@@ -659,7 +659,7 @@ class SumkDFTTools(SumkDFT):
                 broadening=0.001
         else:
             assert 0, "ReFreqMesh input required for calculations without real frequency self-energy"
-        mesh_val = numpy.linspace(mesh.omega_min,mesh.omega_max,len(mesh))
+        mesh_val = numpy.linspace(mesh.w_min,mesh.w_max,len(mesh))
         n_om = len(mesh)
         om_minplot = mesh_val[0] - 0.001
         om_maxplot = mesh_val[-1] + 0.001
@@ -765,7 +765,7 @@ class SumkDFTTools(SumkDFT):
                        resolved into orbital contributions. Empty if proj_type = None
         """
 
-        mesh_val = numpy.linspace(mesh.omega_min,mesh.omega_max,len(mesh))
+        mesh_val = numpy.linspace(mesh.w_min,mesh.w_max,len(mesh))
         n_om = len(mesh)
         om_minplot = mesh_val[0] - 0.001
         om_maxplot = mesh_val[-1] + 0.001