Add functionality to handle a vector of region choices

By using RegionChoice=c('Inner', 'Outer'), regions that are not solely bordered by External, but also by other regions, will now contain all surfaces shared with other regions as well.