Knowledge of the poloidal distribution of zonal flow (ZF) drive has been a missing yet crucial component of both our understanding of turbulence–flow interactions and the validation of gyrokinetic codes. We present the first analysis of the distribution of electrostatic ZF drive due to Reynolds stress. Using gyrokinetic flux-tube simulations we examine the effect of strong axisymmetric shaping, including elongation, triangularity, and aspect ratio, as well as up–down asymmetric equilibria. With increased shaping, the ZF drive develops local maxima near those of poloidal curvature, with a tendency to shift towards the so-called bad curvature side. Thus, depending on the shaping, the nonlinear ZF drive may not peak at the outboard midplane where the turbulent fluctuations are the strongest and where turbulence diagnostics are usually located. Our results therefore suggest that a shaping correction must be taken into account when extrapolating from a poloidally localized measurement to the total electrostatic ZF drive on a flux-surface.