Properly reduce phase into [0, 2pi] range.  The fits no longer complain about variable being set out of bounds.