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

Fix SOIL restart initialization

Browse Source
This commit is contained in:
Tor Harald Sandve 2017-01-05 10:02:50 +01:00
parent edf883e747
commit 5990201ee9
1 changed files with 13 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
opm/autodiff/Compat.hpp
View File

@ -120,6 +120,7 @@ inline void solutionToSim( const data::Solution& sol,
SimulationDataContainer& state ) { SimulationDataContainer& state ) {
const auto stride = phases.num_phases; const auto stride = phases.num_phases;
if( sol.has( "SWAT" ) ) { if( sol.has( "SWAT" ) ) {
stripe( sol.data( "SWAT" ), stripe( sol.data( "SWAT" ),
stride, stride,
@ -134,6 +135,18 @@ inline void solutionToSim( const data::Solution& sol,
state.saturation() ); state.saturation() );
} }
for (size_t c = 0; c < state.numCells(); ++c) {
double& so = state.saturation()[phases.num_phases*c + phases.phase_pos[ BlackoilPhases::Liquid ]];
so = 1.0;
if (phases.phase_used[ BlackoilPhases::Aqua]) {
so -= state.saturation()[phases.num_phases*c + phases.phase_pos[ BlackoilPhases::Aqua ]];
}
if (phases.phase_used[ BlackoilPhases::Vapour]) {
so -= state.saturation()[phases.num_phases*c + phases.phase_pos[ BlackoilPhases::Vapour ]];
}
}
if( sol.has( "PRESSURE" ) ) { if( sol.has( "PRESSURE" ) ) {
state.pressure() = sol.data( "PRESSURE" ); state.pressure() = sol.data( "PRESSURE" );
} }