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

      dimension a(5,5),temp(5), rho(5)
      data a /9.246777e+0, -2.097036,   -1.267003,  5.490678e-2,
     1    2.463979e-3, -1.265387e+1, 1.716563e-1,   1.088986e+0,
     2   -4.576479e-2, -6.966255e-3, 4.539412e+0, -1.427574e-1,
     3     -3.499339e-1,  1.589662e-2, 3.421386e-3, -5.842722e-1,
     4    3.930197e-2,  4.702301e-2,-2.529960e-3, -5.814372e-4,
     5    2.617848e-2, -2.791279e-3,-2.238884e-3, 
     6    1.442103e-4,  3.234929e-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
      hlcnhe = sum 
      return
      end 

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

