77 std::string fileformat = param.getDefault<std::string>(
"fileformat",
"cartesian");
78 if (fileformat ==
"eclipse") {
79 std::string ecl_file = param.get<std::string>(
"filename");
82 auto deck = parser.parseFile(ecl_file);
83 if (param.has(
"z_tolerance")) {
84 std::cerr <<
"****** Warning: z_tolerance parameter is obsolete, use PINCH in deck input instead\n";
86 bool periodic_extension = param.getDefault<
bool>(
"periodic_extension",
false);
87 bool clip_z = param.getDefault<
bool>(
"clip_z",
false);
88 bool turn_normals = param.getDefault<
bool>(
"turn_normals",
false);
90 Opm::EclipseGrid inputGrid(deck);
91 grid.processEclipseFormat(&inputGrid,
nullptr, periodic_extension, turn_normals, clip_z,
true);
94 if (param.getDefault(
"output_ecl",
false)) {
95 OPM_THROW(std::runtime_error,
"Saving to EGRID files is not yet implemented");
103 double perm_threshold_md = param.getDefault(
"perm_threshold_md", 0.0);
104 double perm_threshold = Opm::unit::convert::from(perm_threshold_md, Opm::prefix::milli*Opm::unit::darcy);
105 std::string rock_list = param.getDefault<std::string>(
"rock_list",
"no_list");
106 std::string* rl_ptr = (rock_list ==
"no_list") ? 0 : &rock_list;
107 bool use_j = param.getDefault(
"use_jfunction_scaling",
useJ<ResProp<3> >());
111 sigma = param.getDefault(
"sigma", sigma);
112 theta = param.getDefault(
"theta", theta);
114 if (param.has(
"viscosity1") || param.has(
"viscosity2")) {
115 double v1 = param.getDefault(
"viscosity1", 0.001);
116 double v2 = param.getDefault(
"viscosity2", 0.003);
117 res_prop.setViscosities(v1, v2);
119 res_prop.init(deck, grid.globalCell(), perm_threshold, rl_ptr,
120 use_j, sigma, theta);
121 }
else if (fileformat ==
"cartesian") {
122 std::array<int, 3> dims = {{ param.getDefault<
int>(
"nx", 1),
123 param.getDefault<
int>(
"ny", 1),
124 param.getDefault<
int>(
"nz", 1) }};
125 std::array<double, 3> cellsz = {{ param.getDefault<
double>(
"dx", 1.0),
126 param.getDefault<
double>(
"dy", 1.0),
127 param.getDefault<
double>(
"dz", 1.0) }};
128 grid.createCartesian(dims, cellsz);
129 double default_poro = param.getDefault(
"default_poro", 0.2);
130 double default_perm_md = param.getDefault(
"default_perm_md", 100.0);
131 double default_perm = Opm::unit::convert::from(default_perm_md, Opm::prefix::milli*Opm::unit::darcy);
132 OPM_MESSAGE(
"Warning: For generated cartesian grids, we use uniform reservoir properties.");
133 res_prop.init(grid.size(0), default_poro, default_perm);
135 OPM_THROW(std::runtime_error,
136 "Unknown file format string: " + fileformat);
138 if (param.getDefault(
"use_unique_boundary_ids",
false)) {
139 grid.setUniqueBoundaryIds(
true);
148 bool periodic_extension,
152 double perm_threshold,
153 const std::string& rock_list,
154 bool use_jfunction_scaling,
158 ResProp<3>& res_prop)
160 Opm::EclipseGrid eg(deck);
161 const std::string* rl_ptr = (rock_list ==
"no_list") ? 0 : &rock_list;
162 grid.processEclipseFormat(&eg,
nullptr, periodic_extension, turn_normals, clip_z,
true);
163 res_prop.init(deck, grid.globalCell(), perm_threshold, rl_ptr, use_jfunction_scaling, sigma, theta);
165 grid.setUniqueBoundaryIds(
true);
void setupGridAndPropsEclipse(const Opm::Deck &deck, bool periodic_extension, bool turn_normals, bool clip_z, bool unique_bids, double perm_threshold, const std::string &rock_list, bool use_jfunction_scaling, double sigma, double theta, Dune::CpGrid &grid, ResProp< 3 > &res_prop)
Definition setupGridAndProps.hpp:147
void setupGridAndProps(const Opm::ParameterGroup ¶m, Dune::CpGrid &grid, ResProp< 3 > &res_prop)
Definition setupGridAndProps.hpp:71