With[{b = 0.985}, Plot3D[Re[QPochhammer[x + I y]], {x, -b, b}, {y, -b, b}, RegionFunction -> (#1^2 + #2^2 < b^2 &), Mesh -> None, PlotStyle -> Directive[Yellow, Specularity[White, 30]], PlotPoints -> 35]]