Methods of numerical simulation of MHD processes are discussed for the case of finite electrical conductivity. Three main problems are mentioned. The first is the exact description of thin dissipative layers around resonant toroidal magnetic surfaces – this is connected with the fast growth of conductivity in the external limiting layer, which makes problematic the construction of the difference schemes in its vicinity. The second is that in the presence of finite electric conductivity the differential operators of the corresponding MHD equations are no longer selfadjoint, which leads to complex eigenvalues in the linear approximations and finally to a continuous spectrum instead of a discrete one. The third is the recombination of magnetic surfaces if the conductivity is finite. This work proposes, with the short outline of formerly published works, a way to overcome the mentioned difficulties.