Stress relaxation in single crystalline thin films on substrates subjected to thermal loading is studied using a recently proposed nonlocal continuum crystal plasticity theory. The theory is founded on a statistical-mechanics description of the collective behaviour of dislocations in multiple slip, which is coupled to a small-strain continuum crystal plasticity description. The theory is inherently nonlocal with the length scale being determined by the evolving dislocation density. Symmetric double slip is considered with the film being in plane strain. The predicted stress versus temperature response and the evolution of the dislocation structure are analyzed for different orientations and film thicknesses. The effect of film size is assoc...