[en] The exotic nuclei have attracted more and more attention in nuclear physics. As a candidate for the description of exotic nuclei, the Relativistic Continuum Hartree-Bogoliubov (RCHB) [1] theory provides a fully self-consistent treatment of pairing correlations in the presence of the continuum and has well reproduced the exotic halo phenomena in 1^1Li. But the theory encounters numerical difficulties in the extension to the deformed systems in the coordinate space. Among all the state of arts recipes available in the market, imaginary time step (ITS) [2] is one of the most promising methods for deformed exotic system. However, whether or not the ITS method can be applied for relativistic systems has received wide attention due to the existence of the Dirac sea. Specifically, this question has been raised by several world-leading physicists in recent international conferences and workshops. In this work [3], the discrete single-particle spectra in both the Fermi and Dirac sea have been calculated by the ITS method for the Schrodinger-like equation for the first time after avoiding the 'tsunami' of the Dirac sea, i.e., the diving behavior of the single-particle level into the Dirac sea in the direct application of the ITS method for the Dirac equation. It is found that by the transform from the Dirac equation to the Schrodinger-like equation, the single-particle spectra, which extend from the positive to the negative infinity, are separated into two branches confined in the Fermi sea and in the Dirac sea respectively. Identical results with those in the conventional shooting method have been obtained via the ITS evolution for the equivalent Schrodinger-like equation, which demonstrates the feasibility, practicality and reliability of the present algorithm and dispels the doubts on the ITS method in the relativistic system.(author)