Slightly reshuffle enumerated boundary condition types; Add 'bcstatic' to reflect artificially fixed cells for boundaries.
Set default thermal cutoff for radiative shock to 1.1 T0 rather than 1.0 T0; Current explicit cooling algorithm Eint(t+dt) -> Eint(t) - lambda dt experiences problems otherwise.
Change free radiation call to reflect default cutoff temperature of 1.1T0 rather than T0; Reflect change in dt provided from source.m as well.
Alter source.m to reflect that it is taking a timestep equal to two fluid timesteps as part of F S F operator splitting.
Remove debug array readouts from fluid/relaxingFluid.m
Run the Sod tube slightly longer to use most of the grid and autosets # of timesteps so you don't have to worry about that.
Reshuffle BC modes along side reshuffle in ENUM.m
Const BC now does zero gradient extrapolation, move hold-cells-fixed-forever BC to BCMODE_STATIC.
Prevent frame loader from creating blank magnetic arrays for nonmagnetic frames.
Prevent frame loader from creating blank magnetic arrays for nonmagnetic frames.