Ran comfortably for B = 2 (before it ran only till B = 1, before the code revamp), with some small pattern near the core. Also ran with B = 10, but crashed after around 100 timesteps, and gives completely wrong results. Also given, but the divergence goes to like 100, so could be really a divergence cleaning issue. Results for all runs improved with higher resolution.