If one builds GWP basis functions of order 1, then the shape functions associated with the face are not added to the list functions. Line 150 in gwpglobal.jl, it should skip there only for the lowest order (order<1).