REAL*8 FUNCTION VN02B9(J1,J2,J3,L1,L2,L3)
C
C  CALCULATION OF WIGNER SIXJ SYMBOL WITH THE UPPER ROW SYMBOLS
C  J1,J2,J3  AND THE LOWER ROW SYMBOLS L1,L2,L3
C
      IMPLICIT REAL*8 (A-H,O-Z)
      COMMON/VN06FC/FCT(40)
      CC=0.0
      IF(J1+J2-J3) 20,1,1
  1   IF(IABS( J1-J2)-J3) 2,2,20
   2  IF (J1+J2+J3-2*((J1+J2+J3)/2)) 20,3,20
  3   IF(J1+L2-L3) 20,4,4
  4   IF (IABS(J1-L2)-L3) 5,5,20
  5   IF (J1+L2+L3-2*((J1+L2+L3)/2)) 20,6,20
  6   IF(L1+J2-L3) 20,7,7
  7   IF(IABS(L1-J2)-L3) 8,8,20
  8   IF(L1+J2+L3-2*((L1+J2+L3)/2)) 20,9,20
  9   IF(L1+L2-J3) 20,10,10
  10  IF(IABS(L1-L2)-J3) 11,11,20
  11  IF(L1+L2+J3-2*((L1+L2+J3)/2)) 20,12,20
  12  OMEGA=0.0
      IF(J3) 37,38,37
  37  IF(L3) 40,39,40
  38  VN02B9=(-1.)**((J1+L2+L3)/2)/SQRT((FLOAT(J1)+1.)*(FLOAT(L2)+1.))
      GO TO 41
  39  VN02B9=(-1.)**((J1+J2+J3)/2)/SQRT((FLOAT(J1)+1.)*(FLOAT(J2)+1.))
      GO TO 41
  40  IWMIN=J1+J2+J3
      IF(IWMIN-J1-L2-L3) 13,14,14
  13  IWMIN=J1+L2+L3
  14  IF(IWMIN-L1-J2-L3) 15,16,16
  15  IWMIN=L1+J2+L3
  16  IF(IWMIN-L1-L2-J3) 17,18,18
  17  IWMIN=L1+L2+J3
  18  IWMAX=J1+J2+L1+L2
      IF(IWMAX-J2-J3-L2-L3) 22,22,23
  23  IWMAX= J2+J3+L2+L3
  22  IF (IWMAX-J1-J3-L1-L3) 24,24,25
  25  IWMAX=J1+J3+L1+L3
  24  IF(IWMIN-IWMAX) 26,26,20
  26  DO 701 IW = IWMIN,IWMAX,2
      IW1=IW/2+2
      IW2=(IW-J1-J2-J3)/2+1
      IW3=(IW-J1-L2-L3)/2+1
      IW4=(IW-L1-J2-L3)/2+1
      IW5=(IW-L1-L2-J3)/2+1
      IW6=(J1+J2+L1+L2-IW)/2+1
      IW7=(J1+J3+L1+L3-IW)/2+1
      IW8=(J2+J3+L2+L3-IW)/2+1
      IF(IW-4*(IW/4)) 30,31,30
  31  PH=1.0
      GO TO 35
  30  PH=-1.0
  35  OMEGA=OMEGA+PH*FCT(IW1)/FCT(IW2)/FCT(IW3)/FCT(IW4)/FCT(IW5)
     1/FCT(IW6)/FCT(IW7)/FCT(IW8)
  701 CONTINUE
      CC=OMEGA*VN02BB(J1,J2,J3,FCT)*VN02BB(J1,L2,L3,FCT)*VN02BB(L1,J2,L3
     1,FCT)*VN02BB(L1,L2,J3,FCT)
  20  VN02B9=CC*10.0
 41   RETURN
      END
      REAL*8 FUNCTION VN02BB(J1,J2,J3,FCT)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION FCT(40)
      IW1=(J1+J2-J3)/2+1
      IW2=(J1-J2+J3)/2+1
      IW3=(-J1+J2+J3)/2+1
      IW4=(J1+J2+J3+2)/2+1
      FDELTA=SQRT(FCT(IW1)*FCT(IW2)*FCT(IW3)/FCT(IW4))
      VN02BB=FDELTA/3.16227765
      RETURN
      END
      PROGRAM SIXJ
C
C  TEST PROGRAM TO EVALUATE THE WIGNER SIXJ SYMBOLS
C
      COMMON/VN06FC/FCT(40)
      IMPLICIT REAL*8(A-H,O-Z)
      FCT(1)=1.
      DO 400 I=2,40
 400  FCT(I)=FCT(I-1)*(I-1.)/10.
 4    READ(1,2)J1,J2,J3,L1,L2,L3,ICON
 2    FORMAT(7I2)
      X=VN02B9(J1,J2,J3,L1,L2,L3)
      WRITE(6,3)J1,J2,J3,L1,L2,L3,X
 3    FORMAT(1H ,6I4,D22.12)
      IF(ICON.EQ.0)GOTO 5
      GOTO 4
 5    STOP
      END