fold - fold the ends inside
ierr=fold_lop(adj,add,x,y)
Standard operator interface
Fold portion of y outside x back in on itself
the triangle manpage
geef90