diff --git a/python/converters/vasp_converter.py b/python/converters/vasp_converter.py
index 01b38a19..d2a04ecc 100644
--- a/python/converters/vasp_converter.py
+++ b/python/converters/vasp_converter.py
@@ -162,7 +162,9 @@ class VaspConverter(ConverterTools):
                     for i, ion in enumerate(ion_list):
                         pars = {}
                         pars['atom'] = ion
-                        pars['sort'] = sh['ion_sort']
+# We set all sites inequivalent
+#                        pars['sort'] = sh['ion_sort']
+                        pars['sort'] = ion
                         pars['l'] = sh['lorb']
                         pars['dim'] = sh['ndim']
                         pars['SO'] = SO
@@ -180,6 +182,9 @@ class VaspConverter(ConverterTools):
 #        to define equivalence classes of sites.
             n_inequiv_shells, corr_to_inequiv, inequiv_to_corr = ConverterTools.det_shell_equivalence(self, corr_shells)
 
+            if mpi.is_master_node():
+                print "  No. of inequivalent shells:", n_inequiv_shells
+
 # NB!: these rotation matrices are specific to Wien2K! Set to identity in VASP
             use_rotations = 1
             rot_mat = [numpy.identity(corr_shells[icrsh]['dim'],numpy.complex_) for icrsh in range(n_corr_shells)]