Put in spherical source terms for NR/SR MHD

It's unclear yet whether these are working. Stability is bad near the
inner boundary with MHD switched on, and also depends on the range of
polar and azimuthal angle.