SUBROUTINE wigner(jtot, lam, k, beta, djmkp )
!
!            P U R P O S E   O F    S U B R O U T I N E
! Calculates parity possessing Wigner d functions.

!            I N P U T    A R G U M E N T S
!  nx
!  ny
!  nxdim
!  nydmi
!  jr
!  jrjtmax
!  betqf
!  djmkp
!  ic
!  jtot
!  lam is the aph index.
!  a
!  cnorm
!  parity      Desired parity.

!            O U T P U T    A R G U M E N T S
!
!this routine is called by: AphDel_Old or AphDel_New
!this routine calls
!
!-----------------------------------------------------------------------
! <><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>

!USE FileUnits_Module
USE Numeric_Kinds_Module

IMPLICIT NONE

REAL(dp) betqf, djmkp, beta, a, cnorm, djmk
INTEGER nx, ny, nxdim, nydim, jrjtmax, ic, jtot, k,lam, parity, iy, ix
EXTERNAL djmk
DIMENSION a(200)

!--------------------------------------------------------------------
!  now construct parity possessing Wigner small d functions
!-------------------------------------------------------------------
  cnorm=1.d0/Sqrt(2.d0)
  IF (lam.eq.0) cnorm=0.5d0

  djmkp=djmk(jtot,k,lam,beta,a,200)+(-1)**(jtot+lam+parity)*djmk(jtot,k,-lam,beta,a,200)
  djmkp=cnorm*djmkp

RETURN
END