Inconsistent handling of B0 and B1 #28
Description
In 2d MHD solar wind box with strong background dipole field, simulation stops with negative pressure if entire box is also initialized with B1 = (44,13,0) nT, but runs fine if this field is added to B0 instead. In principle result should be identical in both cases. Problem doesn't seem to be related to boundaries.