Wind flow perturbations, recirculations and turbulence generated by buildings often dominate air pollutant distributions around buildings. This paper describes dispersion of contaminants in the vicinity of a building by solving the concentration equation based on previously simulated wind flow field. Turbulence closure is achieved by using the standard k-ε two-equation model. The paper shows application of the CIP method for solving a species concentration equation of contaminant gas around a rectangular building for two different sources under conditions of neutral atmospheric stratification. Results have been compared to the experimental data and the previous numerical results by hybrid scheme. The computational results of concentration profiles by the CIP method agree well with experimental data.