Looks like norm of the wavefunction is conserved!
Tully I/O appears to be working beautifully.