Mathematical modeling of the melt draw step in film casting is used to understand the impact of viscoelasticity on gage variations. The thin film approximation is used to simplify the governing equations to a set of two-dimensional partial differential equations which are solved using the finite element method. Simulations using the Giesekus constitutive model show that the magnitude of thickness variations increases with increasing flow elasticity or Deborah number.