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 ************************************************************************