gengclqu.f90

Go to the documentation of this file.

! Copyright (C) 2025 Wenhan Chen, J. K. Dewhurst and S. Sharma.
! This file is distributed under the terms of the GNU General Public License.
! See the file COPYING for license details.

subroutine gengclqu
use modmain
use modulr
implicit none
! local variables
integer iq
if (q0cut >= 0.d0) then
! hard cut-off
  do iq=1,nqpt
    gclq(iq)=merge(fourpi/qc(iq)**2,0.d0,qc(iq) > q0cut)
  end do
else
! Yukawa-type screening
  gclq(1)=0.d0
  gclq(2:nqpt)=fourpi/(qc(2:nqpt)**2+q0cut**2)
end if
end subroutine