An additive finite-difference scheme for numerical approximation of initial-boundary value problem for two-dimensional fractional in time diffusion equation is proposed. Its stability is investigated and a convergence rate estimate is obtained.