Abstract. We derive a thin film model for viscoelastic liquids under strong slip which obey the stress tensor dynamics of corotational Jeffreys fluids.PACS. 83.60.Bc Viscoelasticity -47.50.+d Non-newtonian fluid flows -68.15.+e Liquid thin films Dewetting liquid polymer films on nonwetting substrates, such as silicone wafers grafted with a monolayer of brushes, play a prominent role in many nanotechnological applications. It is known that for these situations, polymer films on the scale of a few hundred nanometers typically show large slippage [1]. Furthermore, for highly entangled polymers, the assumption of a Newtonian fluid will seize to be valid. To understand the interplay of viscous and viscoelastic properties of liquid polymers on hydrophobically coated substrates, there is a need for refined theoretical methods that are able to capture and evolve the emerging morphologies and their longtime dynamics. Dimension-reduced thin film models have shown in the past to be extremly successful to enable quantitative predictions that are hardly being attained simply via the underlying free-boundary problem.In this paper we make an important step in that direction by developing a new thin film model that combines large slippage with viscoelastic properties. In [2] a family of thin film models ranging from weak to strong slip regimes could be derived depending on the order of magnitude of the slip length (see also [3]). Modelling the viscoelastic properties of such polymers generalized Maxwell and Jeffreys models have been widely used. In [4] a weak slip model could be combined with the linearized Jeffreys model to discuss effects of viscoelastic relaxation. More recently [5] we have shown that the strong slip limit can also be recovered for the linear Jeffreys model. In this Rapid Note we show that, for the strong slip regime, we are able to fully incorporate the general corotational Jeffreys model into our thin film model.We begin by presenting the underlying free boundary problem for incompressible, viscous flow with velocity u,where we assume that the traceless part of the symmetric stress tensor τ obeys the corotational Jeffreys model [6]Here, D/Dt denotes the Jaumann derivative which for arbitrary tensor fields Λ is given bywhereγ and ω denote the rate of strain tensor and the vorticity tensor, given bẏrespectively; d/dt is the material derivative ∂ t + u · ∇. We assume the viscosity µ as well as the relaxation parameters λ 1 , λ 2 to be constant material parameters.