A three-dimensional wind field model based on the variational technique has been developed for estimating the overall wind patterns over a complex terrain. The three-dimensional elliptic partial differential equations on Cartesian and terrain-following coordinates have been established to obtain the Lagrangian multiplier and the adjusted wind velocity. The simulations were performed to evaluate the variations of the velocity vectors on the hemisphere, half-cylinder, and saddle type obstacles. Also, the wind field model in the terrain-following coordinate has been applied for evaluating the characteristics of wind patterns according to the variations of Gauss precision moduli on the hemispheric topography. The results showed that the horizontal and vertical wind components were strongly governed by the selection of the values of Gauss precision moduli.