Merge pull request #846 from pfebrer/840_bugnc

fixed bug reading geometry from ncSileSiesta

Author: zerothi

src/sisl/io/siesta/siesta_nc.py

@@ -143,9 +143,7 @@ def read_geometry(self) -> Geometry:
         xyz.shape = (-1, 3)
 
         if "BASIS" in self.groups:
-            basis = self.read_basis()
-            species = self.groups["BASIS"].variables["basis"][:] - 1
-            atom = Atoms([basis[i] for i in species])
+            atom = self.read_basis()
         else:
             atom = Atom(1)