Commit 81708430 authored by Matthias Redies's avatar Matthias Redies

bugfix

parent 39553231
......@@ -58,7 +58,7 @@ CONTAINS
!Put the charge on the grid, in GGA case also calculate gradients
CALL pw_to_grid(xcpot,input%jspins,noco%l_noco,stars,cell,den%pw,grad,rho)
xcpot%is_lapl%grid = grad%laplace
if(allocated(grad%laplace)) xcpot%is_lapl%grid = grad%laplace
ALLOCATE(v_xc,mold=rho)
ALLOCATE(v_x,mold=rho)
......
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