The routine handling the center of charges has been more generalized. Now, at the end it is printed a coupled cluster calculation and the substraction needed to have the incremental energy
README.md improved