27#ifndef OPM_VTK_TRACER_MODULE_HPP
28#define OPM_VTK_TRACER_MODULE_HPP
30#include <dune/common/fvector.hh>
55 template <
class TypeTag>
68 static constexpr int vtkFormat = getPropValue<TypeTag, Properties::VtkOutputFormat>();
85 Parameters::Register<Parameters::VtkWriteTracerConcentration>
86 (
"Include the tracer concentration in the VTK output files");
95 if (eclTracerConcentrationOutput_()){
96 const auto& tracerModel = this->
simulator_.problem().tracerModel();
97 eclFreeTracerConcentration_.resize(tracerModel.numTracers());
98 eclSolTracerConcentration_.resize(tracerModel.numTracers());
99 const auto& enableSolTracers = tracerModel.enableSolTracers();
101 for (std::size_t tracerIdx = 0; tracerIdx < eclFreeTracerConcentration_.size(); ++tracerIdx) {
103 if (enableSolTracers[tracerIdx])
116 if (!Parameters::Get<Parameters::EnableVtkOutput>()) {
120 if (eclTracerConcentrationOutput_()) {
121 const auto& tracerModel = elemCtx.problem().tracerModel();
122 const auto& enableSolTracers = tracerModel.enableSolTracers();
124 for (std::size_t tracerIdx = 0; tracerIdx < eclFreeTracerConcentration_.size(); ++tracerIdx) {
126 for (
unsigned dofIdx = 0; dofIdx < elemCtx.numPrimaryDof(0); ++dofIdx) {
127 unsigned globalDofIdx = elemCtx.globalSpaceIndex(dofIdx, 0);
128 eclFreeTracerConcentration_[tracerIdx][globalDofIdx] = tracerModel.freeTracerConcentration(tracerIdx, globalDofIdx);
131 if (enableSolTracers[tracerIdx]) {
132 for (
unsigned dofIdx = 0; dofIdx < elemCtx.numPrimaryDof(0); ++dofIdx) {
133 unsigned globalDofIdx = elemCtx.globalSpaceIndex(dofIdx, 0);
134 eclSolTracerConcentration_[tracerIdx][globalDofIdx] = tracerModel.solTracerConcentration(tracerIdx, globalDofIdx);
150 if (eclTracerConcentrationOutput_()){
151 const auto& tracerModel = this->
simulator_.problem().tracerModel();
152 const auto& enableSolTracers = tracerModel.enableSolTracers();
154 for (std::size_t tracerIdx = 0; tracerIdx < eclFreeTracerConcentration_.size(); ++tracerIdx) {
155 const std::string tmp =
"freeTracerConcentration_" + tracerModel.name(tracerIdx);
156 this->
commitScalarBuffer_(baseWriter, tmp.c_str(), eclFreeTracerConcentration_[tracerIdx]);
157 if (enableSolTracers[tracerIdx]) {
158 const std::string tmp2 =
"solTracerConcentration_" + tracerModel.name(tracerIdx);
159 this->
commitScalarBuffer_(baseWriter, tmp2.c_str(), eclSolTracerConcentration_[tracerIdx]);
166 static bool eclTracerConcentrationOutput_()
168 static bool val = Parameters::Get<Parameters::VtkWriteTracerConcentration>();
172 std::vector<ScalarBuffer> eclFreeTracerConcentration_;
173 std::vector<ScalarBuffer> eclSolTracerConcentration_;
Declares the properties required by the black oil model.
The base class for writer modules.
Definition: baseoutputmodule.hh:67
BaseOutputWriter::ScalarBuffer ScalarBuffer
Definition: baseoutputmodule.hh:85
const Simulator & simulator_
Definition: baseoutputmodule.hh:434
void commitScalarBuffer_(BaseOutputWriter &baseWriter, const char *name, ScalarBuffer &buffer, BufferType bufferType=DofBuffer)
Add a buffer containing scalar quantities to the result file.
Definition: baseoutputmodule.hh:244
void resizeScalarBuffer_(ScalarBuffer &buffer, BufferType bufferType=DofBuffer)
Allocate the space for a buffer storing a scalar quantity.
Definition: baseoutputmodule.hh:156
The base class for all output writers.
Definition: baseoutputwriter.hh:44
Simplifies writing multi-file VTK datasets.
Definition: vtkmultiwriter.hh:66
VTK output module for the tracer model's parameters.
Definition: VtkTracerModule.hpp:57
static void registerParameters()
Register all run-time parameters for the tracer VTK output module.
Definition: VtkTracerModule.hpp:83
void commitBuffers(BaseOutputWriter &baseWriter)
Add all buffers to the VTK output writer.
Definition: VtkTracerModule.hpp:144
void processElement(const ElementContext &elemCtx)
Modify the internal buffers according to the intensive quantities relevant for an element.
Definition: VtkTracerModule.hpp:114
VtkTracerModule(const Simulator &simulator)
Definition: VtkTracerModule.hpp:75
void allocBuffers()
Allocate memory for the scalar fields we would like to write to the VTK file.
Definition: VtkTracerModule.hpp:93
Declare the properties used by the infrastructure code of the finite volume discretizations.
Definition: blackoilnewtonmethodparams.hpp:31
Definition: blackoilboundaryratevector.hh:37
typename Properties::Detail::GetPropImpl< TypeTag, Property >::type::type GetPropType
get the type alias defined in the property (equivalent to old macro GET_PROP_TYPE(....
Definition: propertysystem.hh:235
This file provides the infrastructure to retrieve run-time parameters.
The Opm property system, traits with inheritance.
Definition: VtkTracerModule.hpp:45
static constexpr bool value
Definition: VtkTracerModule.hpp:45