IOporosityModelList pZones(mesh);