Commit 749abd62 authored by Matthias Redies's avatar Matthias Redies

remove debugging mistakes

parent 4e05d9ea
...@@ -49,8 +49,6 @@ CONTAINS ...@@ -49,8 +49,6 @@ CONTAINS
REAL, ALLOCATABLE :: rho(:,:), ED_rs(:,:), vTot_rs(:,:), kinED_rs(:,:) REAL, ALLOCATABLE :: rho(:,:), ED_rs(:,:), vTot_rs(:,:), kinED_rs(:,:)
REAL, ALLOCATABLE :: v_x(:,:),v_xc(:,:),e_xc(:,:) REAL, ALLOCATABLE :: v_x(:,:),v_xc(:,:),e_xc(:,:)
REAL, PARAMETER :: fac = 1.5
CALL init_pw_grid(xcpot,stars,sym,cell) CALL init_pw_grid(xcpot,stars,sym,cell)
!Put the charge on the grid, in GGA case also calculate gradients !Put the charge on the grid, in GGA case also calculate gradients
...@@ -76,7 +74,7 @@ CONTAINS ...@@ -76,7 +74,7 @@ CONTAINS
cell, EnergyDen%pw, tmp_grad, ED_rs) cell, EnergyDen%pw, tmp_grad, ED_rs)
CALL pw_to_grid(xcpot, input%jspins, noco%l_noco, stars, & CALL pw_to_grid(xcpot, input%jspins, noco%l_noco, stars, &
cell, vTot%pw, tmp_grad, vTot_rs) cell, vTot%pw, tmp_grad, vTot_rs)
CALL calc_kinEnergyDen(fac * ED_rs, vTot_rs, rho, kinED_rs) CALL calc_kinEnergyDen(ED_rs, vTot_rs, rho, kinED_rs)
ENDIF ENDIF
!calculate the ex.-cor energy density !calculate the ex.-cor energy density
......
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