In this paper we consider Iwahori Whittaker functions on n-fold metaplectic covers G of G(F ) with G a split reductive group over a non-archimedean local field F . For every element φ of a basis of Iwahori Whittaker functions, and for every g ∈ G, we evaluate φ(g) by recurrence relations over the Weyl group using "vector Demazure-Whittaker operators." Specializing to the case of G = GL r , we exhibit a solvable lattice model whose partition function equals φ(g). These models are of a new type associated with the quantum affine super group U q ( gl(r|n)). The recurrence relations on the representation theory side then correspond to solutions to Yang-Baxter equations for the lattice models.