Mesoscale modeling of organic semiconductors relies on solving an appropriately parametrized master equation. Essential ingredients of the parametrization are site energies (driving forces), which enter the charge transfer rate between pairs of neighboring molecules. Site energies are often Gaussian-distributed and are spatially correlated. Here, we propose an algorithm that generates these energies with a given Gaussian distribution and spatial correlation function. The method is tested on an amorphous organic semiconductor, DPBIC, illustrating that the accurate description of correlations is essential for the quantitative modeling of charge transport in amorphous mesophases.