Context Navigation

scratch5

-              v3
+              v4
 Next, I calculated the energy injection rate using the luminosity of a FHSC:
+L=3.9e+29 s [[br]]
+E(radiated in a diffusion time)=L*tdiff=8e+35 erg
+I convert this into an energy density, in computational units, by dividing by pscale,
+e(cu)=8e+35erg/Pscale=4.74e+43
+From this I can get an energy injection rate (i.e. the amount of energy density injected into the kernel each timestep),
+de/dt=e(cu)/tsim(cu)=4.74e+43/6.13e-9=7.73246e+51.
+Running this high of an de/dt gave problems, so I decreased this to:
+de/dt=e(cu)/tsim(cu)=4.74e+43/6.13e-9=7.73246e+21.
+L=3.9e+29 erg s^-1^ [[br]]
+Lscale=Pscale*(lscale)^3^/timescale[[br]]
+L(cu)=2.32 e-33
 Next came scaling the opacities.
 The code takes specific opacity (this is in cm^2/g), and then asks whether it should be constant, or multiplied by density (or temperature) to some power. I set it to be a power law scaling.
+The code takes specific opacity (this is in cm^2^/g), and then asks whether it should be constant, or multiplied by density (or temperature) to some power. I set it to be a power law scaling.
 tpow=0 [[br]]