ienode(ne,2):: ienode(i,1->2) node numbers at each end of element edge i 

iec(ne,2) :: nbe(i,1->2) cell number of cell(s) connected to edge i


