There is an increasing demand for accurate atomic data due to advancements in experimental techniques and investments in large scale research facilities. In astrophysics the quality and resolution of solar and stellar spectra has so improved that the accuracy of atomic data is frequently a limiting factor in the interpretation. Accurate atomic data are also required in plasma physics and in other emerging areas such as laser spectroscopy on isotope separators, X-ray lithography, and lighting research. The needs include accurate transition energies, fine- and hyperfine structures, mass- and field shifts as well as parameters related to interaction with external magnetic fields. Also there is a constant need for transition rates of different multipolarities between excited states. Data are needed for a wide range of elements and ionization stages. To meet the demands for accurate atomic data the COMPutational Atomic Structure (COMPAS) group has been formed. The group is involved in developing state of the art computer codes for atomic structure calculations in the non-relativistic scheme with relativistic corrections in the Breit-Pauli approximation [1] as well as in the fully relativistic domain. Here we describe new developments of the GRASP2K relativistic atomic structure code [2,3]. We present results for a number of systems and properties to illustrate the potential and restriction of modern computational atomic structure. Among the properties are transition rates, hyperfine- and magnetically induced rates, energy structure, and isotope shifts. We also discuss current code developments and plans for future work. The codes developed by the COMPAS group, along with detailed user manuals, are freely available at http://ddwap.mah.se/tsjoek/compas/ .