The fractional Fokker--Planck equation is an important physical model for simulating anomalous diffusions with external forces. Because of the nonlocal property of the fractional derivative an interesting problem is to explore high accuracy numerical methods for fractional differential equations. In this paper, a space-time spectral method is presented for the numerical solution of the time fractional Fokker--Planck initial-boundary value problem. The proposed method employs the Jacobi polynomials for the temporal discretization and Fourier-like basis functions for the spatial discretization. Due to the diagonalizable trait of the Fourier-like basis functions, this leads to a reduced representation of the inner product in the Galerkin analysis. We prove that the time fractional Fokker--Planck equation attains the same approximation order as the time fractional diffusion equation developed in [X. Li and C. Xu, SIAM J. Numer. Anal., 47 (2009), pp. 2108--2131] by using the present method. That indicates an e...