saturation BH loop