From 47784d4cc44ee0421a916f7dcb0bb3a1c11043f2 Mon Sep 17 00:00:00 2001
From: QuanSheng Wu <wuquansheng@gmail.com>
Date: Thu, 12 Apr 2018 12:20:50 +0200
Subject: [PATCH] fixed a bug in readinput.f90

---
 soc/readinput.f90 | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/soc/readinput.f90 b/soc/readinput.f90
index ae056f83..564d475d 100644
--- a/soc/readinput.f90
+++ b/soc/readinput.f90
@@ -552,6 +552,8 @@ subroutine readinput
         write(stdout, '(a,2f10.4)')' cell_volume vs cell_volume-new', cell_volume, cell_volume2
         write(stdout, '(a)')" However, don't worry, WannierTools will help you to find a suitable rotation matrix."
         write(stdout, '(a)')" I am looking for new unit cell atuomatically: "
+     endif     !> print out the new basis
+     if (abs(abs(cell_volume2)-abs(cell_volume))> 0.001d0) then
         call FindTheThirdLatticeVector()
         R1= Umatrix(1, 1)*Rua+ Umatrix(1, 2)*Rub+ Umatrix(1, 3)*Ruc
         R2= Umatrix(2, 1)*Rua+ Umatrix(2, 2)*Rub+ Umatrix(2, 3)*Ruc