Commit 16bf1f54 authored by Uliana Alekseeva's avatar Uliana Alekseeva

made PGI compiler even happier

parent f1f26cc0
......@@ -32,7 +32,11 @@ CONTAINS
CALL gaussp(atoms%lmaxd, rx,wt)
! generate the lattice harmonics on the angular mesh
ALLOCATE ( ylh(nsp,0:sphhar%nlhd,sphhar%ntypsd))
IF (l_grad) ALLOCATE(ylht,ylhtt,ylhf,ylhff,ylhtf,MOLD=ylh )
IF (l_grad) ALLOCATE(ylht,MOLD=ylh )
IF (l_grad) ALLOCATE(ylhtt,MOLD=ylh )
IF (l_grad) ALLOCATE(ylhf,MOLD=ylh )
IF (l_grad) ALLOCATE(ylhff,MOLD=ylh )
IF (l_grad) ALLOCATE(ylhtf,MOLD=ylh )
IF (l_grad) THEN
CALL lhglptg(sphhar,atoms,rx,nsp,xcpot,sym,&
......
......@@ -52,7 +52,8 @@ 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)
ALLOCATE(v_xc,v_x,mold=rho)
ALLOCATE(v_xc,mold=rho)
ALLOCATE(v_x,mold=rho)
CALL xcpot%get_vxc(input%jspins,rho,v_xc,v_x,grad)
......
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