xmin := 0
xmax := 2*Pi
ymin := 0
ymax := 2*Pi
xxx[x_, y_] := (3 + Cos[x])*Cos[3*y]*y/(Pi*2)
yyy[x_, y_] := (3 + Cos[x])*Sin[3*y]*y/(Pi*2)
zzz[x_, y_] := Sin[x] - 6*y/(Pi*2) + 3
ParametricPlot3D[{yyy[x, y], -xxx[x, y], zzz[x, y]}, {x, xmin, 
        xmax}, {y, ymin, ymax}, PlotPoints -> {15, 51}]