“…These combined fluxes at face centers, f (n ) x,i+1/2,j,k , f (n ) y ,i,j +1/2,k , and f (n ) z ,i,j,k+1/2 , are used to enforce ∇ · B = 0 as well as to update the magnetic field components to the next time step, as the advection only fluxes for the ideal MHD case do. In fact, we used the fluxinterpolated CT (Constraint Transport) scheme developed in Balsara & Spicer(1999).…”