wftsca Subroutine

public subroutine wftsca(drdi, efac, pz, qz, fz, sz, nsra, pzlm, qzlm, pzekdr, qzekdr, ek, loflm, irmind, irmd, irmin, irmax, lmaxd, lmmaxd)

Get wavefunctions of same magnitude by scaling with efac

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in), dimension(*) :: drdi

Derivative dr/di

complex(kind=dp), intent(out), dimension(lmmaxd) :: efac

efac(lm) = sqrt(e)**l/(2l - 1)

complex(kind=dp), intent(in), dimension(irmd, 0:lmaxd) :: pz
complex(kind=dp), intent(in), dimension(irmd, 0:lmaxd) :: qz
complex(kind=dp), intent(in), dimension(irmd, 0:lmaxd) :: fz
complex(kind=dp), intent(in), dimension(irmd, 0:lmaxd) :: sz
integer, intent(in) :: nsra
complex(kind=dp), intent(out), dimension(lmmaxd, irmind:irmd, 2) :: pzlm
complex(kind=dp), intent(out), dimension(lmmaxd, irmind:irmd, 2) :: qzlm
complex(kind=dp), intent(out), dimension(lmmaxd, irmind:irmd, 2) :: pzekdr
complex(kind=dp), intent(out), dimension(lmmaxd, irmind:irmd, 2) :: qzekdr
complex(kind=dp), intent(in) :: ek
integer, intent(in), dimension(*) :: loflm

l of lm=(l,m) (GAUNT)

integer, intent(in) :: irmind

irmd - irnsd

integer, intent(in) :: irmd

Maximum number of radial points

integer, intent(in) :: irmin

Max R for spherical treatment

integer, intent(in) :: irmax
integer, intent(in) :: lmaxd

Maximum l component in wave function expansion

integer, intent(in) :: lmmaxd

(KREL+KORBIT+1)(LMAX+1)*2