Inforce homogeneous Neumann BCs using DG FEM

General Discussion
4

You can just multiply by (nTonEdge-1) to delete the boundary edges,
see About intalledges and internal edges - #2 by fb77