Extreme value analyses of a large number of relatively short time series are in increasing demand in environmental sciences and design. Here, we present an automated procedure for the peaks-over-threshold (POT) approach to extreme value theory and use it to provide a climatology of extreme hourly precipitation in Switzerland. The POT approach fits the generalized Pareto distribution (GPD) to independent exceedances above some high threshold. To guarantee independence, the time series is pruned: exceedances separated by less than a fixed interval called the run parameter are considered a cluster, and all but the cluster maxima are discarded. We propose the automation of an existing graphical method for joint selection of threshold and run parameter. Hourly precipitation is analyzed at 59 stations of the MeteoSwiss observational network over the period 1981-2010. The four seasons are considered separately. When necessary, a simple detrending is applied. Results suggest that unnecessarily large run parameters have adverse effects on the estimation of the GPD parameters. The proposed method yields mean cluster sizes that reflect the seasonal and geographical variation of lag dependence of hourly precipitation. The climatology, as represented by the return level maps and Alpine cross-section, mirror known aspects of the Swiss climate. Unlike for daily precipitation, summer thunderstorm tracks are visible in the seasonal frequency of events, rather than in the amplitude of rare events.