From d2576102d20466205d348d7a31b249fce949924a Mon Sep 17 00:00:00 2001
From: Junfeng Qiao <qiaojunfeng@outlook.com>
Date: Sat, 18 Jan 2025 18:27:55 +0100
Subject: [PATCH] Print phase factor for plotting MLWF

---
 src/plot.F90 | 3 +++
 1 file changed, 3 insertions(+)

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
           !