Commit 703b8e51 authored by Daniel Wortmann's avatar Daniel Wortmann

Bugfix for last commit (c_1 was defined twice), Changed initialization of ccof

parent b1509ec4
......@@ -94,6 +94,10 @@ MODULE m_setabc1locdn
!---> set up enough, nbasf0 and initialize nkvec
natom = 0
nbasf = lapw%nv(jspin)
!---> initialize ccof
IF (iintsp.NE.2) THEN
ccof(:,:,:,:)=cmplx(0.,0.)
ENDIF
if (noco%l_ss) nbasf=lapw%nv(iintsp)
DO ntyp = 1,atoms%ntype
DO nn = 1,atoms%neq(ntyp)
......@@ -111,10 +115,6 @@ MODULE m_setabc1locdn
nbasf0(lo,sym%invsatnr(natom)) = nbasf
nbasf = nbasf + 2* (2*l+1)
END IF
!---> initialize ccof
IF (iintsp.NE.2) THEN
ccof(:,:,lo,natom)=cmplx(0.,0.)
ENDIF
END DO
END DO
END DO
......
......@@ -177,7 +177,7 @@ C***********************************************************************
REAL, INTENT (OUT) :: exc(mgrid) ! xc energy
!
! .. Local Scalars ..
REAL c_1, ex
REAL ex
REAL :: rho, rh1, rh2 ! total, spin up & spin down charge density
REAL :: x, y1, y2, s, thfpi, c_1, rs, beta, bs41, alc
REAL :: xpx, xfx, xlx, xpx0, xfx0, xlx0, ecp, ecf, fs, ec
......
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