We describe a package of programs for the implementation of the collisional-radiative model to complex configurations. The number of levels taken into account may be several hundreds. The heart of the package is a very efficient program for excitation cross sections in the Distorted Wave framework, using the Relativistic Parametric Potential wave functions. The basic jj coupling scheme actually simplified the computations, enabling a useful factorization into radial and angular parts. Intermediate coupling and configuration interactions are accounted for. We computed ratios of intensities of 3d9 − 3d84s (E2) to 3d9 −3d84p (El) transitions as functions of ne and Te in Xe XXVIII and other Co-like spectra. The atomic model involves all the levels of configurations (3p6)3d9; −3d84s, −3d84p, −3d84d, −3d84f, and (3p5) −3d10, −3d94p. (275 levels) and all the transitions between them. Results compare very well with experimental spectra from TFR.