Analyzing or optimizing wind farm layouts often requires reduced-order wake models to estimate turbine wake interactions and wind velocity. We propose a wake model for vertical-axis wind turbines (VAWTs) in streamwise and crosswind directions. Using vorticity data from computational fluid dynamic (CFD) simulations and cross-validated Gaussian distribution fitting, we produced a wake model that can estimate normalized wake velocity deficits of an isolated VAWT using normalized downstream and lateral positions, tip-speed ratio, and solidity. Compared to CFD, taking over a day to run one simulation, our wake model predicts a velocity deficit in under a second with an appropriate accuracy and computational cost necessary for wind farm optimization. The model agreed with two experimental studies producing percent differences of the maximum wake deficit of 6.3% and 14.6%. The wake model includes multiple wake interactions and blade aerodynamics to calculate power, allowing its use in wind farm layout analysis and optimization. Keywords vertical-axis wind turbine, wind farm optimization, wake model, computational fluid dynamics, data parameterization u x-component of velocity (downstream direction) v y-component of velocity (lateral direction) w wake integration lateral extent