OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Parameter trace_start now controls trace direction.

Browse Source 
Allowable values are Injectors or Producers.
This affects both tof and tracer, and properly switches
flux directions for Producers.
This commit is contained in:
Atgeirr Flø Rasmussen 2015-01-21 20:55:23 +01:00
parent 79ad0c1417
commit 2a8beb1abc
1 changed files with 15 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
examples/compute_tof.cpp
View File

@ -202,14 +202,13 @@ try
use_multidim_upwind = param.getDefault("use_multidim_upwind", false);
}
bool compute_tracer = param.getDefault("compute_tracer", false);
bool trace_injectors = true;
if (compute_tracer) {
std::string trace = param.getDefault<std::string>("trace", "Injectors");
if (trace == "Producers") {
trace_injectors = false;
} else if (trace != "Injectors") {
OPM_THROW(std::runtime_error, "Unknown trace specification (allowed is Injectors or Producers): " << trace);
}
bool start_from_injectors = true;
std::string trace_start = param.getDefault<std::string>("trace_start", "Injectors");
if (trace_start == "Producers") {
start_from_injectors = false;
} else if (trace_start != "Injectors") {
OPM_THROW(std::runtime_error, "Unknown trace_start specification (allowed is Injectors or Producers): " << trace_start);
}
// Write parameters used for later reference.
@ -253,6 +252,13 @@ try
std::cout << "Pressure solver took: " << pt << " seconds." << std::endl;
ptime += pt;
// Turn direction of flux if starting from producers.
if (!start_from_injectors) {
for (auto it = flux.begin(); it != flux.end(); ++it) {
(*it) = -(*it);
}
}
// Process transport sources (to include bdy terms and well flows).
std::vector<double> src(num_cells, 0.0);
std::vector<double> transport_src;
@ -265,7 +271,7 @@ try
std::vector<double> tracer;
Opm::SparseTable<int> tracerheads;
if (compute_tracer) {
buildTracerheadsFromWells(wells, trace_injectors, tracerheads);
buildTracerheadsFromWells(wells, start_from_injectors, tracerheads);
}
if (use_dg) {
if (compute_tracer) {