Slightly reorder operations in calculation of apar_m to improve cancellation

For half the calls to set_source we'd expect apar = aparnew and hence apar_m (= dA/dt) == 0. However, with the existing operation order this cancellation would not occur exactly. This is unlikely to have a significant impact on anything but can help with some detailed testing.