Multivariate simulation is necessary for cases where non-Gaussian processes at spatially distributed locations are desired. A simulation algorithm to generate non-Gaussian wind pressure fields is proposed. Gaussian sample fields are generated based on the spectral representation method using wavelet transforms method and then mapped into non-Gaussian sample fields with the aid of a CDF mapping transformation technique. To illustrate the procedure, this approach is applied to experimental results obtained from wind tunnel tests on the domes. A multivariate Gaussian simulation technique is developed and then extended to multivariate non-Gaussian simulation using the CDF mapping technique. It is proposed to develop a new wavelet-based CDF mapping technique for simulation of multivariate non-Gaussian wind pressure process. The efficiency of the proposed methodology for the non-Gaussian nature of pressure fluctuations on separated flow regions of different rise-span ratios of domes is also discussed.