We give a combinatorial evaluation of Iwahori Whittaker functions for unramified genuine principal series representations on metaplectic covers of the general linear group over a non-archimedean local field. To describe the combinatorics, we introduce new combinatorial data that we call colored data: colored Lusztig data, colored Gelfand-Tsetlin patterns, and colored lattice models. We show that all three are equivalent. To achieve the result, we give an explicit Iwahori decomposition for the maximal unipotent subgroup of a split reductive group which gives the parametrization of the generalized Mirković-Vilonen cycles in the affine flag varieties and is of interest in itself. Our result is based and naturally extends Peter McNamara's evaluation of the metaplectic spherical Whittaker function in terms of Lusztig data.