Commit 802af398 authored by Matthias Redies's avatar Matthias Redies

check if lapl is allocated

parent 74296079
......@@ -126,11 +126,13 @@ CONTAINS
REAL :: q(input%jspins), qis(input%jspins), qmt(atoms%ntype,input%jspins),&
qvac(2,input%jspins), qtot, qistot
call integrate_grid(xcpot, stars, atoms, sym, vacuum, input, cell, oneD, sphhar,noco,&
xcpot%is_lapl, xcpot%mt_lapl, &
q, qis, qmt, qvac, qtot, qistot)
if(allocated(xcpot%is_lapl%grid)) then
call integrate_grid(xcpot, stars, atoms, sym, vacuum, input, cell, oneD, sphhar,noco,&
xcpot%is_lapl, xcpot%mt_lapl, &
q, qis, qmt, qvac, qtot, qistot)
call print_qs(q,qis,qmt,qvac,qtot,qistot,"laplace")
call print_qs(q,qis,qmt,qvac,qtot,qistot,"laplace")
endif
END SUBROUTINE integrate_lapl
SUBROUTINE integrate_kED_schr(xcpot, stars, atoms, sym, vacuum, input, cell, oneD, sphhar,noco)
......
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