Add renormalization of wavefunction

[adam-srut]

The norm of the WF is during the run calculated via a simple dot product. While this is efficient, it would be more convenient to have properly normalized WF (with trapezoidal integration) as an output.