[mesa-users] Radiation pressure at zero optical depth

[Bill Paxton]

here's what the code looks like now.   thanks again!

         
         ! P = (tau*g/kap)*(1 + 1.6d-4*kap*(L/Lsun)/(M/Msun))
         ! the factor comes from including nonzero Prad at tau=0
         ! see, e.g., Cox & Giuli, Section 20.1
         P0 = tau*g/kap         
         Pextra = kap*L/(6d0*pi*clight*cgrav*tau*M)
         ! thanks to Konstantin Pavlovskii for correcting an earlier version of this