Multilevel Monte Carlo method for parabolic stochastic partial differential equations