We consider the two-dimensional classical XY model on a square lattice in the thermodynamic limit using tensor renormalization group and precisely determine the critical temperature corresponding to the Berezinskii-Kosterlitz-Thouless (BKT) phase transition to be 0.89290(5) which is an improvement compared to earlier studies using tensor network methods.