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},Shading->False]
