Apply RF pulse in coordinate system rotating with the RF pulse
This makes computing the effect of the RF pulses more robust and allows for larger step sizes even for spins with large off-resonances.
I've also fixed the saturation for RF pulses with a very short T2.