anchor.diff:        implicit none
         integer l1,l2,l3,m1,m2,m3
         integer l_min,l_mid,l_max,m_min,m_mid,m_max,
  -     &    ll_min,ll_mid,ll_max,l_sum
  +     &    ll,ll_min,ll_mid,ll_max,l_sum
         integer i,j,k,l,m,n,phase
  -      double precision result,wa,wb,wc
  +      double precision result,wa,wb
   *
         result = 0.0d00
         m3 = -(m1+m2)
  -C      write (6,'(3I4)') l1,l2,l3
  -C      write (6,'(3I4)') m1,m2,m3
         if (abs(m1) .gt. l1) return
         if (abs(m2) .gt. l2) return
         if (abs(m3) .gt. l3) return
  
positive.code:         else
            wb = (l_max+m_mid)*(l_max+m_mid+1)/2
          end if
  C        write (6,*) wa,wb
          result = sqrt(wb)/wa
  *
  *   General code for l_min > 1
  *
        else
          l_sum = l1+l2+l3
          ll_max = l_sum
          ll_min = 2*l_min
          ll_mid = 2*l_mid
          m = 1
  
negative.code: C$Id: ecp_3j_prod.f,v 1.5 1997-05-22 23:52:50 mg201 Exp $
  ************************************************************************
  *                                                                      *
        subroutine ecp_3j_prod (l1,l2,l3,m1,m2,result)
  