if the subgrid dimension is one, the subgrid method sorts the coordinates of the resulting grid without permuting the CellNodes, BFaceNodes etc. This may lead to unexpected behaviour in situations beyond nodal interpolation and plotting or when other projections than just restricting to one coordinate are applied.
In PR#49 the sorting is switched off if project = false at least.