Indeed, it does look like only the real parts are output. You can also output the imaginary parts by modifying the code near lines 968 in core/Hdf5Interface.f90. The easiest is probably to just output the imaginary parts under a different name, and then put the real and imaginary parts back together in postprocessing.
It seems like the TEBD code produces a symmetrized <adag a> without the imaginary parts? Is there a way I can get the imaginary parts? I want to calculate transient currents and momentum distributions.