Careful modeling of the radio channel characteristics is an important issue for system level simulations. While a number of detailed path loss, shadowing and fast fading models can be found in literature, many of them entail large computational efforts and hence lead to long simulation times, or require huge amounts of memory. In this paper, we focus on shadowing models with two-dimensional correlation properties. We present handover performance evaluations as a use case from which requirements to the modeling of the shadow fading are derived. For different implementations of the shadow fading in system level simulations, we discuss their pros and cons and then focus on the sum of sinusoid model proposed by Cai and Giannakis [1]. Our main contribution lies in a new frequency sampling method to determine the coefficients of the sinusoidal waveforms. We show that this new method, denoted as Power Sampling Method (PSM), allows for a significant reduction of the number of sinusoids while the correlation properties of the shadow fading channel are preserved.