f = Interpolation[Last /@ Sort /@ GatherBy[l[[All, 1 ;; 2]], #[[1]] &]] ListContourPlot[l, RegionFunction -> (#2 < f[#1] &)] ![Mathematica graphics](https://i.sstatic.net/aKL7v.png)