In the TL simulations, the total concentration of the trapping state (N) is chosen to be a fixed value during the different simulation stages (the irradiation and the heating stages). However, in the irradiation stage, this parameter is supposed to be dependent on the dose of excitation (D*t) where D is the dose rate and t is the time of irradiation. Therefore, this parameter N should be given as a function of D and t in the irradiation stage, not a fixed value.