[appl] add average grid bounding box

appl/porescale/channel/main.cc

@@ -190,9 +190,11 @@ int main(int argc, char** argv)
     // Make a grid for the volume averaged elements
     using AveragedGrid = Dune::YaspGrid<2>;
     DimensionPair numVolumes = getParam<DimensionPair>("VolumeAveraging.NumVolumes");
+    std::pair<GlobalPosition, GlobalPosition> averageBox = {getParam<GlobalPosition>("VolumeAveraging.AverageBoxMin"),
+                                                            getParam<GlobalPosition>("VolumeAveraging.AverageBoxMax")};
     auto averageGrid = Dune::StructuredGridFactory<AveragedGrid>::createCubeGrid(
-                                                   {porousMediaBox.first},
-                                                   {porousMediaBox.second},
+                                                   {averageBox.first},
+                                                   {averageBox.second},
                                                    {numVolumes[0], numVolumes[1]});
     using AveragedGridGeometry = CCTpfaFVGridGeometry<typename AveragedGrid::LeafGridView>;
     auto averagedGridGeometry = std::make_shared<AveragedGridGeometry>(averageGrid->leafGridView());

appl/porescale/channel/params.input

@@ -26,6 +26,8 @@ GridInputFileName = gridfiles/PoreGeometry_Random-1.csv # gridfiles/*.csv
 
 [VolumeAveraging]
 NumVolumes = 10 5
+AverageBoxMin = 0.0 0.0
+AverageBoxMax = 1.0 0.5
 NumPoreScaleCellsPerBlock = 20
 
 [Newton]

appl/porescale/filter/main.cc

@@ -193,9 +193,11 @@ int main(int argc, char** argv)
     // Make a grid for the volume averaged elements
     using AveragedGrid = Dune::YaspGrid<2>;
     DimensionPair numVolumes = getParam<DimensionPair>("VolumeAveraging.NumVolumes");
+    std::pair<GlobalPosition, GlobalPosition> averageBox = {getParam<GlobalPosition>("VolumeAveraging.AverageBoxMin"),
+                                                            getParam<GlobalPosition>("VolumeAveraging.AverageBoxMax")};
     auto averageGrid = Dune::StructuredGridFactory<AveragedGrid>::createCubeGrid(
-                                                   {porousMediaBox.first},
-                                                   {porousMediaBox.second},
+                                                   {averageBox.first},
+                                                   {averageBox.second},
                                                    {numVolumes[0], numVolumes[1]});
     using AveragedGridGeometry = CCTpfaFVGridGeometry<typename AveragedGrid::LeafGridView>;
     auto averagedGridGeometry = std::make_shared<AveragedGridGeometry>(averageGrid->leafGridView());

appl/porescale/filter/params.input

@@ -26,6 +26,8 @@ GridInputFileName = gridfiles/PoreGeometry_Random-1.csv # gridfiles/*.csv
 
 [VolumeAveraging]
 NumVolumes = 10 5
+AverageBoxMin = 0.0 0.0
+AverageBoxMax = 1.0 0.5
 NumPoreScaleCellsPerBlock = 20
 
 [Newton]