Dear Khani,
I did not check the units of the magnetic field, but I can tell you that as it
is, it has some effects on the density of states,
and you can even see the Landau levels.
The only tweak to do is to increase the energy resolution of
`kwant.kpm.SpectralDensity`. You can either modify the
argument `num_moments` or `energy_resolution`.
I tried with
```
spectrum = kwant.kpm.SpectralDensity(fsyst, num_vectors=None,
energy_resolution=0.01,
vector_factory=s_factory,
rng=0)
```
and seems to work fine.
Regards,
Pablo
