Commit fb019632 authored by Matthias Redies's avatar Matthias Redies

use new nsp() function

parent 19de5dee
......@@ -105,9 +105,8 @@ CONTAINS
REAL :: theta,phi
REAL,ALLOCATABLE :: ch(:,:)
nsp=(atoms%lmaxd+1+MOD(atoms%lmaxd+1,2))*(2*atoms%lmaxd+1)*atoms%jmtd
ALLOCATE(ch(nsp,4))
nsp=nsp/atoms%jmtd
nsp=atoms%nsp()
ALLOCATE(ch(nsp*atoms%jmtd,4))
CALL xcpot%init("vwn",.FALSE.,1)
CALL init_mt_grid(nsp,4,atoms,sphhar,xcpot,sym)
......
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