Prism coupling of radiation inside a linear film over a nonlinear substrate is studied with a numerical method, showing that the change of coupling as a function of the input power in the case of a quasi-continuous operation (thermal nonlinearity) is influenced by a thermal gradient generated in the linear film. A comparison with experimental data confirms these calculations.