subroutine plasmon(flam,rho,t2,mue,cv,cvp,ca,cap,qplas) * * Plasmon Neutrinos * Have temp dependent coefficients when in low T regime. * I'm including coefficients valid both above and below * log T = 7.8 The entire T range is covered. * implicit double precision (a-h,o-z) dimension a(0:3),b(3) double precision mue xi=(((rho/mue)/1d9)**(0.33333333333d0))/flam if(t2.lt.7.8)then a(0) = 2.787d-7*t2*t2 - 3.936d-6*t2 + 1.408d-5 a(1) = 5.408d-7*t2*t2 - 7.715d-6*t2 + 2.751d-5 a(2) = -1.812d-7*t2*t2 + 2.731d-6*t2 - 1.024d-5 a(3) = 1.639d-8*t2*t2 - 2.445d-7*t2 + 9.137d-7 b(1) = 1.148d-2 b(2) = 1.209d-1 b(3) = 2.432d-4 c = -0.7683 + 0.1798*t2 else a(0) = 2.320d-7 a(1) = 8.449d-8 a(2) = 1.787d-8 a(3) = 0.0 b(1) = 2.581d-2 b(2) = 1.734d-2 b(3) = 6.990d-4 c = 0.56457d0 endif call qcalc2(a,b,c,xi,flam,f) n = 2 qplas=(cv**2+n*cvp**2)*((rho/mue)**3)*f return end ************************************************************************