bode_plot phase wrapping incorrect for multiple systems.

Minimum example:

# two RHP's
G = zpk([],[1,1], gain=1)  
bode_plot(G)
bode_plot((G,G))

The second call shows a figure with different phases for each (identical) system.

bode1
bode2