Commit 8010cf4a authored by Daniel Wortmann's avatar Daniel Wortmann

Added deletion of off-diagonal part of potential in spin-spiral force-theorem calculation

parent 1ead9491
......@@ -50,12 +50,17 @@ CONTAINS
!Average out magnetic part of potential/charge in INT+Vacuum
potden%pw(:,1)=(potden%pw(:,1)+potden%pw(:,2))/2.0
potden%pw(:,2)=potden%pw(:,1)
IF (SIZE(potden%pw,2)==3) potden%pw(:,3)=0.0
potden%vacz(:,:,1)=(potden%vacz(:,:,1)+potden%vacz(:,:,2))/2.0
potden%vacxy(:,:,:,1)=(potden%vacxy(:,:,:,1)+potden%vacxy(:,:,:,2))/2.0
potden%vacz(:,:,2)=potden%vacz(:,:,1)
potden%vacxy(:,:,:,2)=potden%vacxy(:,:,:,1)
!Off diagonal part
IF (SIZE(potden%pw,2)==3) THEN
potden%pw(:,3)=0.0
potden%vacz(:,:,3:)=0.0
potden%vacxy(:,:,:,3)=0.0
END IF
END SUBROUTINE ssdisp_start
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment