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

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