Thermally induced phase mismatching (TIPM) has been proven to be an influential issue in nonlinear phenomena. It occurs when refractive indices of crystal are changed due to temperature rise. In this work, the authors report on a modeling of spatiotemporal dependence of TIPM in a repetitively pulsed pumping KTP crystal. Gaussian profiles for both spatial and temporal dependences of pump beam were used to generate second-harmonic waves in a type II configuration. This modeling is of importance in predicting the nonlinear conversion efficiency of crystals when heat is loaded in the system. To this end, at first, an approach to solve the heat equation in a repetitively pulsed pumping system with consideration of the temperature dependence of thermal conductivity and realistic cooling mechanisms such as conduction, convection, and radiation, is presented. The TIPM is then calculated through the use of experimental thermal dispersion relations of KTP crystal. The results show how accumulative behaviors of temperature and TIPM (with its reverse sign) happen when the number of pulses is increased. Fluctuations accompanying temperature and TIPM were observed which were attributed to the off-time between successive pulses. Moreover, in this work, a numerical procedure for solving a repetitively pulsed pumped crystal is introduced. This procedure enables us to solve the problem with home-used computing machines.