Had a problem with domain area:
Finding the area of the bounding box of pore coordinates was calculating the rectangle of outermost pore centers. We miss one Lc on a cubic.
I wanted to add pore radius to this, but boundary pores have r=0.
So, I'm finding which plane the face is on: x-normal, y-normal, or z-normal. Then I'm finding the area of the bounding box of the entire network.