We consider products of independent square non-Hermitian random matrices. More precisely, let X1, . . . , Xn be independent N × N random matrices with independent entries (real or complex with independent real and imaginary parts) with zero mean and variance 1 N . Soshnikov-O'Rourke [15] and Götze-Tikhomirov [11] showed that the empirical spectral distribution of the product of n random matrices with iid entries converges toWe prove that if the entries of the matrices X1, . . . , Xn satisfy uniform subexponential decay condition, then in the bulk the convergence of the ESD of X1 · · · Xn to (1) holds up to the scale N −1/2+ε .