Stuck trying to track down small discrepancies from paper Fig 2 and 3. Created a whole new version with an advection step solver that mimics the layout of
ATHENA. No matter what I try, the modifications cause the NR solve to fail by nt=6.