From 5401baf5332a32980854df8b20ac4f56e543d409 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Wed, 24 Apr 2013 10:40:56 +0200
Subject: [PATCH] Add call to DG tracer computations if user requests.

Used to throw with an 'unimplemented' message.
---
 examples/compute_tof_from_files.cpp | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/examples/compute_tof_from_files.cpp b/examples/compute_tof_from_files.cpp
index 723df3783..a6580c2f8 100644
--- a/examples/compute_tof_from_files.cpp
+++ b/examples/compute_tof_from_files.cpp
@@ -148,9 +148,6 @@ main(int argc, char** argv)
     } else {
         use_multidim_upwind = param.getDefault("use_multidim_upwind", false);
     }
-    if (use_dg && compute_tracer) {
-        THROW("DG for tracer not yet implemented.");
-    }
 
     // Write parameters used for later reference.
     bool output = param.getDefault("output", true);
@@ -178,7 +175,11 @@ main(int argc, char** argv)
     std::vector<double> tof;
     std::vector<double> tracer;
     if (use_dg) {
-        dg_solver->solveTof(&flux[0], &porevol[0], &src[0], tof);
+        if (compute_tracer) {
+            dg_solver->solveTofTracer(&flux[0], &porevol[0], &src[0], tracerheads, tof, tracer);
+        } else {
+            dg_solver->solveTof(&flux[0], &porevol[0], &src[0], tof);
+        }
     } else {
         Opm::TofReorder tofsolver(grid, use_multidim_upwind);
         if (compute_tracer) {