saturation M[H] loop