Meshing for non-contiguous phases

[fussegli]

Currently, if the domain to be meshed has non-contiguous phases (i.e., isolated clusters), the meshed region are assigned to the label 0. The post-processing used to re-assign such cells is perfectible. I was not aware this issue was related to connectivity, now that i know what produce the region label 0, I can make a better post-processing.