Set the following scales in physics.data:

lscale	3.08 e+18 cm
rscale	2.02 e-16 g/cm3
tempscale	1 K

From that initialized an ambient object with the following density and pressure,

ambient%rho=2.03d-16/rscale
ambient%pressure=(rhoOut/rscale)*(tempOut/tempscale)

Have my grid running from, gxbounds= -.0068 -.0068 0 .0068 .0068

I ran the sim to frame 0 to get the timescale from scales.data.

I used this to get the tsim=tdiff in computational units,

tsim=2.074e+6s/!338402526430554s = 6.13e-9

Next, I calculated the energy injection rate using the luminosity of a FHSC:

L=3.9e+29 erg s^-1
Lscale=Pscale*(lscale)³/timescale
L(cu)=2.32 e-33

Next came scaling the opacities.

The code takes specific opacity (this is in cm²/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
dpow=1

so that,

opacity=kappa*rho

I scaled the planck and rosseland opacities into computational units by:

kappa_r=0.23*mscale/(lscale)²=144.07
kappa_p=.4*mscale/(lscale)²=250.56

Radiation was set to 'diffusion and energy exchange only', rad limiter was set to 'diffusion limit'.

This all seems right to me, so then here are the results.

Run 1

de/dt=7.73246e+51

The code failed to run past frame zero.