We present, for the first time, three dimensional (3D) radiative transfer (RT) results with realistic atomic models (multilevel) and without using the local thermodinamical equilibrium approximation (non-LTE). We have developed a new code based on efficient iterative methods (Trujillo Bueno, and Fabiani Bendicho 1995; Fabiani Bendicho, Trujillo Bueno and Auer 1997) characterized by a very high convergence rate. With this 3D multilevel code and using a schematic atmospheric model we are able to demonstrate that one may need self-consistent multidimensional RT calculations in order to interpret high spatial resolution solar spectroscopic observations.