Commit 45f7a51b authored by Frank Freimuth's avatar Frank Freimuth

Update wann_uHu_vac.F

parent 6470ca44
......@@ -28,7 +28,7 @@ c***************************************c
implicit none
TYPE(t_zmat), INTENT(IN) :: zMat, zMat_b
TYPE(t_mat), INTENT(IN) :: zMat, zMat_b
c .. scalar Arguments..
logical, intent (in) :: l_noco,l_soc,zrfs
......@@ -276,13 +276,13 @@ c-----> loop over basis functions
IF(zMat%l_real) THEN
do n = 1,nslibd
ac(l,n) = ac(l,n) + zMat%z_r(k+addnoco,n)*av
bc(l,n) = bc(l,n) + zMat%z_r(k+addnoco,n)*bv
ac(l,n) = ac(l,n) + zMat%data_r(k+addnoco,n)*av
bc(l,n) = bc(l,n) + zMat%data_r(k+addnoco,n)*bv
enddo
ELSE
do n = 1,nslibd
ac(l,n) = ac(l,n) + zMat%z_c(k+addnoco,n)*av
bc(l,n) = bc(l,n) + zMat%z_c(k+addnoco,n)*bv
ac(l,n) = ac(l,n) + zMat%data_c(k+addnoco,n)*av
bc(l,n) = bc(l,n) + zMat%data_c(k+addnoco,n)*bv
enddo
END IF
enddo
......@@ -323,13 +323,13 @@ c-----> loop over basis functions
IF(zMat_b%l_real) THEN
do n = 1,nslibd_b
ac_b(l,n) = ac_b(l,n) + zMat_b%z_r(k+addnoco2,n)*av
bc_b(l,n) = bc_b(l,n) + zMat_b%z_r(k+addnoco2,n)*bv
ac_b(l,n) = ac_b(l,n) + zMat_b%data_r(k+addnoco2,n)*av
bc_b(l,n) = bc_b(l,n) + zMat_b%data_r(k+addnoco2,n)*bv
enddo
ELSE
do n = 1,nslibd_b
ac_b(l,n) = ac_b(l,n) + zMat_b%z_c(k+addnoco2,n)*av
bc_b(l,n) = bc_b(l,n) + zMat_b%z_c(k+addnoco2,n)*bv
ac_b(l,n) = ac_b(l,n) + zMat_b%data_c(k+addnoco2,n)*av
bc_b(l,n) = bc_b(l,n) + zMat_b%data_c(k+addnoco2,n)*bv
enddo
END IF
enddo
......
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