-- read the output of FPLO
print("--Read FPLO output--\n")
FPLOOut = FileReadDresdenFPLO("DFT/out.wan")


-- from the DFT output we can create a tight binding Hamiltonian
print("--Create the Tight Binding Hamiltonian--\n")
print("Create the tight binding Hamiltonian for the crystal\n")
TB = TightBindingDefFromDresdenFPLO(FPLOOut)
print("Basis")
print(TB.Atoms)
print("Cell")
print(TB.Cell)

print("Define a cluster")

NewCluster = FindAllAtomsInsideSphere(TB.Atoms,TB.Cell,{0,0,0},5)
print("Quanty detected the following cluster:")
print("Cluster")
print(NewCluster)