Integral cohomology of certain Picard modular surfaces

Let Gamma be a congruence subgroup of the Picard modular group of an imaginary number field k, and let D be the associated symmetric space. We describe a method to compute the integral cohomology of the locally symmetric space Gamma\D. The method is implemented for the case k=Q(i) and k=Q(sqrt(-3)), and the cohomology is computed for various Gamma.