You have insufficient data to plot a 3D Figure, as it would require a measurement for every observation angle (combination of azimuth and polar angle). It seems like you have the data of all angles θ for a given ϕ, and all angles ϕ for a given θ, but not every combination (this would require a very long measurement and special equipmentto do, and is not very practical).
However, a 2d polar plot is possible, but simply using polarplot.
You will have to unwrap your phase a bit, since it ranges from -90 to 90 and goes back, but it should go from -180 to 180 or 0 to 360;
Once done, you'll get the plot below, using your data from the last two columns (phi vs S21). You can change the axis, ticks, and all the other things to your liking. In fact, it looks almost identical to the picture you already posted, but mirrored, did you already use the data to plot the picture?: