OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-25 20:16:25 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #4309 from atgeirr/write-global-cell-index-to-damaris

Browse Source 
Write global cell index to damaris
This commit is contained in:
Atgeirr Flø Rasmussen 2022-12-07 15:38:40 +01:00 committed by GitHub
commit e6a8686c77
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
4 changed files with 17 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
ebos/collecttoiorank.hh
View File

@ -110,7 +110,12 @@ public:
{ return toIORankComm_.size() > 1; }
int localIdxToGlobalIdx(unsigned localIdx) const;
const std::vector<int>& localIdxToGlobalIdxMapping() const
{
return localIdxToGlobalIdx_;
}
bool doesNeedReordering() const
{ return needsReordering;}

5
ebos/eclwriter.hh
View File

@ -309,8 +309,11 @@ public:
// N.B. damarisUpdate_ should be set to true if at any time the model geometry changes
if (this->damarisUpdate_) {
const auto& gridView = simulator_.gridView();
const int numElements = gridView.size(/*codim=*/0); // This it is the local ranks model size.
const auto& interior_elements = elements(gridView, Dune::Partitions::interior);
const int numElements = std::distance(interior_elements.begin(), interior_elements.end());
Opm::DamarisOutput::setupDamarisWritingPars(simulator_.vanguard().grid().comm(), numElements);
const std::vector<int>& local_to_global = this->collectToIORank_.localIdxToGlobalIdxMapping();
damaris_write("GLOBAL_CELL_INDEX", local_to_global.data());
// By default we assume static grid
this->damarisUpdate_ = false;
}

5
opm/simulators/utils/DamarisOutputModule.cpp
View File

@ -137,5 +137,10 @@ setupDamarisWritingPars(Parallel::Communication comm, const int n_elements_local
OpmLog::error("Damaris library produced an error result for "
"damaris_set_position(\"PRESSURE\", temp_int64_t);");
}
damaris_err = damaris_set_position("GLOBAL_CELL_INDEX", temp_int64_t);
if (damaris_err != DAMARIS_OK && rank == 0) {
OpmLog::error("Damaris library produced an error result for "
"damaris_set_position(\"GLOBAL_CELL_INDEX\", temp_int64_t);");
}
}
} // namespace Opm::DamarisOutput

2
opm/simulators/utils/initDamarisXmlFile.cpp
View File

@ -47,6 +47,8 @@ std::string initDamarisXmlFile()
<parameter name="n_elements_local" type="int" value="1" />
<parameter name="n" type="int" value="1" />
<layout name="zonal_layout_usmesh_integer" type="int" dimensions="n_elements_local" global="n_elements_total" comment="For the field data e.g. Pressure" />
<variable name="GLOBAL_CELL_INDEX" layout="zonal_layout_usmesh_integer" type="scalar" visualizable="false" time-varying="false" centering="zonal" />
<layout name="zonal_layout_usmesh" type="double" dimensions="n_elements_local" global="n_elements_total" comment="For the field data e.g. Pressure" />
<variable name="PRESSURE" layout="zonal_layout_usmesh" type="scalar" visualizable="false" unit="Pa" centering="zonal" store="_MYSTORE_OR_EMPTY_REGEX_" />
_MORE_VARIABLES_REGEX_