[Kwant] superconducting leads

Hello, I am trying to how understand how to compute the wavefunctions of a system with superconducting leads as a function of the phase difference, just as in the case of normal leads, using fsys.wave_function() where I would do it as a function of the energy. However, the result is always an empty array. Is there a way to make it work ? Thanks in advance, Anil