In this work we address the issue of nonlinear modes in a two dimensional waveguide array, spatially distributed in the Lieb lattice geometry, modeled by the discrete nonlinear Schrodinger equation. In particular, we analyzed the existence and stability of vortex-type solutions in this system and we found two main kind of vortex modes, namely the on-site and off-site, ranging from S = 1 to S = 3. We study their stability in function of coupling anisotropy effect, finding different behaviours according to the topological charge of solutions.