New parameter 'injected_porevolumes_per_day' implemented.

- Only used if no wells present.
 - Defaults to 0.0 if gravity is present.
 - Otherwise defaults to 0.1.
These defaults are intended to preserve behaviour from before the
initialisation code changed away from using the 'scenario' parameter.
This commit is contained in:
Atgeirr Flø Rasmussen
2012-03-29 21:16:45 +02:00
parent 5b8d060b1a
commit c53fb7eb15

View File

@@ -391,10 +391,9 @@ main(int argc, char** argv)
if (wells->c_wells()) {
Opm::wellsToSrc(*wells->c_wells(), num_cells, src);
} else {
double flow_per_sec = 0.1*tot_porevol_init/Opm::unit::day;
if (param.has("injection_rate_per_day")) {
flow_per_sec = param.get<double>("injection_rate_per_day")/Opm::unit::day;
}
const double default_injection = use_gravity ? 0.0 : 0.1;
const double flow_per_sec = param.getDefault<double>("injected_porevolumes_per_day", default_injection)
*tot_porevol_init/Opm::unit::day;
src[0] = flow_per_sec;
src[num_cells - 1] = -flow_per_sec;
}