Adding routines to deal with atomic orbitals. For the moment this is not yet fully used.