diff --git a/opm/autodiff/GeoProps.hpp b/opm/autodiff/GeoProps.hpp index e47f497b9..dbb00277d 100644 --- a/opm/autodiff/GeoProps.hpp +++ b/opm/autodiff/GeoProps.hpp @@ -27,7 +27,7 @@ #include #include - +#include #include #include @@ -90,6 +90,10 @@ namespace Opm ntg = eclState->getDoubleGridProperty("NTG")->getData(); } + // get grid from parser. + + // Get original grid cell volume. + EclipseGridConstPtr eclgrid = eclState->getEclipseGrid(); // Pore volume for (int cellIdx = 0; cellIdx < numCells; ++cellIdx) { int cartesianCellIdx = AutoDiffGrid::globalCell(grid)[cellIdx]; @@ -97,7 +101,8 @@ namespace Opm props.porosity()[cellIdx] * multpv[cartesianCellIdx] * ntg[cartesianCellIdx] - * AutoDiffGrid::cellVolume(grid, cellIdx); + * eclgrid->getCellVolume(cartesianCellIdx); + // * AutoDiffGrid::cellVolume(grid, cellIdx); } // Transmissibility