Saturate Integral at Hypotenuse of 2D Signal

Hi,

I am integrating a 2D signal, which I want to saturate with the integrator block to avoid overflow. However, the saturation limits for each dimension are dependent on one another.

Essentially what I am doing is integrating a velocity vector to get position, but I want to put a limit on the hypotenuse of the position vector. You can think of this as a radial limit on cartesian coordinates. Hopefully that makes sense.

The built-in integrator block only allows me to put fixed values as saturation limits, so this won't work for me. I tried creating my own saturation loop but I can't do this without creating an algebraic loop.

I want to believe this is possible to do, since it makes sense physically. But mathematically, I'm stuck. Is it possible to do or am I going to need to resort to a discrete solution? Maybe if I do something with polar coordinates?

Thanks!