      function hlcnc(d,t)
*
* Don Lamb's fit to the Hubbard-Lampe carbon conductive opacities
*
      implicit double precision(a-h,o-z)
      dimension a(25),w(25)

      data a/-5.18321e0,2.88323e-2,-1.68700e0,-3.80388e-1,-1.42808e0, 
     1       1.42606e0,-6.77026e-2,2.99633e-1,3.02198e-1,-1.94700e-1, 
     2       -2.72927e-3,4.24382e-2,-9.77295e-2,-1.50302e-2,8.16693e-3,
     3       4.80447e-4,1.04512e-3,-8.95778e-3,1.27411e-2,-3.37754e-4,
     4       4.17955e-5,-7.52181e-5,-1.38498e-4,5.92069e-4,-5.46719e-4/

      w(1)=1.
      m=0 
      n=1 
      do 2 jmax=1,4 
        n=n+1
        do 1 j=1,jmax 
          m=m+1
          w(n)=d*w(m)
          n=n+1
1       continue
        w(n)=t*w(m)
2     continue
      do 3 i=16,25
        w(i)=d*w(i-5) 
3     continue
      hlcnc=0
      do 4 i=1,25
        hlcnc=hlcnc+a(i)*w(i) 
4     continue
      return
      end 

************************************************************************

