Let u be a positive harmonic function on the unit disk. Bourgain showed that the radial variation of u is finite for many points θ, and moreover, that the set is dense in the unit circle T and its Hausdorff dimension equals one. In the paper, this result is generalized to a class of smooth domains in Rd, d ≥ 3.