From 5aede81b99498c2bb77e66ebe186d1edc96ba70e Mon Sep 17 00:00:00 2001
From: Joakim Hove <joakim.hove@gmail.com>
Date: Tue, 23 Nov 2021 12:56:56 +0100
Subject: [PATCH] Add function to set tracer concentration in a cell

---
 ebos/eclgenerictracermodel.cc | 9 +++++++++
 ebos/eclgenerictracermodel.hh | 1 +
 2 files changed, 10 insertions(+)

diff --git a/ebos/eclgenerictracermodel.cc b/ebos/eclgenerictracermodel.cc
index 61ae6fcc3..4081e981d 100644
--- a/ebos/eclgenerictracermodel.cc
+++ b/ebos/eclgenerictracermodel.cc
@@ -126,6 +126,15 @@ tracerConcentration(int tracerIdx, int globalDofIdx) const
     return tracerConcentration_[tracerIdx][globalDofIdx];
 }
 
+
+template<class Grid,class GridView, class DofMapper, class Stencil, class Scalar>
+void EclGenericTracerModel<Grid,GridView,DofMapper,Stencil,Scalar>::
+setTracerConcentration(int tracerIdx, int globalDofIdx, Scalar value)
+{
+    this->tracerConcentration_[tracerIdx][globalDofIdx] = value;
+}
+
+
 template<class Grid,class GridView, class DofMapper, class Stencil, class Scalar>
 void EclGenericTracerModel<Grid,GridView,DofMapper,Stencil,Scalar>::
 doInit(bool enabled, size_t numGridDof,
diff --git a/ebos/eclgenerictracermodel.hh b/ebos/eclgenerictracermodel.hh
index aa02163e5..42dc4c6bf 100644
--- a/ebos/eclgenerictracermodel.hh
+++ b/ebos/eclgenerictracermodel.hh
@@ -67,6 +67,7 @@ public:
      * \brief Return the tracer concentration for tracer index and global DofIdx
      */
     Scalar tracerConcentration(int tracerIdx, int globalDofIdx) const;
+    void setTracerConcentration(int tracerIdx, int globalDofIdx, Scalar value);
 
     /*!
     * \brief Return well tracer rates