In the paper a new, fractional order, discrete model of a two-dimensional temperature field is addressed. The proposed model uses Grünwald-Letnikov definition of the fractional operator. Such a model has not been proposed yet. Elementary properties of the model: practical stability, accuracy and convergence are analysed. Analytical conditions of stability and convergence are proposed and they allow to estimate the orders of the model. Theoretical considerations are validated using exprimental data obtained with the use of a thermal imaging camera. Results of analysis supported by experiments point that the proposed model assures good accuracy and convergence for low order and relatively short memory length.