We constructed Wannier functions to represent the lower energy bands
as good as possible. There was a bit of experimenting involved to find
the best energy window. Since the bands are all entangled we need to
take a real space cutoff of 35 Bohr radii. This leads to replica at
large distances. This is seen in the +WFstat... files which show the
orbital contribution weights of each Wannier function (WF). If it
falls of at larger distance and starts to grow again for very large
distances you see a replica. This is an artefact of the finite
k-integration. Consequently, we used ``nkx=24`` instead of the
default ``12``. This removes the replicas.
The allover quality of the tight-binding WF representation is not
perfect but as  good as it gets.

The :file:`=.wandef` file was created via `makewandef.py`. This has the
advantage that the Wannier function names (labels) have a defined convenient
format.

To complete the example first run ``fplo15.00-50-x86_64`` where you
insert the correct fplo version (>14).
This should produce :file:`+hamdata`, :file:`+wanband` and :file:`+wanbandtb`

The resulting wannier bands and tight binding model bands are compared
via ``xfplo wband.cmd``
