making sure the aquifer cells have the correct cell volume

because AQUNUM can specify different cells volume, which is indepedent
of the grid coordinates.

ebos/eclcpgridvanguard.hh

@@ -341,6 +341,11 @@ protected:
     void createGrids_()
     {
         grid_.reset(new Dune::CpGrid());
+
+        std::unordered_map<size_t, double> aquifer_cell_volumes;
+        if (this->eclState().aquifer().hasNumericalAquifer()) {
+            aquifer_cell_volumes = this->eclState().aquifer().numericalAquifers().cellVolumes();
+        }
         grid_->processEclipseFormat(mpiRank == 0 ? &this->eclState().getInputGrid()
                                                  : nullptr,
                                     /*isPeriodic=*/false,
@@ -348,7 +353,8 @@ protected:
                                     /*clipZ=*/false,
                                     mpiRank == 0 ? this->eclState().fieldProps().porv(true)
                                                  : std::vector<double>(),
-                                    this->eclState().getInputNNC());
+                                    this->eclState().getInputNNC(),
+                                    aquifer_cell_volumes);
 
         // we use separate grid objects: one for the calculation of the initial condition
         // via EQUIL and one for the actual simulation. The reason is that the EQUIL code