This article considers an inverse problem of time fractional parabolic partial differential equations with the nonlocal boundary condition. Dirichlet-measured output data are used to distinguish the unknown coefficient. A finite difference scheme is constructed and a numerical approximation is made. Examples and numerical experiments, such as man-made noise, are provided to show the stability and efficiency of this numerical method.