In many of the plasmas generated in many fields of the high energy density physics, the radiation can significantly alter the material response. A proper microscopic description of these plasmas entails integrated computer codes that self-consistently combines large-scale atomic kinetics and radiation transport. Due to the inherent complexity of this type of codes and its interest in the area of high energy density physics, new developments in this field are welcomed. In this work, we present MIXKIP/RAPCAL, an integrated computational package to perform 1D and 2D large-scale non-local thermodynamic equilibrium atomic kinetics and radiation transfer coupled simulations for high energy density plasmas. This package includes different modules that allow simulations of non-equilibrium plasmas under different degrees of detail and accuracy, depending on the requirements of the situations to analyze. Comparisons with experimental results of homogeneous optically thick plasmas of gold and xenon are presented in order to check its accuracy. The influence of different approaches of the spatial discretization of the plasma in the radiation dependent atomic kinetics simulations is also analyzed. Finally, this study is also made for a non-homogeneous optically thick aluminum plasma.