Now fleur.F90 calls a subroutine that does most of the hybrid stuff.
In eigen.F90 only the non-local potential is added to the Hamiltonian