Commit d07b18cf authored by Daniel Wortmann's avatar Daniel Wortmann

Fixed problem in vacuum part in case only one vacuum exists

parent 536795a2
......@@ -855,6 +855,14 @@ CONTAINS
v%pw(i,js)=vpw_w(i,js)/stars%nstr(i)
ENDDO
ENDDO
IF (vacuum%nvac==1) THEN
v%vacz(:,2,:)=v%vacz(:,1,:)
IF (sym%invs) THEN
v%vacxy(:,:,2,:)= cmplx(v%vacxy(:,:,1,:))
ELSE
v%vacxy(:,:,2,:)=v%vacxy(:,:,1,:)
ENDIF
ENDIF
CALL writePotential(stars,vacuum,atoms,cell,sphhar,input,sym,oneD,POT_ARCHIVE_TYPE_TOT_const,&
v%iter,v%mt,v%pw,v%vacz,v%vacxy)
......
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