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

output temperature when TEMP is set

Browse Source
This commit is contained in:
Tor Harald Sandve
2020-09-21 13:58:09 +02:00
parent 6a592a8c55
commit ae95570693
1 changed files with 10 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

17
ebos/ecloutputblackoilmodule.hh
View File

@@ -348,9 +348,6 @@ public:
} }
} }
// always allocate memory for temperature
temperature_.resize(bufferSize, 0.0);
// field data should be allocated // field data should be allocated
// 1) when we want to restart // 1) when we want to restart
// 2) when it is ask for by the user via restartConfig // 2) when it is ask for by the user via restartConfig
@@ -370,6 +367,12 @@ public:
rstKeywords["PRES"] = 0; rstKeywords["PRES"] = 0;
rstKeywords["PRESSURE"] = 0; rstKeywords["PRESSURE"] = 0;
// allocate memory for temperature
if (enableEnergy || rstKeywords["TEMP"]) {
temperature_.resize(bufferSize, 0.0);
rstKeywords["TEMP"] = 0;
}
if (FluidSystem::phaseIsActive(oilPhaseIdx)) if (FluidSystem::phaseIsActive(oilPhaseIdx))
rstKeywords["SOIL"] = 0; rstKeywords["SOIL"] = 0;
if (FluidSystem::phaseIsActive(gasPhaseIdx)) if (FluidSystem::phaseIsActive(gasPhaseIdx))
@@ -560,7 +563,7 @@ public:
Opm::Valgrind::CheckDefined(oilPressure_[globalDofIdx]); Opm::Valgrind::CheckDefined(oilPressure_[globalDofIdx]);
} }
if (enableEnergy) { if (temperature_.size() > 0) {
temperature_[globalDofIdx] = Opm::getValue(fs.temperature(oilPhaseIdx)); temperature_[globalDofIdx] = Opm::getValue(fs.temperature(oilPhaseIdx));
Opm::Valgrind::CheckDefined(temperature_[globalDofIdx]); Opm::Valgrind::CheckDefined(temperature_[globalDofIdx]);
} }
@@ -958,7 +961,7 @@ public:
sol.insert("PRESSURE", Opm::UnitSystem::measure::pressure, std::move(oilPressure_), Opm::data::TargetType::RESTART_SOLUTION); sol.insert("PRESSURE", Opm::UnitSystem::measure::pressure, std::move(oilPressure_), Opm::data::TargetType::RESTART_SOLUTION);
} }
if (enableEnergy) { if (temperature_.size() > 0) {
sol.insert("TEMP", Opm::UnitSystem::measure::temperature, std::move(temperature_), Opm::data::TargetType::RESTART_SOLUTION); sol.insert("TEMP", Opm::UnitSystem::measure::temperature, std::move(temperature_), Opm::data::TargetType::RESTART_SOLUTION);
} }
@@ -1702,7 +1705,7 @@ public:
if (oilPressure_.size() > 0 && sol.has("PRESSURE")) if (oilPressure_.size() > 0 && sol.has("PRESSURE"))
oilPressure_[elemIdx] = sol.data("PRESSURE")[globalDofIndex]; oilPressure_[elemIdx] = sol.data("PRESSURE")[globalDofIndex];
if (enableEnergy && sol.has("TEMP")) if (temperature_.size() > 0 && sol.has("TEMP"))
temperature_[elemIdx] = sol.data("TEMP")[globalDofIndex]; temperature_[elemIdx] = sol.data("TEMP")[globalDofIndex];
if (rs_.size() > 0 && sol.has("RS")) if (rs_.size() > 0 && sol.has("RS"))
rs_[elemIdx] = sol.data("RS")[globalDofIndex]; rs_[elemIdx] = sol.data("RS")[globalDofIndex];
@@ -1756,7 +1759,7 @@ public:
} }
} }
if (enableEnergy) if (temperature_.size() > 0)
fs.setTemperature(temperature_[elemIdx]); fs.setTemperature(temperature_[elemIdx]);
if (rs_.size() > 0) if (rs_.size() > 0)
fs.setRs(rs_[elemIdx]); fs.setRs(rs_[elemIdx]);