Skip to content
Snippets Groups Projects
Commit 4b1c7666 authored by Philipp Rüssmann's avatar Philipp Rüssmann
Browse files

Fix VCONST+IVSHIFT mode where constant potential is only added to one atom

parent e6f5b4a1
No related branches found
No related tags found
No related merge requests found
......@@ -496,7 +496,7 @@ contains
do j = 1, irmin(ih)
vm2z(j, i) = vm2z(j, i) + vconst
end do
else if (abs(vconst)>eps) then
else if (abs(vconst)>eps .and. ivshift==0) then
write (1337, *) 'shifting potential by VCONST=', vconst
do j = 1, nr
vm2z(j, i) = vm2z(j, i) + vconst
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment