Found critical error in advection step: not dividing IV by C in the first substep
Also, fixed so that advection substep IV only is on for nonzero interaction terms
Successfully reproduced fig 3b. How do I find the analytic solution in this case?