This paper presents a three dimensional numerical model of flow and movement of gas condensate spills based on Navier-Stokes and continuity equations with Boussinesq approximation involving various surface wind forcings. The model simulates the surface movement of gas condensate slick from spill accidents in Assalouyeh Marine Region. For the advection term an upwind weighted, multidimensional positive definite advection transport algorithm (MPDATA) was used. This algorithm uses an explicit finite difference scheme with an antidiffusive velocity for equilibrium diffusion. It also uses a generalized-conjugate-residual (GCR) method for the solution. The model is run for gas condensate spill accidents in Assalouyeh Marine Region in summer and winter of 2005. Numerical results show that gas condensate particles spread torwards the shore in summer, while in winter it mostly spreads towards east. The spreadings follow the flow fields that are in good agreement with flow field observations. Diffusion of gas condensate particles in the water due to more turbulence in winter is larger, while gas condensate particles are observed on the water surface due to more stability and buoyancy force in summer.