lo = 2*(j-1) + (i+1)/2
IF ( facet_dims(lo).EQ.0 .AND. facet_dims(ll).GT.0 ) THEN
facet_dims(lo) = facet_dims(ll)
& ' facet',j,'.',edge(i), ' set dim', lo, ' = dim', ll,
INTEGER i, j, ii, jj, lo, ll
lo = facet_dims( lo )
IF ( lo.NE.ll ) THEN
& edge(i), '.Edge facet #', j, ' (=',lo,')'
facet_oi(i,j) = lo+1
facet_oj(i,j) = lo+1
INTEGER i, j, ii, jj, i1, j1, k, lo, ll