potcut Subroutine

public subroutine potcut(imt1, irc1, ins, lmpot, r, vm2z, vspsme, vins, z1, irmd, irmind)

Set potential equal zero between muffin-tin and outer sphere

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: imt1

R point at MT radius

integer, intent(in) :: irc1

R point for potential cutting

integer, intent(in) :: ins

(LPOT+1)**2

integer, intent(in) :: lmpot

(LPOT+1)**2

real(kind=dp), intent(in), dimension(*) :: r

Radial mesh ( in units a Bohr)

real(kind=dp), intent(inout), dimension(*) :: vm2z
real(kind=dp), intent(inout), dimension(*) :: vspsme
real(kind=dp), intent(inout), dimension(irmind:irmd, *) :: vins

Non-spherical part of the potential

real(kind=dp), intent(in) :: z1

Nuclear charge

integer, intent(in) :: irmd

Maximum number of radial points

integer, intent(in) :: irmind

irmd - irnsd