# Phase wrapping in Bode plot

I'm trying to model the transfer function of a linear system (a MOEMS mirror) based on a datasheet plot of its response. I'm developing the model in the python-control model and I think that I've got the model pretty close, except that the Bode plot appears to wrap the phase response to keep it within $\pm 180^\circ$. My question is threefold:

1. Is there a material difference between my model and what's in the datasheet with respect to the frequency wrapping?

2. Is there a simple way in the python-control library to unwrap the plot

3. Are there any other suggestions to improve my model?

Here is the response plots from the datasheet:

Here is my model:

In order to improve you model I would take a slightly better at the pole location, since their phase does not seem to drop instantly, so a small real part. There also might be a small delay ($e^{-Ts}$), since the phase of the data drops further than -270°.