Support using multiple finite-area regions per volume region

The goal of this work is to have different finite-area shells on the same volume mesh region. These might correspond to different physics etc (eg, a thermal shell on one part, a film elsewhere).