1
2! Copyright (C) 2020 J. K. Dewhurst and S. Sharma.
3! This file is distributed under the terms of the GNU General Public License.
4! See the file COPYING for license details.
5
6subroutine zvcldisp(zvclmt)
7use modmain
8use modtddft
9implicit none
10! arguments
11complex(8), intent(inout) :: zvclmt(npmtmax,natmtot)
12! local variables
13integer is,ia,ias,np,i
14real(8) t1
15do is=1,nspecies
16  np=npmt(is)
17  do ia=1,natoms(is)
18    ias=idxas(ia,is)
19    do i=1,3
20      t1=-datposc(i,0,ia,is)
21! add the gradient of the Coulomb potential of the nucleus plus static density
22      zvclmt(1:np,ias)=zvclmt(1:np,ias)+t1*gvclns(1:np,i,ias)
23    end do
24  end do
25end do
26end subroutine
27
28