Fix rho calculation in vdw-DF.F for parallel threads

Author: ebylaska

src/nwpw/pspw/lib/exchange-correlation/vdw-DF/vdw-DF.F

@@ -393,7 +393,7 @@ subroutine vdw_DF_Generate_rho(ispin,n2ft3d,dn,rho)
       nthr = Parallel_nthreads()
 
       do i=tid+1,n2ft3d,nthr
-         rho(i) = dn(i,1)+dn(1,ispin)+dncut
+         rho(i) = dn(i,1)+dn(i,ispin)+dncut
       end do
 !$OMP BARRIER