      function hlcnh(d,t)
*
* Gilles Fontaine's fit to the Hubbard-Lampe Hydrogen conductive opacities
*
      implicit double precision(a-h,o-z)

      dimension a(5,5),temp(5), rho(5)
      data a/12.40928,     -2.545420,   -3.019045,    6.746942e-2,
     1    4.686838e-2,-15.29732,      .8154491,   2.387849, 
     2   -6.939705e-2, -4.129124e-2, 5.162517,   -3.839103e-1,
     3   -6.899530e-1,  2.568522e-2, 1.281765e-2, -.6347533,
     4    7.203151e-2,  8.405360e-2,-3.965926e-3,-1.652045e-3,
     5    2.699066e-2, -4.243979e-3,-3.648879e-3, 
     6    2.127303e-4,  7.509063e-5/

      sum = 0.
      temp(1) = 1.
      rho(1)  = 1.
      do 10 i = 2,5 
        temp(i) = temp(i-1) * t
        rho(i)  = rho(i-1) * d
10    continue
      do 2 i=1,5
        do 1 j=1,5
           sum = sum + a(i,j)*temp(j)*rho(i)
1       continue
2     continue
      hlcnh = sum
      return
      end 

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

