Commit 92b2de19 authored by Daniel Wortmann's avatar Daniel Wortmann

Modified scaling factor in laplacian

parent f24144ac
......@@ -82,7 +82,8 @@ CONTAINS
END IF
IF (ALLOCATED(grad%laplace)) THEN
!Lapacian of density
fac=MERGE(2,1,jspins==1)
!fac=MERGE(2,1,jspins==1)
fac=1
DO js=1,jspins
DO i=1,nsp
grad%laplace(kt+i,js) = (rhdrr(i,js) + 2.e0*rhdr(i,js)/rv1 +&
......
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