Implemented bracketed Newton-Raphson from Numerical Recipes (rtsafe). Saved previous attempts at multi-D NR in leftover_nr.m. Restructured layout of mono2D. Fixed incorrect CFL with dimensional speed of light c
debugging failure of NR convergence in 0.3C beam test. Switched advection flux terms to be added after absorption, but this causes a faster failure of the simulaiton
Added funciton for updating rad moments and eddington variables
Clean up code formatting. Still debugging moving medium equilibration problem