flank.region function

Get adjacent flanks from regions