Merge pull request wannier-developers#337 from hjunlee/develop

fix the issue in plot_bvec
manxkim · Jul 17, 2020 · 461b96a · 461b96a
2 parents 359c628 + 9c00bc2
commit 461b96a
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 3 deletions.
diff --git a/src/plot.F90 b/src/plot.F90
@@ -1702,7 +1702,6 @@ subroutine plot_bvec()
     character(len=9)  :: cdate, ctime
     !
     file_unit = io_file_unit()
-    open (file_unit, file=trim(seedname)//'.bvec', form='formatted', status='unknown', err=101)
     call io_date(cdate, ctime)
     header = 'written on '//cdate//' at '//ctime
     !
@@ -1711,7 +1710,7 @@ subroutine plot_bvec()
     write (file_unit, *) num_kpts, nntot
     do nkp = 1, num_kpts
       do nn = 1, nntot
-        write (file_unit, '(4F12.6)') bk(:, nn, nkp), wb(nn)
+        write (file_unit, '(4F14.8)') bk(:, nn, nkp), wb(nn)
       enddo
     enddo
     close (file_unit)

diff --git a/src/wannier_lib.F90 b/src/wannier_lib.F90
@@ -387,7 +387,7 @@ subroutine wannier_run(seed__name, mp_grid_loc, num_kpts_loc, &
   time2 = io_time()
   write (stdout, '(1x,a25,f11.3,a)') 'Time for wannierise      ', time2 - time1, ' (sec)'
 
-  if (wannier_plot .or. bands_plot .or. fermi_surface_plot .or. write_hr) then
+  if (wannier_plot .or. bands_plot .or. fermi_surface_plot .or. write_hr .or. write_bvec) then
     call plot_main()
     time1 = io_time()
     write (stdout, '(1x,a25,f11.3,a)') 'Time for plotting        ', time1 - time2, ' (sec)'