We address the issue of nonlinear modes in a two-dimensional waveguide array, spatially distributed in the Lieb lattice geometry, and modeled by a saturable non-linear Schrodinger equation. In particular, we analyze the ...