diff --git a/examples/spu_2p.cpp b/examples/spu_2p.cpp index 1a584134..f08aa20a 100644 --- a/examples/spu_2p.cpp +++ b/examples/spu_2p.cpp @@ -50,6 +50,8 @@ #include #include +#include +#include #include #include #include @@ -324,6 +326,7 @@ main(int argc, char** argv) bool use_deck = param.has("deck_filename"); boost::scoped_ptr grid; boost::scoped_ptr props; + boost::scoped_ptr wells; if (use_deck) { std::string deck_filename = param.get("deck_filename"); Opm::EclipseGridParser deck(deck_filename); @@ -333,6 +336,8 @@ main(int argc, char** argv) const int* gc = grid->c_grid()->global_cell; std::vector global_cell(gc, gc + grid->c_grid()->number_of_cells); props.reset(new Opm::IncompPropertiesFromDeck(deck, global_cell)); + // Wells init. + wells.reset(new Opm::WellsManager(deck, *grid->c_grid(), props->permeability())); } else { // Grid init. const int nx = param.getDefault("nx", 100); @@ -344,6 +349,8 @@ main(int argc, char** argv) grid.reset(new Opm::GridManager(nx, ny, nz, dx, dy, dz)); // Rock and fluid init. props.reset(new Opm::IncompPropertiesBasic(param, grid->c_grid()->dimensions, grid->c_grid()->number_of_cells)); + // Wells init. + wells.reset(new Opm::WellsManager()); } // Extra rock init.