A parallel method for time-discretization of parabolic problems based on contour integral representation and quadrature