diff --git a/src/plot.F90 b/src/plot.F90
index b7e4b32f..cfcbb565 100644
--- a/src/plot.F90
+++ b/src/plot.F90
@@ -1938,7 +1938,10 @@ subroutine plot_wannier(wannier_plot, wvfn_read, wannier_data, print_output, u_m
             end do
             wmod = wmod/sqrt(real(wmod)**2 + aimag(wmod)**2)
             wann_func(:, :, :, loop_w) = wann_func(:, :, :, loop_w)/wmod
+            write (stdout, '(6x,a,i4,7x,a,f11.6,SP,f11.6,"i")') "Wannier Function Num: ", &
+              wannier_plot%list(loop_w), "Phase Factor = ", 1/wmod
           end do
+          write (stdout, *) ''
           !
           ! Check the 'reality' of the WF
           !