Indices to the one body interactions have been added (in the case that π and σ orbitals are split in the localization). Now two body interaction input file is almost ready for production.