Natural ventilation and the use of fans are recognized as sustainable design strategies to reduce energy use while reaching thermal comfort. A big challenge for designers is to predict ventilation rates of buildings in dense urban areas. One significant factor for calculating the ventilation rate is the wind pressure coefficient (Cp). Cp values can be obtained at a high cost, via real measurements, wind tunnel experiments, or high computational effort via computational fluid dynamic (CFD) simulation. A fast surrogate model to predict Cp for a schematic urban environment is required for the integration in building performance simulations. There are well-known surrogate models for Cp. The average surface pressure coefficient model integrated in EnergyPlus considers only a box-shaped building, without surrounding buildings. CpCalc, a surrogate model for Cp, considers only one height of neighbouring buildings. The Toegepast Natuurwetenschappelijk Onderzoek (TNO) Cp Generator model was available via web interface, and could include several box-shaped buildings in the surrounding area. These models are complex for fast integration in a natural ventilation simulation. For optimization processes, with thousands of simulation runs, speed is even more essential. Our study proposes a new surrogate model for Cp estimation based on data obtained from the TNO CP Generator model. The new model considers the effect of different neighbouring buildings in a simplified urban configuration, with an orthogonal street pattern, box-shaped buildings, and repetitive dimensions. The developed surrogate model is fast, and can easily be integrated in a dynamic energy simulation tool like EnergyPlus for optimization of natural ventilation in the urban areas.