Add explicit calculation of M Neumann condition

(hoping that it might fix |m| conservation)